Characterization of phosphate binding in the active site of barnase by site-directed mutagenesis and NMR

Phosphate is a competitive inhibitor of transesterification of GpC by the ribonuclease barnase. Barnase is significantly stabilized in the presence of phosphate against urea denaturation. The data are consistent with the existence of a single phosphate binding site in barnase with a dissociation constant, Kd, of 1.3 mM. The 2D 1H NMR spectrum of wild-type barnase with bound phosphate is assigned. Changes in chemical shifts and NOEs for wild type with bound phosphate compared with free wild type indicate that phosphate binds in the active site and that only small conformational changes occur on binding. Site-directed mutagenesis of the active site residues His-102, Lys-27, and Arg-87 to Ala increases the magnitude of Kd for phosphate by more than 20-fold. The 2D 1H NMR spectra of the mutants His-102----Ala, Lys-27----Ala, and Arg-87----Ala are assigned. Comparison with the spectra of wild-type barnase reveals that His-102----Ala and Lys-27----Ala have essentially the same structure as weild type, while some structural changes occur in Arg-87----Ala. It appears that phosphate binding by barnase is effected mainly by positively charge residues including His-102, Lys-27, and Arg-87. This may have applications for the design of phosphate binding sites in other proteins.     

Interaction energy decomposition in protein-protein association: a quantum mechanical study of barnase-barstar complex

Protein-protein interactions are very important in the function of a cell. Computational studies of these interactions have been of interest, but often they have utilized classical modelling techniques. In recent years, quantum mechanical (QM) treatment of entire proteins has emerged as a powerful approach to study biomolecular systems. Herein, we apply a semi-empirical divide and conquer (DC) methodology coupled with a dielectric continuum model for the solvent, to explore the contribution of electrostatics, polarization and charge transfer to the interaction energy between barnase and barstar in their complex form. Molecular dynamic (MD) simulation was performed to account for the dynamic behavior of the complex. The results show that electrostatics, charge transfer and polarization favor the formation of the complex. Our study shows that electrostatics dominates the interaction between barnase and barstar ( approximately 73%), while charge transfer and polarization are approximately 21% and approximately 6%, respectively. Close inspection of the polarization and charge-transfer effects on the charge distribution of the complex reveals the existence of two, well localized, regions in barstar. The first region includes the residues between P27 and Y47 and the second region is between N65 and D83. Since no such regions could be detected in barnase clearly suggests that barstar is well optimized for efficiently binding barnase. Furthermore, using our interaction energy decomposition scheme, we were able to identify all residues that have been experimentally determined to be important for the complex formation and to suggest other residues never have been investigated. This suggests that our approach will be useful as an aid in further understanding protein-protein contacts for the ultimate goal to produce successful inhibitors for protein complexes.     

Downhill binding energy surface of the barnase–barstar complex

We have employed biased molecular dynamics simulations in explicit solvent to characterize the one‐dimensional potential of mean force for the dissociation process of the barnase–barstar protein–protein complex. Unbinding of barstar from wild‐type barnase was compared with dissociation from four charge‐deletion mutants of barnase. Interestingly, we find in all cases that unbinding of barnase and barstar is an uphill process on a smooth, tilted energy landscape. The total free energy difference between the dissociated and bound state was similar for wild‐type barnase–barstar and for the R87A mutant of barnase. The values for the three other mutant barnase mutants K27A, R59A, and R83Q were only about half as much. Besides, we have analyzed the conformational dynamics of important residues at the barnase–barstar interface. In the bound state, their conformational fluctuations are reduced relatively to the free state because of the formation of intermolecular contacts. Interestingly, we find that some residues also show decreased mobility at intermediate stages of the unbinding process suggesting that these residues may be involved in the first contacts being formed on binding. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 977–985, 2010.     

Barstar is electrostatically optimized for tight binding to barnase

We used a novel charge optimization technique to study the small ribonuclease barnase and to analyze its interaction with a natural tight binding inhibitor, the protein barstar. The approach uses a continuum model to explicitly determine the charge distributions that lead to the most favorable electrostatic contribution to binding when competing desolvation and interaction effects are included. Given its backbone fold, barstar is electrostatically optimized for tight binding to barnase when compared with mutants where residues have been substituted with one of the 20 common amino acids. Natural proteins thus appear to use optimization of electrostatic interactions as one strategy for achieving tight binding.     

