Hot‐spot analysis to dissect the functional protein–protein interface of a tRNA‐modifying enzyme

Interference with protein–protein interactions of interfaces larger than 1500 Ų by small drug‐like molecules is notoriously difficult, particularly if targeting homodimers. The tRNA modifying enzyme Tgt is only functionally active as a homodimer. Thus, blocking Tgt dimerization is a promising strategy for drug therapy as this protein is key to the development of Shigellosis. Our goal was to identify hot‐spot residues which, upon mutation, result in a predominantly monomeric state of Tgt. The detailed understanding of the spatial location and stability contribution of the individual interaction hot‐spot residues and the plasticity of motifs involved in the interface formation is a crucial prerequisite for the rational identification of drug‐like inhibitors addressing the respective dimerization interface. Using computational analyses, we identified hot‐spot residues that contribute particularly to dimer stability: a cluster of hydrophobic and aromatic residues as well as several salt bridges. This in silico prediction led to the identification of a promising double mutant, which was validated experimentally. Native nano‐ESI mass spectrometry showed that the dimerization of the suggested mutant is largely prevented resulting in a predominantly monomeric state. Crystal structure analysis and enzyme kinetics of the mutant variant further support the evidence for enhanced monomerization and provide first insights into the structural consequences of the dimer destabilization. Proteins 2014; 82:2713–2732. © 2014 Wiley Periodicals, Inc.     

Monte Carlo simulations of the peptide recognition at the consensus binding site of the constant fragment of human immunoglobulin G: the energy landscape analysis of a hot spot at the intermolecular interface

Monte Carlo simulations of molecular recognition at the consensus binding site of the constant fragment (Fc) of human immunoglobulin G (Ig) protein have been performed to analyze structural and thermodynamic aspects of binding for the 13-residue cyclic peptide DCAWHLGELVWCT. The energy landscape analysis of a hot spot at the intermolecular interface using alanine scanning and equilibrium-simulated tempering dynamics with the simplified, knowledge-based energy function has enabled the role of the protein hot spot residues in providing the thermodynamic stability of the native structure to be determined. We have found that hydrophobic interactions between the peptide and the Met-252, Ile-253, His-433, and His-435 protein residues are critical to guarantee the thermodynamic stability of the crystallographic binding mode of the complex. Binding free energy calculations, using a molecular mechanics force field and a solvation energy model, combined with alanine scanning have been conducted to determine the energetic contribution of the protein hot spot residues in binding affinity. The conserved Asn-434, Ser-254, and Tyr-436 protein residues contribute significantly to the binding affinity of the peptide-protein complex, serving as an energetic hot spot at the intermolecular interface. The results suggest that evolutionary conserved hot spot protein residues at the intermolecular interface may be partitioned in fulfilling thermodynamic stability of the native binding mode and contributing to the binding affinity of the complex.     

Lethal mutations in the major homology region and their suppressors act by modulating the dimerization of the rous sarcoma virus capsid protein C‐terminal domain

An infective retrovirus requires a mature capsid shell around the viral replication complex. This shell is formed by about 1500 capsid protein monomers, organized into hexamer and pentamer rings that are linked to each other by the dimerization of the C‐terminal domain (CTD). The major homology region (MHR), the most highly conserved protein sequence across retroviral genomes, is part of the CTD. Several mutations in the MHR appear to block infectivity by preventing capsid formation. Suppressor mutations have been identified that are distant in sequence and structure from the MHR and restore capsid formation. The effects of two lethal and two suppressor mutations on the stability and function of the CTD were examined. No correlation with infectivity was found for the stability of the lethal mutations (D155Y‐CTD, F167Y‐CTD) and suppressor mutations (R185W‐CTD, I190V‐CTD). The stabilities of three double mutant proteins (D155Y/R185W‐CTD, F167Y/R185W‐CTD, and F167Y/I190V‐CTD) were additive. However, the dimerization affinity of the mutant proteins correlated strongly with biological function. The CTD proteins with lethal mutations did not dimerize, while those with suppressor mutations had greater dimerization affinity than WT‐CTD. The suppressor mutations were able to partially correct the dimerization defect caused by the lethal MHR mutations in double mutant proteins. Despite their dramatic effects on dimerization, none of these residues participate directly in the proposed dimerization interface in a mature capsid. These findings suggest that the conserved sequence of the MHR has critical roles in the conformation(s) of the CTD that are required for dimerization and correct capsid maturation. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Biophysical Analysis of the MHR Motif in Folding and Domain Swapping of the HIV Capsid Protein C-Terminal Domain

