Action-at-a-distance interactions enhance protein binding affinity

Energy landscape reshaped by strain-specific mutations underlies epistasis in NS1 evolution of influenza A virus

Elucidating how individual mutations affect the protein energy landscape is crucial for understanding how proteins evolve. However, predicting mutational effects remains challenging because of epistasis—the nonadditive interactions between mutations. Here, we investigate the biophysical mechanism of strain-specific epistasis in the nonstructural protein 1 (NS1) of influenza A viruses (IAVs). We integrate structural, kinetic, thermodynamic, and conformational dynamics analyses of four NS1s of influenza strains that emerged between 1918 and 2004. Although functionally near-neutral, strain-specific NS1 mutations exhibit long-range epistatic interactions with residues at the p85β-binding interface. We reveal that strain-specific mutations reshaped the NS1 energy landscape during evolution. Using NMR spin dynamics, we find that the strain-specific mutations altered the conformational dynamics of the hidden network of tightly packed residues, underlying the evolution of long-range epistasis. This work shows how near-neutral mutations silently alter the biophysical energy landscapes, resulting in diverse background effects during molecular evolution.

Influenza A virus (IAV) nonstructural protein 1 (NS1) is a multifunctional virulence factor that interacts with several host factors such as phosphatidylinositol-3-kinase (PI3K). NS1 binds specifically to the p85β regulatory subunit of PI3K and subsequently activates PI3K signaling. Here, Kim et al. show that functionally near-neutral, strain-specific NS1 mutations lead to variations in binding kinetics to p85β exhibit long-range epistatic interactions. Applying NMR they provide evidence that the structural dynamics of the NS1 hydrophobic core have evolved over time and contributed to epistasis.

Structural basis of human ACE2 higher binding affinity to currently circulating Omicron SARS-CoV-2 sub-variants BA.2 and BA.1.1

The currently circulating Omicron sub-variants are the SARS-CoV-2 strains with the highest number of known mutations. Herein, we found that human angiotensin-converting enzyme 2 (hACE2) binding affinity to the receptor-binding domains (RBDs) of the four early Omicron sub-variants (BA.1, BA.1.1, BA.2, and BA.3) follows the order BA.1.1 > BA.2 > BA.3 ≈ BA.1. The complex structures of hACE2 with RBDs of BA.1.1, BA.2, and BA.3 reveal that the higher hACE2 binding affinity of BA.2 than BA.1 is related to the absence of the G496S mutation in BA.2. The R346K mutation in BA.1.1 majorly affects the interaction network in the BA.1.1 RBD/hACE2 interface through long-range alterations and contributes to the higher hACE2 affinity of the BA.1.1 RBD than the BA.1 RBD. These results reveal the structural basis for the distinct hACE2 binding patterns among BA.1.1, BA.2, and BA.3 RBDs.

Insight into Structure-Function Relationships of β-Lactamase and BLIPs Interface Plasticity using Protein-Protein Interactions

BACKGROUND

Mostly BLIPs are identified in soil bacteria Streptomyces and originally isolated from Streptomyces clavuligerus and can be utilized as a model system for biophysical, structural, mutagenic and computational studies. BLIP possess homology with two proteins viz., BLIP-I (Streptomyces exofoliatus) and BLP (beta-lactamase inhibitory protein like protein from S. clavuligerus). BLIP consists of 165 amino acid, possessing two homologues domains comprising helix-loop-helix motif packed against four stranded beta-sheet resulting into solvent exposed concave surface with extended four stranded beta-sheet. BLIP-I is a 157 amino acid long protein obtained from S. exofoliatus having 37% sequence identity to BLIP and inhibits beta-lactamase.

METHODS

This review is intended to briefly illustrate the beta-lactamase inhibitory activity of BLIP via proteinprotein interaction and aims to open up a new avenue to combat antimicrobial resistance using peptide based inhibition.

RESULTS