A novel cell ablation strategy blocks tobacco anther dehiscence.

We utilized a new cell ablation strategy to ablate specific anther cell types involved in the dehiscence process. The tobacco TA56 gene promoter is active within the circular cell cluster, stomium, and connective regions of the anther at different developmental stages. We introduced a cytotoxic TA56/barnase gene into tobacco plants together with three different anticytotoxic barstar genes. The anticytotoxic barstar genes were used to protect subsets of anther cell types from the cytotoxic effects of the TA56/barnase gene. The chimeric barstar genes were fused with (1) the tobacco TP12 gene promoter that is active at high levels in most anther cell types; (2) the soybean lectin gene promoter that is active earlier in the connective, and at lower levels in the circular cell cluster and stomium, than is the TA56 promoter; and (3) the tobacco TA20 gene promoter that is active at high levels in most anther cell types but has a different developmental profile than does the TP12 promoter. Normal anther development and dehiscence occurred in plants containing the TA56/barnase and TP12/barstar genes, indicating that barstar protects diverse anther cell types from the cytotoxic effects of barnase. Anthers containing the TA56/barnase and lectin/barstar genes also developed normally but failed to dehisce because of extensive ablation of the circular cell cluster, stomium, and contiguous connective regions. Anthers containing the TA56/barnase and TA20/barstar genes failed to dehisce as well. However, only the stomium region was ablated in these anthers. The connective, circular cell cluster, and adjacent wall regions were protected from ablation by the formation of barnase/barstar complexes. We conclude that anther dehiscence at flower opening depends on the presence of a functional stomium region and that chimeric barnase and barstar genes containing promoters that are active in several overlapping cell types can be used for targeted cell ablation experiments.     

How optimal are the binding energetics of barnase and barstar?

The extracellular ribonuclease barnase and its intracellular inhibitor barstar bind fast and with high affinity. Although extensive experimental and theoretical studies have been carried out on this system, it is unclear what the relative importance of different contributions to the high affinity is and whether binding can be improved through point mutations. In this work, we first applied Poisson-Boltzmann electrostatic calculations to 65 barnase-barstar complexes with mutations in both barnase and barstar. The continuum electrostatic calculations with a van der Waals surface dielectric boundary definition result in the electrostatic interaction free energy providing the dominant contribution favoring barnase-barstar binding. The results show that the computed electrostatic binding free energy can be improved through mutations at W44/barstar and E73/barnase. Furthermore, the determinants of binding affinity were quantified by applying COMparative BINding Energy (COMBINE) analysis to derive quantitative structure-activity relationships (QSARs) for the 65 complexes. The COMBINE QSAR model highlights approximately 20 interfacial residue pairs as responsible for most of the differences in binding affinity between the mutant complexes, mainly due to electrostatic interactions. Based on the COMBINE model, together with Brownian dynamics simulations to compute diffusional association rate constants, several mutants were designed to have higher binding affinities than the wild-type proteins.     

Barnase and barstar. Expression of its cloned inhibitor permits expression of a cloned ribonuclease

Barnase is the extracellular ribonuclease of Bacillus amyloliquefaciens and barstar its specific intracellular inhibitor. The gene for barstar has now been cloned and sequenced. When the wild-type gene for barnase is reconstructed from its previously cloned parts on the same plasmid as the barstar gene, the lethal effect of its expression is suppressed. A plasmid has been devised which directs the secretion of 100 mg per active barnase liter by Escherichia coli and another which provides large (500 to 1000 mg/l) yields of barstar. The structure of these plasmids and the derived 89 amino acid sequence of barstar are reported.     

Diffusional encounter of barnase and barstar

