Electrical potentials in trypsin isozymes

Electrostatic free energies carry structural information on nucleic acid molecules in solution

Over the last several decades, a range of experimental techniques from x-ray crystallography and atomic force microscopy to nuclear magnetic resonance and small angle x-ray scattering have probed nucleic acid structure and conformation with high resolution both in the condensed state and in solution. We present a computational study that examines the prospect of using electrostatic free energy measurements to detect 3D conformational properties of nucleic acid molecules in solution. As an example, we consider the conformational difference between A- and B-form double helices whose structures differ in the values of two key parameters—the helical radius and rise per basepair. Mapping the double helix onto a smooth charged cylinder reveals that electrostatic free energies for molecular helices can, indeed, be described by two parameters: the axial charge spacing and the radius of a corresponding equivalent cylinder. We show that electrostatic free energies are also sensitive to the local structure of the molecular interface with the surrounding electrolyte. A free energy measurement accuracy of 1%, achievable using the escape time electrometry (ET e) technique, could be expected to offer a measurement precision on the radius of the double helix of approximately 1 Å. Electrostatic free energy measurements may, therefore, not only provide information on the structure and conformation of biomolecules but could also shed light on the interfacial hydration layer and the size and arrangement of counterions at the molecular interface in solution.

Isolating influential regions of electrostatic focusing in protein and DNA structure

Electrostatic focusing is a general phenomenon that occurs in cavities and grooves on the molecular surface of biomolecules. Narrow surface features can partially shield charged atoms from the high-dielectric solvent, enhancing electrostatic potentials inside the cavity and projecting electric field lines outward into the solvent. This effect has been observed in many instances and is widely considered in the human examination of molecular structure, but it is rarely integrated into the digital representations used in protein structure comparison software. To create a computational representation of electrostatic focusing, that is compatible with structure comparison algorithms, this paper presents an approach that generates three-dimensional solids that approximate regions where focusing occurs. We verify the accuracy of this representation against instances of focusing in proteins and DNA. Noting that this representation also identifies thin focusing regions on the molecular surface that are unlikely to affect binding, we describe a second algorithm that conservatively isolates larger focusing regions. The resulting 3D solids can be compared with Boolean set operations, permitting a new range of analyses on the regions where electrostatic focusing occurs. They also represent a novel integration of molecular shape and electrostatic focusing into the same structure comparison framework.

Assessment of two theoretical methods to estimate potentiometric titration curves of peptides: comparison with experiment

We compared the ability of two theoretical methods of pH-dependent conformational calculations to reproduce experimental potentiometric titration curves of two models of peptides: Ac-K5-NHMe in 95% methanol (MeOH)/5% water mixture and Ac-XX(A)7OO-NH2 (XAO) (where X is diaminobutyric acid, A is alanine, and O is ornithine) in water, methanol (MeOH), and dimethyl sulfoxide (DMSO), respectively. The titration curve of the former was taken from the literature, and the curve of the latter was determined in this work. The first theoretical method involves a conformational search using the electrostatically driven Monte Carlo (EDMC) method with a low-cost energy function (ECEPP/3 plus the SRFOPT surface-solvation model, assumming that all titratable groups are uncharged) and subsequent reevaluation of the free energy at a given pH with the Poisson-Boltzmann equation, considering variable protonation states. In the second procedure, molecular dynamics (MD) simulations are run with the AMBER force field and the generalized Born model of electrostatic solvation, and the protonation states are sampled during constant-pH MD runs. In all three solvents, the first pKa of XAO is strongly downshifted compared to the value for the reference compounds (ethylamine and propylamine, respectively); the water and methanol curves have one, and the DMSO curve has two jumps characteristic of remarkable differences in the dissociation constants of acidic groups. The predicted titration curves of Ac-K5-NHMe are in good agreement with the experimental ones; better agreement is achieved with the MD-based method. The titration curves of XAO in methanol and DMSO, calculated using the MD-based approach, trace the shape of the experimental curves, reproducing the pH jump, while those calculated with the EDMC-based approach and the titration curve in water calculated using the MD-based approach have smooth shapes characteristic of the titration of weak multifunctional acids with small differences between the dissociation constants. Nevertheless, quantitative agreement between theoretically predicted and experimental titration curves is not achieved in all three solvents even with the MD-based approach, which is manifested by a smaller pH range of the calculated titration curves with respect to the experimental curves. The poorer agreement obtained for water than for the nonaqueous solvents suggests a significant role of specific solvation in water, which cannot be accounted for by the mean-field solvation models.

Monitoring the DNA Binding Kinetics of a Binuclear Ruthenium Complex by Energy Transfer: Evidence for Slow Shuffling

