Acetylcholinesterase: Role of the enzyme's charge distribution in steering charged ligands toward the active site

The electrostatic steering of charged ligands toward the active site of Torpedo californica acetylcholinesterase is investigated by Brownian dynamics simulations of wild type enzyme and several mutated forms, in which some normally charged residues are neutralized. The simulations reveal that the total ligand influx through a surface of 42 Å radius centered in the enzyme monomer and separated from the protein surface by I-14 Å is not significantly influenced by electrostatic interactions. Electrostatic effects are visible for encounters with a surface of 32 Å radius, which is partially hidden inside the protein, but mostly within the solvent. A clear accumulation of encounter events for that sphere is observed in the area directly above the entrance to the active site gorge. In this area, the encounter events are increased by 40% compared to the case of a neutral ligand. However, the differences among the encounter rates for the various mutants considered here are not pronounced, all rate constants being within ±10% of the average value. The enzyme charge distribution becomes more important as the charged ligand moves toward the bottom of the gorge, where the active site is located. We show that neither the enzyme's total charge, nor its dipole moment, fully account for the electrostatic steering of ligand to the active site. Higher moments of the enzyme's charge distribution are also important. However, for a series of mutations for which the direction of the enzyme dipole moment is constant within a few degrees, one observes a gradual decrease in the diffusional encounter rate constant with the number of neutralized residues. On the other hand, for other mutants that change the direction of the dipole moment from that of the wild type, the calculated encounter rate constants can be very close to that of the wild type. The present work yields two new insights to the kinetics of acetylcholinesterase. First, evolution appears to have built a redundant electrostatic steering capability into this important enzyme through the overall distribution of its thousands of partially charged atoms. And second, roughly half of the rate enhancement due to electrostatics arises from steering of the substrate outside the enzyme; the other half of the rate enhancement arises from improved trapping of the substrate after it has entered the gorge. The computational results reproduce qualitatively, and help to rationalize, many surprising experimental results obtained recently for human acetylcholinesterase. © 1996 John Wiley & Sons, Inc.     

Multiple drugs and multiple targets: An analysis of the electrostatic determinants of binding between non‐nucleoside HIV‐1 reverse transcriptase inhibitors and variants of HIV‐1 RT

We present a systematic, computational analysis of the electrostatic component of binding of three HIV‐1 RT inhibitors—nevirapine (NVP), efavirenz (EFV), and the recently approved rilpivirine (RPV)—to wild‐type (WT) and mutant variants of RT. Electrostatic charge optimization was applied to determine how suited each molecule's charge distribution is for binding WT and individual mutants of HIV‐1 RT. Although the charge distributions of NVP and EFV are rather far from being optimal for tight binding, RPVs charge distribution is close to the theoretical, optimal charge distribution for binding WT HIV‐1 RT, although slight changes in charge can dramatically impact binding energetics. Moreover, toward the L100I/K103N double mutant, RPVs charge distribution is quite far from optimal. We also determine the contributions of chemical moieties on each molecule toward the electrostatic component of binding and show that different regions of a drug molecule may be used for recognition by different RT variants. The electrostatic contributions of certain RT residues toward drug binding are also computed to highlight critical residues for each interaction. Finally, the charge distribution of RPV is optimized to promiscuously bind to three RT variants rather than to each one in turn, with the resulting charge distribution being a compromise between the optimal charge distributions to each individual variant. Taken together, this work demonstrates that even in a binding site considered quite hydrophobic, electrostatics play a subtle yet varying role that must be considered in designing next‐generation molecules that recognize rapidly mutating targets. Proteins 2012. © 2011 Wiley Periodicals, Inc.     

Neuroprotective effect of 2-(4-methoxyphenyl)ethyl-2-acetamido-2-deoxy-β-d-pyranoside against sodium nitroprusside-induced neurotoxicity in HT22 cells