Infection by human immunodeficiency virus (HIV) depends on the function, in virion morphogenesis and other stages of the viral cycle, of a highly conserved structural element, the major homology region (MHR), within the carboxyterminal domain (CTD) of the capsid protein. In a modified CTD dimer, MHR is swapped between monomers. While no evidence for MHR swapping has been provided by structural models of retroviral capsids, it is unknown whether it may occur transiently along the virus assembly pathway. Whatever the case, the MHR-swapped dimer does provide a novel target for the development of anti-HIV drugs based on the concept of trapping a nonnative capsid protein conformation. We have carried out a thermodynamic and kinetic characterization of the domain-swapped CTD dimer in solution. The analysis includes a dissection of the role of conserved MHR residues and other amino acids at the dimerization interface in CTD folding, stability, and dimerization by domain swapping. The results revealed some energetic hotspots at the domain-swapped interface. In addition, many MHR residues that are not in the protein hydrophobic core were nevertheless found to be critical for folding and stability of the CTD monomer, which may dramatically slow down the swapping reaction. Conservation of MHR residues in retroviruses did not correlate with their contribution to domain swapping, but it did correlate with their importance for stable CTD folding. Because folding is required for capsid protein function, this remarkable MHR-mediated conformational stabilization of CTD may help to explain the functional roles of MHR not only during immature capsid assembly but in other processes associated with retrovirus infection. This energetic dissection of the dimerization interface in MHR-swapped CTD may also facilitate the design of anti-HIV compounds that inhibit capsid assembly by conformational trapping of swapped CTD dimers.     

Functional Roles of the Dimer-Interface Residues in Human Ornithine Decarboxylase

Ornithine decarboxylase (ODC) catalyzes the decarboxylation of ornithine to putrescine and is the rate-limiting enzyme in the polyamine biosynthesis pathway. ODC is a dimeric enzyme, and the active sites of this enzyme reside at the dimer interface. Once the enzyme dissociates, the enzyme activity is lost. In this paper, we investigated the roles of amino acid residues at the dimer interface regarding the dimerization, protein stability and/or enzyme activity of ODC. A multiple sequence alignment of ODC and its homologous protein antizyme inhibitor revealed that 5 of 9 residues (residues 165, 277, 331, 332 and 389) are divergent, whereas 4 (134, 169, 294 and 322) are conserved. Analytical ultracentrifugation analysis suggested that some dimer-interface amino acid residues contribute to formation of the dimer of ODC and that this dimerization results from the cooperativity of these interface residues. The quaternary structure of the sextuple mutant Y331S/Y389D/R277S/D332E/V322D/D134A was changed to a monomer rather than a dimer, and the K _d value of the mutant was 52.8 µM, which is over 500-fold greater than that of the wild-type ODC (ODC_WT). In addition, most interface mutants showed low but detectable or negligible enzyme activity. Therefore, the protein stability of these interface mutants was measured by differential scanning calorimetry. These results indicate that these dimer-interface residues are important for dimer formation and, as a consequence, are critical for enzyme catalysis.     

The C-Terminal Domain of ERp29 Mediates Polyomavirus Binding, Unfolding, and Infection

Penetration of the endoplasmic reticulum (ER) membrane by polyomavirus (PyV) is a decisive step in virus entry. We showed previously that the ER-resident factor ERp29 induces the local unfolding of PyV to initiate the ER membrane penetration process. ERp29 contains an N-terminal thioredoxin domain (NTD) that mediates its dimerization and a novel C-terminal all-helical domain (CTD) whose function is unclear. The NTD-mediated dimerization of ERp29 is critical for its unfolding activity; whether the CTD plays any role in PyV unfolding is unknown. We now show that three hydrophobic residues within the last helix of the ERp29 CTD that were individually mutated to either lysine or alanine abolished ERp29's ability to stimulate PyV unfolding and infection. This effect was not due to global misfolding of the mutant proteins, as they dimerize and do not form aggregates or display increased protease sensitivity. Moreover, the mutant proteins stimulated secretion of the secretory protein thyroglobulin with an efficiency similar to that of wild-type ERp29. Using a cross-linking coimmunoprecipitation assay, we found that the physical interaction of the ERp29 CTD mutants with PyV is inefficient. Our data thus demonstrate that the ERp29 CTD plays a crucial role in PyV unfolding and infection, likely by serving as part of a substrate-binding domain.     