The semirigid binuclear ruthenium complex Delta,Delta-[mu-(11,11'-bidppz)(phen)(4)Ru(2)](4+) has been shown to rearrange slowly from an initial groove-bound nonluminescent state to a final intercalated emissive state by threading one of its bulky Ru(phen)(2) moieties through the DNA base stack. When this complex binds to poly[d(A-T)(2)], a further increase in emission from the complex is observed after completion of the intercalation, assigned to reorganization of the intercalated complex. We here report a study of the threading process in poly[d(A-T)(2)], in which the minor groove binding dye DAPI is used as an energy transfer probe molecule to assess the distribution of ruthenium complex during and also after the actual threading phase. The emission from DAPI is found to change with the same rate as the emission from the ruthenium complex, and furthermore, DAPI does not disturb the binding kinetics of the latter, justifying it as a good probe of both the threading and the reorganization processes. We conclude from the change in the emission from both DAPI and the ruthenium complex with time that DAPI-ruthenium interactions are most pronounced during the process of threading of the complex, suggesting that the complexes are initially threaded slightly anticooperatively and thereafter redistribute along the DNA to reach their thermodynamically most favorable distribution. The final distribution is characterized by a small but significant binding cooperativity, probably as a result of hydrophobic interactions between the complex ions despite their tetravalent positive charges. The mechanism of "shuffling" the complex along the DNA chain is discussed, i.e., whether the ruthenium complex remains threaded (requiring sequential base-pair openings) or if unthreading followed by lateral diffusion within the ionic atmosphere of the DNA and rethreading occurs.

Interplay of charge distribution and conformation in peptides: Comparison of theory and experiment

We assessed the correlation between charge distribution and conformation of flexible peptides by comparing the theoretically calculated potentiometric-titration curves of two model peptides, Ac-Lys5-NHMe (a model of poly-L-lysine) and Ac-Lys-Ala11-Lys-Gly2-Tyr-NH2 (P1) in water and methanol, with the experimental curves. The calculation procedure consisted of three steps: (i) global conformational search of the peptide under study using the electrostatically driven Monte Carlo (EDMC) method with the empirical conformational energy program for peptides (ECEPP)/3 force field plus the surface-hydration (SRFOPT) or the generalized Born surface area (GBSA) solvation model as well as a molecular dynamics method with the assisted model building and energy refinement (AMBER)99/GBSA force field; (ii) reevaluation of the energy in the pH range considered by using the modified Poisson-Boltzmann approach and taking into account all possible protonation microstates of each conformation, and (iii) calculation of the average degree of protonation of the peptide at a given pH value by Boltzmann averaging over conformations. For Ac-Lys5-NHMe, the computed titration curve agrees qualitatively with the experimental curve of poly-L-lysine in 95% methanol. The experimental titration curves of peptide P1 in water and methanol indicate a remarkable downshift of the first pK(a) value compared to the values for reference compounds (n-butylamine and phenol, respectively), suggesting the presence of a hydrogen bond between the tyrosine hydroxyl oxygen and the H(epsilon) proton of a protonated lysine side chain. The theoretical titration curves agree well with the experimental curves, if conformations with such hydrogen bonds constitute a significant part of the ensemble; otherwise, the theory predicts too small a downward pH shift.

Recognition of duplex RNA by helix-threading peptides

Many important biological processes, from the interferon antiviral response to the generation of microRNA regulators of translation, involve duplex RNA. Small molecules capable of binding duplex RNA structures with high affinity and selectivity will be useful in regulating these processes and, as such, are valuable research tools and potentially therapeutic. In this paper, the synthesis and duplex RNA-binding properties of EDTA.Fe-modified peptide-intercalator conjugates (PICs) are described. Peptide appendages at the 4- and 9-positions of the planar acridine ring system render these PICs threading intercalators, directing the substituents into both grooves of double helical RNA simultaneously. Directed hydroxyl radical cleavage experiments conducted with varying RNA stem-loop structures indicate a preferred binding polarity with the N- and C-termini of the PIC in the minor and major grooves, respectively. However, this binding polarity is shown to be dependent on both the structure of the PIC and the RNA secondary structure adjacent to the intercalation site. Definition of the minimal RNA structure required for binding to one of these PICs led to the identification of an intercalation site in a pre-microRNA from Caenorhabditis elegans. Results presented will guide both rational design and combinatorial approaches for the generation of new RNA binding PICs and will continue to facilitate the identification of naturally occurring RNA targets for these small molecules.

Implementation and testing of stable, fast implicit solvation in molecular dynamics using the smooth-permittivity finite difference Poisson-Boltzmann method

A fast stable finite difference Poisson-Boltzmann (FDPB) model for implicit solvation in molecular dynamics simulations was developed using the smooth permittivity FDPB method implemented in the OpenEye ZAP libraries. This was interfaced with two widely used molecular dynamics packages, AMBER and CHARMM. Using the CHARMM-ZAP software combination, the implicit solvent model was tested on eight proteins differing in size, structure, and cofactors: calmodulin, horseradish peroxidase (with and without substrate analogue bound), lipid carrier protein, flavodoxin, ubiquitin, cytochrome c, and a de novo designed 3-helix bundle. The stability and accuracy of the implicit solvent simulations was assessed by examining root-mean-squared deviations from crystal structure. This measure was compared with that of a standard explicit water solvent model. In addition we compared experimental and calculated NMR order parameters to obtain a residue level assessment of the accuracy of MD-ZAP for simulating dynamic quantities. Overall, the agreement of the implicit solvent model with experiment was as good as that of explicit water simulations. The implicit solvent method was up to eight times faster than the explicit water simulations, and approximately four times slower than a vacuum simulation (i.e., with no solvent treatment).

Peptide quinoline conjugates: a new class of RNA-binding molecules

[reaction: see text] A synthesis of 4,8-disubstituted 2-phenylquinoline amino acids is reported with the incorporation of one example into a peptide by solid-phase synthesis. The phenylquinoline-containing peptide binds an RNA target with nanomolar affinity (K(D) = 208 nM). The strategy can be used to prepare a variety of 2-substituted quinoline amino acids for alteration of affinity in intercalator peptides. Since quinolones represent an important class of antibacterials, these compounds may be useful in the discovery of new antibacterial agents.

Selective intercalation of charge neutral intercalators into GG and CG steps: implication of HOMO-LUMO interaction for sequence-selective drug intercalation into DNA

We have synthesized naphthopyranone epoxide 4 from D-isoascorbic acid together with its three diastereoisomers. DNA alkylation of ODNs containing 5'XGT3' and 5'TGY3' by 4 (11R, 13R), where X and Y are any nucleotide bases, occurred at all G residues except at G of the 5'TGC3' sequence. In contrast, the three other diastereoisomers of 4 showed only weak G alkylation activity. Differential (1)H NMR NOE of the 4-G adduct confirmed the G-N7 alkylation at the epoxide carbon of 4 with concomitant S(N)2 ring opening of the epoxide. Quantitative HPLC analysis of G alkylation efficiency for 4 showed the order of G alkylation susceptibility as TGGT approximately CGT >> TGA > AGT > TGT >> TGC. The order was fully consistent with those reported for aflatoxin B(1) oxide and kapurimycin A(3), suggesting that the sequence selectivity observed for these DNA alkylating agents is not structure dependent but most likely due to the intrinsic property of DNA sequences. We found that the order of G alkylation susceptibility obtained for 4 completely matched the calculated HOMO energy level of G-containing sequences. These results underscore that 4 is a unique molecular probe for ranking the HOMO level of G-containing sequences by well-known G alkylation chemistry and suggests that the intercalation of charge neutral intercalators is a HOMO-controlled process.

Engineering the pH-optimum of a triglyceride lipase: from predictions based on electrostatic computations to experimental results

The optimisation of enzymes for particular purposes or conditions remains an important target in virtually all protein engineering endeavours. Here, we present a successful strategy for altering the pH-optimum of the triglyceride lipase cutinase from Fusarium solani pisi. The computed electrostatic pH-dependent potentials in the active site environment are correlated with the experimentally observed enzymatic activities. At pH-optimum a distinct negative potential is present in all the lipases and esterases that we studied so far. This has prompted us to propose the "The Electrostatic Catapult Model" as a model for product release after cleavage of the ester bond. The origin of the negative potential is associated with the titration status of specific residues in the vicinity of the active site cleft. In the case of cutinase, the role of Glu44 was systematically investigated by mutations into Ala and Lys. Also, the neighbouring Thr45 was mutated into Proline, with the aim of shifting the spatial location of Glu44. All the charge mutants displayed altered titration behaviour of active site electrostatic potentials. Typically, the substitution of the residue Glu44 pushes the onset of the active site negative potential towards more alkaline conditions. We, therefore, predicted more alkaline pH optima, and this was indeed the experimentally observed. Finally, it was found that the pH-dependent computed Coulombic energy displayed a strong correlation with the observed melting temperatures of native cutinase.

Enantioselective DNA threading dynamics by phenazine-linked

The interactions between the stereoisomers of the chiral bis-intercalator [mu-C4(cpdppz)(2)-(phen)(4)Ru(2)](4+) and DNA reveal interesting dynamic discrimination properties. The two enantiomers Delta-Delta and Lambda-Lambda both form very strong complexes with calf thymus DNA with similar thermodynamic affinities. By contrast, they display considerable variations in their binding kinetics. The Delta-Delta enantiomer has higher affinity for calf thymus DNA than for [poly(dA-dT)](2), and the association kinetics of the dimer to DNA, as well as to polynucleotides, requires a multiexponential fitting function. The dissociation reaction, on the other hand, could be described by a single exponential for [poly(dA-dT)](2), whereas two exponentials were required for mixed-sequence DNA. To understand the key mechanistic steps of the reaction, the kinetics was studied at varied salt concentration for different choices of DNA and chirality of the threading complex. The enantiomers were found to have markedly different dissociation rates, the Lambda-Lambda enantiomer dissociating about an order of magnitude faster than the Delta-Delta enantiomer. Also, the salt dependence of the dissociation rate constants differed between the enantiomers, being stronger for the Lambda-Lambda enantiomer than for the Delta-Delta enantiomer. Since the dissociation reaction requires unthreading of bulky parts of the bis-intercalator through the DNA helix, a considerable conformational change of the DNA must be involved, possibly defining the rate-limiting step.

Alternate modes of binding in two crystal structures of alkaline phosphatase-inhibitor complexes

Two high resolution crystal structures of Escherichia coli alkaline phosphatase (AP) in the presence of phosphonate inhibitors are reported. The phosphonate compounds, phosphonoacetic acid (PAA) and mercaptomethylphosphonic acid (MMP), bind competitively to AP with dissociation constants of 5.5 and 0.6 mM, respectively. The structures of the complexes of AP with PAA and MMP were refined at high resolution to crystallographic R-values of 19.0 and 17.5%, respectively. Refinement of the AP-inhibitor complexes was carried out using X-PLOR. The final round of refinement was done using SHELXL-97. Crystallographic analyses of the inhibitor complexes reveal different binding modes for the two phosphonate compounds. The significant difference in binding constants can be attributed to these alternative binding modes observed in the high resolution X-ray structures. The phosphinyl group of PAA coordinates to the active site zinc ions in a manner similar to the competitive inhibitor and product inorganic phosphate. In contrast, MMP binds with its phosphonate moiety directed toward solvent. Both enzyme-inhibitor complexes exhibit close contacts, one of which has the chemical and geometrical potential to be considered an unconventional hydrogen bond of the type C-H...X.

Comparison of anionic and cationic trypsinogens: the anionic activation domain is more flexible in solution and differs in its mode of BPTI binding in the crystal structure

Unlike bovine cationic trypsin, rat anionic trypsin retains activity at high pH. This alkaline stability has been attributed to stabilization of the salt bridge between the N-terminal Ile16 and Asp194 by the surface negative charge (Soman K, Yang A-S, Honig B, Fletterick R., 1989, Biochemistry 28:9918-9926). The formation of this salt bridge controls the conformation of the activation domain in trypsin. In this work we probe the structure of rat trypsinogen to determine the effects of the surface negative charge on the activation domain in the absence of the Ile16-Asp194 salt bridge. We determined the crystal structures of the rat trypsin-BPTI complex and the rat trypsinogen-BPTI complex at 1.8 and 2.2 A, respectively. The BPTI complex of rat trypsinogen resembles that of rat trypsin. Surprisingly, the side chain of Ile16 is found in a similar position in both the rat trypsin and trypsinogen complexes, although it is not the N-terminal residue and cannot form the salt bridge in trypsinogen. The resulting position of the activation peptide alters the conformation of the adjacent autolysis loop (residues 142-153). While bovine trypsinogen and trypsin have similar CD spectra, the CD spectrum of rat trypsinogen has only 60% of the intensity of rat trypsin. This lower intensity most likely results from increased flexibility around two conserved tryptophans, which are adjacent to the activation domain. The NMR spectrum of rat trypsinogen contains high field methyl signals as observed in bovine trypsinogen. It is concluded that the activation domain of rat trypsinogen is more flexible than that of bovine trypsinogen, but does not extend further into the protein core.

Enzyme-inhibitor association thermodynamics: explicit and continuum solvent studies

Studying the thermodynamics of biochemical association reactions at the microscopic level requires efficient sampling of the configurations of the reactants and solvent as a function of the reaction pathways. In most cases, the associating ligand and receptor have complementary interlocking shapes. Upon association, loosely connected or disconnected solvent cavities at and around the binding site are formed. Disconnected solvent regions lead to severe statistical sampling problems when simulations are performed with explicit solvent. It was recently proposed that, when such limitations are encountered, they might be overcome by the use of the grand canonical ensemble. Here we investigate one such case and report the association free energy profile (potential of mean force) between trypsin and benzamidine along a chosen reaction coordinate as calculated using the grand canonical Monte Carlo method. The free energy profile is also calculated for a continuum solvent model using the Poisson equation, and the results are compared to the explicit water simulations. The comparison shows that the continuum solvent approach is surprisingly successful in reproducing the explicit solvent simulation results. The Monte Carlo results are analyzed in detail with respect to solvation structure. In the binding site channel there are waters bridging the carbonyl oxygen groups of Asp189 with the NH2 groups of benzamidine, which are displaced upon inhibitor binding. A similar solvent-bridging configuration has been seen in the crystal structure of trypsin complexed with bovine pancreatic trypsin inhibitor. The predicted locations of other internal waters are in very good agreement with the positions found in the crystal structures, which supports the accuracy of the simulations.

Demonstration that 1-trans-epoxysuccinyl-L-leucylamido-(4-guanidino) butane (E-64) is one of the most effective low Mr inhibitors of trypsin-catalysed hydrolysis. Characterization by kinetic analysis and by energy minimization and molecular dynamics simulation of the E-64-beta-trypsin complex

1-trans-Epoxysuccinyl-L-leucylamido(4-guanidino)butane (E-64) was shown to inhibit beta-trypsin by a reversible competitive mechanism; this contrasts with the widely held view that E-64 is a class-specific inhibitor of the cysteine proteinases and reports in the literature that it does not inhibit a number of other enzymes including, notably, trypsin. The K1, value (3 x 10(-5) M) determined by kinetic analysis of the hydrolysis of N alpha-benzoyl-L-arginine 4-nitroanilide in Tris/HCl buffer, pH 7.4, at 25 degrees C, I = 0.1, catalysed by beta-trypsin is comparable with those for the inhibition of trypsin by benzamidine and 4-aminobenzamidine, which are widely regarded as the most effective low Mr inhibitors of this enzyme. Computer modelling of the beta-trypsin-E64 adsorptive complex, by energy minimization, molecular dynamics simulation and Poisson-Boltzmann electrostatic-potential calculations, was used to define the probable binding mode of E-64; the ligand lies parallel to the active-centre cleft, anchored principally by the dominant electrostatic interaction of the guanidinium cation at one end of the E-64 molecule with the carboxylate anion of Asp-171 (beta-trypsin numbering from Ile-1) in the S1-subsite, and by the interaction of the carboxylate substituent on C-2 of the epoxide ring at the other end of the molecule with Lys-43; the epoxide ring of E-64 is remote from the catalytic site serine hydroxy group. The possibility that E-64 might bind to the cysteine proteinases clostripain (from Clostridium histolyticum) and alpha-gingivain (one of the extracellular enzymes from phyromonas gingivalis) in a manner analogous to that deduced for the beta-trypsin-E-64 complex is discussed.

The structural puzzle of how serpin serine proteinase inhibitors work

Serine proteinase cleavage of proteins is essential to a wide variety of biological processes and is primarily regulated by protein inhibitors. Many inhibitors are conformationally rigid simulations of optimal serine proteinase substrates, which makes them highly efficient competitive inhibitors of target proteinases. In contrast, members of the serpin family of serine proteinase inhibitors display extensive flexibility and polymorphism, particularly in their reactive site segments and in beta-sheet secondary structure, which can take up and expel strands. Reactive site and beta-sheet polymorphism appear to be coupled in the serpins and may account for the extreme stability of serpin-proteinase complexes through the insertion of the reactive site strand into a beta-sheet. These unusual properties may have opened an adaptive pathway of proteinase regulation that was unavailable to the conformationally rigid proteinase inhibitors.

Classical electrostatics in biology and chemistry

A major revival in the use of classical electrostatics as an approach to the study of charged and polar molecules in aqueous solution has been made possible through the development of fast numerical and computational methods to solve the Poisson-Boltzmann equation for solute molecules that have complex shapes and charge distributions. Graphical visualization of the calculated electrostatic potentials generated by proteins and nucleic acids has revealed insights into the role of electrostatic interactions in a wide range of biological phenomena. Classical electrostatics has also proved to be successful quantitative tool yielding accurate descriptions of electrical potentials, diffusion limited processes, pH-dependent properties of proteins, ionic strength-dependent phenomena, and the solvation free energies of organic molecules.

A structural model for the prostate disease marker, human prostate-specific antigen

Prostate-specific antigen (PSA) provides an excellent serum marker for prostate cancer, the most frequent form of cancer in American males. PSA is a 237-residue protease based on sequence homology to kallikrein-like enzymes. To predict the 3-dimensional structure of PSA, homology modeling studies were performed based on sequence and structural alignments with tonin, pancreatic kallikrein, chymotrypsin, and trypsin. The structurally conserved regions of the 4 reference X-ray proteins provided the core structure of PSA, whereas the loop structures were modeled on the loops of tonin and kallikrein. The unique "kallikrein loop" insert, between Ser 95b and Pro 95k of kallikrein, was constructed using molecular mechanics, dynamics, and electrostatics calculations. In the resulting PSA structure, the catalytic triad, involving residues His 57, Asp 102, and Ser 195, and hydrophobic and electrostatic interactions typical of serine proteases were extremely well conserved. Similarly, the 5-disulfide bonds of kallikrein were also conserved in PSA. These results, together with the fact that no major steric clashes arose during the modeling process, provide strong evidence for the validity of the PSA model. Calculation of the electrostatic potential contours of kallikrein and PSA was carried out using the finite difference Poisson-Boltzmann method. The calculations revealed matching areas of negative potential near the catalytic triad, but differences in the positive potential surrounding the active site. The PSA glycosylation site, Asn 61, is fully accessible to the solvent and is enclosed in a positive region of the isopotential map. The bottom of the substrate specificity pocket, residue S1, is a serine (Ser 189) as in chymotrypsin, rather than aspartate (Asp 189) as in tonin, kallikrein, and trypsin. This fact, plus other features of the S1 binding-pocket region, suggest that PSA would prefer substrates with hydrophobic residues at the P1 position. The location of a potential zinc ion binding site involving the side chain of histidines 91, 101, and 233 is also suggested. This PSA model should facilitate the understanding and prediction of structural and functional properties of this important cancer marker.

Identification and analysis of a template-primer (ds-DNA) binding cleft in E. coli DNA polymerase I: an electrostatic potential contour pattern of the modeled structure

In the modeled structure of the Klenow fragment of E. coli DNA polymerase I, we have identified a distinct region that exhibits a strong electropositive potential contour. The examination of the distribution of the electropositive and negative potential across the two-dimensional slices of the modeled structure revealed that the positive potential was concentrated around the cleft. The approximate size and shape of the region appears well suited to accommodate eight base pairs of duplex DNA and is consistent with the position of the dsDNA binding cleft reported in the crystal structure [Beese et al., Science (1993) 260, 352-355].

The blind watchmaker and rational protein engineering

In the present review some scientific areas of key importance for protein engineering are discussed, such as problems involved in deducting protein sequence from DNA sequence (due to posttranscriptional editing, splicing and posttranslational modifications), modelling of protein structures by homology, NMR of large proteins (including probing the molecular surface with relaxation agents), simulation of protein structures by molecular dynamics and simulation of electrostatic effects in proteins (including pH-dependent effects). It is argued that all of these areas could be of key importance in most protein engineering projects, because they give access to increased and often unique information. In the last part of the review some potential areas for future applications of protein engineering approaches are discussed, such as non-conventional media, de novo design and nanotechnology.

Electrostatic effects in trypsin reactions. Influence of salts

The influence of inorganic salts on trypsin-catalyzed reactions has been studied. It is shown that: (a) monovalent cations are reversible competitive inhibitors of tryptic hydrolysis of cationic substrates, whereas their binding has no effect on the reaction of neutral substrates; (b) a nonelectrostatic salt effect on the binding of both cationic and non-ionic substrates is caused by changes in the thermodynamic activity coefficient of the substrate; (c) the rate of trypsin active-site acylation is not affected by inorganic salts with monovalent cations. The data suggest that low-molecular-mass substrates are extracted into the enzyme microphase during substrate binding and further chemical transformations proceed without an access from surrounding medium. It is proposed that formation of a properly oriented dipole in the trypsin binding pocket by the cationic group of the substrate and Asp189 carboxyl is responsible for the elevated acylation rate of trypsin active site by substrates containing lysine and arginine. Introduction of additional negative charges into the enzyme molecule by chemical modification of lysyl residues by pyromellitic anhydride increased the specificity of trypsin towards cationic substrates and inhibitors. Lysine residues are therefore considered as suitable targets for site-directed mutagenesis aimed at the improvement of selectivity and catalytic properties of trypsin.

Contribution of long-range electrostatic interactions to the stabilization of the catalytic transition state of the serine protease subtilisin BPN'

The possible role of long-range electrostatic interactions on the catalytic activity of the serine protease subtilisin BPN' is investigated using protein engineering techniques. Charged residues on the surface of the enzyme some 13-15 A from the active site were mutated to either neutral or oppositely charged residues. The effect of these mutations on the stability of a complex formed between subtilisin BPN' and Z-Ala-Ala-Pro-Phe-trifluoromethyl ketone, a transition-state inhibitor of the enzyme, was measured. The values of Ki for the complex between the trifluoromethyl ketone and wild-type and mutant subtilisins were used to study the possible contribution of long-range electrostatics in stabilizing the charge distribution in the complex and thus, by analogy, on the transition state of hydrolysis for subtilisin BPN'. Measurement of kon, koff, and Ki for the inhibition of wild-type and mutant subtilisins showed that charged mutations distant from the active site can affect koff and Ki but have little effect on kon. The experimental results show that there is a small, 0.10-0.46 kcal mol-1, but significant contribution to the binding energy from distant surface charges, at low ionic strength. The experimental results were compared to theoretical results, calculated using the DelPhi program for different charge distributions in the complex. The experimental results were found to be most consistent with a complex in which an ion pair is formed between the protonated active site histidine and the ionized oxyanion. Both experimental and theoretical results suggest that long-range electrostatic interactions do play a role in stabilizing the transition-state complex formed between enzyme and inhibitor.(ABSTRACT TRUNCATED AT 250 WORDS)

Structural prototypes for an extended family of flavoprotein reductases: comparison of phthalate dioxygenase reductase with ferredoxin reductase and ferredoxin

The structure of phthalate dioxygenase reductase (PDR), a monomeric iron-sulfur flavoprotein that delivers electrons from NADH to phthalate dioxygenase, is compared to ferredoxin-NADP+ reductase (FNR) and ferredoxin, the proteins that reduce NADP+ in the final reaction of photosystem I. The folding patterns of the domains that bind flavin, NAD(P), and [2Fe-2S] are very similar in the two systems. Alignment of the X-ray structures of PDR and FNR substantiates the assignment of features that characterize a family of flavoprotein reductases whose members include cytochrome P-450 reductase, sulfite and nitrate reductases, and nitric oxide synthase. Hallmarks of this subfamily of flavoproteins, here termed the FNR family, are an antiparallel beta-barrel that binds the flavin prosthetic group, and a characteristic variant of the classic pyridine nucleotide-binding fold. Despite the similarities between FNR and PDR, attempts to model the structure of a dissociable FNR:ferredoxin complex by analogy with PDR reveal features that are at odds with chemical crosslinking studies (Zanetti, G., Morelli, D., Ronchi, S., Negri, A., Aliverti, A., & Curti, B., 1988, Biochemistry 27, 3753-3759). Differences in the binding sites for flavin and pyridine nucleotides determine the nucleotide specificities of FNR and PDR. The specificity of FNR for NADP+ arises primarily from substitutions in FNR that favor interactions with the 2' phosphate of NADP+. Variations in the conformation and sequences of the loop adjoining the flavin phosphate affect the selectivity for FAD versus FMN. The midpoint potentials for reduction of the flavin and [2Fe-2S] groups in PDR are higher than their counterparts in FNR and spinach ferredoxin, by about 120 mV and 260 mV, respectively. Comparisons of the structure of PDR with spinach FNR and with ferredoxin from Anabaena 7120, along with calculations of electrostatic potentials, suggest that local interactions, including hydrogen bonds, are the dominant contributors to these differences in potential.

On the calculation of pKas in proteins

This paper describes a general method to calculate the pKas of ionizable groups in proteins. Electrostatic calculations are carried out using the finite difference Poisson-Boltzmann (FDPB) method. A formal treatment of the calculation of pKas within the framework of the FDPB method is presented. The major change with respect to previous work is the specific incorporation of the complete charge distribution of both the neutral and charged forms of each ionizable group into the formalism. This is extremely important for the treatment of salt bridges. A hybrid statistical mechanical/Tanford-Roxby method, which is found to be significantly faster than previous treatments, is also introduced. This simplifies the problem of summing over the large number of possible ionization states for a complex polyion. Applications to BPTI and serine proteases suggest that the calculations can be quite reliable. However, the necessity of including bound waters in the treatment of the Asp-70... His-31 salt bridge in T4 lysozyme and experience with other proteins suggest that additional factors ultimately need to be considered in a comprehensive treatment of pKas in proteins.

The crystal structure of the Bacillus lentus alkaline protease, subtilisin BL, at 1.4 A resolution

The crystal structure of subtilisin BL, an alkaline protease from Bacillus lentus with activity at pH 11, has been determined to 1.4 A resolution. The structure was solved by molecular replacement starting with the 2.1 A structure of subtilisin BPN' followed by molecular dynamics refinement using X-PLOR. A final crystallographic R-factor of 19% overall was obtained. The enzyme possesses stability at high pH, which is a result of the high pI of the protein. Almost all of the acidic side-chains are involved in some type of electrostatic interaction (ion pairs, calcium binding, etc.). Furthermore, three of seven tyrosine residues have potential partners for forming salt bridges. All of the potential partners are arginine with a pK around 12. Lysine would not function well in a salt bridge with tyrosine as it deprotonates at around the same pH as tyrosine ionizes. Stability at high pH is acquired in part from the pI of the protein, but also from the formation of salt bridges (which would affect the pI). The overall structure of the enzyme is very similar to other subtilisins and shows that the subtilisin fold is more highly conserved than would be expected from the differences in amino acid sequence. The amino acid side-chains in the hydrophobic core are not conserved, though the inter-residue interactions are. Finally, one third of the serine side-chains in the protein have multiple conformations. This presents an opportunity to correlate computer simulations with observed occupancies in the crystal structure.

Electrostatic forces in two lysozymes: calculations and measurements of histidine pKa values

In order to examine the electrostatic forces in globular proteins, pKa values and their ionic strength dependence of His residues of hen egg white lysozyme (HEWL) and human lysozyme (HUML) were measured, and they were compared with those calculated numerically. pKa values of His residues in HEWL, HUML, and short oligopeptides were determined from chemical shift changes of His side chains by 1H-nmr measurements. The associated changes in pKa values in HEWL and HUML were calculated by solving the Poisson-Boltzmann equations numerically for macroscopic dielectric models. The calculated pKa changes and their ionic strength dependence agreed fairly well with the observed ones. The contribution from each residue of each alpha-helix dipole to the pKa values and their ionic strength dependence was analyzed using Green's reciprocity theorem. The results indicate that (1) the pKa of His residues are largely affected by surrounding ionized and polar groups; (2) the ionic strength dependence of the pKa values is determined by the overall charge distributions and their accessibilities to solvent; and (3) alpha-helix dipoles make a significant contribution to the pKa, when the His residue is close to the helix terminus and not fully exposed to the solvent.

Electrostatic control of midpoint potentials in the cytochrome subunit of the Rhodopseudomonas viridis reaction center

The photosynthetic reaction center of Rhodopseudomonas viridis has four hemes in a linear arrangement with alternating high- and low-potential sites. Their midpoints are -60, 20, 310, and 380 mV [Dracheva, S. M., Drachev, L. A., Konstantinov, A. A., Semenov, A. Y., Skulachev, V. P., Arutjunjan, A. M., Shuvalov, V. A. & Zaberezhnaya, S. M. (1988) Eur. J. Biochem. 171, 253-264]. Electrostatic calculations reproduce the 440-mV midpoint spread and assignments of high- and low-potential hemes. When calculations on model compounds to connect the theoretical midpoints to the standard hydrogen electrode are used, the absolute electrochemical midpoints for the reaction center hemes are also in good agreement with experiment. The free energy of oxidation is found to be dependent on pairwise interactions with charged amino acids, heme propionic acids, previously oxidized hemes, and axial ligands and on the reaction field induced by heme oxidation.

The electrostatic potential of Escherichia coli dihydrofolate reductase

Escherichia coli dihydrofolate reductase (DHFR) carries a net charge of -10 electrons yet it binds ligands with net charges of -4 (NADPH) and -2 (folate or dihydrofolate). Evaluation and analysis of the electrostatic potential of the enzyme give insight as to how this is accomplished. The results show that the enzyme is covered by an overall negative potential (as expected) except for the ligand binding sites, which are located inside "pockets" of positive potential that enable the enzyme to bind the negatively charged ligands. The electrostatic potential can be related to the asymmetric distribution of charged residues in the enzyme. The asymmetric charge distribution, along with the dielectric boundary that occurs at the solvent-protein interface, is analogous to the situation occurring in superoxide dismutase. Thus DHFR is another case where the shape of the active site focuses electric fields out into solution. The positive electrostatic potential at the entrance of the ligand binding site in E. coli DHFR is shown to be a direct consequence of the presence of three positively charged residues at positions 32, 52, and 57--residues which have also been shown recently to contribute significantly to electronic polarization of the ligand folate. The latter has been postulated to be involved in the catalytic process. A similar structural motif of three positively charged amino acids that gives rise to a positive potential at the entrance to the active site is also found in DHFR from chicken liver, and is suggested to be a common feature in DHFRs from many species. It is noted that, although the net charges of DHFRs from different species vary from +3 to -10, the enzymes are able to bind the same negatively charged ligands, and perform the same catalytic function.

Linkage of thioredoxin stability to titration of ionizable groups with perturbed pKa

The highly conserved, buried, Asp 26 in Escherichia coli thioredoxin has a pKa = 7.5, and its titration is associated with a sizable destabilization of the protein [Langsetmo, K., Fuchs, J., & Woodward, C. (1991) Biochemistry (preceding paper in this issue)]. A fit of the experimental pH dependence of thioredoxin stability to a theoretical expression for the pH/stability relation in proteins agrees closely with a pKa value of 7.5 for Asp 26. The agreement between the experimental and theoretical changes in protein stability due to substitution of Asp 26 by alanine is also good. The local structure in the vicinity of Asp 26 in the low-pH crystal structure (with uncharged Asp 26) is hydrophobic, indicating that the aspartate would be highly destabilized. In theoretical calculations, the desolvation penalty for deprotonating Asp 26 in this environment is similar to the total protein folding energy. As a consequence, the Asp 26 pKa would be much greater than 7.5, and/or the protein might not fold. This suggests that a compensating process partially stabilizes the Asp 26 carboxyl group when it is charged. A simple model for this proposed, whereby the Lys 57 side chain rotates to form a salt bridge with Asp 26 when it is deprotonated.

Crystal structure of rat trypsin-S195C at -150 degrees C. Analysis of low activity of recombinant and semisynthetic thiol proteases

The X-ray crystal structure of trypsin-S195C, a rat anionic trypsin mutant in which the active site serine has been replaced by cysteine, was determined at -150 degrees C and room temperature to 1.6 A resolution, R = 15.4% and 1.8 A resolution, R = 15.0%, respectively. Cryo-crystallography was employed to improve the quality of the diffraction data and the resulting structure by eliminating radiation damage and decreasing atomic thermal motion. The average temperature factor decreased by 10 A2 relative to that of the room temperature structure. No radiation-induced decay of the data was detected. The side-chains of the catalytic cysteine and histidine of trypsin-S195C are found with 25% occupancy in secondary orientations rotated 104 degrees and 90 degrees out of the active site, respectively. These alterations, as well as more subtle changes in the active site may be caused by the oxidation of the catalytic sulfur to sulfenic acid. The position of the carbonyl carbon of the tetrahedral intermediate analog, p-amidinophenylpyruvic acid, modeled into trypsin-S195C, is 1.1 A from the catalytic sulfur. The large size and altered approach of the catalytic sulfur to substrates could account for the observed low catalytic activity relative to wild-type trypsin. In addition to the benzamidine in the specificity pocket, two additional binding sites for benzamidine are characterized. One of these mediates an intermolecular contact that appears to maintain the crystal lattice.

Mutational remodeling of enzyme specificity

With the advent of genetic engineering techniques has come the ability to modify proteins as desired. Given this stunning capability, the question remains what residues should be altered, and how should they be changed to achieve a particular specificity pattern. The goals of such modifications are likely to fall into either of two categories: probing the function of a protein or attempting to alter its properties. In either case, our understanding of the consequences of a mutation, as ascertained by our ability to predict the results, is currently quite limited. The problem is extraordinarily complex; our understanding of how to calculate the energetics involved is still incomplete, and we are just beginning to accumulate experimental data which may help guide us. On the positive side, theoretical methods are now being developed and refined that should prove useful in the drive to engineer enzyme specificity. What may be most important at this juncture is to expand the experimental database interrelating sequence, function, and structure. That is, there should be a concerted effort to combine functional analysis of mutant proteins with structural analysis. Only from this combined examination of the effects of mutations can sufficient data be accumulated to test and improve both qualitative and quantitative approaches or methods for remodeling enzyme specificity.