In the present investigation, we used directed evolution approach to engineer a lipase from metagenomic origin. A variant S311C, was generated, characterized in detail and compared with wild type. Wild type and variant lipases were overexpressed and purified to homogeneity. The temperature optima of the purified lipases (Variant and wild type) were almost same, and found to be 45 and 50 °C, respectively. The variant protein was highly thermostable (54 times) as compared with the wild type at 60 °C. The variant displayed very high kinetic efficiency over the wild type protein. Analysis of the homology models of wild type and variant lipase showed that the substitution is on the surface of the protein. This substitution, along with hydrophobic residues in near vicinity may be involved in formation of strong hydrophobic channel leading to active site. This study identifies the role of hydrophobic interactions in protein stability along with enhancement of enzyme activity.     

Rapid binding of a cationic active site inhibitor to wild type and mutant mouse acetylcholinesterase: Brownian dynamics simulation including diffusion in the active site gorge

It is known that anionic surface residues play a role in the long-range electrostatic attraction between acetylcholinesterase and cationic ligands. In our current investigation, we show that anionic residues also play an important role in the behavior of the ligand within the active site gorge of acetylcholinesterase. Negatively charged residues near the gorge opening not only attract positively charged ligands from solution to the enzyme, but can also restrict the motion of the ligand once it is inside of the gorge. We use Brownian dynamics techniques to calculate the rate constant k_on for wild type and mutant acetylcholinesterase with a positively charged ligand. These calculations are performed by allowing the ligand to diffuse within the active site gorge. This is an extension of previously reported work in which a ligand was allowed to diffuse only to the enzyme surface. By setting the reaction criteria for the ligand closer to the active site, better agreement with experimental data is obtained. Although a number of residues influence the movement of the ligand within the gorge, Asp74 is shown to play a particularly important role in this function. Asp74 traps the ligand within the gorge, and in this way helps to ensure a reaction. © 1998 John Wiley & Sons, Inc. Biopoly 46: 465–474, 1998     

Synthesis and biochemical characterization of obligatory dimers of the sugar non-specific nuclease from Serratia marcescens using specifically designed bismaleimidoalkanes as SH-specific crosslinking reagents

The genetically engineered S140C variant of the homodimeric nuclease from Serratia marcescens was crosslinked across the dimer interface at the Cys 140 residues using bifunctional SH-specific 1,1-alkanediyl-bis-pyrrole-2,5-diones of different lengths. These bismaleimidoalkanes were synthesized by the condensation of n-alkyldiamines with maleic anhydride and subsequent cyclization with acetic anhydride and sodium acetate. Bismaleimidohexane (BMH) which gave the best crosslinking yield was used to produce in preparative amounts crosslinked Serratia nuclease. The crosslinked protein has the same secondary structure and exhibits the same guanidinium chloride unfolding behavior as the wild type enzyme or the non-covalently linked S140C variant. In contrast, in thermal unfolding experiments the crosslinked dimer behaves differently from the wild type enzyme or the non-covalently linked S140C variant. CD-spectra recorded during temperature rise showed only minor changes of the secondary structure composition for the wild type enzyme and the non-covalently linked S140C variant, whereas in the case of the crosslinked S140C dimer a distinct increase of the CD effect was observed corresponding to an increase in -helix. Our results demonstrate that bismaleimidoalkanes are very well suited to covalently link subunits of proteins, provided suitably located cysteine residues are present.     

Charge-directed targeting of antimicrobial protein-nanoparticle conjugates

Use of antimicrobial enzymes covalently attached to nanoparticles is of great interest as an antibiotic-free approach to treat microbial infections. Intrinsic properties of nanoparticles can also be used to add functionality to their conjugates with biomolecules. Here, we show in a model system that nanoparticle charge can be used to enhance delivery and increase bactericidal activity of an antimicrobial enzyme, lysozyme. Hen egg lysozyme was covalently attached to two types of polystyrene latex nanoparticles: positively charged, containing aliphatic amine surface groups, and negatively charged, containing sulfate and chloromethyl surface groups. In the case of bacterial lysis assay with a Gram-positive bacteria Micrococcus lysodeikticus, activity of lysozyme conjugated to positively charged nanoparticles was approximately twice as large as that of free lysozyme, while lysozyme conjugated to negatively charged nanoparticles showed little detectable activity. At the same time, when assayed using a low-molecular weight oligosaccharide substrate, lysozyme attached to both positively and negatively charged nanoparticles showed slightly lower activity than free enzyme. A possible explanation of these results is that lysozyme attached to negatively charged nanoparticles cannot be effectively targeted to the bacteria because of the electrostatic Coulombic repulsion from the negatively charged bacterial cell walls, whereas lysozyme conjugated to positively charged nanoparticles was targeted better than free enzyme due to stronger electrostatic attraction to bacteria. Zeta potential measurements confirmed the validity of this hypothesis. Thus, nanoparticle charge is an important factor that can be used to control targeting and activity of protein-nanoparticle conjugates.     