Extraction of molecular features for the drug discovery targeting protein-protein interaction of Helicobacter pylori CagA and tumor suppressor protein ASSP2

Half of the world population is infected by the Gram-negative bacterium Helicobacter pylori (H. pylori). It colonizes in the stomach and is associated with severe gastric pathologies including gastric cancer and peptic ulceration. The most virulent factor of H. pylori is the cytotoxin-associated gene A (CagA) that is injected into the host cell. CagA interacts with several host proteins and alters their function, thereby causing several diseases. The most well-known target of CagA is the tumor suppressor protein ASPP2. The subdomain I at the N-terminus of CagA interacts with the proline-rich motif of ASPP2. Here, in this study, we carried out alanine scanning mutagenesis and an extensive molecular dynamics simulation summing up to 3.8 μs to find out hot spot residues and discovered some new protein-protein interaction (PPI)-modulating molecules. Our findings are in line with previous biochemical studies and further suggested new residues that are crucial for binding. The alanine scanning showed that mutation of Y207 and T211 residues to alanine decreased the binding affinity. Likewise, dynamics simulation and molecular mechanics with generalized Born surface area (MMGBSA) analysis also showed the importance of these two residues at the interface. A four-feature pharmacophore model was developed based on these two residues, and top 10 molecules were filtered from ZINC, NCI, and ChEMBL databases. The good binding affinity of the CHEMBL17319 and CHEMBL1183979 molecules shows the reliability of our adopted protocol for binding hot spot residues. We believe that our study provides a new insight for using CagA as the therapeutic target for gastric cancer treatment and provides a platform for a future experimental study.     

Kinetic computational alanine scanning: application to p53 oligomerization

We have developed a novel computational alanine scanning approach that involves analysis of ensemble unfolding kinetics at high temperature to identify residues that are critical for the stability of a given protein. This approach has been applied to dimerization of the oligomerization domain (residues 326-355) of tumor suppressor p53. As validated by experimental results, our approach has reasonable success in identifying deleterious mutations, including mutations that have been linked to cancer. We discuss a method for determining the effect of mutations on the location of the dimerization transition state.     

A binding free energy hot spot in the ankyrin repeat protein GABPbeta mediated protein-protein interaction

The frequently observed ankyrin repeat motif represents a structural scaffold evolved for mediating protein-protein interactions. As such, these repeats modulate a diverse range of cellular functions. We thermodynamically characterized the heterodimeric GA-binding protein (GABP) alphabeta complex and focused specifically on the interaction mediated by the ankyrin repeat domain of the GABPbeta. Our isothermal titration calorimetric analysis of the interaction between the GABP subunits determined an association constant (K(A)) of 6.0 x 10(8) M(-1) and that the association is favorably driven by a significant change in enthalpy (DeltaH) and a minor change in entropy (-TDeltaS). A total of 16 GABPbeta interface residues were chosen for alanine scanning mutagenesis. The calorimetrically measured differences in the free energy of binding were compared to computationally calculated values resulting in a correlation coefficient r = 0.71. We identified three spatially contiguous hydrophobic and aromatic residues that form a binding free energy hot spot (DeltaDeltaG > 2.0 kcal/mol). One residue provides structural support to the hot spot residues. Three non-hot spot residues are intermediate contributors (DeltaDeltaG approximately 1.0 kcal/mol) and create a canopy-like structure over the hot spot residues to possibly occlude solvent and orientate the subunits. The remaining interface residues are located peripherally and have weak contributions. Finally, our mutational analysis revealed a significant entropy-enthalpy compensation for this interaction.     

APIS: accurate prediction of hot spots in protein interfaces by combining protrusion index with solvent accessibility

Background  

It is well known that most of the binding free energy of protein interaction is contributed by a few key hot spot residues. These residues are crucial for understanding the function of proteins and studying their interactions. Experimental hot spots detection methods such as alanine scanning mutagenesis are not applicable on a large scale since they are time consuming and expensive. Therefore, reliable and efficient computational methods for identifying hot spots are greatly desired and urgently required.     

Computational detection of the binding-site hot spot at the remodeled human growth hormone-receptor interface