We present an analysis of trajectories from Brownian dynamics simulations of diffusional protein-protein encounter for the well-studied system of barnase and barstar. This analysis reveals details about the optimal association pathways, the regions of the encounter complex, possible differences of the pathways for dissociation and association, the coupling of translational and rotation motion, and the effect of mutations on the trajectories. We found that a small free-energy barrier divides the energetically most favorable region into a region of the encounter complex above the barnase binding interface and a region around a second energy minimum near the RNA binding loop. When entering the region of the encounter complex from the region near the RNA binding loop, barstar has to change its orientation to increase the electrostatic attraction between the proteins. By concentrating the analysis on the successful binding trajectories, we found that the region of the second minimum is not essential for the binding of barstar to barnase. Nevertheless, this region may be helpful to steer barstar into the region of the encounter complex. When applying the same analysis to several barnase mutants, we found that single mutations may drastically change the free-energy landscape and may significantly alter the population of the two minima. Therefore, certain protein-protein pairs may require careful adaptation of the positions of encounter and transition states when interpreting mutation effects on kinetic rates of association and/or dissociation.     

Structural response to mutation at a protein-protein interface

We have crystallised three mutants of the barnase-barstar complex in which interactions across the interface have been deleted by simultaneous mutation of both residues involved in the interaction. Each mutant deletes a different type of interaction at the interface: the first complex bnHis102-->Ala-bsTyr29-->Phe (bn, barnase; bs, barstar), deletes a van der Waals packing interaction; the second complex, bnLys27-->Ala-bsThr42-->Ala, deletes a hydrogen bond; the third, bnLys27-->Ala-bsAsp35-->Ala, deletes a long-range charge-charge interaction. The contribution of each of these side-chains to the stability of the complex is known; the coupling energy between the deleted side-chains is also known. Despite each of the double mutants being significantly destabilised compared with the wild-type, the effects of mutation are local. Only small movements in the main-chain surrounding the sites of mutation and some larger movements of neighbouring side-chains are observed in the mutant complexes. The exact response to mutation is context-dependent and for the same mutant can vary depending upon the environment within the crystal. In some double mutant complexes, interfacial pockets, which are accessible to bulk solvent are formed, whereas interfacial cavities which are isolated from bulk solvent, are formed in others. In all double mutants, water molecules fill the created pockets and cavities. These water molecules mimic the deleted side-chains by occupying positions close to the non-carbon atoms of truncated side-chains and re-making many hydrogen bonds made by the truncated side-chains in the wild-type. It remains extremely difficult, however, to correlate energetic and structural responses to mutation because of unknown changes in entropy and entropy-enthalpy compensation.     

Foldability, enzymatic activity, and interacting ability of barnase mutants obtained by permutation of secondary structure units

Barnase, a well-characterized ribonuclease, has been decomposed into six modules (M1-M6) or secondary structure units (S1-S6). We have studied the foldability and activity of the barnase mutants obtained by permutation of the four internal modules (M2-M5) or secondary structure units (S2-S5) to investigate whether permutation of these building blocks is a useful way to create foldable and/or functional proteins. In this study, we found that one of the secondary structure unit mutants was expressed in Escherichia coli only when His102 was substituted by alanine, which is a catalytic residue of wild-type barnase. This mutant (S2354H102A) had ordered conformations, which unfolded cooperatively during urea-induced unfolding experiments. S2354H102A interacted with other barnase mutants to show a distinct RNase activity, although its own activity was quite weak. This interaction was specific, because S2354H102A interacted with only barnase mutants having His 102 and certain orders of the secondary structure units giving a distinct RNase activity. These results suggest that secondary structure units permuted in barnase mutants maintain their intrinsic "interacting ability" that is used for the folding of wild-type barnase, and the units can form certain conformations that complement those of the appropriate counterparts. Seven of 23 secondary structure unit mutants and only 2 of 23 module mutants had RNase activity. On the basis of the results of analyses of foldability and RNase activity of the mutants performed in this and previous studies, we conclude that secondary structure units are more suitable than modules as building blocks to create novel foldable and/or functional proteins in the case of barnase.     

Application of fusion protein 4D5 scFv-dibarnase:barstar-gold complex for studying P185HER2 receptor distribution in human cancer cells

Overexpression of the P185^HER2 protein determines the malignancy and unfavorable prognosis of ovarian and breast tumors. In this work, the distribution of P185^HER2 in human cancer cells was studied by electron microscopy, using a novel approach. It is based on the interaction between barnase (a ribonuclease from Bacillus amyloliquefaciens) and its specific inhibitor barstar. The monoclonal antibody 4D5 scFv to extracellular P185^HER2 domain fused with two molecules of barnase was used as a recognizing agent, and the conjugate of colloidal gold with barstar, as an electron dense label for electron microscopic visualization. For labeling, we used supramolecular complexes 4D5 scFv-dibarnase:barstar-Au.     