Structure of a stabilizing disulfide bridge mutant that closes the active-site cleft of T4 lysozyme.

The engineered disulfide bridge between residues 21 and 142 of phage T4 lysozyme spans the active-site cleft and can be used as a switch to control the activity of the enzyme (Matsumura, M. & Matthews, B.W., 1989, Science 243, 792-794). In the oxidized form the disulfide increases the melting temperature of the protein by 11 degrees C at pH 2. The crystal structure of this mutant lysozyme has been determined in both the reduced and oxidized forms. In the reduced form, the crystal structure of the mutant is shown to be extremely similar to that of wild type. In the oxidized form, however, the formation of the disulfide bridge causes the alpha-carbons of Cys 21 and Cys 142, on opposite sides of the active-site cleft, to move toward each other by 2.5 A. In association with this movement, the amino-terminal domain of the protein undergoes a rigid-body rotation of 5.1 degrees relative to the carboxy-terminal domain. This rotation occurs about an axis passing through the junction of the amino-terminal and carboxy-terminal domains and is also close to the axis that best fits the apparent thermal motion of the amino-terminal domain seen previously in crystals of wild-type lysozyme. Even though the engineered Cys 21-Cys 142 disulfide links together the amino-terminal and carboxy-terminal domains of T4 lysozyme, it does not reduce the apparent mobility of the one domain relative to the other. The pronounced "hinge-bending" mobility of the amino-terminal domain that is suggested by the crystallographic thermal parameters of wild-type lysozyme persists in the oxidized (and reduced) mutant structures. In the immediate vicinity of the introduced disulfide bridge the mutant structure is more mobile (or disordered) than wild type, so much so that the exact conformation of Cys 21 remains obscure. As with the previously described disulfide bridge between residues 9 and 164 of T4 lysozyme (Pjura, P.E., Matsumura, M., Wozniak, J.A., & Matthews, B.W., 1990, Biochemistry 29, 2592-2598), the engineered cross-link substantially enhances the stability of the protein without making the folded structure more rigid.     

Structural analysis of the temperature-sensitive mutant of bacteriophage T4 lysozyme, glycine 156----aspartic acid

The structure of the mutant of bacteriophage T4 lysozyme in which Gly-156 is replaced by aspartic acid is described. The lysozyme was isolated by screening for temperature-sensitive mutants and has a melting temperature at pH 6.5 that is 6.1 degrees C lower than wild type. The mutant structure is destabilized, in part, because Gly-156 has conformational angles (phi, psi) that are not optimal for a residue with a beta-carbon. High resolution crystallographic refinement of the mutant structure (R = 17.7% at 1.7 A resolution) shows that the Gly----Asp substitution does not significantly alter the configurational angles (phi, psi) but forces the backbone to move, as a whole, approximately 0.6 A away from its position in wild-type lysozyme. This induced strain weakens a hydrogen bond network that exists in the wild-type structure and also contributes to the reduced stability of the mutant lysozyme. The introduction of an acidic side chain reduces the overall charge on the molecule and thereby tends to increase the stability of the mutant structure relative to wild type. However, at neutral pH this generalized electrostatic stabilization is offset by specific electrostatic repulsion between Asp-156 and Asp-92. The activity of the mutant lysozyme is approximately 50% that of wild-type lysozyme. This reduction in activity might be due to introduction of a negative charge and/or perturbation of the surface of the molecule in the region that is assumed to interact with peptidoglycan substrates.     

Lysozyme dimerization: brownian dynamics simulation

