THz frequency spectrum of protein–solvent interaction energy using a recurrence plot‐based Wiener–Khinchin method

The dynamics of a protein and the water surrounding it are coupled via nonbonded energy interactions. This coupling can exhibit a complex, nonlinear, and nonstationary nature. The THz frequency spectrum for this interaction energy characterizes both the vibration spectrum of the water hydrogen bond network, and the frequency range of large amplitude modes of proteins. We use a Recurrence Plot based Wiener–Khinchin method RPWK to calculate this spectrum, and the results are compared to those determined using the classical auto‐covariance‐based Wiener–Khinchin method WK. The frequency spectra for the total nonbonded interaction energy extracted from molecular dynamics simulations between the β‐Lactamase Inhibitory Protein BLIP, and water molecules within a 10 Å distance from the protein surface, are calculated at 150, 200, 250, and 310 K, respectively. Similar calculations are also performed for the nonbonded interaction energy between the residues 49ASP, 53TYR, and 142PHE in BLIP, with water molecules within 10 Å from each residue respectively at 150, 200, 250, and 310 K. A comparison of the results shows that RPWK performs better than WK, and is able to detect some frequency data points that WK fails to detect. This points to the importance of using methods capable of taking the complex nature of the protein–solvent energy landscape into consideration, and not to rely on standard linear methods. In general, RPWK can be a valuable addition to the analysis tools for protein molecular dynamics simulations. Proteins 2016; 84:1549–1557. © 2016 Wiley Periodicals, Inc.     

The solution structure of the outer membrane lipoprotein OmlA from Xanthomonas axonopodis pv. citri reveals a protein fold implicated in protein-protein interaction

The outer membrane lipoprotein A (OmlA) belongs to a family of bacterial small lipoproteins widely distributed across the beta and gamma proteobacteria. Although the role of numerous bacterial lipoproteins is known, the biological function of OmlA remains elusive. We found that in the citrus canker pathogen, Xanthomonas axonopodis pv. citri (X. citri), OmlA is coregulated with the ferric uptake regulator (Fur) and their expression is enhanced when X. citri is grown on citrus leaves, suggesting that these proteins are involved in plant-pathogen interaction. To gain insights into the function of OmlA, its conformational and dynamic features were determined by nuclear magnetic resonance. The protein has highly flexible N- and C- termini and a structurally well defined core composed of three beta-strands and two small alpha-helices, which pack against each other forming a two-layer alpha/beta scaffold. This protein fold resembles the domains of the beta-lactamase inhibitory protein BLIP, involved in protein-protein binding. In conclusion, the structure of OmlA does suggest that this protein may be implicated in protein-protein interactions required during X. citri infection.     

Insights into Positive and Negative Requirements for Protein-Protein Interactions by Crystallographic Analysis of the β-Lactamase Inhibitory Proteins BLIP, BLIP-I, and BLP

β-Lactamase inhibitory protein (BLIP) binds a variety of β-lactamase enzymes with wide-ranging specificity. Its binding mechanism and interface interactions are a well-established model system for the characterization of protein-protein interactions. Published studies have examined the binding of BLIP to diverse target β-lactamases (e.g., TEM-1, SME-1, and SHV-1). However, apart from point mutations of amino acid residues, variability on the inhibitor side of this enzyme-inhibitor interface has remained unexplored. Thus, we present crystal structures of two likely BLIP relatives: (1) BLIP-I (solved alone and in complex with TEM-1), which has β-lactamase inhibitory activity very similar to that of BLIP; and (2) β-lactamase-inhibitory-protein-like protein (BLP) (in two apo forms, including an ultra-high-resolution structure), which is unable to inhibit any tested β-lactamase. Despite categorical differences in species of origin and function, BLIP-I and BLP share nearly identical backbone conformations, even at loop regions differing in BLIP.We describe interacting residues and provide a comparative structural analysis of the interactions formed at the interface of BLIP-I·TEM-1 versus those formed at the interface of BLIP·TEM-1. Along with initial attempts to functionally characterize BLP, we examine its amino acid residues that structurally correspond to BLIP/BLIP-I binding hotspots to explain its inability to bind and inhibit TEM-1. We conclude that the BLIP family fold is a robust and flexible scaffold that permits the formation of high-affinity protein-protein interactions while remaining highly selective. Comparison of the two naturally occurring, distinct binding interfaces built upon this scaffold (BLIP and BLIP-I) shows that there is substantial variation possible in the subnanomolar binding interaction with TEM-1. The corresponding (non-TEM-1-binding) BLP surface shows that numerous favorable backbone-backbone/backbone-side-chain interactions with a protein partner can be negated by the presence of a few, strongly unfavorable interactions, especially electrostatic repulsions.     