Design of multivalent complexes using the barnase*barstar module

The ribonuclease barnase (12 kDa) and its inhibitor barstar (10 kDa) form a very tight complex in which all N and C termini are accessible for fusion. Here we exploit this system to create modular targeting molecules based on antibody scFv fragment fusions to barnase, to two barnase molecules in series and to barstar. We describe the construction, production and purification of defined dimeric and trimeric complexes. Immobilized barnase fusions are used to capture barstar fusions from crude extracts to yield homogeneous, heterodimeric fusion proteins. These proteins are stable, soluble and resistant to proteolysis. Using fusions with anti-p185(HER2-ECD) 4D5 scFv, we show that the anticipated gain in avidity from monomer to dimer to trimer is obtained and that favorable tumor targeting properties are achieved. Many permutations of engineered multispecific fusion proteins become accessible with this technology of quasi-covalent heterodimers.     

A comparative study of the interaction of Bacillus amyloliquefaciens ribonuclease (barnase) and Bacillus intermedius ribonuclease (binase) with barstar by brownian dynamics simulation

A comparative study of the interaction of two RNases (binase and barnase) with the polypeptide inhibitor barstar was performed by Brownian dynamics simulation. It was demonstrated that this method adequately reproduced the dependence of the association rate on the pH of solution as well as the effect of mutations at individual amino acid residues on the inhibition of barnase by barstar. Two types of energy-favorable binase-barstar encounter complexes were found. In type I complex, the amino acid residues of the binase active center are involved in formation of the complex; in type II complex, the active center remains free. It is suggested that temporary binding of free barstar into type II complex competes with the inhibition reaction. Presumably, this explains the decrease in the rate of binase inhibition by barstar as compared with the analogous reaction of barnase.     

Recognition between a bacterial ribonuclease, barnase, and its natural inhibitor, barstar

BACKGROUND: Protein-protein recognition is fundamental to most biological processes. The information we have so far on the interfaces between proteins comes largely from several protease-inhibitor and antigen-antibody complexes. Barnase, a bacterial ribonuclease, and barstar, its natural inhibitor, form a tight complex which provides a good model for the study and design of protein-protein non-covalent interactions. RESULTS: Here we report the structure of a complex between barnase and a fully functional mutant of barstar determined by X-ray analysis. Barstar is composed of three parallel alpha-helices stacked against a three-stranded parallel, beta-sheet, and sterically blocks the active site of the enzyme with an alpha-helix and adjacent loop. The buried surface in the interface between the two molecules totals 1630 A2. The barnase-barstar complex is predominantly stabilized by charge interactions involving positive charges in the active site of the enzyme. Asp39 of barstar binds to the phosphate-binding site of barnase, mimicking enzyme-substrate interactions. CONCLUSION: The phosphate-binding site of the enzyme is the anchor point for inhibitor binding. We propose that this is also likely to be the case for other ribonuclease inhibitors.     

Propagation of dynamic changes in barnase upon binding of barstar: an NMR and computational study

NMR spectroscopy and computer simulations were used to examine changes in chemical shifts and in dynamics of the ribonuclease barnase that result upon binding to its natural inhibitor barstar. Although the spatial structures of free and bound barnase are very similar, binding results in changes of the dynamics of both fast side-chains, as revealed by (2)H relaxation measurements, and NMR chemical shifts in an extended beta-sheet that is located far from the binding interface. Both side-chain dynamics and chemical shifts are sensitive to variations in the ensemble populations of the inter-converting molecular states, which can escape direct structural observation. Molecular dynamics simulations of free barnase and barnase in complex with barstar, as well as a normal mode analysis of barnase using a Gaussian network model, reveal relatively rigid domains that are separated by the extended beta-sheet mentioned above. The observed changes in NMR parameters upon ligation can thus be rationalized in terms of changes in inter-domain dynamics and in populations of exchanging states, without measurable structural changes. This provides an alternative model for the propagation of a molecular response to ligand binding across a protein that is based exclusively on changes in dynamics.     

Ribonuclease-charged vector for facile direct cloning with positive selection