D49A mutation in BLIP-I results in a decrease in affinity for TEM-1 from 0.5 nM to 10 nM (Ki). It is capable of inhibiting TEM-1 and bactopenemase and differs from BLIP only in modulating cell wall synthesis enzyme. Whereas, BLP is a 154 amino acid long protein isolated from S. clavuligerus via DNA sequencing analysis of Cephamycin-Clavulanate gene bunch. It shares 32% sequence similarity with BLIP and 42% with BLIP-I. Its biological function is unclear and lacks beta-lactamase inhibitory activity.

CONCLUSION

Protein-protein interactions mediate a significant role in regulation and modulation of cellular developments and processes. Specific biological markers and geometric characteristics are manifested by active site binding clefts of protein surfaces which determines the specificity and affinity for their targets. TEM1.BLIP is a classical model to study protein-protein interaction. β-Lactamase inhibitory proteins (BLIPs) interacts and inhibits various β-lactamases with extensive range of affinities.

Gangliosides interact with synaptotagmin to form the high-affinity receptor complex for botulinum neurotoxin B

Botulinum neurotoxin type B (BoNT/B) recognizes nerve terminals by binding to 2 receptor components: a polysialoganglioside, predominantly GT1b, and synaptotagmin 1/2. It is widely thought that BoNT/B initially binds to GT1b then diffuses in the plane of the membrane to interact with synaptotagmin. We have addressed the hypothesis that a GT1b-synaptotagmin cis complex forms the BoNT/B receptor. We identified a consensus glycosphingolipid-binding motif in the extracellular juxtamembrane domain of synaptotagmins 1/2 and confirmed by Langmuir monolayer, surface plasmon resonance, and circular dichroism that GT1b interacts with synaptotagmin peptides containing this sequence, inducing α-helical structure. Molecular modeling and tryptophan fluorescence spectroscopy were consistent with the intertwining of GT1b and synaptotagmin, involving cis interactions between the oligosaccharide and ceramide moieties of GT1b and the juxtamembrane and transmembrane domains of synaptotagmin, respectively. Furthermore, a point mutation on synaptotagmin, located outside of the BoNT/B-binding segment, inhibited GT1b binding and blocked GT1b-induced potentiation of BoNT/B binding to synaptotagmin-expressing cells. Our findings are consistent with a model in which a preassembled GT1b-synaptotagmin complex constitutes the high-affinity BoNT/B receptor.

Tackling the Antibiotic Resistance Caused by Class A β-Lactamases through the Use of β-Lactamase Inhibitory Protein

β-Lactams are the most widely used and effective antibiotics for the treatment of infectious diseases. Unfortunately, bacteria have developed several mechanisms to combat these therapeutic agents. One of the major resistance mechanisms involves the production of β-lactamase that hydrolyzes the β-lactam ring thereby inactivating the drug. To overcome this threat, the small molecule β-lactamase inhibitors (e.g., clavulanic acid, sulbactam and tazobactam) have been used in combination with β-lactams for treatment. However, the bacterial resistance to this kind of combination therapy has evolved recently. Therefore, multiple attempts have been made to discover and develop novel broad-spectrum β-lactamase inhibitors that sufficiently work against β-lactamase producing bacteria. β-lactamase inhibitory proteins (BLIPs) (e.g., BLIP, BLIP-I and BLIP-II) are potential inhibitors that have been found from soil bacterium Streptomyces spp. BLIPs bind and inhibit a wide range of class A β-lactamases from a diverse set of Gram-positive and Gram-negative bacteria, including TEM-1, PC1, SME-1, SHV-1 and KPC-2. To the best of our knowledge, this article represents the first systematic review on β-lactamase inhibitors with a particular focus on BLIPs and their inherent properties that favorably position them as a source of biologically-inspired drugs to combat antimicrobial resistance. Furthermore, an extensive compilation of binding data from β-lactamase–BLIP interaction studies is presented herein. Such information help to provide key insights into the origin of interaction that may be useful for rationally guiding future drug design efforts.