The lysozyme dimerization reaction has been studied within the framework of encounter-complex (EC) formation theory using the MacroDox software package. Two types of energetically favorite ECs were determined. In the first of them, active-center amino acids of lysozyme take part in the complex formation or the second molecule blocks accessibility to active center sterically. Epitope amino-acid residues are involved in the complex of type II. The existence of both types of complexes does not contradict experimental data. Dimer-formation rate constants for different kinds of EC were calculated. Increasing the pH from 2.0 to 10.0 decreases the total positive lysozyme charge and eliminates the unfavorable repulsive electrostatic interaction. The rate constant of EC formation is inversely proportional to the protein total charge. The association rate constant was also enhanced by an increase of ionic strength that screened repulsive electrostatic interaction between positively charged proteins.     

Electrostatic attraction by surface charge does not contribute to the catalytic efficiency of acetylcholinesterase.

Acetylcholinesterases (AChEs) are characterized by a high net negative charge and by an uneven surface charge distribution, giving rise to a negative electrostatic potential extending over most of the molecular surface. To evaluate the contribution of these electrostatic properties to the catalytic efficiency, 20 single- and multiple-site mutants of human AChE were generated by replacing up to seven acidic residues, vicinal to the rim of the active-center gorge (Glu84, Glu285, Glu292, Asp349, Glu358, Glu389 and Asp390), by neutral amino acids. Progressive simulated replacement of these charged residues results in a gradual decrease of the negative electrostatic potential which is essentially eliminated by neutralizing six or seven charges. In marked contrast to the shrinking of the electrostatic potential, the corresponding mutations had no significant effect on the apparent bimolecular rate constants of hydrolysis for charged and non-charged substrates, or on the Ki value for a charged active center inhibitor. Moreover, the kcat values for all 20 mutants are essentially identical to that of the wild type enzyme, and the apparent bimolecular rate constants show a moderate dependence on the ionic strength, which is invariant for all the enzymes examined. These findings suggest that the surface electrostatic properties of AChE do not contribute to the catalytic rate, that this rate is probably not diffusion-controlled and that long-range electrostatic interactions play no role in stabilization of the transition states of the catalytic process.     

Genetically engineered alginate lyase-PEG conjugates exhibit enhanced catalytic function and reduced immunoreactivity

Alginate lyase enzymes represent prospective biotherapeutic agents for treating bacterial infections, particularly in the cystic fibrosis airway. To effectively deimmunize one therapeutic candidate while maintaining high level catalytic proficiency, a combined genetic engineering-PEGylation strategy was implemented. Rationally designed, site-specific PEGylation variants were constructed by orthogonal maleimide-thiol coupling chemistry. In contrast to random PEGylation of the enzyme by NHS-ester mediated chemistry, controlled mono-PEGylation of A1-III alginate lyase produced a conjugate that maintained wild type levels of activity towards a model substrate. Significantly, the PEGylated variant exhibited enhanced solution phase kinetics with bacterial alginate, the ultimate therapeutic target. The immunoreactivity of the PEGylated enzyme was compared to a wild type control using in vitro binding studies with both enzyme-specific antibodies, from immunized New Zealand white rabbits, and a single chain antibody library, derived from a human volunteer. In both cases, the PEGylated enzyme was found to be substantially less immunoreactive. Underscoring the enzyme's potential for practical utility, >90% of adherent, mucoid, Pseudomonas aeruginosa biofilms were removed from abiotic surfaces following a one hour treatment with the PEGylated variant, whereas the wild type enzyme removed only 75% of biofilms in parallel studies. In aggregate, these results demonstrate that site-specific mono-PEGylation of genetically engineered A1-III alginate lyase yielded an enzyme with enhanced performance relative to therapeutically relevant metrics.     

Anisotropic behavior of the pH dependence part of Young's modulus in a lysozyme triclinic crystal