A hierarchical computational approach is used to identify the engineered binding-site cavity at the remodeled intermolecular interface between the mutants of human growth hormone (hGH) and the extracellular domain of its receptor (hGHbp). Multiple docking simulations are conducted with the remodeled hGH-hGHbp complex for a panel of potent benzimidazole-containing inhibitors that can restore the binding affinity of the wild-type complex, and for a set of known nonactive small molecules that contain different heterocyclic motifs. Structural clustering of ligand-bound conformations and binding free-energy calculations, using the AMBER force field and a continuum solvation model, can rapidly locate and screen numerous ligand-binding modes on the protein surface and detect the binding-site hot spot at the intermolecular interface. Structural orientation of the benzimidazole motif in the binding-site cavity closely mimics the position of the hot spot residue W104 in the crystal structure of the wild-type complex, which is recognized as an important structural requirement for restoring binding affinity. Despite numerous pockets on the protein surface of the mutant hGH-hGHbp complex, the binding-site cavity presents the energetically favorable hot spot for the benzimidazole-containing inhibitors, whereas for a set of nonactive molecules, the lowest energy ligand conformations do not necessarily bind in the engineered cavity. The results reveal a dominant role of the intermolecular van der Waals interactions in providing favorable ligand-protein energetics in the redesigned interface, in agreement with the experimental and computational alanine scanning of the hGH-hGHbp complex.     

The Extended Transmembrane Orai1 N-terminal (ETON) Region Combines Binding Interface and Gate for Orai1 Activation by STIM1

STIM1 and Orai1 represent the two molecular key components of the Ca²⁺ release-activated Ca²⁺ channels. Their activation involves STIM1 C terminus coupling to both the N terminus and the C terminus of Orai. Here we focused on the extended transmembrane Orai1 N-terminal (ETON, aa73–90) region, conserved among the Orai family forming an elongated helix of TM1 as recently shown by x-ray crystallography. To identify “hot spot” residues in the ETON binding interface for STIM1 interaction, numerous Orai1 constructs with N-terminal truncations or point mutations within the ETON region were generated. N-terminal truncations of the first four residues of the ETON region or beyond completely abolished STIM1-dependent Orai1 function. Loss of Orai1 function resulted from neither an impairment of plasma membrane targeting nor pore damage, but from a disruption of STIM1 interaction. In a complementary approach, we monitored STIM1-Orai interaction via Orai1 V102A by determining restored Ca²⁺ selectivity as a consequence of STIM1 coupling. Orai1 N-terminal truncations that led to a loss of function consistently failed to restore Ca²⁺ selectivity of Orai1 V102A in the presence of STIM1, demonstrating impairment of STIM1 binding. Hence, the major portion of the ETON region (aa76–90) is essential for STIM1 binding and Orai1 activation. Mutagenesis within the ETON region revealed several hydrophobic and basic hot spot residues that appear to control STIM1 coupling to Orai1 in a concerted manner. Moreover, we identified two basic residues, which protrude into the elongated pore to redound to Orai1 gating. We suggest that several hot spot residues in the ETON region contribute in aggregate to the binding of STIM1, which in turn is coupled to a conformational reorientation of the gate.     

Connecting Mutations of the RNA Polymerase II C-Terminal Domain to Complex Phenotypic Changes Using Combined Gene Expression and Network Analyses

The C-terminal domain (CTD) of the largest subunit in DNA-dependent RNA polymerase II (RNAP II) is essential for mRNA synthesis and processing, through coordination of an astounding array of protein-protein interactions. Not surprisingly, CTD mutations can have complex, pleiotropic impacts on phenotype. For example, insertions of five alanine residues between CTD diheptads in yeast, which alter the CTD''s overall tandem structure and physically separate core functional units, dramatically reduce growth rate and result in abnormally large cells that accumulate increased DNA content over time. Patterns by which specific CTD-protein interactions are disrupted by changes in CTD structure, as well as how downstream metabolic pathways are impacted, are difficult to target for direct experimental analyses. In an effort to connect an altered CTD to complex but quantifiable phenotypic changes, we applied network analyses of genes that are differentially expressed in our five alanine CTD mutant, combined with established genetic interactions from the Saccharomyces cerevisiae Genome Database (SGD). We were able to identify candidate genetic pathways, and several key genes, that could explain how this change in CTD structure leads to the specific phenotypic changes observed. These hypothetical networks identify links between CTD-associated proteins and mitotic function, control of cell cycle checkpoint mechanisms, and expression of cell wall and membrane components. Such results can help to direct future genetic and biochemical investigations that tie together the complex impacts of the CTD on global cellular metabolism.     