Structure-guided design of a potent peptide inhibitor targeting the interaction between CRK and ABL kinase

CT-10 regulator of kinase (CRK) proteins play important roles in human cancer metastasis and invasion. Moreover, CRK proteins are the major phosphorylation substrates of ABL kinase and its oncogenic mutant BCR-ABL kinase. The interaction between CRK and BCR-ABL plays important roles in chronic myeloid leukemia. Hence, inhibiting the interaction of CRK with BCR-ABL is an attractive way to attenuate cancer metastasis. Herein, we report the development of a peptide inhibitor, PRM-3, targeting the interaction between CRK-II and ABL kinase. PRM-3 binds to the N-terminal SH3 (nSH3) domain in CRK-II with a 10 nM affinity and prevents the interaction between CRK-II and ABL kinase. An in vitro biochemical assay demonstrated that PRM-3 inhibits the ABL-dependent phosphorylation of CRK-II more effectively than imatinib. Remarkably, PRM-3 also inhibited the CRK phosphorylation by T315I-ABL kinase, which is resistant to all first- and second-generation tyrosine kinase inhibitors. Our study provides a promising alternative approach to overcome the drug resistance of ABL kinase.

Using thermal scanning assays to test protein-protein interactions of inner-ear cadherins

Protein-protein interactions play a crucial role in biological processes such as cell-cell adhesion, immune system-pathogen interactions, and sensory perception. Understanding the structural determinants of protein-protein complex formation and obtaining quantitative estimates of their dissociation constant (KD) are essential for the study of these interactions and for the discovery of new therapeutics. At the same time, it is equally important to characterize protein-protein interactions in a high-throughput fashion. Here, we use a modified thermal scanning assay to test interactions of wild type (WT) and mutant variants of N-terminal fragments (EC1+2) of cadherin-23 and protocadherin-15, two proteins essential for inner-ear mechanotransduction. An environmentally sensitive fluorescent dye (SYPRO orange) is used to monitor melting temperature (Tm) shifts of protocadherin-15 EC1+2 (pcdh15) in the presence of increasing concentrations of cadherin-23 EC1+2 (cdh23). These Tm shifts are absent when we use proteins containing deafness-related missense mutations known to disrupt cdh23 binding to pcdh15, and are increased for some rationally designed mutants expected to enhance binding. In addition, surface plasmon resonance binding experiments were used to test if the Tm shifts correlated with changes in binding affinity. We used this approach to find a double mutation (cdh23(T15E)- pcdh15(G16D)) that enhances binding affinity of the cadherin complex by 1.98 kJ/mol, roughly two-fold that of the WT complex. We suggest that the thermal scanning methodology can be used in high-throughput format to quickly compare binding affinities (KD from nM up to 100 μM) for some heterodimeric protein complexes and to screen small molecule libraries to find protein-protein interaction inhibitors and enhancers.

The Molecular Mechanisms Underlying the Hijack of Host Proteins by the 1918 Spanish Influenza Virus

The 1918 Spanish influenza A virus (IAV) caused one of the most serious pandemics in history. The nonstructural protein 1 (NS1) of the 1918 IAV hijacks the interaction between human CrkII and JNK1. Little is, however, known about its molecular mechanism. Here, we performed X-ray crystallography, NMR relaxation dispersion experiment, and fluorescence spectroscopy to determine the structural, kinetic, and thermodynamic mechanisms underlying the hijacking of CrkII by 1918 IAV NS1. We observed that the interaction between a proline-rich motif in NS1 and the N-terminal SH3 domain of CrkII displays strikingly rapid kinetics and exceptionally high affinity with 100-fold faster k_on and 3300-fold lower K_d compared to those for the CrkII-JNK1 interaction. These results provide molecular insight into the mechanism by which 1918 IAV NS1 hijacks CrkII and disrupts its interactions with critical cellular signaling proteins.

Structural insight into the calcium ion modulated interdomain electron transfer in cellobiose dehydrogenase