A numerical simulation of the Young's modulus pH dependence in a triclinic crystal of hen egg white lysozyme is performed. On the basis of calculation of electrostatic interactions between the molecules in the crystal at different pH values, we obtain the Young's modulus quantity as a function of pH. The analysis is carried out along the main crystal axes as well as along the [011] direction. In the last case, a good agreement between our calculations of Young's modulus pH dependence and the experimental data [V. N. Morosov and T. Ya. Morozova (1981) Biopolymers, Vol. 20, pp. 451–467] is obtained. The calculations show a strong anisotropy of both the absolute value of the electrostatic part of Young's modulus and its pH dependence along different directions. A detailed analysis of the electrostatic potential inside the crystal, as well as the charged group distribution in the contact areas, permit us to find the reason for the anisotropy. It is caused by a nonuniform distribution of the charged residues in the lysozyme molecule and especially by the titration behavior of certain acidic groups. It is also shown that Young's modulus depends on the ionic strength, being anisotropic. The reason is that the cavities are not uniformly distributed within the crystal. They are located predominantly along the c axis. © 1995 John Wiley & Sons, Inc.     

The bifunctional active site of S-adenosylmethionine synthetase. Roles of the basic residues

S-adenosylmethionine (AdoMet) synthetase catalyzes a unique two-step enzymatic reaction leading to formation of the primary biological alkylating agent. The crystal structure of Escherichia coli AdoMet synthetase shows that the active site, which lies between two subunits, contains four lysines and one histidine as basic residues. In order to test the proposed charge and hydrogen bonding roles in catalytic function, each lysine has been changed to an uncharged methionine or alanine, and the histidine has been altered to asparagine. The resultant enzyme variants are all tetramers like the wild type enzyme; however, circular dichroism spectra show reductions in helix content for the K245*M and K269M mutants. (The asterisk denotes that the residue is in the second subunit.) Four mutants have k(cat) reductions of approximately 10(3)-10(4)-fold in AdoMet synthesis; however, the k(cat) of K165*M variant is only reduced 2-fold. In each mutant, there is a smaller catalytic impairment in the partial reaction of tripolyphosphate hydrolysis. The K165*A enzyme has a 100-fold greater k(cat) for tripolyphosphate hydrolysis than the wild type enzyme, but this mutant is not activated by AdoMet in contrast to the wild type enzyme. The properties of these mutants require reassessment of the catalytic roles of these residues.     

Electrostatics of Cytochrome-c assemblies

Electrostatic potentials along with computational mutagenesis are used to obtain atomic level insights into Cytochrome-c in order to design efficient bionanosensors. The electrostatic properties of wild type and mutant Cytochrome-c are examined in the context of their assembly, i.e. are examined in the absence and presence of neighboring molecules from the assembly. An intense increase in the positive potential ensues when the neighboring molecules are taken into account. This suggests that in the extrapolation of electric field effects upon the design of assemblies, considering the properties of only the central molecule may not be sufficient. Additionally, the influence of the uncharged residues becomes quite diminished when the molecule is considered in an assembly. This could pave the way for making mutants that might be more soluble in different media used in the construction of devices.FigureThe electrostatic potential, calculated using the program DELPHI [20] mapped on to the surface of Cytochrome-c when it is considered by itself (in the left column) and in the presence of the electrostatic field generated by the presence of the surrounding 4 molecules on the right (see Fig. 3). The potentials range from –10kT in red to +10kT in blue. The central figure shows the regions that have been mutated to positively charged residues by placing a unit positive charge at the terminal atom of the respective side chain. The figures range from the wild type in the first row, followed by the Gln12, Asn70, Asp50, Glu90 and Ala83 mutants.     

Surface charge engineering of PQQ glucose dehydrogenase for downstream processing

The ion-exchange chromatography behavior of recombinant glucose dehydrogenase harboring pyrroloquinoline quinone (PQQGDH) was modified to greatly simplify its purification. The surface charge of PQQGDH was engineered by either fusing a three-arginine tail to the C-terminus of PQQGDH (PQQGDH+Arg3) or by substituting three residues exposed on the surface of the enzyme to Arg by site-directed mutagenesis (3RPQQGDH). During cation exchange chromatography, both surface charge-engineered enzymes eluted at much higher salt concentrations than the wild-type enzyme. After the chromatography purification step, both PQQGDH+Arg3 and 3RPQQGDH appeared as single bands on SDS-PAGE, while extra bands appeared with the wild-type protein sample. Although all tested kinetic parameters of both engineered enzymes are similar to those of wild type, both modifications resulted in enzymes with increased thermal stability. Our achievements have resulted in the greater production of an improved quality PQQGDH by a simplified process.     