Plasmid vectors for positive selection of cloned inserts in Escherichia coli were devised, based on an expression plasmid (pMT416) for the bacterial ribonuclease barnase. In addition to the barnase gene under control of a synthetic tac promoter, these plasmids carry the gene for the barnase inhibitor, barstar, the constitutive expression of which protects the bacterium from the detrimental effects of moderate barnase production. Full expression of the barnase gene overcomes protection by barstar and becomes lethal. Having a unique SmaI/XmaI site in the barnase structural gene, pMT416 itself can be used as a selective vector: uncut or religated pMT416 will preclude growth while plasmids with inserts in the barnase gene will allow the cells to survive. The entire pUC polylinker was inserted into the barnase gene in place of the Val-36 codon. This insert of nineteen largely hydrophilic amino acids does not prevent the lethal effect of full expression of the gene. The resulting plasmid, pMT440, is a generally useful selective cloning vector representing the “kill-the-rest” approach. Received: 25 April 1998 / Accepted: 26 May 1998     

I87E mutation prevents barstar dimerization

The C40,82A;I87E mutant of barstar, an intracellular inhibitor of the ribonuclease barnase from Bacillus amyloliquefaciens, was obtained, and its physicochemical properties were studied. It was produced as a fusion protein with thioredoxin and then cleaved from this by EKmax enterokinase. The mutant was shown by NMR to retain the spatial structure of the wild-type protein but, in contrast to barstar, does not form the homodimers characteristic of barstar in aqueous solution. The mutant protein binds barnase with the dissociation constant (6.6 +/- 1.1) x 10(-11) M and exhibits other physicochemical properties similar to those of the wild-type barstar. This allows the use of C40,82A;I87E mutant instead of wild-type barstar in investigations where the protein dimerization is undesirable. The English version of the paper: Russian Journal of Bioorganic Chemistry, 2004, vol. 30, no. 6; see also http://www.maik.ru.     

Simulation of the diffusional association of barnase and barstar.

The rate of protein association places an upper limit on the response time due to protein interactions, which, under certain circumstances, can be diffusion-controlled. Simulations of model proteins show that diffusion-limited association rates are approximately 10(6)-10(7) M-1 s-1 in the absence of long-range forces (Northrup, S. H., and H. P. Erickson. 1992. Kinetics of protein-protein association explained by Brownian dynamics computer simulations. Proc. Natl. Acad. Sci. U.S.A. 89:3338-3342). The measured association rates of barnase and barstar are 10(8)-10(9) M-1 s-1 at 50 mM ionic strength, and depend on ionic strength (Schreiber, G., and A. R. Fersht. 1996. Rapid, electrostatically assisted association of proteins. Nat. Struct. Biol. 3:427-431), implying that their association is electrostatically facilitated. We report Brownian dynamics simulations of the diffusional association of barnase and barstar to compute association rates and their dependence on ionic strength and protein mutation. Crucial to the ability to reproduce experimental rates is the definition of encounter complex formation at the endpoint of diffusional motion. Simple definitions, such as a required root mean square (RMS) distance to the fully bound position, fail to explain the large influence of some mutations on association rates. Good agreement with experiments could be obtained if satisfaction of two intermolecular residue contacts was required for encounter complex formation. In the encounter complexes, barstar tends to be shifted from its position in the bound complex toward the guanine-binding loop on barnase.     

Bisexual sterility conferred by the differential expression of Barnase and Barstar: a simple and efficient method of transgene containment

To establish a simple and an efficient system to minimize the environmental risk of genetically modified plants, we tested the applicability of the barnase/barstar system in conferring bisexual sterility; that is, in preventing plants setting seeds by self-fertilization and out-crossing. Transgenic tobacco plants were generated to express barnase, a cell death inducing ribonuclease, under the control of the gamete-specific AtDMC1 promoter, and barstar, a specific inhibitor of barnase, under the control of the ACT2 promoter, which is constitutively active in almost all tissues except gametes. In contrast to control plants harboring the barstar expression unit only, which set seeds normally with self-pollination, all transformants harboring both barnase and barstar were bisexually sterile. They produced aberrant anthers containing no detectable pollen and failed to set seeds even after pollination with wild-type tobacco pollen.An erratum to this article can be found at