Engineering trypsin for inhibitor resistance

The development of effective protease therapeutics requires that the proteases be more resistant to naturally occurring inhibitors while maintaining catalytic activity. A key step in developing inhibitor resistance is the identification of key residues in protease-inhibitor interaction. Given that majority of the protease therapeutics currently in use are trypsin-fold, trypsin itself serves as an ideal model for studying protease-inhibitor interaction. To test the importance of several trypsin-inhibitor interactions on the prime-side binding interface, we created four trypsin single variants Y39A, Y39F, K60A, and K60V and report biochemical sensitivity against bovine pancreatic trypsin inhibitor (BPTI) and M84R ecotin. All variants retained catalytic activity against small, commercially available peptide substrates [k_cat/K_M = (1.2 ± 0.3) × 10⁷ M⁻¹ s⁻¹. Compared with wild-type, the K60A and K60V variants showed increased sensitivity to BPTI but less sensitivity to ecotin. The Y39A variant was less sensitive to BPTI and ecotin while the Y39F variant was more sensitive to both. The relative binding free energies between BPTI complexes with WT, Y39F, and Y39A were calculated based on 3.5 µs combined explicit solvent molecular dynamics simulations. The BPTI:Y39F complex resulted in the lowest binding energy, while BPTI:Y39A resulted in the highest. Simulations of Y39F revealed increased conformational rearrangement of F39, which allowed formation of a new hydrogen bond between BPTI R17 and H40 of the variant. All together, these data suggest that positions 39 and 60 are key for inhibitor binding to trypsin, and likely more trypsin-fold proteases.     

Native-like interactions favored in the unfolded bovine pancreatic trypsin inhibitor have different roles in folding

Folding kinetics of a series of bovine pancreatic trypsin inhibitor (BPTI) variants with similar stabilities and structures have been measured. All are strongly destabilized relative to WT. In Y21A, F22A, Y23A, G37A, and F45A, the three native disulfide bonds are retained. In RM(14-38), Cys14 and Cys38 thiols are methylated while C30-C51 and C5-C55 disulfides remain intact. At pH 2 and 20 degrees C, relaxation rate constants of the major kinetic phase range from approximately 10 ms to 0.71 s in the absence of denaturant. All mutants except G37A exhibit standard two-state behavior. Y21A, F22A, and Y23A fold much more slowly than other mutants. The experiments were designed to test the hypothesis that native-like structure detected in the unfolded BPTI is important in folding. Two native-like contacts are implied by NOEs in reduced and unfolded BPTI, between residues Tyr23 and Ala25, and between Gly37 NH and the Tyr35 ring. The results support an earlier hypothesis that formation of the central beta-hairpin, monitored by a local native interaction between Tyr23 and Ala25, is crucial to initiation of BPTI folding. The second native-like contact is important, not in folding initiation, but in preventing a kinetic trap later in the process. Evidence for this comes from mutant G37A, which behaves very differently from the others in displaying a phenomenon called rollover. G37A is, to our knowledge, the first reported case in which a single-site replacement causes rollover, while the wild type and all other known mutants of the same protein show typical two-state chevron plots. The best explanation is that the G37A mutation introduces a kinetic trap of the type described by Chan and Dill [(1998) Proteins 30, 2-33]. In native BPTI, there is an unusual polar interaction between the ring of Tyr35 and the backbone NH of Gly37. Our results suggest that the NH-aromatic interaction between residues 37 and 35 is important throughout folding in stabilizing native-like loop conformations and in preventing the flexible loops from being trapped in nonfunctional conformations during later stages of folding.     

A highly destabilizing mutation, G37A, of the bovine pancreatic trypsin inhibitor retains the average native conformation but greatly increases local flexibility

A point mutation, G37A, on the surface of bovine pancreatic trypsin inhibitor (BPTI) destabilizes the protein by approximately 5 kcal/mol, which is very high for addition of one methyl group. In wild-type (WT) BPTI, Gly 37 HN is in an unusual NH-aromatic-NH network of interactions with the ring of Tyr 35 and the side chain HN of Asn 44. G37A was designed to disrupt this interaction, since the phi and psi backbone angles of G37 are not favorable for an amino acid containing a beta-carbon. Investigations of the structure and dynamics by NMR methods show that G37A retains the average WT structure. The NH-aromatic-NH interactions remain intact, as indicated by NOEs and the large upfield ring current shift (approximately 4 ppm) of A37 HN. The NMR structure, confirmed by molecular modeling calculations, requires phi and psi backbone angles that are highly destabilizing when alanine is in position 37. Although the average structure is essentially unchanged, the dynamics are altered dramatically. Many residues in the region of the mutation have increased flexibility, as probed by aromatic ring flip rates and native state hydrogen exchange. We conclude that a large fraction of the destabilization arises from maintaining A37 in a high-energy conformation. This suggests that disruption of the NH-aromatic-NH network is energetically very costly, and may involve other cooperatively linked interactions. The results illustrate the importance of the Gly-Gly sequence at positions 36 and 37 and the 37 HN-35 aromatic interaction to the stability, folding, and dynamics of the BPTI.     

The response of T4 lysozyme to large-to-small substitutions within the core and its relation to the hydrophobic effect.

To further examine the structural and thermodynamic basis of hydrophobic stabilization in proteins, all of the bulky non-polar residues that are buried or largely buried within the core of T4 lysozyme were substituted with alanine. In 25 cases, including eight reported previously, it was possible to determine the crystal structures of the variants. The structures of four variants with double substitutions were also determined. In the majority of cases the "large-to-small" substitutions lead to internal cavities. In other cases declivities or channels open to the surface were formed. In some cases the structural changes were minimal (mainchain shifts < or = 0.3 A); in other cases mainchain atoms moved up to 2 A. In the case of Ile 29 --> Ala the structure collapsed to such a degree that the volume of the putative cavity was zero. Crystallographic analysis suggests that the occupancy of the engineered cavities by solvent is usually low. The mutants Val 149 --> Ala (V149A) and Met 6 --> Ala (M6A), however, are exceptions and have, respectively, one and two well-ordered water molecules within the cavity. The Val 149 --> Ala substitution allows the solvent molecule to hydrogen bond to polar atoms that are occluded in the wild-type molecule. Similarly, the replacement of Met 6 with alanine allows the two solvent molecules to hydrogen bond to each other and to polar atoms on the protein. Except for Val 149 --> Ala the loss of stability of all the cavity mutants can be rationalized as a combination of two terms. The first is a constant for a given class of substitution (e.g., -2.1 kcal/mol for all Leu --> Ala substitutions) and can be considered as the difference between the free energy of transfer of leucine and alanine from solvent to the core of the protein. The second term can be considered as the energy cost of forming the cavity and is consistent with a numerical value of 22 cal mol(-1) A(-3). Physically, this term is due to the loss of van der Waal''s interactions between the bulky sidechain that is removed and the atoms that form the wall of the cavity. The overall results are consistent with the prior rationalization of Leu --> Ala mutants in T4 lysozyme by Eriksson et al. (Eriksson et al., 1992, Science 255:178-183).     

Altered conductance and permeability of Cx40 mutations associated with atrial fibrillation

Gap junctions ensure the rapid propagation of the action potential throughout the myocardium. Three mutant forms of connexin40 (Cx40; A96S, M163V, and G38D), the primary component of the atrial gap junction channel, are associated with atrial fibrillation and retain the ability to form functional channels. We determined the biophysical properties of these mutant gap junctions in transiently transfected HeLa and N2A cells. All three mutants showed macroscopic junctional conductances over the range of 0.5 to 40 nS, and voltage dependences comparable to those of wild-type (WT) Cx40. However, the unitary conductance of G38D channels was ∼1.6-fold higher than that of WT Cx40 channels (∼220 vs. ∼135 pS), whereas the unitary conductances of the A96S and M163V mutants were similar to that of WT Cx40. Furthermore, the M163V and G38D channels exhibited approximately two- and approximately fivefold higher permeability to the anionic dye Lucifer yellow (LY) relative to K⁺ (LY/K⁺) compared with that of WT Cx40, whereas A96S LY transfer was similar to that of WT (G38D > M163V > A96S ≈ Cx40WT). In contrast, G38D channels were almost impermeable to cationic ethidium bromide (EtBr), suggesting that G38D alters channel selectivity. Conversely, A96S and M163V channels showed enhanced EtBr permeability relative to WT Cx40, with the following permeability order: M163V > A96S > Cx40WT > G38D. Altered conductive and permeability properties of mutant channels suggest an essential role for Cx40-mediated biochemical and electrical coupling in cardiac tissues. The altered properties of the three single-base substitution mutants may play a role in mechanisms of reentry arrhythmias.     

Binding of the bovine basic pancreatic trypsin inhibitor (Kunitz) to the 33,000 Mr and 54,000 Mr species of human urokinase: thermodynamic study

The effect of pH and temperature on the apparent association equilibrium constant (Ka) for the binding of the bovine basic pancreatic trypsin inhibitor (BPTI, Kunitz inhibitor) to the 33,000 Mr and 54,000 Mr species of human urokinase (EC 3.4.21.31) has been investigated. Under all the experimental conditions, values of Ka for BPTI binding to the 33,000 Mr and 54,000 Mr species of human urokinase are identical. On lowering the pH from 9.5 to 4.5, values of Ka (at 21.0 degrees C) for BPTI binding to human urokinase (33,000 Mr and 54,000 Mr species) decrease thus reflecting the acidic pK-shift of the His-57 catalytic residue from 6.9, in the free enzyme, to 5.1, in the proteinase:inhibitor complex. At pH 8.0, values of the apparent thermodynamic parameters for BPTI binding to human urokinase (33,000 Mr and 54,000 Mr species) are: Ka = 4.9 x 10(4) M-1, delta G degree = -6.3 kcal/mol, and delta S degree = -37 entropy units (all at 21.0 degrees C); and delta H degree = +4.6 kcal/mol (temperature independent over the explored range, from 5.0 degrees C to 45.0 degrees C). Thermodynamics of BPTI binding to human urokinase (33,000 Mr and 54,000 Mr species) have been analyzed in parallel with those of related serine (pro)enzyme Kazal- and /Kunitz-type inhibitor systems. Considering the known molecular models, the observed binding behaviour of BPTI to human urokinase (33,000 Mr and 54,000 Mr species) was related to the inferred stereochemistry of the proteinase/inhibitor contact region.     

Pairwise coupling analysis of helical junction hydrogen bonding interactions in luteoviral RNA pseudoknots

A 28-nucleotide mRNA pseudoknot that overlaps the P1 and P2 genes of sugarcane yellow leaf virus (ScYLV) stimulates -1 ribosomal frameshifting. The in vitro frameshifting efficiency is decreased >or=8-fold upon substitution of the 3'-most loop 2 nucleotide (C27) with adenosine, which accepts a hydrogen bond from the 2'-OH group of C14 in stem S1. The solution structures of the wild-type (WT) and C27A ScYLV RNA pseudoknots show that while the RNAs adopt virtually identical overall structures, there are significant structural differences at the helical junctions of the two RNAs. Specifically, C8(+) in loop L1 in the C8(+).(G12.C28) L1-S2 major groove base triple is displaced by approximately 2.3 A relative to the accepting stem 2 base pair (G12.C28) in the C27A RNA. Here, we use a double mutant cycle approach to analyze the pairwise coupling of the C8(+).(G12.C28)...C27.(C14-G7) and ...A27.(C14-G7) hydrogen bonds in the WT and C27A ScYLV RNAs, respectively, and compare these findings with previous results from the beet western yellows virus (BWYV) RNA. We find that the pairwise coupling free energy (delta(AB)(i)) is favorable for the WT RNA (-0.7 +/- 0.1 kcal/mol), thus revealing that formation of these two hydrogen bonds is positively cooperative. In contrast, delta(AB)(i) is 0.9 +/- 0.4 kcal/mol for the poorly functional C27A ScYLV RNA, indicative of nonadditive hydrogen bond formation. These results reveal that cooperative hydrogen bond formation across the helical stem junction in H-type pseudoknots correlates with enhanced frameshift stimulation by luteoviral mRNA pseudoknots.     

Homonyms, synonyms and mutations of the sequence/structure vocabulary

The effect of pH and temperature on the association equilibrium constant (Ka) for bovine basic pancreatic trypsin inhibitor (BPTI, Kunitz inhibitor) binding to human urinary kallikrein and porcine pancreatic beta-kallikreins A and B has been investigated. Ka values decrease with decreasing pH, reflecting the acid-midpoint and pK shifts, upon BPTI binding, of a three-proton co-operative transition, between pH 3 and 5, and of a single ionizable group, between pH 5 and 9. At pH 8, the values of delta H degree (between 7 degrees C and 42 degrees C) and delta S degree (at 21 degrees C) for BPTI binding to the glandular kallikreins considered were determined. In particular, the delta H degree values have been found to be independent of temperature and the following values have been obtained by van't Hoff plots: +1.8 kcal/mol, +2.3 kcal/mol and +2.4 kcal/mol (1 kcal = 4184 J) for the inhibitor binding to human urinary kallikrein and porcine pancreatic beta-kallikreins A and B, respectively. Considering the known molecular structures of free porcine pancreatic beta-kallikrein A and BPTI, and of their complex, the stereochemistry of the enzyme : inhibitor contact regions was analysed for the three serine proteinases, in relation to their respective types of behaviour.     

Binding of the bovine basic pancreatic trypsin inhibitor (Kunitz) to human Lys77-plasmin

The effect of pH and temperature on the association equilibrium constant (Ka) for the binding of the bovine basic pancreatic trypsin inhibitor (BPTI Kunitz inhibitor) to human Lys77-plasmin has been investigated. Ka values decrease with decreasing pH, reflecting the acid-pK and -midpoint shifts, upon BPTI binding, of a single ionizable group, between pH 5 and 9, and of a three-proton transition, between pH 3 and 5. At pH 8.0, values of thermodynamic parameters for BPTI binding to human Lys77-plasmin are: Ka = 1.2 X 10(9) M-1, delta G degree = -12.2 kcal/mol, and delta S degree = +49 entropy units (at 21 degrees C); and delta H degree = +2.3 kcal/mol (temperature independent between 5 degrees C and 45 degrees C; 1 kcal = 4184 J). BPTI binding properties of human Lys77-plasmin have been analysed in parallel with those of serine (pro)enzymes acting on cationic and non-cationic substrates. Considering the known molecular structures of homologous serine (pro)enzymes, or Kunitz and Kazal-type inhibitors and of their complexes, the observed binding behaviour of BPTI to human Lys77-plasmin was related to the inferred stereochemistry of the enzyme-inhibitor contact region.     

Differential scanning calorimetric study of the thermal unfolding of mutant forms of phage T4 lysozyme.

In two recent papers, we reported the effects of several point mutations on the thermodynamics of the thermal unfolding of the lysozyme of phage T4 as determined by differential scanning calorimetry. The mutants studied were R96H [Kitamura, S., & Sturtevant, J.M. (1989) Biochemistry 28, 3788-3792] and T157 replaced by A, E, I, L, N, R, and V [Connelly, P., Ghosaini, L., Hu, C.-Q., Kitamura, S., Tanaka, A., & Sturtevant, J.M. (1991) Biochemistry 30, 1887-1891]. Here we report the results of a similar study of the single mutations A82P, A93P, and G113A and the double mutation C54T:C97A. The three single mutants all show small apparent stabilization at pH 2.5 and 46.2 degrees C (the denaturational temperature of the wild-type protein), amounting to -0.5 +/- 0.4 kcal mol-1 in free energy, whereas the double mutant shows a weak apparent destabilization, +0.8 +/- 0.4 kcal mol-1. As in all our previous studies of mutant proteins, the enthalpy changes produced by these mutations are in general of much larger magnitude than the corresponding free energy changes and frequently of opposite sign.     

Comparison of the amino acid residues in the sixth transmembrane domains accessible in the binding-site crevices of mu, delta, and kappa opioid receptors.

We have mapped the residues in the sixth transmembrane domains (TMs 6) of the mu, delta, and kappa opioid receptors that are accessible in the binding-site crevices by the substituted cysteine accessibility method (SCAM). We previously showed that ligand binding to the C7.38S mutants of the mu and kappa receptors and the wild-type delta receptor was relatively insensitive to methanethiosulfonate ethylammonium (MTSEA), a positively charged sulfhydryl-specific reagent. These MTSEA-insensitive constructs were used as the templates, and 22 consecutive residues in TM6 (excluding C6.47) of each receptor were mutated to cysteine, 1 at a time. Most mutants retained binding affinities for [3H]diprenorphine, a nonselective opioid antagonist, similar to that of the template receptors. Treatment with MTSEA significantly inhibited [3H]diprenorphine binding to 11 of 22 mutants of the rat mu receptor and 9 of 22 mutants of the human delta receptor and 10 of 22 mutants of the human kappa receptor. Naloxone or diprenorphine protected all sensitive mutants, except the A6.42(287)C mu mutant. Thus, V6.40, F6.44, W6.48, I6.51, Y6.54, V6.55, I6.56, I6.57, K6.58, and A6.59 of the mu receptor; F6.44, I6.51, F6.54, V6.55, I6.56, V6.57, W6.58, T6.59, and L6.60 of the delta receptor; and F6.44, W6.48, I6.51, F6.54, I6.55, L6.56, V6.57, E6.58, A6.59, and L6.60 of the kappa receptor are on the water-accessible surface of the binding-site crevices. The accessibility patterns of residues in the TMs 6 of the mu, delta, and kappa opioid receptors are consistent with the notion that the TMs 6 are in alpha-helical conformations with a narrow strip of accessibility on the intracellular side of 6.54 and a wider area of accessibility on the extracellular side of 6.54, likely due to a proline kink at 6.50 that bends the helix in toward the binding pocket and enables considerable motion in this region. The wider exposure of residues 6.55-6.60 to the binding-site crevice, combined with the divergent amino acid sequences, is consistent with the inferred role of residues in this region in determining ligand binding selectivity. The conservation of the accessibility pattern on the cytoplasmic side of 6.54 suggests that this region may be important for receptor activation. This accessibility pattern is similar to that of the D2 dopamine receptor, the only other GPCR in which TM6 has been mapped by SCAM. That opioid receptors and the remotely related D2 dopamine receptor have similar accessibility patterns in TM6 suggest that these segments of GPCRs in the rhodopsin-like subfamily not only share secondary structure but also are packed similarly into the transmembrane bundle and thus have similar tertiary structure.     

Destabilizing effect of proline substitutions in two helical regions of T4 lysozyme: leucine 66 to proline and leucine 91 to proline.

A class of temperature-sensitive (ts) mutants of T4 lysozyme with reduced activity at 30 degrees C and no activity at 43 degrees C has been selected. These mutants, designated "tight" ts mutants, differ from most other T4 lysozyme mutants that are active at 43 degrees C, but only manifest their ts lesion by a reduced halo size around phage plaques after exposure of the growth plates to chloroform vapors. For example, in the series of T4 lysozyme mutants at position 157, the original randomly selected mutant, T1571, is the least stable of the series, yet, apart from the halo assay and subsequent in vitro protein stability measurements, this mutant is indistinguishable from wild type (WT) even at 43 degrees C. Two mutants were identified: L91P and L66P. Both insert proline residues into alpha-helical regions of the WT protein structure. The stabilities (delta delta G) as determined by urea denaturation are 8.2 kcal/mol for L91P and 7.1 kcal/mol for L66P. CD spectra indicate that no major conformational changes have occurred in the mutant structures. The structures of the mutants were modeled with a 40-ps molecular dynamics simulation using explicit solvent. For L91P, the reduction of stability appears to be due to an unsatisfied hydrogen bond in the alpha-helix and to a new buried cavity. For L66P, the reduction of stability appears to be due to a disruption of the interdomain alpha-helix, at least two unsatisfied hydrogen bonds, and a newly formed solvent-filled pocket that protrudes into the hydrophobic core, possibly reducing the stabilizing contribution of a partially buried intrachain salt bridge.     

A thermodynamic scale for leucine zipper stability and dimerization specificity: e and g interhelical interactions.

The leucine zipper is a dimeric coiled-coil protein structure composed of two amphipathic alpha-helices with the hydrophobic surfaces interacting to create the dimer interface. This structure has been found to mediate the dimerization of two abundant classes of DNA binding proteins: the bZIP and bHLH-Zip proteins. Several workers have reported that amino acids in the e and g positions of the coiled coil can modulate dimerization stability and specificity. Using the bZIP protein VBP as a host molecule, we report a thermodynamic scale (delta delta G) for 27 interhelical interactions in 35 proteins between amino acids in the g and the following e positions (g<==>e') of a leucine zipper coiled coil. We have examined the four commonly occurring amino acids in the e and g positions of bZIP proteins, lysine (K), arginine (R), glutamine (Q), glutamic acid (E), as well as the only other remaining charged amino acid aspartic acid (D), and finally alanine (A) as a reference amino acid. These results indicate that E<==>R is the most stable interhelical pair, being 0.35 kcal/mol more stable than E<==>K. A thermodynamic cycle analysis shows that the E<==>R pair is 1.33 kcal/mol more stable than A<==>A with -1.14 kcal/mol of coupling energy (delta delta Gint) coming from the interaction of E with R. The E<==>K coupling energy is only -0.14 kcal/mol. E interacts with more specificity than Q. The R<==>R pair is less stable than the K<==>K by 0.24 kcal/mol. R interacts with more specificity than K. Q forms more stable pairs with the basic amino acids K and R rather than with E. Changing amino acids in the e position to A creates bZIP proteins that form tetramers.     

Thermodynamics of single peptide bond cleavage in bovine pancreatic trypsin inhibitor (BPTI)

A major goal of this paper was to estimate a dynamic range of equilibrium constant for the opening of a single peptide bond in a model protein, bovine pancreatic trypsin inhibitor (BPTI). Ten mutants of BPTI containing a single Xaa-->Met substitution introduced in different parts of the molecule were expressed in Escherichia coli. The mutants were folded, purified to homogeneity, and cleaved with cyanogen bromide to respective cleaved forms. Conformation of the intact mutants was similar to the wildtype, as judged from their circular dichroism spectra. Substantial conformational changes were observed on the chemical cleavage of three single peptide bonds--Met46-Ser, Met49-Cys, and Met53-Thr--located within the C-terminal helix. Cleavage of those peptide bonds caused a significant destabilization of the molecule, with a drop of the denaturation temperature by 56.4 degrees C to 68 degrees C at pH 4.3. Opening of the remaining seven peptide bonds was related to a 10.8 degrees C to 39.4 degrees C decrease in T(den). Free energies of the opening of 10 single peptide bonds in native mutants (Delta G(op,N)) were estimated from the thermodynamic cycle that links denaturation and cleavage free energies. To calculate those values, we assumed that the free energy of opening of a single peptide bond in the denatured state (Delta G(op,D)) was equal to -2.7 kcal/mole, as reported previously. Calculated Delta G(op,N) values in BPTI were in the range from 0.2 to 10 kcal/mole, which was equivalent to a >1 million-fold difference in equilibrium constants. The values of Delta G(op,N) were the largest for peptide bonds located in the C-terminal helix and significantly lower for peptide bonds in the beta-structure or loop regions. It appears that opening constants for single peptide bonds in various proteins span across 33 orders of magnitude. Typical equilibrium values for a single peptide bond opening in a protein containing secondary structure elements fall into negligibly low values, from 10(-3) to 10(-8), and are efficient to ensure stability against proteolysis.     

Thermodynamic effects of formamide on DNA stability.

Formamide lowers melting temperatures (Tm) of DNAs linearly by 2.4-2.9 degrees C/mole of formamide (C(F)) depending on the (G+C) composition, helix conformation and state of hydration. The inherent cooperativity of melting is unaffected by the denaturant. dTm/dC(F)for 11 plasmid domains of 0.23 < (G+C)<0.71 generally fit to a linear dependence on (G+C)-content, which, however, is consistent with a (G+C)-independent alteration in the apparent equilibrium constant for thermally induced helix <--> coil transitions. Results indicate that formamide has a destabilizing effect on the helical state, and that sequence-dependent variations in hydration patterns are primarily responsible for small variations in sensitivity to the denaturant. The average unit transition enthalpy delta H(m)[see text for complete expression], exhibits a biphasic dependence on formamide concentration. The initial drop of -0.8 kcal/mol bp at low formamide concentrations is attributable to a delta delta H(m)[see text for complete expression], for exchange of solvent in the vicinity of the helix: displacement by formamide of weakly bound hydrate or counterion. The phenomenological effects are equivalent to lowering the bulk counterion concentration. Poly(dA.dT) exhibits a much lower sensitivity to formamide, due to the specific pattern of tightly bound, immobilized water bridges that buttress the helix from within the narrow minor groove. Tracts of three (A.T)-pairs behave normally, but tracts of six exhibit the same level of reduced sensitivity as the polymer, suggesting a conformational shift as tracts are elongated beyond some critical length [McCarthy J.G. and Rich,A. (1991) Nucleic Acids Res. 19, 3421-3429].     

Crevice-forming mutants in the rigid core of bovine pancreatic trypsin inhibitor: crystal structures of F22A, Y23A, N43G, and F45A.

Crystal structures of four mutants of bovine pancreatic trypsin inhibitor (F22A, Y23A, N43G, and F45A), engineered to alter their stability properties, have been determined. The mutated residues, which are highly conserved among Kunitz-type inhibitors, are located in the rigid core of the molecule. Replacement of the partially buried bulky residues of the wild-type protein with smaller residues resulted in crevices open to the exterior of the molecule. The overall three-dimensional structure of these mutants is very similar to that of the wild-type protein and only small rearrangements are observed among the atoms lining the crevices.     

On the mechanism of preferential incorporation of dAMP at abasic sites in translesional DNA synthesis. Role of proofreading activity of DNA polymerase and thermodynamic characterization of model template-primers containing an abasic site.

DNA polymerase preferentially incorporate dAMP opposite abasic sites (A-rule). The mechanism of the A-rule can be studied by analyzing three dissected stages of the reaction including (i) initial nucleotide insertion, (ii) proofreading excision of the inserted nucleotide and (iii) extension of the nascent primer terminus. To assess the role of the stage (ii) in the A-rule, kinetic parameters of the proofreading excision of primer terminus nucleotides opposite abasic sites were determined using E.coli DNA polymerase I Klenow fragment. The relative efficiency of the excision (Vmax/Km) revealed that removal of A was the least favored of the four nucleotides, but the differences in the efficiencies between excision of A and the other nucleotides was less than 2-fold. In addition, in an attempt to reconcile kinetic data associated with the stage (i) or (ii), the differences in free energy changes (delta delta G degrees) for the formation of model template-primer termini containing XN pairs (X = abasic site, N = A, G, C or T) were determined by temperature dependent UV-melting measurements. The order of delta delta G degrees was XG > XA = XC > or = XT, with delta delta G degrees being 0.5 kcal/mol for the most stable XG and the least stable XT. Based on these data, the role of the stage (ii) and energetic aspects of the A-rule are discussed.     

A molecular model for lipid-protein interaction in membranes: the role of hydrophobic mismatch.

The interaction free energy between a hydrophobic, transmembrane, protein and the surrounding lipid environment is calculated based on a microscopic model for lipid organization. The protein is treated as a rigid hydrophobic solute of thickness dP, embedded in a lipid bilayer of unperturbed thickness doL. The lipid chains in the immediate vicinity of the protein are assumed to adjust their length to that of the protein (e.g., they are stretched when dP > doL) in order to bridge over the lipid-protein hydrophobic mismatch (dP-doL). The bilayer's hydrophobic thickness is assumed to decay exponentially to its asymptotic, unperturbed, value. The lipid deformation free energy is represented as a sum of chain (hydrophobic core) and interfacial (head-group region) contributions. The chain contribution is calculated using a detailed molecular theory of chain packing statistics, which allows the calculation of conformational properties and thermodynamic functions (in a mean-field approximation) of the lipid tails. The tails are treated as single chain amphiphiles, modeled using the rotational isometric state scheme. The interfacial free energy is represented by a phenomenological expression, accounting for the opposing effects of head-group repulsions and hydrocarbon-water surface tension. The lipid deformation free energy delta F is calculated as a function of dP-doL. Most calculations are for C14 amphiphiles which, in the absence of a protein, pack at an average area per head-group ao approximately equal to 32 A2 (doL approximately 24.5 A), corresponding to the fluid state of the membrane. When dP = doL, delta F > 0 and is due entirely to the loss of conformational entropy experienced by the chains around the protein. When dP > doL, the interaction free energy is further increased due to the enhanced stretching of the tails. When dP < doL, chain flexibility (entropy) increases, but this contribution to delta F is overcounted by the increase in the interfacial free energy. Thus, delta F obtains a minimum at dP-doL approximately 0. These qualitative interpretations are supported by detailed numerical calculations of the various contributions to the interaction free energy, and of chain conformational properties. The range of the perturbation of lipid order extends typically over few molecular diameters. A rather detailed comparison of our approach to other models is provided in the discussion.     

Contribution of a tyrosine side chain to ribonuclease A catalysis and stability.

An intricate architecture of covalent bonds and noncovalent interactions appear to position the side chain of Lys 41 properly within the active site of bovine pancreatic ribonuclease A (RNase A). One of these interactions arises from Tyr 97, which is conserved in all 41 RNase A homologues of known sequence. Tyr 97 has a solvent-inaccessible side chain that donates a hydrogen bond to the main-chain oxygen of Lys 41. Here, the role of Tyr 97 was examined by replacing Tyr 97 with a phenylalanine, alanine, or glycine residue. All three mutant proteins have diminished catalytic activity, with the value of Kcat being perturbed more significantly than that of Km. The free energies with which Y97F, Y97A, and Y97G RNase A bind to the rate-limiting transition state during the cleavage of poly(cytidylic acid) are diminished by 0.74, 3.3, and 3.8 kcal/mol, respectively. These results show that even though Tyr 97 is remote from the active site, its side chain contributes to catalysis. The role of Tyr 97 in the thermal stability of RNase A is large. The conformational free energies of native Y97F, Y97A, and Y97G RNase A are decreased by 3.54, 12.0, and 11.7 kcal/mol, respectively. The unusually large decrease in stability caused by the Tyr-->Phe mutation could result from a decrease in the barrier to isomerization of the Lys 41-Pro 42 peptide bond.