Binding hot spots in the TEM1-BLIP interface in light of its modular architecture

Proteins bind one another in aqua's solution to form tight and specific complexes. Previously we have shown that this is achieved through the modular architecture of the interaction network formed by the interface residues, where tight cooperative interactions are found within modules but not between them. Here we extend this study to cover the entire interface of TEM1 beta-lactamase and its protein inhibitor BLIP using an improved method for deriving interaction maps based on REDUCE to add hydrogen atoms and then by evaluating the interactions using modifications of the programs PROBE, NCI and PARE. An extensive mutagenesis study of the interface residues indeed showed that each module is energetically independent on other modules, and that cooperativity is found only within a module. By solving the X-ray structure of two interface mutations affecting two different modules, we demonstrated that protein-protein binding occur via the structural reorganization of the binding modules, either by a "lock and key" or an induced fit mechanism. To explain the cooperativity within a module, we performed multiple-mutant cycle analysis of cluster 2 resulting in a high-resolution energy map of this module. Mutant studies are usually done in reference to alanine, which can be regarded as a deletion of a side-chain. However, from a biological perspective, there is a major interest to understand non-Ala substitutions, as they are most common. Using X-ray crystallography and multiple-mutant cycle analysis we demonstrated the added complexity in understanding non-Ala mutations. Here, a double mutation replacing the wild-type Glu,Tyr to Tyr,Asn on TEM1 (res id 104,105) caused a major backbone structural rearrangement of BLIP, changing the composition of two modules but not of other modules within the interface. This shows the robustness of the modular approach, yet demonstrates the complexity of in silico protein design.     

Specificity and cooperativity at β-lactamase position 104 in TEM-1/BLIP and SHV-1/BLIP interactions

Establishing a quantitative understanding of the determinants of affinity in protein–protein interactions remains challenging. For example, TEM‐1/β‐lactamase inhibitor protein (BLIP) and SHV‐1/BLIP are homologous β‐lactamase/β‐lactamase inhibitor protein complexes with disparate K_d values (3 nM and 2 μM, respectively), and a single substitution, D104E in SHV‐1, results in a 1000‐fold enhancement in binding affinity. In TEM‐1, E104 participates in a salt bridge with BLIP K74, whereas the corresponding SHV‐1 D104 does not in the wild type SHV‐1/BLIP co‐structure. Here, we present a 1.6 Å crystal structure of the SHV‐1 D104E/BLIP complex that demonstrates that this point mutation restores this salt bridge. Additionally, mutation of a neighboring residue, BLIP E73M, results in salt bridge formation between SHV‐1 D104 and BLIP K74 and a 400‐fold increase in binding affinity. To understand how this salt bridge contributes to complex affinity, the cooperativity between the E/K or D/K salt bridge pair and a neighboring hot spot residue (BLIP F142) was investigated using double mutant cycle analyses in the background of the E73M mutation. We find that BLIP F142 cooperatively stabilizes both interactions, illustrating how a single mutation at a hot spot position can drive large perturbations in interface stability and specificity through a cooperative interaction network. Proteins 2011. © 2011 Wiley‐Liss, Inc.     

Detecting protein atom correlations using correlation of probability of recurrence