Cellobiose dehydrogenase (CDH) from wood degrading fungi represents a subclass of oxidoreductases with unique properties. Consisting of two domains exhibiting interdomain electron transfer, this is the only known flavocytochrome involved in wood degradation. High resolution structures of the separated domains were solved, but the overall architecture of the intact protein and the exact interface of the two domains is unknown. Recently, it was shown that divalent cations modulate the activity of CDH and its pH optimum and a possible mechanism involving bridging of negative charges by calcium ions was proposed. Here we provide a structural explanation of this phenomenon confirming the interaction between negatively charged surface patches and calcium ions at the domain interface.

Rigid-body motions of interacting proteins dominate multispecific binding of ubiquitin in a shape-dependent manner

To understand the dynamic aspects of multispecificity of ubiquitin, we studied nine ubiquitin-ligand (partner protein) complexes by normal mode analysis based on an elastic network model. The coupling between ubiquitin and ligand motions was analyzed by decomposing it into rigid-body (external) and vibrational (internal) motions of each subunit. We observed that in total the external motions in one of the subunits largely dominated the coupling. The combination of external motions of ubiquitin and the ligands showed general trends of rotations and translations. Moreover, we observed that the rotational motions of ubiquitin were correlated to the ligand orientations. We also identified ubiquitin atomic vibrations that differentiated the orientation of the ligand molecule. We observed that the extents of coupling were correlated to the shapes of the ligands, and this trend was more pronounced when the coupling involved vibrational motions of the ligand. In conclusion, an intricate interplay between internal and external motions of ubiquitin and the ligands help understand the dynamics of multispecificity, which is mostly guided by the shapes of the ligands and the complex.

Optimization of affinity, specificity and function of designed influenza inhibitors using deep sequencing

We show that comprehensive sequence-function maps obtained by deep sequencing can be used to reprogram interaction specificity and to leapfrog over bottlenecks in affinity maturation by combining many individually small contributions not detectable in conventional approaches. We use this approach to optimize two computationally designed inhibitors against H1N1 influenza hemagglutinin and, in both cases, obtain variants with subnanomolar binding affinity. The most potent of these, a 51-residue protein, is broadly cross-reactive against all influenza group 1 hemagglutinins, including human H2, and neutralizes H1N1 viruses with a potency that rivals that of several human monoclonal antibodies, demonstrating that computational design followed by comprehensive energy landscape mapping can generate proteins with potential therapeutic utility.

A survey of aspartate-phenylalanine and glutamate-phenylalanine interactions in the protein data bank: searching for anion-π pairs

Protein structures are stabilized using noncovalent interactions. In addition to the traditional noncovalent interactions, newer types of interactions are thought to be present in proteins. One such interaction, an anion-π pair, in which the positively charged edge of an aromatic ring interacts with an anion, forming a favorable anion-quadrupole interaction, has been previously proposed [Jackson, M. R., et al. (2007) J. Phys. Chem. B111, 8242-8249]. To study the role of anion-π interactions in stabilizing protein structure, we analyzed pairwise interactions between phenylalanine (Phe) and the anionic amino acids, aspartate (Asp) and glutamate (Glu). Particular emphasis was focused on identification of Phe-Asp or -Glu pairs separated by less than 7 Å in the high-resolution, nonredundant Protein Data Bank. Simplifying Phe to benzene and Asp or Glu to formate molecules facilitated in silico analysis of the pairs. Kitaura-Morokuma energy calculations were performed on roughly 19000 benzene-formate pairs and the resulting energies analyzed as a function of distance and angle. Edgewise interactions typically produced strongly stabilizing interaction energies (-2 to -7.3 kcal/mol), while interactions involving the ring face resulted in weakly stabilizing to repulsive interaction energies. The strongest, most stabilizing interactions were identified as preferentially occurring in buried residues. Anion-π pairs are found throughout protein structures, in helices as well as β strands. Numerous pairs also had nearby cation-π interactions as well as potential π-π stacking. While more than 1000 structures did not contain an anion-π pair, the 3134 remaining structures contained approximately 2.6 anion-π pairs per protein, suggesting it is a reasonably common motif that could contribute to the overall structural stability of a protein.