Structural and functional analysis of the MutS C-terminal tetramerization domain

The Escherichia coli DNA mismatch repair (MMR) protein MutS is essential for the correction of DNA replication errors. In vitro, MutS exists in a dimer/tetramer equilibrium that is converted into a monomer/dimer equilibrium upon deletion of the C-terminal 53 amino acids. In vivo and in vitro data have shown that this C-terminal domain (CTD, residues 801–853) is critical for tetramerization and the function of MutS in MMR and anti-recombination. We report the expression, purification and analysis of the E.coli MutS-CTD. Secondary structure prediction and circular dichroism suggest that the CTD is folded, with an α-helical content of 30%. Based on sedimentation equilibrium and velocity analyses, MutS-CTD forms a tetramer of asymmetric shape. A single point mutation (D835R) abolishes tetramerization but not dimerization of both MutS-CTD and full-length MutS. Interestingly, the in vivo and in vitro MMR activity of MutS^CF/D835R is diminished to a similar extent as a truncated MutS variant (MutS800, residues 1–800), which lacks the CTD. Moreover, the dimer-forming MutS^CF/D835R has comparable DNA binding affinity with the tetramer-forming MutS, but is impaired in mismatch-dependent activation of MutH. Our data support the hypothesis that tetramerization of MutS is important but not essential for MutS function in MMR.     

Key dimer interface residues impact the catalytic activity of 3CLpro,the main protease of SARS-CoV-2

3C-like protease (3CLpro) processes and liberates functional viral proteins essential for the maturation and infectivity of severe acute respiratory syndrome coronavirus 2, the virus responsible for COVID-19. It has been suggested that 3CLpro is catalytically active as a dimer, making the dimerization interface a target for antiviral development. Guided by structural analysis, here we introduced single amino acid substitutions at nine residues at three key sites of the dimer interface to assess their impact on dimerization and activity. We show that at site 1, alanine substitution of S1 or E166 increased by twofold or reduced relative activity, respectively. At site 2, alanine substitution of S10 or E14 eliminated activity, whereas K12A exhibited ∼60% relative activity. At site 3, alanine substitution of R4, E290, or Q299 eliminated activity, whereas S139A exhibited 46% relative activity. We further found that the oligomerization states of the dimer interface mutants varied; the inactive mutants R4A, R4Q, S10A/C, E14A/D/Q/S, E290A, and Q299A/E were present as dimers, demonstrating that dimerization is not an indication of catalytically active 3CLpro. In addition, present mostly as monomers, K12A displayed residual activity, which could be attributed to the conspicuous amount of dimer present. Finally, differential scanning calorimetry did not reveal a direct relationship between the thermodynamic stability of mutants with oligomerization or catalytic activity. These results provide insights on two allosteric sites, R4/E290 and S10/E14, that may promote the design of antiviral compounds that target the dimer interface rather than the active site of severe acute respiratory syndrome coronavirus 2 3CLpro.     

Rational disruption of the oligomerization of the mini-ferritin E. coli DPS through protein-protein interface mutation

DNA-binding protein from starved cells (DPS), a mini-ferritin capable of self-assembling into a 12-meric nano-cage, was chosen as the basis for an alanine-shaving mutagenesis study to investigate the importance of key amino acid residues, located at symmetry-related protein-protein interfaces, in controlling protein stability and self-assembly. Nine mutants were designed through simple inspection, synthesized, and subjected to transmission electron microscopy, circular dichroism, size exclusion chromatography, and "virtual alanine scanning" computational analysis. The data indicate that many of these residues may be hot spot residues. Most remarkably, two residues, R83 and R133, were observed to shift the oligomerization state to ~50% dimer. Based on the hypothesis that these two residues constitute a "hot strip," located at the ferritin-like threefold axis, the double mutant was generated which completely shuts down detectable formation of 12-mer in solution, favoring a cooperatively folded dimer. The fact that this effect logically builds upon the single mutants emphasizes that complex self-assembly has the potential to be manipulated rationally. This study should have an impact on the fundamental understanding of the assembly of DPS protein cages specifically and protein quaternary structure in general. In addition, as there is much interest in applying these and similar systems to the templation of nano-materials and drug delivery, the ability to control this ferritin's oligomerization state and stability could prove especially valuable.