The dynamic cross‐correlation Map(DCCM) technique has been used extensively to study protein dynamics. In this work, we introduce the use of the method of correlation of probability of recurrence (CPR) as a complementary method to detect correlations between protein residue atoms. Time series of the distances of the C_α atoms of the β‐lactamase inhibitory protein (BLIP) from a reference position are analyzed using CPR and mutual information (MI). The results are compared to those provided by DCCM. In comparison to MI, CPR is found to detect more of the correlations present in DCCM. It is also able to detect a small number of significant correlations between distant residues that are not detected by DCCM. Proteins 2014; 82:2180–2189. © 2014 Wiley Periodicals, Inc.     

Protein�Cprotein binding affinities by pulse proteolysis: Application to TEM-1/BLIP protein complexes

Efficient methods for quantifying dissociation constants have become increasingly important for high‐throughput mutagenesis studies in the postgenomic era. However, experimentally determining binding affinity is often laborious, requires large amounts of purified protein, and utilizes specialized equipment. Recently, pulse proteolysis has been shown to be a robust and simple method to determine the dissociation constants for a protein–ligand pair based on the increase in thermodynamic stability upon ligand binding. Here, we extend this technique to determine binding affinities for a protein–protein complex involving the β‐lactamase TEM‐1 and various β‐lactamase inhibitor protein (BLIP) mutants. Interaction with BLIP results in an increase in the denaturation curve midpoint, C_m, of TEM‐1, which correlates with the rank order of binding affinities for several BLIP mutants. Hence, pulse proteolysis is a simple, effective method to assay for mutations that modulate binding affinity in protein–protein complexes. From a small set (n = 4) of TEM‐1/BLIP mutant complexes, a linear relationship between energy of stabilization (dissociation constant) and ΔC_m was observed. From this “calibration curve,” accurate dissociation constants for two additional BLIP mutants were calculated directly from proteolysis‐derived ΔC_m values. Therefore, in addition to qualitative information, armed with knowledge of the dissociation constants from the WT protein and a limited number of mutants, accurate quantitation of binding affinities can be determined for additional mutants from pulse proteolysis. Minimal sample requirements and the suitability of impure protein preparations are important advantages that make pulse proteolysis a powerful tool for high‐throughput mutagenesis binding studies.     

Interdomain flexibility and interfacial integrity of β-lactamase inhibitory protein (BLIP) modulate its binding to class A β-lactamases

β-Lactamase inhibitory protein (BLIP) consists of a tandem repeat of αβ domains conjugated by an interdomain loop and can effectively bind and inactivate class A β-lactamases, which are responsible for resistance of bacteria to β-lactam antibiotics. The varied ability of BLIP to bind different β-lactamases and the structural determinants for significant enhancement of BLIP variants with a point mutation are poorly understood. Here, we investigated the conformational dynamics of BLIP upon binding to three clinically prevalent class A β-lactamases (TEM1, SHV1, and PC1) with dissociation constants between subnanomolar and micromolar. Hydrogen deuterium exchange mass spectrometry revealed that the flexibility of the interdomain region was significantly suppressed upon strong binding to TEM1, but was not significantly changed upon weak binding to SHV1 or PC1. E73M and K74G mutations in the interdomain region improved binding affinity toward SHV1 and PC1, respectively, showing significantly increased flexibility of the interdomain region compared to the wild-type and favorable conformational changes upon binding. In contrast, more rigidity of the interfacial loop 135–145 was observed in these BLIP mutants in both free and bound states. Consistently, molecular dynamics simulations of BLIP exhibited drastic changes in the flexibility of the loop 135–145 in all complexes. Our results indicated for the first time that higher flexibility of the interdomain linker, as well as more rigidity of the interfacial loop 135–145, could be desirable determinants for enhancing inhibition of BLIP to class A β-lactamases. Together, these findings provide unique insights into the design of enhanced inhibitors.     

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

β-lactamases are enzymes that catalyze the hydrolysis of β-lactam antibiotics. β-lactamase/β-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 β-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 β-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 Å 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.