Understanding the molecular activity of alkaline sphingomyelinase (NPP7) by computer modeling

The enzymes in the nucleotide pyrophosphatase/phosphodiesterase (NPP) family have various substrates such as nucleotides, phospholipids, and sphingolipids. The substrate specificity in relation to their structures is largely unknown because no mammalian NPP complex has been crystallized. NPP7, also called alkaline sphingomyelinase (alk-SMase), is a NPP family member that may have important implications in carcinogenesis and cholesterol absorption. The sequence of NPP7 is 36% similar to that of the closest NPP member, but NPP7 has no activity against nucleotides. In this work, we predict the three-dimensional structure of NPP7 by homology modeling using a recently crystallized NPP from bacteria. Using the model, we studied the substrate specificity of the enzyme by docking. The model generated explains the functional changes in previous mutagenesis studies and rationalizes the structural basis for the lack of activity toward nucleotides. An effort to shift the substrate specificity from sphingomyelin (SM) to nucleotide was not successful but revealed a site-directed mutation that increased activity toward SM. In conclusion, this is the first study to predict the structure of a mammalian NPP and its substrate specificity by molecular modeling. The information may be helpful in understanding the functional differences of NPP members.

Energetic decomposition with the generalized-born and Poisson-Boltzmann solvent models: lessons from association of G-protein components

Continuum electrostatic models have been shown to be powerful tools in providing insight into the energetics of biomolecular processes. While the Poisson-Boltzmann (PB) equation provides a theoretically rigorous approach to computing electrostatic free energies of solution in such a model, computational cost makes its use for large ensembles of states impractical. The generalized-Born (GB) approximation provides a much faster alternative, although with a weaker theoretical framework. While much attention has been given to how GB recapitulates PB energetics for the overall stability of a biomolecule or the affinity of a complex, little attention has been given to how the contributions of individual functional groups are captured by the two methods. Accurately capturing these individual electrostatic components is essential both for the development of a mechanistic understanding of biomolecular processes and for the design of variant sequences and structures with desired properties. Here, we present a detailed comparison of the group-wise decomposition of both PB and GB electrostatic free energies of binding, using association of various components of the heterotrimeric-G-protein complex as a model. We find that, while net binding free energies are strongly correlated in the two models, the correlations of individual group contributions are highly variable; in some cases, strong correlation is seen, while in others, there is essentially none. Structurally, the GB model seems to capture the magnitude of direct, short-range electrostatic interactions quite well but performs more poorly with moderate-range "action-at-a-distance" interactions--GB has a tendency to overestimate solvent screening over moderate distances, and to underestimate the costs of desolvating charged groups somewhat removed from the binding interface. Despite this, however, GB does seem to be quite effective as a predictor of those groups that will be computed to be most significant in a PB-based model.

Hub promiscuity in protein-protein interaction networks

Hubs are proteins with a large number of interactions in a protein-protein interaction network. They are the principal agents in the interaction network and affect its function and stability. Their specific recognition of many different protein partners is of great interest from the structural viewpoint. Over the last few years, the structural properties of hubs have been extensively studied. We review the currently known features that are particular to hubs, possibly affecting their binding ability. Specifically, we look at the levels of intrinsic disorder, surface charge and domain distribution in hubs, as compared to non-hubs, along with differences in their functional domains.

Fine mapping of the sequence requirements for binding of beta-lactamase inhibitory protein (BLIP) to TEM-1 beta-lactamase using a genetic screen for BLIP function