Inhibition of human pancreatic ribonuclease by the human ribonuclease inhibitor protein

The ribonuclease inhibitor protein (RI) binds to members of the bovine pancreatic ribonuclease (RNase A) superfamily with an affinity in the femtomolar range. Here, we report on structural and energetic aspects of the interaction between human RI (hRI) and human pancreatic ribonuclease (RNase 1). The structure of the crystalline hRI x RNase 1 complex was determined at a resolution of 1.95 A, revealing the formation of 19 intermolecular hydrogen bonds involving 13 residues of RNase 1. In contrast, only nine such hydrogen bonds are apparent in the structure of the complex between porcine RI and RNase A. hRI, which is anionic, also appears to use its horseshoe-shaped structure to engender long-range Coulombic interactions with RNase 1, which is cationic. In accordance with the structural data, the hRI.RNase 1 complex was found to be extremely stable (t(1/2)=81 days; K(d)=2.9 x 10(-16) M). Site-directed mutagenesis experiments enabled the identification of two cationic residues in RNase 1, Arg39 and Arg91, that are especially important for both the formation and stability of the complex, and are thus termed "electrostatic targeting residues". Disturbing the electrostatic attraction between hRI and RNase 1 yielded a variant of RNase 1 that maintained ribonucleolytic activity and conformational stability but had a 2.8 x 10(3)-fold lower association rate for complex formation and 5.9 x 10(9)-fold lower affinity for hRI. This variant of RNase 1, which exhibits the largest decrease in RI affinity of any engineered ribonuclease, is also toxic to human erythroleukemia cells. Together, these results provide new insight into an unusual and important protein-protein interaction, and could expedite the development of human ribonucleases as chemotherapeutic agents.     

Stabilization of substilisin E in organic solvents by site-directed mutagenesis

Subtilisin E was rationally engineered to improve its stability in polar organic solvents such as dimethylformamide (DMF). A charged surface residue, Asp248, was substituted by three amino acids of increasing hydrophobicity, Asn, Ala, and Leu; all three variants were stabilized with respect to wild type in 80% DMF. This stabilization was only observed in the presence of high concentrations of the organic solvent: no stability enhancements were observed in 40% DMF. In contrast, the mutation Asn218 --> Ser alters internal hydrogen bonding interactions and stabilizes subtilisin E in both 40% and 80% DMF. This study provides additional evidence that substitution of surface-charged residues is a generally useful mechanism for stabilizing enzymes in organic media and that the stabilizing effects of such substitutions are unique to highly altered solvent environments. The effects of the single amino acid substitutions on free energies of stabilization are additive in the Asp248 --> Asn + Asn218 --> Ser combination variant, yielding an enzyme that is 3.4 times more stable than wild type in 80% DMF.     

Increasing the hydrolysis constant of the reactive site upon introduction of an engineered Cys14?Cys39 bond into the ovomucoid third domain from silver pheasant

P14C/N39C is the disulfide variant of the ovomucoid third domain from silver pheasant (OMSVP3) introducing an engineered Cys¹⁴? Cys³⁹ bond near the reactive site on the basis of the sequence homology between OMSVP3 and ascidian trypsin inhibitor. This variant exhibits a narrower inhibitory specificity. We have examined the effects of introducing a Cys¹⁴? Cys³⁹ bond into the flexible N‐terminal loop of OMSVP3 on the thermodynamics of the reactive site peptide bond hydrolysis, as well as the thermal stability of reactive site intact inhibitors. P14C/N39C can be selectively cleaved by Streptomyces griseus protease B at the reactive site of OMSVP3 to form a reactive site modified inhibitor. The conversion rate of intact to modified P14C/N39C is much faster than that for wild type under any pH condition. The pH‐independent hydrolysis constant (K_hyd°) is estimated to be approximately 5.5 for P14C/N39C, which is higher than the value of 1.6 for natural OMSVP3. The reactive site modified form of P14C/N39C is thermodynamically more stable than the intact one. Thermal denaturation experiments using intact inhibitors show that the temperature at the midpoint of unfolding at pH 2.0 is 59 °C for P14C/N39C and 58 °C for wild type. There have been no examples, except P14C/N39C, where introducing an engineered disulfide causes a significant increase in K_hyd°, but has no effect on the thermal stability. The site‐specific disulfide introduction into the flexible N‐terminal loop of natural Kazal‐type inhibitors would be useful to further characterize the thermodynamics of the reactive site peptide bond hydrolysis. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

