1
|
Mutational signatures in GATA3 transcription factor and its DNA binding domain that stimulate breast cancer and HDR syndrome. Sci Rep 2021; 11:22762. [PMID: 34815386 PMCID: PMC8611019 DOI: 10.1038/s41598-021-01832-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2021] [Accepted: 10/14/2021] [Indexed: 11/08/2022] Open
Abstract
Transcription factors (TFs) play important roles in many biochemical processes. Many human genetic disorders have been associated with mutations in the genes encoding these transcription factors, and so those mutations became targets for medications and drug design. In parallel, since many transcription factors act either as tumor suppressors or oncogenes, their mutations are mostly associated with cancer. In this perspective, we studied the GATA3 transcription factor when bound to DNA in a crystal structure and assessed the effect of different mutations encountered in patients with different diseases and phenotypes. We generated all missense mutants of GATA3 protein and DNA within the adjacent and the opposite GATA3:DNA complex models. We mutated every amino acid and studied the new binding of the complex after each mutation. Similarly, we did for every DNA base. We applied Poisson-Boltzmann electrostatic calculations feeding into free energy calculations. After analyzing our data, we identified amino acids and DNA bases keys for binding. Furthermore, we validated those findings against experimental genetic data. Our results are the first to propose in silico modeling for GATA:DNA bound complexes that could be used to score effects of missense mutations in other classes of transcription factors involved in common and genetic diseases.
Collapse
|
2
|
Sequeiros-Borja CE, Surpeta B, Brezovsky J. Recent advances in user-friendly computational tools to engineer protein function. Brief Bioinform 2021; 22:bbaa150. [PMID: 32743637 PMCID: PMC8138880 DOI: 10.1093/bib/bbaa150] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2020] [Revised: 06/03/2020] [Accepted: 06/16/2020] [Indexed: 12/14/2022] Open
Abstract
Progress in technology and algorithms throughout the past decade has transformed the field of protein design and engineering. Computational approaches have become well-engrained in the processes of tailoring proteins for various biotechnological applications. Many tools and methods are developed and upgraded each year to satisfy the increasing demands and challenges of protein engineering. To help protein engineers and bioinformaticians navigate this emerging wave of dedicated software, we have critically evaluated recent additions to the toolbox regarding their application for semi-rational and rational protein engineering. These newly developed tools identify and prioritize hotspots and analyze the effects of mutations for a variety of properties, comprising ligand binding, protein-protein and protein-nucleic acid interactions, and electrostatic potential. We also discuss notable progress to target elusive protein dynamics and associated properties like ligand-transport processes and allosteric communication. Finally, we discuss several challenges these tools face and provide our perspectives on the further development of readily applicable methods to guide protein engineering efforts.
Collapse
Affiliation(s)
- Carlos Eduardo Sequeiros-Borja
- Laboratory of Biomolecular Interactions and Transport, Department of Gene Expression, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University and the International Institute of Molecular and Cell Biology in Warsaw, Warsaw, Poland
| | - Bartłomiej Surpeta
- Laboratory of Biomolecular Interactions and Transport, Department of Gene Expression, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University and the International Institute of Molecular and Cell Biology in Warsaw, Warsaw, Poland
| | - Jan Brezovsky
- Laboratory of Biomolecular Interactions and Transport, Department of Gene Expression, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University and the International Institute of Molecular and Cell Biology in Warsaw
| |
Collapse
|
3
|
Narkhede YB, Gautam AK, Hsu RV, Rodriguez W, Zewde NT, Harrison RES, Arantes PR, Gaieb Z, Gorham RD, Kieslich C, Morikis D, Sahu A, Palermo G. Role of Electrostatic Hotspots in the Selectivity of Complement Control Proteins Toward Human and Bovine Complement Inhibition. Front Mol Biosci 2021; 8:618068. [PMID: 33829039 PMCID: PMC8020814 DOI: 10.3389/fmolb.2021.618068] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2020] [Accepted: 01/08/2021] [Indexed: 11/13/2022] Open
Abstract
Poxviruses are dangerous pathogens, which can cause fatal infection in unvaccinated individuals. The causative agent of smallpox in humans, variola virus, is closely related to the bovine vaccinia virus, yet the molecular basis of their selectivity is currently incompletely understood. Here, we examine the role of the electrostatics in the selectivity of the smallpox protein SPICE and vaccinia protein VCP toward the human and bovine complement protein C3b, a key component of the complement immune response. Electrostatic calculations, in-silico alanine-scan and electrostatic hotspot analysis, as introduced by Kieslich and Morikis (PLoS Comput. Biol. 2012), are used to assess the electrostatic complementarity and to identify sites resistant to local perturbation where the electrostatic potential is likely to be evolutionary conserved. The calculations suggest that the bovine C3b is electrostatically prone to selectively bind its VCP ligand. On the other hand, the human isoform of C3b exhibits a lower electrostatic complementarity toward its SPICE ligand. Yet, the human C3b displays a highly preserved electrostatic core, which suggests that this isoform could be less selective in binding different ligands like SPICE and the human Factor H. This is supported by experimental cofactor activity assays revealing that the human C3b is prone to bind both SPICE and Factor H, which exhibit diverse electrostatic properties. Additional investigations considering mutants of SPICE and VCP that revert their selectivity reveal an “electrostatic switch” into the central modules of the ligands, supporting the critical role of the electrostatics in the selectivity. Taken together, these evidences provide insights into the selectivity mechanism of the complement regulator proteins encoded by the variola and vaccinia viruses to circumvent the complement immunity and exert their pathogenic action. These fundamental aspects are valuable for the development of novel vaccines and therapeutic strategies.
Collapse
Affiliation(s)
- Yogesh B Narkhede
- Department of Bioengineering, University of California, Riverside, CA
| | - Avneesh K Gautam
- National Centre for Cell Science, Pune University Campus, Ganeshkhind, India
| | - Rohaine V Hsu
- Department of Bioengineering, University of California, Riverside, CA
| | - Wilson Rodriguez
- Department of Bioengineering, University of California, Riverside, CA
| | - Nehemiah T Zewde
- Department of Bioengineering, University of California, Riverside, CA
| | - Reed E S Harrison
- Department of Bioengineering, University of California, Riverside, CA
| | - Pablo R Arantes
- Department of Bioengineering, University of California, Riverside, CA
| | - Zied Gaieb
- Department of Bioengineering, University of California, Riverside, CA
| | - Ronald D Gorham
- Department of Bioengineering, University of California, Riverside, CA
| | - Chris Kieslich
- Department of Bioengineering, University of California, Riverside, CA.,Department of Chemical Engineering, Auburn University, Auburn, AL
| | - Dimitrios Morikis
- Department of Bioengineering, University of California, Riverside, CA
| | - Arvind Sahu
- National Centre for Cell Science, Pune University Campus, Ganeshkhind, India
| | - Giulia Palermo
- Department of Bioengineering, University of California, Riverside, CA.,Department of Chemistry, University of California, Riverside, CA
| |
Collapse
|
4
|
AESOP: A Python Library for Investigating Electrostatics in Protein Interactions. Biophys J 2017; 112:1761-1766. [PMID: 28494947 DOI: 10.1016/j.bpj.2017.04.005] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2016] [Revised: 02/07/2017] [Accepted: 04/06/2017] [Indexed: 11/24/2022] Open
Abstract
Electric fields often play a role in guiding the association of protein complexes. Such interactions can be further engineered to accelerate complex association, resulting in protein systems with increased productivity. This is especially true for enzymes where reaction rates are typically diffusion limited. To facilitate quantitative comparisons of electrostatics in protein families and to describe electrostatic contributions of individual amino acids, we previously developed a computational framework called AESOP. We now implement this computational tool in Python with increased usability and the capability of performing calculations in parallel. AESOP utilizes PDB2PQR and Adaptive Poisson-Boltzmann Solver to generate grid-based electrostatic potential files for protein structures provided by the end user. There are methods within AESOP for quantitatively comparing sets of grid-based electrostatic potentials in terms of similarity or generating ensembles of electrostatic potential files for a library of mutants to quantify the effects of perturbations in protein structure and protein-protein association.
Collapse
|
5
|
López-Perrote A, Harrison RES, Subías M, Alcorlo M, Rodríguez de Córdoba S, Morikis D, Llorca O. Ionic tethering contributes to the conformational stability and function of complement C3b. Mol Immunol 2017; 85:137-147. [PMID: 28254726 DOI: 10.1016/j.molimm.2016.12.015] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2016] [Revised: 12/13/2016] [Accepted: 12/15/2016] [Indexed: 11/28/2022]
Abstract
C3b, the central component of the alternative pathway (AP) of the complement system, coexists as a mixture of conformations in solution. These conformational changes can affect interactions with other proteins and complement regulators. Here we combine a computational model for electrostatic interactions within C3b with molecular imaging to study the conformation of C3b. The computational analysis shows that the TED domain in C3b is tethered ionically to the macroglobulin (MG) ring. Monovalent counterion concentration affects the magnitude of electrostatic forces anchoring the TED domain to the rest of the C3b molecule in a thermodynamic model. This is confirmed by observing NaCl concentration dependent conformational changes using single molecule electron microscopy (EM). We show that the displacement of the TED domain is compatible with C3b binding to Factor B (FB), suggesting that the regulation of the C3bBb convertase could be affected by conditions that promote movement in the TED domain. Our molecular model also predicts mutations that could alter the positioning of the TED domain, including the common R102G polymorphism, a risk variant for developing age-related macular degeneration. The common C3b isoform, C3bS, and the risk isoform, C3bF, show distinct energetic barriers to displacement in the TED that are related to a network of electrostatic interactions at the interface of the TED and MG-ring domains of C3b. These computational predictions agree with experimental evidence that shows differences in conformation observed in C3b isoforms purified from homozygous donors. Altogether, we reveal an ionic, reversible attachment of the TED domain to the MG ring that may influence complement regulation in some mutations and polymorphisms of C3b.
Collapse
Affiliation(s)
- Andrés López-Perrote
- Centro de Investigaciones Biológicas, Spanish National Research Council (CSIC), Madrid, Spain
| | - Reed E S Harrison
- Department of Bioengineering, University of California, Riverside, CA 92521, USA
| | - Marta Subías
- Centro de Investigaciones Biológicas, Spanish National Research Council (CSIC), Madrid, Spain; Centro de Investigación Biomédica en Enfermedades Raras, Madrid, Spain
| | - Martín Alcorlo
- Centro de Investigaciones Biológicas, Spanish National Research Council (CSIC), Madrid, Spain
| | - Santiago Rodríguez de Córdoba
- Centro de Investigaciones Biológicas, Spanish National Research Council (CSIC), Madrid, Spain; Centro de Investigación Biomédica en Enfermedades Raras, Madrid, Spain.
| | - Dimitrios Morikis
- Department of Bioengineering, University of California, Riverside, CA 92521, USA.
| | - Oscar Llorca
- Centro de Investigaciones Biológicas, Spanish National Research Council (CSIC), Madrid, Spain.
| |
Collapse
|
6
|
El-Assaad A, Dawy Z, Nemer G. Electrostatic study of Alanine mutational effects on transcription: application to GATA-3:DNA interaction complex. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2016; 2015:4005-8. [PMID: 26737172 DOI: 10.1109/embc.2015.7319272] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Protein-DNA interaction is of fundamental importance in molecular biology, playing roles in functions as diverse as DNA transcription, DNA structure formation, and DNA repair. Protein-DNA association is also important in medicine; understanding Protein-DNA binding kinetics can assist in identifying disease root causes which can contribute to drug development. In this perspective, this work focuses on the transcription process by the GATA Transcription Factor (TF). GATA TF binds to DNA promoter region represented by `G,A,T,A' nucleotides sequence, and initiates transcription of target genes. When proper regulation fails due to some mutations on the GATA TF protein sequence or on the DNA promoter sequence (weak promoter), deregulation of the target genes might lead to various disorders. In this study, we aim to understand the electrostatic mechanism behind GATA TF and DNA promoter interactions, in order to predict Protein-DNA binding in the presence of mutations, while elaborating on non-covalent binding kinetics. To generate a family of mutants for the GATA:DNA complex, we replaced every charged amino acid, one at a time, with a neutral amino acid like Alanine (Ala). We then applied Poisson-Boltzmann electrostatic calculations feeding into free energy calculations, for each mutation. These calculations delineate the contribution to binding from each Ala-replaced amino acid in the GATA:DNA interaction. After analyzing the obtained data in view of a two-step model, we are able to identify potential key amino acids in binding. Finally, we applied the model to GATA-3:DNA (crystal structure with PDB-ID: 3DFV) binding complex and validated it against experimental results from the literature.
Collapse
|
7
|
Abstract
![]()
Electrostatic effects
are ubiquitous in protein interactions and
are found to be pervasive in the complement system as well. The interaction
between complement fragment C3d and complement receptor 2 (CR2) has
evolved to become a link between innate and adaptive immunity. Electrostatic
interactions have been suggested to be the driving factor for the
association of the C3d:CR2 complex. In this study, we investigate
the effects of ionic strength and mutagenesis on the association of
C3d:CR2 through Brownian dynamics simulations. We demonstrate that
the formation of the C3d:CR2 complex is ionic strength-dependent,
suggesting the presence of long-range electrostatic steering that
accelerates the complex formation. Electrostatic steering occurs through
the interaction of an acidic surface patch in C3d and the positively
charged CR2 and is supported by the effects of mutations within the
acidic patch of C3d that slow or diminish association. Our data are
in agreement with previous experimental mutagenesis and binding studies
and computational studies. Although the C3d acidic patch may be locally
destabilizing because of unfavorable Coulombic interactions of like
charges, it contributes to the acceleration of association. Therefore,
acceleration of function through electrostatic steering takes precedence
to stability. The site of interaction between C3d and CR2 has been
the target for delivery of CR2-bound nanoparticle, antibody, and small
molecule biomarkers, as well as potential therapeutics. A detailed
knowledge of the physicochemical basis of C3d:CR2 association may
be necessary to accelerate biomarker and drug discovery efforts.
Collapse
Affiliation(s)
- Rohith R Mohan
- Department of Bioengineering, University of California , Riverside, California 92521, United States
| | - Gary A Huber
- Department of Chemistry and Biochemistry, University of California , San Diego, California 92093, United States
| | - Dimitrios Morikis
- Department of Bioengineering, University of California , Riverside, California 92521, United States
| |
Collapse
|
8
|
Electrostatic Interactions between Complement Regulator CD46(SCR1-2) and Adenovirus Ad11/Ad21 Fiber Protein Knob. Mol Biol Int 2015; 2015:967465. [PMID: 26357573 PMCID: PMC4556874 DOI: 10.1155/2015/967465] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2015] [Revised: 07/12/2015] [Accepted: 07/16/2015] [Indexed: 11/24/2022] Open
Abstract
Adenoviruses bind to a variety of human cells to cause infection. Both the B2 adenovirus 11 and B1 adenovirus 21 use protein knobs to bind to complement regulator CD46(SCR1-2) in order to gain entry into host cells. In each complex, the two proteins are highly negatively charged but bind to each other at an interface with oppositely charged surface patches. We computationally generated single-alanine mutants of charged residues in the complexes CD46(SCR1-2)-Ad11k and CD46(SCR1-2)-Ad21k. We used electrostatic clustering and Poisson-Boltzmann free energy calculations to propose a hypothesis on the role of electrostatics in association. Our results delineate specific interfacial electrostatic interactions that are critical for association in both CD46(SCR1-2)-Ad11k and CD46(SCR1-2)-Ad21k. These results will serve as a predictive tool in the selection of mutants with desired binding affinity in experimental mutagenesis studies. This study will also serve as a foundation for the design of inhibitors to treat adenovirus infections.
Collapse
|
9
|
Liu Y, Kieslich CA, Morikis D, Liao J. Engineering pre-SUMO4 as efficient substrate of SENP2. Protein Eng Des Sel 2015; 27:117-26. [PMID: 24671712 DOI: 10.1093/protein/gzu004] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023] Open
Abstract
SUMOylation, one of the most important protein post-translational modifications, plays critical roles in a variety of physiological and pathological processes. SENP (Sentrin/SUMO-specific protease), a family of SUMO-specific proteases, is responsible for the processing of pre-SUMO and removal of SUMO from conjugated substrates. SUMO4, the latest discovered member in the SUMO family, has been found as a type 1 diabetes susceptibility gene and its maturation is not understood so far. Despite the 14 amino acid differences between pre-SUMO4 and SUMO2, pre-SUMO4 is not processed by SENP2 but pre-SUMO2 does. A novel interdisciplinary approach involving computational modeling and a FRET-based protease assay was taken to engineer pre-SUMO4 as a substrate of SENP2. Given the difference in net charge between pre-SUMO4 and pre-SUMO2, the computational framework analysis of electrostatic similarities of proteins was applied to determine the contribution of each ionizable amino acid in a model of SENP2-(pre-SUMO4) binding, and to propose pre-SUMO4 mutations. The specificities of the SENP2 toward different pre-SUMO4 mutants were determined using a quantitative FRET assay by characterizing the catalytic efficiencies (kcat/KM). A single amino acid mutation made pre-SUMO4 amenable to SENP2 processing and a combination of two amino acid mutations made it highly accessible as SENP2 substrate. The combination of the two approaches provides a powerful protein engineering tool for future SUMOylation studies.
Collapse
Affiliation(s)
- Yan Liu
- Department of Bioengineering, University of California at Riverside, 900 University Avenue, Riverside, CA 92521, USA
| | | | | | | |
Collapse
|
10
|
E S Harrison R, Gorham RD, Morikis D. Energetic evaluation of binding modes in the C3d and Factor H (CCP 19-20) complex. Protein Sci 2015; 24:789-802. [PMID: 25628052 DOI: 10.1002/pro.2650] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2014] [Revised: 12/27/2014] [Accepted: 01/25/2015] [Indexed: 01/01/2023]
Abstract
As a part of innate immunity, the complement system relies on activation of the alternative pathway (AP). While feed-forward amplification generates an immune response towards foreign surfaces, the process requires regulation to prevent an immune response on the surface of host cells. Factor H (FH) is a complement protein secreted by native cells to negatively regulate the AP. In terms of structure, FH is composed of 20 complement-control protein (CCP) modules that are structurally homologous but vary in composition and function. Mutations in these CCPs have been linked to states of autoimmunity. In particular, several mutations in CCP 19-20 are correlated to atypical hemolytic uremic syndrome (aHUS). From crystallographic structures there are three putative binding sites of CCP 19-20 on C3d. Since there has been some controversy over the primary mode of binding from experimental studies, we approach characterization of binding using computational methods. Specifically, we compare each binding mode in terms of electrostatic character, structural stability, dissociative and associative properties, and predicted free energy of binding. After a detailed investigation, we found two of the three binding sites to be similarly stable while varying in the number of contacts to C3d and in the energetic barrier to complex dissociation. These sites are likely physiologically relevant and may facilitate multivalent binding of FH CCP 19-20 to C3b and either C3d or host glycosaminoglycans. We propose thermodynamically stable binding with modules 19 and 20, the latter driven by electrostatics, acting synergistically to increase the apparent affinity of FH for host surfaces.
Collapse
Affiliation(s)
- Reed E S Harrison
- Department of Bioengineering, Bourns College of Engineering, University of California, Riverside, California
| | | | | |
Collapse
|
11
|
A theoretical view of the C3d:CR2 binding controversy. Mol Immunol 2014; 64:112-22. [PMID: 25433434 DOI: 10.1016/j.molimm.2014.11.006] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2014] [Accepted: 11/06/2014] [Indexed: 11/23/2022]
Abstract
The C3d:CR2(SCR1-2) interaction plays an important role in bridging innate and adaptive immunity, leading to enhanced antibody production at sites of complement activation. Over the past decade, there has been much debate over the binding mode of this interaction. An initial cocrystal structure (PDB: 1GHQ) was published in 2001, in which the only interactions observed were between the SCR2 domain of CR2 and a side-face of C3d whereas a cocrystal structure (PDB: 3OED) published in 2011 showed both the SCR1 and SCR2 domains of CR2 interacting with an acidic patch on the concave surface of C3d. The initial 1GHQ structure is at odds with the majority of existing biochemical data and the publication of the 3OED structure renewed uncertainty regarding the physiological relevance of 1GHQ, suggesting that crystallization may have been influenced by the presence of zinc acetate in the crystallization process. In our study, we used a variety of computational approaches to gain insight into the binding mode between C3d and CR2 and demonstrate that the binding site at the acidic patch (3OED) is electrostatically more favorable, exhibits better structural and dissociative stability, specifically at the SCR1 domain, and has higher binding affinity than the 1GHQ binding mode. We also observe that nonphysiological zinc ions enhance the formation of the C3d:CR2 complex at the side face of C3d (1GHQ) through increases in electrostatic favorability, intermolecular interactions, dissociative character and overall energetic favorability. These results provide a theoretical basis for the association of C3d:CR2 at the acidic cavity of C3d and provide an explanation for binding of CR2 at the side face of C3d in the presence of nonphysiological zinc ions.
Collapse
|
12
|
Gorham RD, Rodriguez W, Morikis D. Molecular analysis of the interaction between staphylococcal virulence factor Sbi-IV and complement C3d. Biophys J 2014; 106:1164-73. [PMID: 24606940 DOI: 10.1016/j.bpj.2014.01.033] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2013] [Revised: 12/29/2013] [Accepted: 01/23/2014] [Indexed: 02/07/2023] Open
Abstract
Staphylococcus aureus expresses numerous virulence factors that aid in immune evasion. The four-domain staphylococcal immunoglobulin binding (Sbi) protein interacts with complement component 3 (C3) and its thioester domain (C3d)-containing fragments. Recent structural data suggested two possible modes of binding of Sbi domain IV (Sbi-IV) to C3d, but the physiological binding mode remains unclear. We used a computational approach to provide insight into the C3d-Sbi-IV interaction. Molecular dynamics (MD) simulations showed that the first binding mode (PDB code 2WY8) is more robust than the second (PDB code 2WY7), with more persistent polar and nonpolar interactions, as well as conserved interfacial solvent accessible surface area. Brownian dynamics and steered MD simulations revealed that the first binding mode has faster association kinetics and maintains more stable intermolecular interactions compared to the second binding mode. In light of available experimental and structural data, our data confirm that the first binding mode represents Sbi-IV interaction with C3d (and C3) in a physiological context. Although the second binding mode is inherently less stable, we suggest a possible physiological role. Both binding sites may serve as a template for structure-based design of novel complement therapeutics.
Collapse
Affiliation(s)
- Ronald D Gorham
- Department of Bioengineering, University of California, Riverside, California
| | - Wilson Rodriguez
- Department of Bioengineering, University of California, Riverside, California
| | - Dimitrios Morikis
- Department of Bioengineering, University of California, Riverside, California.
| |
Collapse
|
13
|
Ojha H, Panwar HS, Gorham RD, Morikis D, Sahu A. Viral regulators of complement activation: structure, function and evolution. Mol Immunol 2014; 61:89-99. [PMID: 24976595 DOI: 10.1016/j.molimm.2014.06.004] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2014] [Revised: 05/30/2014] [Accepted: 06/01/2014] [Indexed: 11/25/2022]
Abstract
The complement system surveillance in the host is effective in controlling viral propagation. Consequently, to subvert this effector mechanism, viruses have developed a series of adaptations. One among these is encoding mimics of host regulators of complement activation (RCA) which help viruses to avoid being labeled as 'foreign' and protect them from complement-mediated neutralization and complement-enhanced antiviral adaptive immunity. In this review, we provide an overview on the structure, function and evolution of viral RCA proteins.
Collapse
Affiliation(s)
- Hina Ojha
- National Centre for Cell Science, Pune University Campus, Ganeshkhind, Pune 411007, India
| | - Hemendra Singh Panwar
- National Centre for Cell Science, Pune University Campus, Ganeshkhind, Pune 411007, India
| | - Ronald D Gorham
- Department of Bioengineering, University of California, Riverside, CA 92521, USA
| | - Dimitrios Morikis
- Department of Bioengineering, University of California, Riverside, CA 92521, USA.
| | - Arvind Sahu
- National Centre for Cell Science, Pune University Campus, Ganeshkhind, Pune 411007, India.
| |
Collapse
|
14
|
Botello-Smith WM, Cai Q, Luo R. Biological applications of classical electrostatics methods. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2014. [DOI: 10.1142/s0219633614400082] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Continuum electrostatics modeling of solvation based on the Poisson–Boltzmann (PB) equation has gained wide acceptance in biomolecular applications such as energetic analysis and structural visualization. Successful application of the PB solvent models requires careful calibration of the solvation parameters. Extensive testing and validation is also important to ensure accuracy in their applications. Limitation in the continuum modeling of solvation is also a known issue in certain biomolecular applications. Growing interest in membrane systems has further spurred developmental efforts to allow inclusion of membrane in the PB solvent models. Despite their past successes due to careful parameterization, algorithm development and parallel implementation, there is still much to be done to improve their transferability from the small molecular systems upon which they were developed and validated to complex macromolecular systems as advances in technology continue to push forward, providing ever greater computational resources to researchers to study more interesting biological systems of higher complexity.
Collapse
Affiliation(s)
- Wesley M. Botello-Smith
- Chemical Physics and Material Physics Graduate Program, University of California, Irvine, CA 92697, USA
- Department of Chemistry, University of California, Irvine, CA 92697, USA
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA 92697, USA
| | - Qin Cai
- Department of Biomedical Engineering, University of California, Irvine, CA 92697, USA
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA 92697, USA
| | - Ray Luo
- Department of Biomedical Engineering, University of California, Irvine, CA 92697, USA
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA 92697, USA
| |
Collapse
|
15
|
Electrostatic similarities between protein and small molecule ligands facilitate the design of protein-protein interaction inhibitors. PLoS One 2013; 8:e75762. [PMID: 24130741 PMCID: PMC3794991 DOI: 10.1371/journal.pone.0075762] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2013] [Accepted: 08/19/2013] [Indexed: 12/26/2022] Open
Abstract
One of the underlying principles in drug discovery is that a biologically active compound is complimentary in shape and molecular recognition features to its receptor. This principle infers that molecules binding to the same receptor may share some common features. Here, we have investigated whether the electrostatic similarity can be used for the discovery of small molecule protein-protein interaction inhibitors (SMPPIIs). We have developed a method that can be used to evaluate the similarity of electrostatic potentials between small molecules and known protein ligands. This method was implemented in a software called EleKit. Analyses of all available (at the time of research) SMPPII structures indicate that SMPPIIs bear some similarities of electrostatic potential with the ligand proteins of the same receptor. This is especially true for the more polar SMPPIIs. Retrospective analysis of several successful SMPPIIs has shown the applicability of EleKit in the design of new SMPPIIs.
Collapse
|
16
|
López de Victoria A, Tamamis P, Kieslich CA, Morikis D. Insights into the structure, correlated motions, and electrostatic properties of two HIV-1 gp120 V3 loops. PLoS One 2012; 7:e49925. [PMID: 23185486 PMCID: PMC3501474 DOI: 10.1371/journal.pone.0049925] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2012] [Accepted: 10/15/2012] [Indexed: 12/11/2022] Open
Abstract
The V3 loop of the glycoprotein 120 (gp120) is a contact point for cell entry of HIV-1 leading to infection. Despite sequence variability and lack of specific structure, the highly flexible V3 loop possesses a well-defined role in recognizing and selecting cell-bound coreceptors CCR5 and CXCR4 through a mechanism of charge complementarity. We have performed two independent molecular dynamics (MD) simulations to gain insights into the dynamic character of two V3 loops with slightly different sequences, but significantly different starting crystallographic structures. We have identified highly populated trajectory-specific salt bridges between oppositely charged stem residues Arg9 and Glu25 or Asp29. The two trajectories share nearly identical correlated motions within the simulations, despite their different overall structures. High occupancy salt bridges play a key role in the major cross-correlated motions in both trajectories, and may be responsible for transient structural stability in preparation for coreceptor binding. In addition, the two V3 loops visit conformations with similarities in spatial distributions of electrostatic potentials, despite their inherent flexibility, which may play a role in coreceptor recognition. It is plausible that cooperativity between overall electrostatic potential, charged residue interactions, and correlated motions could be associated with a coreceptor selection and binding.
Collapse
Affiliation(s)
- Aliana López de Victoria
- Department of Bioengineering, University of California Riverside, Riverside, California, United States of America
| | - Phanourios Tamamis
- Department of Bioengineering, University of California Riverside, Riverside, California, United States of America
| | - Chris A. Kieslich
- Department of Bioengineering, University of California Riverside, Riverside, California, United States of America
| | - Dimitrios Morikis
- Department of Bioengineering, University of California Riverside, Riverside, California, United States of America
| |
Collapse
|
17
|
Bellows-Peterson ML, Fung HK, Floudas CA, Kieslich CA, Zhang L, Morikis D, Wareham KJ, Monk PN, Hawksworth OA, Woodruff T. De novo peptide design with C3a receptor agonist and antagonist activities: theoretical predictions and experimental validation. J Med Chem 2012; 55:4159-68. [PMID: 22500977 PMCID: PMC3349770 DOI: 10.1021/jm201609k] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Targeting the complement component 3a receptor (C3aR) with selective agonists or antagonists is believed to be a viable therapeutic option for several diseases such as stroke, heart attack, reperfusion injuries, and rheumatoid arthritis. We designed a number of agonists, partial agonists, and antagonists of C3aR using our two-stage de novo protein design framework. Of the peptides tested using a degranulation assay in C3aR-transfected rat basophilic leukemia cells, two were prominent agonists (EC(50) values of 25.3 and 66.2 nM) and two others were partial agonists (IC(50) values of 15.4 and 26.1 nM). Further testing of these lead compounds in a calcium flux assay in U937 cells yielded similar results although with reduced potencies compared to transfected cells. The partial agonists also displayed full antagonist activity when tested in a C3aR inhibition assay. In addition, the electrostatic potential profile was shown to potentially discriminate between full agonists and partial agonists.
Collapse
Affiliation(s)
| | - Ho Ki Fung
- Department of Chemical and Biological Engineering, Princeton University, Princeton, NJ
| | | | - Chris A. Kieslich
- Department of Bioengineering, University of California at Riverside, Riverside, CA
| | - Li Zhang
- Department of Bioengineering, University of California at Riverside, Riverside, CA
| | - Dimitrios Morikis
- Department of Bioengineering, University of California at Riverside, Riverside, CA
| | - Kathryn J. Wareham
- Department of Infection and Immunity, The University of Sheffield Medical School, Sheffield, UK
| | - Peter N. Monk
- Department of Infection and Immunity, The University of Sheffield Medical School, Sheffield, UK
| | - Owen A. Hawksworth
- School of Biomedical Sciences, The University of Queensland, Brisbane, Australia
| | - Trent Woodruff
- School of Biomedical Sciences, The University of Queensland, Brisbane, Australia
| |
Collapse
|
18
|
Stein M, Pilli M, Bernauer S, Habermann BH, Zerial M, Wade RC. The interaction properties of the human Rab GTPase family--comparative analysis reveals determinants of molecular binding selectivity. PLoS One 2012; 7:e34870. [PMID: 22523562 PMCID: PMC3327705 DOI: 10.1371/journal.pone.0034870] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2011] [Accepted: 03/06/2012] [Indexed: 01/07/2023] Open
Abstract
Background Rab GTPases constitute the largest subfamily of the Ras protein superfamily. Rab proteins regulate organelle biogenesis and transport, and display distinct binding preferences for effector and activator proteins, many of which have not been elucidated yet. The underlying molecular recognition motifs, binding partner preferences and selectivities are not well understood. Methodology/Principal Findings Comparative analysis of the amino acid sequences and the three-dimensional electrostatic and hydrophobic molecular interaction fields of 62 human Rab proteins revealed a wide range of binding properties with large differences between some Rab proteins. This analysis assists the functional annotation of Rab proteins 12, 14, 26, 37 and 41 and provided an explanation for the shared function of Rab3 and 27. Rab7a and 7b have very different electrostatic potentials, indicating that they may bind to different effector proteins and thus, exert different functions. The subfamily V Rab GTPases which are associated with endosome differ subtly in the interaction properties of their switch regions, and this may explain exchange factor specificity and exchange kinetics. Conclusions/Significance We have analysed conservation of sequence and of molecular interaction fields to cluster and annotate the human Rab proteins. The analysis of three dimensional molecular interaction fields provides detailed insight that is not available from a sequence-based approach alone. Based on our results, we predict novel functions for some Rab proteins and provide insights into their divergent functions and the determinants of their binding partner selectivity.
Collapse
Affiliation(s)
- Matthias Stein
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Heidelberg, Germany
- Max Planck Institute for Dynamics of Complex Technical Systems, Magdeburg, Germany
- * E-mail: (MS); (RW)
| | - Manohar Pilli
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Heidelberg, Germany
| | - Sabine Bernauer
- Max Planck Institute of Molecular Cell Biology and Genetics, Dresden, Germany
| | - Bianca H. Habermann
- Max Planck Institute of Molecular Cell Biology and Genetics, Dresden, Germany
- Max Planck Institute for Biology of Ageing, Cologne, Germany
| | - Marino Zerial
- Max Planck Institute of Molecular Cell Biology and Genetics, Dresden, Germany
| | - Rebecca C. Wade
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Heidelberg, Germany
- * E-mail: (MS); (RW)
| |
Collapse
|
19
|
López de Victoria A, Kieslich CA, Rizos AK, Krambovitis E, Morikis D. Clustering of HIV-1 Subtypes Based on gp120 V3 Loop electrostatic properties. BMC BIOPHYSICS 2012; 5:3. [PMID: 22313935 PMCID: PMC3295656 DOI: 10.1186/2046-1682-5-3] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/13/2011] [Accepted: 02/07/2012] [Indexed: 11/10/2022]
Abstract
BACKGROUND The V3 loop of the glycoprotein gp120 of HIV-1 plays an important role in viral entry into cells by utilizing as coreceptor CCR5 or CXCR4, and is implicated in the phenotypic tropisms of HIV viruses. It has been hypothesized that the interaction between the V3 loop and CCR5 or CXCR4 is mediated by electrostatics. We have performed hierarchical clustering analysis of the spatial distributions of electrostatic potentials and charges of V3 loop structures containing consensus sequences of HIV-1 subtypes. RESULTS Although the majority of consensus sequences have a net charge of +3, the spatial distribution of their electrostatic potentials and charges may be a discriminating factor for binding and infectivity. This is demonstrated by the formation of several small subclusters, within major clusters, which indicates common origin but distinct spatial details of electrostatic properties. Some of this information may be present, in a coarse manner, in clustering of sequences, but the spatial details are largely lost. We show the effect of ionic strength on clustering of electrostatic potentials, information that is not present in clustering of charges or sequences. We also make correlations between clustering of electrostatic potentials and net charge, coreceptor selectivity, global prevalence, and geographic distribution. Finally, we interpret coreceptor selectivity based on the N6X7T8|S8X9 sequence glycosylation motif, the specific positive charge location according to the 11/24/25 rule, and the overall charge and electrostatic potential distribution. CONCLUSIONS We propose that in addition to the sequence and the net charge of the V3 loop of each subtype, the spatial distributions of electrostatic potentials and charges may also be important factors for receptor recognition and binding and subsequent viral entry into cells. This implies that the overall electrostatic potential is responsible for long-range recognition of the V3 loop with coreceptors CCR5/CXCR4, whereas the charge distribution contributes to the specific short-range interactions responsible for the formation of the bound complex. We also propose a scheme for coreceptor selectivity based on the sequence glycosylation motif, the 11/24/25 rule, and net charge.
Collapse
Affiliation(s)
| | - Chris A Kieslich
- Department of Bioengineering, University of California, Riverside 92521, USA
| | - Apostolos K Rizos
- Department of Chemistry, University of Crete and Foundation for Research and Technology-Hellas, FORTH-IESL, GR-71003, Heraklion, Crete, Greece
| | - Elias Krambovitis
- Department of Veterinary Medicine, University of Thessaly, Karditsa, Greece
| | - Dimitrios Morikis
- Department of Bioengineering, University of California, Riverside 92521, USA
| |
Collapse
|
20
|
Poisson–Boltzmann Implicit Solvation Models. ACTA ACUST UNITED AC 2012. [DOI: 10.1016/b978-0-444-59440-2.00006-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
|
21
|
Complement Inhibition by Staphylococcus aureus: Electrostatics of C3d–EfbC and C3d–Ehp Association. Cell Mol Bioeng 2011. [DOI: 10.1007/s12195-011-0195-6] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
|
22
|
Electrostatic exploration of the C3d-FH4 interaction using a computational alanine scan. Mol Immunol 2011; 48:1844-50. [PMID: 21683447 DOI: 10.1016/j.molimm.2011.05.007] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2011] [Accepted: 05/09/2011] [Indexed: 11/20/2022]
Abstract
The complement system is a component of innate immunity and is activated by a cascade of protein interactions whose function is vital to our ability to fight infection. When proper regulation fails, the complement system is unable to recognize "self" from "nonself" and, therefore, attacks own tissues leading to autoimmune diseases. The central protein of the complement system is C3, which is the convergence point of three independently activated but communicating pathways. Regulation of C3 occurs through modular proteins which consist of many repeats of complement control protein (CCP) modules. CCP modules have diverse sequences, similar structures, and diverse physicochemical compositions, with excess of charge being a predominant characteristic. The goal of our study is to understand the electrostatic mechanism that underlies the interaction between the C3d domain of C3 and the fourth module of the complement regulator Factor H (FH4). We have performed a computational alanine scan in which we have replaced every ionizable amino acid, one at a time, with an alanine to generate a family of mutants for the C3d-FH4 complex. We have used Poisson-Boltzmann electrostatic calculations in combination with clustering of spatial distributions of electrostatic potentials and free energy calculations to delineate the contribution of each replaced amino acid to the C3d-FH4 interaction. We have analyzed our data in view of a two-step model which separates association into long-range recognition and short-range binding and we have identified key amino acids that contribute to association. We discuss the complex role of C3d in binding FH4 and the bacterial proteins Efb/Ehp from Staphylococcus aureus.
Collapse
|
23
|
The effect of electrostatics on factor H function and related pathologies. J Mol Graph Model 2011; 29:1047-55. [PMID: 21605993 DOI: 10.1016/j.jmgm.2011.04.010] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2011] [Revised: 04/26/2011] [Accepted: 04/28/2011] [Indexed: 12/15/2022]
Abstract
Factor H (FH) contributes to the regulation of the complement system by binding to polyanionic surfaces and the proteins C3b/C3c/C3d. This implicates charge and electrostatic interactions in recognition and binding of FH. Despite the large amount of experimental and pathology data the exact mechanism at molecular level is not yet known. We have implemented a computational framework for comparative analysis of the charge and electrostatic diversity of FH modules and C3b domains to identify electrostatic hotspots and predict potential binding sites. Our electrostatic potential clustering analysis shows that charge distributions and electrostatic potential distributions are more useful in understanding C3b-FH interactions than net charges alone. We present a model of non-specific electrostatic interactions of FH with polyanion-rich surfaces and specific interactions with C3b, using our computational data and existing experimental data. We discuss the electrostatic contributions to the formation of the C3b-FH complex and the competition between FH and Factor Bb (Bb) for binding to C3b. We also discuss the significance of mutations of charged amino acids in the pathobiology of FH-mediated disease, such as age-related macular degeneration, atypical hemolytic uremic syndrome, and dense deposit disease. Our data can be used to guide future experimental studies.
Collapse
|
24
|
Electrostatic clustering and free energy calculations provide a foundation for protein design and optimization. Ann Biomed Eng 2010; 39:1252-63. [PMID: 21140293 PMCID: PMC3069318 DOI: 10.1007/s10439-010-0226-9] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2010] [Accepted: 11/24/2010] [Indexed: 01/14/2023]
Abstract
Electrostatic interactions are ubiquitous in proteins and dictate stability and function. In this review, we discuss several methods for the analysis of electrostatics in protein–protein interactions. We discuss alanine-scanning mutagenesis, Poisson–Boltzmann electrostatics, free energy calculations, electrostatic similarity distances, and hierarchical clustering of electrostatic potentials. Our recently developed computational framework, known as Analysis of Electrostatic Similarities Of Proteins (AESOP), incorporates these tools to efficiently elucidate the role of electrostatic potentials in protein interactions. We present the application of AESOP to several proteins and protein complexes, for which charge is purported to facilitate protein association. Specifically, we illustrate how recent work has shaped the formulation of electrostatic calculations, the correlation of electrostatic free energies and electrostatic potential clustering results with experimental binding and activity data, the pH dependence of protein stability and association, the design of mutant proteins with enhanced immunological activity, and how AESOP can expose deficiencies in structural models and experimental data. This integrative approach can be utilized to develop mechanistic models and to guide experimental studies by predicting mutations with desired physicochemical properties and function. Alteration of the electrostatic properties of proteins offers a basis for the design of proteins with optimized binding and activity.
Collapse
|