Beta-lactamase inhibitory protein (BLIP) binds and inhibits a diverse collection of class A beta-lactamases with a wide range of affinities. Alanine-scanning mutagenesis was previously performed to identify the amino acid sequence requirements of BLIP for binding the TEM-1, SME-1, SHV-1, and Bla1 beta-lactamases. Twenty-three BLIP residues that contact TEM-1 beta-lactamase in the structure of the complex were mutated to alanine and assayed for inhibition (K(i)) of beta-lactamase to identify two hotspots of binding energy. These studies have been extended by the development of a genetic screen for BLIP function in Escherichia coli. The bla(TEM-1) gene encoding TEM-1 beta-lactamase was inserted into the E. coli pyrF chromosomal locus. Expression of wild-type BLIP from a plasmid in this strain resulted in a large decrease in ampicillin resistance, while introduction of the same plasmid lacking BLIP had no effect on ampicillin resistance. In addition, it was found that when the BLIP alanine-scanning mutants were tested in the strain, the level of ampicillin resistance was proportional to the K(i) of the BLIP mutant. These results indicate that BLIP function can be monitored by the level of ampicillin resistance of the genetic test strain. Each of the 23 BLIP positions examined by alanine scanning was randomized to create libraries containing all possible substitutions at each position. The genetic screen for BLIP function was used to sort the libraries for active mutants, and DNA sequence analysis of functional BLIP mutants identified the sequences required for binding TEM-1 beta-lactamase. The results indicate the BLIP surface is tolerant of substitutions in that many contact positions can be substituted with other amino acid types and retain wild-type levels of function.

Computational redesign of the SHV-1 beta-lactamase/beta-lactamase inhibitor protein interface

Beta-lactamases are enzymes that catalyze the hydrolysis of beta-lactam antibiotics. beta-lactamase/beta-lactamase inhibitor protein (BLIP) complexes are emerging as a well characterized experimental model system for studying protein-protein interactions. BLIP is a 165 amino acid protein that inhibits several class A beta-lactamases with a wide range of affinities: picomolar affinity for K1; nanomolar affinity for TEM-1, SME-1, and BlaI; but only micromolar affinity for SHV-1 beta-lactamase. The large differences in affinity coupled with the availability of extensive mutagenesis data and high-resolution crystal structures for the TEM-1/BLIP and SHV-1/BLIP complexes make them attractive systems for the further development of computational design methodology. We used EGAD, a physics-based computational design program, to redesign BLIP in an attempt to increase affinity for SHV-1. Characterization of several of designs and point mutants revealed that in all cases, the mutations stabilize the interface by 10- to 1000-fold relative to wild type BLIP. The calculated changes in binding affinity for the mutants were within a mean absolute error of 0.87 kcal/mol from the experimental values, and comparison of the calculated and experimental values for a set of 30 SHV-1/BLIP complexes yielded a correlation coefficient of 0.77. Structures of the two complexes with the highest affinity, SHV-1/BLIP (E73M) and SHV-1/BLIP (E73M, S130K, S146M), are presented at 1.7 A resolution. While the predicted structures have much in common with the experimentally determined structures, they do not coincide perfectly; in particular a salt bridge between SHV-1 D104 and BLIP K74 is observed in the experimental structures, but not in the predicted design conformations. This discrepancy highlights the difficulty of modeling salt bridge interactions with a protein design algorithm that approximates side chains as discrete rotamers. Nevertheless, while local structural features of the interface were sometimes miscalculated, EGAD is globally successful in designing complexes with increased affinity.

Modeling the structure of mAb 14B7 bound to the anthrax protective antigen