Second-site revertants of an inactive T4 lysozyme mutant restore activity by restructuring the active site cleft

Substitution of Thr26 by Gln in the lysozyme of bacteriophage T4 produces an enzyme with greatly reduced activity but essentially unaltered stability relative to wild type. Spontaneous second-site revertants of the mutant were selected genetically; two of them were chosen for structural and biochemical characterization. One revertant bears (in addition to the primary mutation) the substitution Tyr18----His, the other, Tyr18----Asp. The primary mutant and both revertant lysozyme genes were reconstructed in a plasmid-based expression system, and the proteins were produced and purified. The two revertant lysozymes exhibit enzymatic activities intermediate between wild type and the primary mutant; both also exhibit melting temperatures approximately 3 degrees C lower than either the wild type or the primary mutant. Crystals suitable for X-ray diffraction analysis were obtained from both revertant lysozymes, but not the primary mutant. Structures of the double mutant lysozymes were refined at 1.8-A resolution to crystallographic residuals of 15.1% (Tyr18----His) and 15.2% (Tyr18----Asp). Model building suggests that the side chain of Gln26 in the primary mutant is forced to protrude into the active site cleft, resulting in low catalytic activity. In contrast, the crystal structures of the revertants reveal that the double substitutions (Gln26 and His18, or Gln26 and Asp18) fit into the same space that is occupied by Thr26 and Tyr18 in the wild-type enzyme; the effect is a restructuring of the surface of the active site cleft, with essentially no perturbation of the polypeptide backbone. This restructuring is effected by a novel series of hydrogen bonds and electrostatic interactions that apparently stabilize the revertant structures.     

Cumulative site-directed charge-change replacements in bacteriophage T4 lysozyme suggest that long-range electrostatic interactions contribute little to protein stability

Bacteriophage T4 lysozyme is a basic molecule with an isoelectric point above 9.0, and an excess of nine positive charges at neutral pH. It might be expected that it would be energetically costly to bring these out-of-balance charges from the extended, unfolded, form of the protein into the compact folded state. To determine the contribution of such long-range electrostatic interactions to the stability of the protein, five positively charged surface residues, Lys16, Arg119, Lys135, Lys147 and Arg154, were individually replaced with glutamic acid. Eight selected double, triple and quadruple mutants were also constructed so as to sequentially reduce the out-of-balance formal charge on the molecule from +9 to +1 units. Each of the five single variant proteins was crystallized and high-resolution X-ray analysis confirmed that each mutant structure was, in general, very similar to the wild-type. In the case of R154E, however, the Arg154 to Glu replacement caused a rearrangement in which Asp127 replaced Glu128 as the capping residue of a nearby alpha-helix. The thermal stabilities of all 13 variant proteins were found to be fairly similar, ranging from 0.5 kcal/mol more stable than wild-type to 1.7 kcal/mol less stable than wild-type. In the case of the five single charge-change variants, for which the structures were determined, the changes in stability can be rationalized in terms of changes in local interactions at the site of the replacement. There is no evidence that the reduction in the out-of-balance charge on the molecule increases the stability of the folded relative to the unfolded form, either at pH 2.8 or at pH 5.3. This indicates that long-range electrostatic interactions between the substituted amino acid residues and other charged groups on the surface of the molecule are weak or non-existent. Furthermore, the relative stabilities of the multiple charge replacement mutant proteins were found to be almost exactly equal to the sums of the relative stabilities of the constituent single mutant proteins. This also clearly indicates that the electrostatic interactions between the replaced charges are negligibly small. The activities of the charge-change mutant lysozymes, as measured by the rate of hydrolysis of cell wall suspensions, are essentially equal to that of the wild-type lysozyme, but on a lysoplate assay the mutant enzymes appear to have higher activity.(ABSTRACT TRUNCATED AT 400 WORDS)