The anthrax protective antigen (PA) is a key component of the tripartite anthrax toxin. Monoclonal antibody (mAb) 14B7 and its engineered, affinity-matured variants have been shown to be effective in blocking PA binding to cellular receptors and mitigating anthrax toxicity. Here, we perform computational structural modeling of the mAb 14B7-PA interaction. Our objectives are to determine the structure of the 14B7-PA complex, to deduce a structural explanation for the affinity maturation from the docking models, and to study the effect of inaccuracies in the antibody homology model on docking. We used the RosettaDock program to dock PA with the mAb 14B7 crystal structure or homology model. Our simulations generate two distinct binding orientations consistent with experimental residue mutations that diminish 14B7-PA binding. Furthermore, the models suggest new site-directed mutations to positively identify one of these two solutions as the correct 14B7-PA docking orientation. The models indicate that PA regions 648-660 and 712-720 may be important for 14B7 binding in addition to the known PA epitope, and the binding interfaces are similar to that seen in the PA complex with cellular receptor CMG2. Antibody residues involved in affinity maturation do not contact the antigen in the docking models, suggesting that affinity maturation in the 14B7 family does not result from direct enhancements of antibody-antigen contacts. Docking the homology model produces low-resolution representations of the crystal structure docking orientations, but homology model docking is frustrated by antibody H3 loop conformation errors. This work demonstrates the usefulness and limitations of computational structure prediction for the development of antibody therapeutics, and reemphasizes the need for flexible backbone docking algorithms to achieve high-resolution docking using homology models.

Computational design of antibody-affinity improvement beyond in vivo maturation

Antibodies are used extensively in diagnostics and as therapeutic agents. Achieving high-affinity binding is important for expanding detection limits, extending dissociation half-times, decreasing drug dosages and increasing drug efficacy. However, antibody-affinity maturation in vivo often fails to produce antibody drugs of the targeted potency, making further affinity maturation in vitro by directed evolution or computational design necessary. Here we present an iterative computational design procedure that focuses on electrostatic binding contributions and single mutants. By combining multiple designed mutations, a tenfold affinity improvement to 52 pM was engineered into the anti-epidermal growth factor receptor drug cetuximab (Erbitux), and a 140-fold improvement in affinity to 30 pM was obtained for the anti-lysozyme model antibody D44.1. The generality of the methods was further demonstrated through identification of known affinity-enhancing mutations in the therapeutic antibody bevacizumab (Avastin) and the model anti-fluorescein antibody 4-4-20. These results demonstrate computational capabilities for enhancing and accelerating the development of protein reagents and therapeutics.

The role of charged surface residues in the binding ability of small hubs in protein-protein interaction networks

Hubs are highly connected proteins in a protein-protein interaction network. Previous work has implicated disordered domains and high surface charge as the properties significant in the ability of hubs to bind multiple proteins. While conformational flexibility of disordered domains plays an important role in the binding ability of large hubs, high surface charge is the dominant property in small hubs. In this study, we further investigate the role of the high surface charge in the binding ability of small hubs in the absence of disordered domains. Using multipole expansion, we find that the charges are highly distributed over the hub surfaces. Residue enrichment studies show that the charged residues in hubs are more prevalent on the exposed surface, with the exception of Arg, which is predominantly found at the interface, as compared to non-hubs. This suggests that the charged residues act primarily from the exposed surface rather than the interface to affect the binding ability of small hubs. They do this through (i) enhanced intra-molecular electrostatic interactions to lower the desolvation penalty, (ii) indirect long - range intermolecular interactions with charged residues on the partner proteins for better complementarity and electrostatic steering, and (iii) increased solubility for enhanced diffusion-controlled rate of binding. Along with Arg, we also find a high prevalence of polar residues Tyr, Gln and His and the hydrophobic residue Met at the interfaces of hubs, all of which have the ability to form multiple types of interactions, indicating that the interfaces of hubs are optimized to participate in multiple interactions.

Progress in computational protein design

Current progress in computational structure-based protein design is reviewed in the areas of methodology and applications. Foundational advances include new potential functions, more efficient ways of computing energetics, flexible treatments of solvent, and useful energy function approximations, as well as ensemble-based approaches to scoring designs for inclusion of entropic effects, improvements to guaranteed and to stochastic search techniques, and methods to design combinatorial libraries for screening and selection. Applications include new approaches and successes in the design of specificity for protein folding, binding, and catalysis, in the redesign of proteins for enhanced binding affinity, and in the application of design technology to study and alter enzyme catalysis. Computational protein design continues to mature and advance.

A Preference for Edgewise Interactions between Aromatic Rings and Carboxylate Anions: The Biological Relevance of Anion−Quadrupole Interactions

Noncovalent interactions are quite important in biological structure-function relationships. To study the pairwise interaction of aromatic amino acids (phenylalanine, tyrosine, tryptophan) with anionic amino acids (aspartic and glutamic acids), small molecule mimics (benzene, phenol or indole interacting with formate) were used at the MP2 level of theory. The overall energy associated with an anion-quadrupole interaction is substantial (-9.5 kcal/mol for a benzene-formate planar dimer at van der Waals contact distance), indicating the electropositive ring edge of an aromatic group can interact with an anion. Deconvolution of the long-range coplanar interaction energy into fractional contributions from charge-quadrupole interactions, higher-order electrostatic interactions, and polarization terms was achieved. The charge-quadrupole term contributes between 30 to 45% of the total MP2 benzene-formate interaction; most of the rest of the interaction arises from polarization contributions. Additional studies of the Protein Data Bank (PDB Select) show that nearly planar aromatic-anionic amino acid pairs occur more often than expected from a random angular distribution, while axial aromatic-anionic pairs occur less often than expected; this demonstrates the biological relevance of the anion-quadrupole interaction. While water may mitigate the strength of these interactions, they may be numerous in a typical protein structure, so their cumulative effect could be substantial.

Rapid Estimation of Relative Protein−Ligand Binding Affinities Using a High-Throughput Version of MM-PBSA

By employing a modified protocol of the Molecular Mechanics with Poisson-Boltzmann Surface Area (MM-PBSA) methodology we substantially decrease the required computation time for calculating relative estimates of protein-ligand binding affinities. The modified method uses a generalized Born implicit solvation model during molecular dynamics to enhance conformational sampling as well as a very efficient Poisson-Boltzmann solver and a computational design based on a distributed-computing paradigm. This construction allows for reduction of the computational cost of the calculations by roughly 2 orders of magnitude compared to the traditional formulation of MM-PBSA. With this high-throughput version of MM-PBSA we show that one can produce efficient physics-based estimates of relative binding free energies with reasonable correlation to experimental data and a total computation time that is sufficiently low such that an industrially relevant throughput can be realized given currently accessible computing resources. We demonstrate this approach by performing a comparison of different MM-PBSA implementations on a set of 18 ligands for the protein target urokinase.

Structural and computational characterization of the SHV-1 beta-lactamase-beta-lactamase inhibitor protein interface

Beta-lactamase inhibitor protein (BLIP) binds a variety of class A beta-lactamases with affinities ranging from micromolar to picomolar. Whereas the TEM-1 and SHV-1 beta-lactamases are almost structurally identical, BLIP binds TEM-1 approximately 1000-fold tighter than SHV-1. Determining the underlying source of this affinity difference is important for understanding the molecular basis of beta-lactamase inhibition and mechanisms of protein-protein interface specificity and affinity. Here we present the 1.6A resolution crystal structure of SHV-1.BLIP. In addition, a point mutation was identified, SHV D104E, that increases SHV.BLIP binding affinity from micromolar to nanomolar. Comparison of the SHV-1.BLIP structure with the published TEM-1.BLIP structure suggests that the increased volume of Glu-104 stabilizes a key binding loop in the interface. Solution of the 1.8A SHV D104K.BLIP crystal structure identifies a novel conformation in which this binding loop is removed from the interface. Using these structural data, we evaluated the ability of EGAD, a program developed for computational protein design, to calculate changes in the stability of mutant beta-lactamase.BLIP complexes. Changes in binding affinity were calculated within an error of 1.6 kcal/mol of the experimental values for 112 mutations at the TEM-1.BLIP interface and within an error of 2.2 kcal/mol for 24 mutations at the SHV-1.BLIP interface. The reasonable success of EGAD in predicting changes in interface stability is a promising step toward understanding the stability of the beta-lactamase.BLIP complexes and computationally assisted design of tight binding BLIP variants.