Refined 1.2 A crystal structure of the complex formed between subtilisin Carlsberg and the inhibitor eglin c. Molecular structure of eglin and its detailed interaction with subtilisin.

The crystal structure of the complex formed between eglin c, an elastase inhibitor from the medical leech, and subtilisin Carlsberg has been determined at 1.2 A resolution by a combination of Patterson search methods and isomorphous replacement techniques. The structure has been refined to a crystallographic R-value of 0.18 (8-1.2 A). Eglin consists of a four-stranded beta-sheet with an alpha-helical segment and the protease-binding loop fixed on opposite sides. This loop, which contains the reactive site Leu45I--Asp46I, is mainly held in its conformation by unique electrostatic/hydrogen bond interactions of Thr44I and Asp46I with the side chains of Arg53I and Arg51I which protrude from the hydrophobic core of the molecule. The conformation around the reactive site is similar to that found in other proteinase inhibitors. The nine residues of the binding loop Gly40I--Arg48I are involved in direct contacts with subtilisin. In this interaction, eglin segment Pro42I--Thr44I forms a three-stranded anti-parallel beta-sheet with subtilisin segments Gly100--Gly102 and Ser125--Gly127. The reactive site peptide bond of eglin is intact, and Ser221 OG of the enzyme is 2.81 A apart from the carbonyl carbon.     

Changing the inhibitory specificity and function of the proteinase inhibitor eglin c by site-directed mutagenesis: functional and structural investigation.

Amino acids in the serine proteinase inhibitor eglin c important for its inhibitory specificity and activity have been investigated by site-directed mutagenesis. The specificity of eglin c could be changed from elastase to trypsin inhibition by the point mutation Leu45----Arg (L45R) in position P1 [nomenclature according to Schechter and Berger (1967) Biochem. Biophys. Res. Commun. 27, 157-162]. Model building studies based on the crystal structure of mutant L45R [Heinz et al. (1991) J. Mol. Biol. 217, 353-371] were used to rationalize this specificity change. Surprisingly, the double mutant L45R/D46S was found to be a substrate of trypsin and various other serine proteinases. Multidimensional NMR studies show that wild-type eglin c and the double mutant have virtually identical conformations. In the double mutant L45R/D46S, however, the N-H bond vector of the scissile peptide bond shows a much higher mobility, indicating that the internal rigidity of the binding loop is significantly weakened due to the loss or destabilization of the internal hydrogen bond of the P1' residue. Mutant T44P was constructed to examine the role of a proline in position P2, which is frequently found in serine proteinase inhibitors [Laskowski and Kato (1980) Annu. Rev. Biochem. 49, 593-626]. The mutant remains a potent elastase inhibitor but no longer inhibits subtilisin, which could be explained by model building. Both Arg51 and Arg53, located in the core of the molecule and participating in the hydrogen bonding network with residues in the binding loop to maintain rigidity around the scissile bond, were individually replaced with the shorter but equally charged amino acid lysine. Both mutants showed a decrease in their inhibitory potential. The crystal structure of mutant R53K revealed the loss of two hydrogen bonds between the core and the binding loop of the inhibitor, which are partially restored by a solvent molecule, leading to a decrease in inhibition of elastase by 2 orders of magnitude.     

Engineered eglin c variants inhibit yeast and human proprotein processing proteases, Kex2 and furin

We engineered eglin c, a potent subtilisin inhibitor, to create inhibitors for enzymes of the Kex2/furin family of proprotein processing proteases. A structural gene was synthesized that encoded "R(1)-eglin", having Arg at P(1) in the reactive site loop in place of Leu(45). Ten additional variants were created by cassette mutagenesis of R(1)-eglin. These polypeptides were expressed in Escherichia coli, purified to homogeneity, and their interactions with secreted, soluble Kex2 and furin were examined. R(1)-eglin itself was a modest inhibitor of Kex2, with a K(a) of approximately 10(7) M(-)(1). Substituting Arg (in R(4)R(1)-eglin) or Met (in M(4)R(1)-eglin) for Pro(42) at P(4) created potent Kex2 inhibitors exhibiting K(a) values of approximately 10(9) M(-)(1). R(4)R(1)-eglin inhibited furin with a K(a) of 4.0 x 10(8) M(-)(1). Introduction of Lys at P(1), in place of Arg in R(4)R(1)-eglin reduced affinity only approximately 3-fold for Kex2 but 15-fold for furin. The stabilities of enzyme-inhibitor complexes were characterized by association and dissociation rate constants and visualized by polyacrylamide gel electrophoresis. R(4)R(1)-eglin formed stable 1:1 complexes with both Kex2 and furin. However, substitution of Lys at P(2) in place of Thr(44) resulted in eglin variants that inhibited both Kex2 and furin but which were eventually cleaved (temporary inhibition). Surprisingly, R(6)R(4)R(1)-eglin, in which Arg was substituted for Gly(40) in R(4)R(1)-eglin, exhibited stable, high-affinity complex formation with Kex2 (K(a) of 3.5 x 10(9) M(-)(1)) but temporary inhibition of furin. This suggests that enzyme-specific interactions can alter the conformation of the reactive site loop, converting a permanent inhibitor into a substrate. Eglin variants offer possible avenues for affinity purification, crystallization, and regulation of proprotein processing proteases.     

Binding of phosphorus-containing inhibitors to thermolysin studied by the Poisson-Boltzmann method.

Zinc endopeptidase thermolysin can be inhibited by a series of phosphorus-containing peptide analogues, Cbz-Gly-psi (PO2)-X-Leu-Y-R (ZGp(X)L(y)R), where X = NH, O, or CH2; Y = NH or O; R = Leu, Ala, Gly, Phe, H, or CH3. The affinity correlation as well as an X-ray crystallography study suggest that these inhibitors bind to thermolysin in an identical mode. In this work, we calculate the electrostatic binding free energies for a series of 13 phosphorus-containing inhibitors with modifications at X, Y, and R moieties using finite difference solution to the Poisson-Boltzmann equation. A method has been developed to include the solvation entropy changes due to binding different ligands to a macromolecule. We demonstrate that the electrostatic energy and empirically derived solvation entropy can account for most of the binding energy differences in this series. By analyzing the binding contribution from individual residues, we show that the energy of a hydrogen bond is not confined to the donor and acceptor. In particular, the positive charges on Zn and Arg 203, which are not the acceptors, contribute significantly to the hydrogen bonds between two amides of ZGpLL and the thermolysin.     

Examination of the reduced affinity of the thymidylate synthase G52S mutation for FdUMP by ab initio and semi-empirical studies

BACKGROUND: The G52S mutation in the Arg50 loop of thymidylate synthase leads to decreased binding of FdUMP. It has been suggested that the mutation affects the Arg50 residue (within the Arg50 loop) responsible for binding the phosphate of FdUMP. The binding of the methylguanidinium moiety as a model for Arg50 to a methylphosphate entity as a model for FdUMP was investigated with theoretical calculations, as well as the structure of the Arg50-Thr51-Gly52 tripeptide in comparison with the Arg50-Thr51-Ser52 tripeptide. METHODS: Gaussian-98 and PC Spartan programs were used to perform Hartree-Fock and Post-Hartree-Fock quantum chemical calculations as well as MNDO (semi-empirical calculations). RESULTS: It was found that the strongest binding occurs between the negative methylphosphate ion and methylguanidine. The replacement of Gly52 by Ser52 leads to a significant displacement of Arg50, which may be responsible for the decreased binding to FdUMP. CONCLUSION: The arginine-phosphate binding appears to be geometry dependent. Thus, the displacement of the Arg50 residue, as observed in these calculated models, upon mutation of Gly52 to Ser may contribute to decreased binding of FdUMP to mTS (G52S).     

Molecular dynamics simulation,binding free energy calculation and molecular docking of human D-amino acid oxidase (DAAO) with its inhibitors

Schizophrenia is a mental illness; most affected people live in developing countries, and neither appropriate treatment nor commercial drugs are currently available. One possibility is to inhibit human-d-amino acid oxidase (h-DAAO). In this study, molecular dynamic simulations of the monomer, dimer and tetramer forms of h-DAAO complexed with the inhibitor 3-hydroxyquinolin-2(1H)-one(2) were performed. Seven residues, Leu51, Gln53, Leu215, Tyr228, Ile230, Arg283 and Gly313, were identified as essential for interacting with the inhibitor. Molecular docking of h-DAAO with pyrrole, quinoline and kojic acid derivatives, representing 69 known or potential h-DAAO inhibitors, was also performed. The results indicated that the activity of the inhibitor can be improved by modifying the compounds to have a substituent group capable of interacting with the side chain of Tyr228. Van der Waals interactions of the inhibitor with the hydrophobic pocket of h-DAAO and electrostatic interactions or H-bonds with Arg283 and Gly313 were important elements in determining the efficiency of the inhibitor. These results provide information on the interaction between h-DAAO and its inhibitors at the molecular level and can aid in the design of novel inhibitors against h-DAAO for new drug development in the treatment of schizophrenia.     

Steered molecular dynamics simulation of the binding of the β2 and β3 regions in domain-swapped human cystatin C dimer

The crystal structure of the human cystatin C (hCC) dimer revealed that a stable twofold-symmetric dimer was formed via 3D domain swapping. Domain swapping with the need for near-complete unfolding has been proposed as a possible route for amyloid fibril initiation. Thus, the interesting interactions that occur between the two molecules may be important for the further aggregation of the protein. In this work, we performed steered molecular dynamics (SMD) simulations to investigate the dissociation of the β2 and β3 strands in the hCC dimer. The energy changes observed during the SMD simulations showed that electrostatic interactions were the dominant interactions involved in stabilizing the two parts of the dimer during the early stages of SMD simulation, whereas van der Waals (VDW) interactions and electrostatic interactions were equally matched during the latter stages. Furthermore, our data indicated that the two parts of the dimer are stabilized by intermolecular hydrogen bonds among the residues Arg51 (β2), Gln48 (β2), Asp65 (β3), and Glu67 (β3), salt bridges among the residues Arg53 (β2), Arg51 (β2), and Asp65 (β3), and VDW interactions among the residues Gln48 (β2), Arg51 (β2), Glu67 (β3), Asp65 (β3), Phe63 (β3), and Asn61 (β3). The residues Gln48 (β2), Arg51 (β2), Asp65 (β3) and Glu67 (β3) appear to be crucial, as they play important roles in both electrostatic and VDW interactions. Thus, the present study determined the key residues involved in the stabilization of the domain-swapped dimer structure, and also provided molecular-level insights into the dissociation process of the hCC dimer.     

Anatomy of specific interactions between lambda repressor and operator DNA

Recognition of specific DNA sequences by proteins is essential for regulation of gene expression. To fully understand the recognition mechanism, it is necessary to understand not only the structure of the specific protein-DNA interactions but also the energetics. We therefore performed a computer analysis in which a phage DNA-binding protein, lambda repressor, was used to examine the changes in binding free energy (DeltaDeltaG) and its energy components caused by single base mutations. We then determined which of the calculated energy components best correlated with the experimental data. The experimental DeltaDeltaG values were well reproduced by the calculations. Component analysis revealed that the electrostatic and hydrogen bond energies were most strongly correlated with the experimental data. Among the 51 single base-substitution mutants examined, positive DeltaDeltaG values, corresponding to weakened binding, were caused by the loss of favorable electrostatic interactions and hydrogen bonds, the introduction of steric collisions and electrostatic repulsion, the loss of favorable interactions with a thymine methyl group, and the increase of unfavorable hydration energy from isolated DNA. This analysis also showed distinct patterns of recognition at A-T and G-C positions, as different combinations of energy components were involved in DeltaDeltaG caused by the two substitution types. We have thus been able to identify the energy components that most strongly correlate with sequence-dependent DeltaDeltaG and determine their contribution to the specificity of DNA sequence recognition by the lambda repressor. Application of this method to other systems should provide additional insight into the molecular mechanism of protein-DNA recognition.     

Metal-tetracycline/H+ antiporter of Escherichia coli encoded by transposon Tn10. The role of a conserved sequence motif, GXXXXRXGRR, in a putative cytoplasmic loop between helices 2 and 3.

The region including the conserved Ser65-Asp66 dipeptide in the tetracycline/H+ antiporter (TET) encoded by transposon Tn10 is thought to play a gating role (Yamaguchi, A., Ono, N., Akasaka, T., Noumi, T., and Sawai, T. (1990) J. Biol. Chem. 265, 15525-15530). The dipeptide is in putative interhelix loop2-3, which also includes the conserved sequence motif, GXXXXRXGRR, found in all TET proteins and sugar/H+ symporters. Through the combination of localized random and site-directed mutagenesis, each residue in loop2-3 was replaced. Among 10 residues in putative loop2-3, the important residues, of which substitution resulted in significant reduction or complete loss of the transport activity, were Gly62, Asp66, Gly69, and Arg70. The defect in the transport activity of the Gly62 and Gly69 substitution mutants corresponded to the steric hindrance by the substituents as to the putative beta-turn structure of the peptide backbone containing these glycines. Of 3 conserved Arg residues, the replacement of only Arg70 caused complete loss of the activity except for replacement with Lys, indicating the importance of a positive charge at this position, which is similar to the essentiality of a negative charge at Asp66. A "charge-neutralizing" intra-loop salt bridge between Asp66 and Arg70 was not likely because the double mutant in which Asp66 and Arg70 were replaced with asparagine and leucine, respectively, showed no transport activity. A triple mutant with only one positive charge at Arg70 in this loop showed about half the wild-type activity, indicating that the polycationic nature of the loop was not critical for the activity. Cys mutants as to the unessential residues in the loop were modifiable with N-ethylmaleimide, except for the Met64----Cys and Arg71----Cys mutants; however, the modification of only the Ser65----Cys mutant caused significant inhibition of the transport activity, indicating that position 65 is a unique position in the structure of loop2-3.     

Crystal and molecular structure of the bovine alpha-chymotrypsin-eglin c complex at 2.0 A resolution.

The crystal structure of the complex between bovine alpha-chymotrypsin and the leech (Hirudo medicinalis) protein proteinase inhibitor eglin c has been refined at 2.0 A resolution to a crystallographic R-factor of 0.167. The structure of the complex includes 2290 protein and 143 solvent atoms. Eglin c is bound to the cognate enzyme through interactions involving 11 residues of the inhibitor (sites P5-P4' in the reactive site loop, P10' and P23') and 17 residues from chymotrypsin. Binding of eglin c to the enzyme causes a contained hinge-bending movement around residues P4 and P4' of the inhibitor. The tertiary structure of chymotrypsin is little affected, with the exception of the 10-13 region, where an ordered structure for the polypeptide chain is observed. The overall binding mode is consistent with those found in other serine proteinase-protein-inhibitor complexes, including those from different inhibition families. Contained, but significant differences are observed in the establishment of intramolecular hydrogen bonds and polar interactions stabilizing the structure of the intact inhibitor, if the structure of eglin c in its complex with chymotrypsin is compared with that of other eglin c-serine proteinase complexes.     

Thermodynamic analysis of the binding of aromatic hydroxamic acid analogues to ferric horseradish peroxidase.

Peroxidases typically bind their reducing substrates weakly, with K(d) values in the millimolar range. The binding of benzhydroxamic acid (BHA) to ferric horseradish peroxidase isoenzyme C (HRPC) [K(d) = 2.4 microM; Schonbaum, G. R. (1973) J. Biol. Chem. 248, 502-511] is a notable exception and has provided a useful tool for probing the environment of the peroxidase aromatic-donor-binding site and the distal heme cavity. Knowledge of the underlying thermodynamic driving forces is key to understanding the roles of the various H-bonding and hydrophobic interactions in substrate binding. The isothermal titration calorimetry results of this study on the binding of aromatic hydroxamic acid analogues to ferric HRPC under nonturnover conditions (no H(2)O(2) present) confirm the significance of H-bonding interactions in the distal heme cavity in complex stabilization. For example, the binding of BHA to HRPC is enthalpically driven at pH 7.0, with the H-bond to the distal Arg38 providing the largest contribution (6.74 kcal/mol) to the binding energy. The overall relatively weak binding of the hydroxamic acid analogues to HRPC is due to large entropic barriers (-11.3 to -37.9 eu) around neutral pH, with the distal Arg38 acting as an "entropic gate keeper". Dramatic enthalpy-entropy compensation is observed for BHA and 2-naphthohydroxamic acid binding to HRPC at pH 4.0. The enthalpic loss and entropic gain are likely due to increased flexibility of Arg38 in the complexes at low pH and greater access by water to the active site. Since the Soret absorption band of HRPC is a sensitive probe of the binding of hydroxamic acids and their analogues, it was used to investigate the binding of six donor substrates over the pH range of 4-12. The negligible pH dependence of the K(d) values corrected for substrate ionization suggests that enthalpy-entropy compensation is operative over a wide pH range. Examination of the thermodynamics of binding of ring-substituted hyrazides to HRPC reveals that the binding affinities of aromatic donors are highly sensitive to the position and nature of the ring substituent.     

Design and synthesis of novel celecoxib analogues as selective cyclooxygenase-2 (COX-2) inhibitors: replacement of the sulfonamide pharmacophore by a sulfonylazide bioisostere

A group of celecoxib analogues in which the para-SO(2)NH(2) substituent on the N(1)-phenyl ring was replaced by a para-sulfonylazido (SO(2)N(3)) 4, or a meta-SO(2)N(3) 8, substituent were designed for evaluation as selective cyclooxygenase-2 (COX-2) inhibitors. In vitro COX-1 and COX-2 inhibition studies showed that 4-[5-(4-methylphenyl)-3-trifluoromethyl-1H-pyrazol-1-yl]benzenesulfonyl azide (4) with a para-SO(2)N(3) substituent was a selective COX-1 inhibitor. In contrast, 3-[5-(4-methylphenyl)-3-trifluoromethylpyrazol-1-yl]benzenesulfonyl azide (8a) having a meta-SO(2)N(3) substituent (COX-1 IC(50) >100microM; COX-2 IC(50)=5.16microM; COX-2 selectivity index >19.3) is a selective COX-2 inhibitor. A molecular modeling (docking) study showed that the SO(2)N(3) group of 8a inserts deep inside the secondary pocket of the COX-2 binding site. The SO(2)N(3) moiety of 8a can undergo a dual H-bonding interaction via one of its SO(2) oxygen-atoms, and an electrostatic (ion-ion) interaction via the terminal azido (N(3)) nitrogen-atom, to the guanidino NH(2) of Arg(513) in the secondary pocket of COX-2. These observations indicate that an appropriately positioned SO(2)N(3) moiety is a novel alternative bioisostere to the traditional SO(2)NH(2) and SO(2)Me pharmacophores present in selective COX-2 inhibitors, that are only capable of H-bonding interactions with the COX-2 isozyme, for use in drug design.     

X-ray crystal structure of the serine proteinase inhibitor eglin c at 1.95 A resolution.

The crystal structure of eglin c, naturally occurring in the leech Hirudo medicinalis, is known from its complexes with various serine proteinases, but the crystallization of free eglin c has not yet been reported. A method is described for growing well-diffracting crystals of free eglin c from highly concentrated protein solutions (approximately 200 mg/ml). The space group of the orthorhombic crystals was determined to be P2(1)2(1)2(1) with unit cell parameters a = 32.6, b = 42.0, c = 44.1 A. The structure of free eglin c was resolved at 1.95 A resolution by Patterson search methods. The final model contains all 70 amino acids of eglin c and 125 water molecules. In comparison to the eglin structure known from its complexes with proteinases, only small differences have been observed in free eglin c. However, the reactive site-binding loop and a few residues on the surface of eglin have been found in different conformations due to crystal contacts. In contrast to the complex structures, the first seven amino acids of the highly flexible amino terminus can be located. Crystallographic refinement comprised molecular dynamics refinement, classical restrained least-squares refinement and individual isotropic atomic temperature refinement. The final R-factor is 15.8%.     

Inhibition of HIV-1 proteinase by non-peptide carboxylates

Some simple dicarboxylates are among the first reported non-peptide inhibitors of HIV-1 proteinase. Only weak inhibition (IC50 greater than or equal to 10 microM) was observed but this may be significant since only two potential enzyme-binding groups are present. Dixon plots and preliminary kinetic data are reported and a possible mechanism for the inhibition is discussed. The dicarboxylates are long enough to engage the carboxylate side chains of Arg 8 and Arg 108 at either end of the 24A long substrate-binding groove. This mode of binding has not been proven but other molecules with similarly separated charged ends are equally effective inhibitors, perhaps indicating a common mechanism of inhibition. There is evidence that placing other functional groups on the inhibitor enables alternative interactions with the enzyme which can reduce inhibitor potency. We propose that incorporation of ionic binding groups in more elaborate and selective non-peptides may potentiate inhibition of HIV-1 proteinase.     

Constraining specificity in the N-domain of tissue inhibitor of metalloproteinases-1; gelatinase-selective inhibitors

The tissue inhibitors of metalloproteinases (TIMPs) are endogenous inhibitors of the matrix metalloproteinases (MMPs). Since unregulated MMP activities are linked to arthritis, cancer, and atherosclerosis, TIMP variants that are selective inhibitors of disease-related MMPs have potential therapeutic value. The structures of TIMP/MMP complexes reveal that most interactions with the MMP involve the N-terminal pentapeptide of TIMP and the C-D beta-strand connector which occupy the primed and unprimed regions of the active site. The loop between beta-strands A and B forms a secondary interaction site for some MMPs, ranging from multiple contacts in the TIMP-2/membrane type-1 (MT1)-MMP complex to none in the TIMP-1/MMP-1 complex. TIMP-1 and its inhibitory domain, N-TIMP-1, are weak inhibitors of MT1-MMP; inhibition is not improved by grafting the longer AB loop from TIMP-2 into N-TIMP-1, but this change impairs binding to MMP-3 and MMP-7. Mutational studies with N-TIMP-1 suggest that its weak inhibition of MT1-MMP, as compared to other N-TIMPs, arises from multiple (>3) sequence differences in the interaction site. Substitutions for Thr2 of N-TIMP-1 strongly influence MMP selectivity; Arg and Gly, that generally reduce MMP affinity, have less effect on binding to MMP-9. When the Arg mutation is added to the N-TIMP-1(AB2) mutant, it produces a gelatinase-specific inhibitor with Ki values of 2.8 and 0.4 nM for MMP-2 and -9, respectively. Interestingly, the Gly mutant has a Ki of 2.1 nM for MMP-9 and >40 muM for MMP-2, indicating that engineered TIMPs can discriminate between MMPs in the same subfamily.     

What can the structures of enzyme-inhibitor complexes tell us about the structures of enzyme substrate complexes?

Proteinases perform many beneficial functions that are essential to life, but they are also dangerous and must be controlled. Here we focus on one of the control mechanisms: the ubiquitous presence of protein proteinase inhibitors. We deal only with a subset of these: the standard mechanism, canonical protein inhibitors of serine proteinases. Each of the inhibitory domains of such inhibitors has one reactive site peptide bond, which serves all the cognate enzymes as a substrate. The reactive site peptide bond is in a combining loop which has an identical conformation in all inhibitors and in all enzyme-inhibitor complexes. There are at least 18 families of such inhibitors. They all share the conformation of the combining loops but each has its own global three-dimensional structure. Many three-dimensional structures of enzyme-inhibitor complexes were determined. They are frequently used to predict the conformation of substrates in very short-lived enzyme-substrate transition state complexes. Turkey ovomucoid third domain and eglin c have a Leu residue at P(1). In complexes with chymotrypsin, these P(1) Leu residues assume the same conformation. The relative free energies of binding of P(1) Leu (relative to either P(1) Gly or P(1) Ala) are within experimental error, the same for complexes of turkey ovomucoid third domain, eglin c, P(1) Leu variant of bovine pancreatic trypsin inhibitor and of a substrate with chymotrypsin. Therefore, the P(1) Leu conformation in transition state complexes is predictable. In contrast, the conformation of P(1) Lys(+) is strikingly different in the complexes of Lys(18) turkey ovomucoid third domain and of bovine pancreatic trypsin inhibitor with chymotrypsin. The relative free energies of binding are also quite different. Yet, the relative free energies of binding are nearly identical for Lys(+) in turkey ovomucoid third domain and in a substrate, thus allowing us to know the structure of the latter. Similar reasoning is applied to a few other systems.     

Synthesis and AT2 receptor-binding properties of angiotensin II analogues.

The present study investigates the importance of the amino acid side chains in the octapeptide angiotensin II (Ang II) for binding to the AT2 receptor. A Gly scan was performed where each amino acid in Ang II was substituted one-by-one with glycine. The resulting set of peptides was tested for affinity to the AT2 receptor (porcine myometrial membranes). For a comparison, the peptides were also tested for affinity to the AT1 receptor (rat liver membranes). Only the substitution of Arg2 reduced affinity to the AT2 receptor considerably (92-fold when compared with Ang II). For the other Gly-substituted analogues the affinity to the AT2 receptor was only moderately affected. To further investigate the role of the Arg2 side chain for receptor binding, we synthesized some N-terminally modified Ang II analogues. According to these studies a positive charge in the N-terminal end of angiotensin III [Ang II (2-8)] is not required for high AT2 receptor affinity but seems to be more important in Ang II. With respect to the AT1 receptor, [Gly2]Ang II and [Gly8]Ang II lacked binding affinity (Ki > 10 microM). Replacement of the Val3 or Ile5 residues with Gly produced only a slight decrease in affinity. Interestingly, substitution of Tyr4 or His6, which are known to be very important for AT1 receptor binding, resulted in only 48 and 14 times reduction in affinity, respectively.     

The high-resolution X-ray crystal structure of the complex formed between subtilisin Carlsberg and eglin c, an elastase inhibitor from the leech Hirudo medicinalis. Structural analysis, subtilisin structure and interface geometry

Triclinic crystals of the complex formed by eglin with subtilisin Carlsberg were analyzed by X-ray diffraction. The crystal and molecular structure of this complex was determined with data that extended to 0.12-nm resolution by a combination of Patterson search methods and isomorphous replacement techniques. Its structure was refined to a crystallographic R value of 0.178 (1.0-0.12 nm) using an energy-restraint least-squares procedure. The complete subtilisin molecule could be traced without ambiguity in the refined electron density. The eglin component, from which an amino-terminal segment is cleaved off, is only defined from Lys8I (i.e. the lysine residue 8 of the inhibitor) onwards. Per unit cell, 436 fixed solvent molecules and 2 calcium ions were located. In spite of 84 amino acid replacements and one deletion, subtilisin Carlsberg exhibits a very similar polypeptide fold to subtilisin BPN'. The root-mean-square deviations of all alpha-carbon atoms (excluding those at the deletion site) from models of subtilisin BPN' [Alden, R. A., Birktoft, J. J., Kraut, J., Robertus, J. D. & Wright, C. S. (1971) Biochem. Biophys. Res. Commun. 45, 337-344] and subtilisin Novo [Drenth, J., Hol, W. G. J., Jansonius, J. N. & Kockoek, R. (1972) Eur. J. Biochem. 25, 177-181] are 0.077 nm and 0.103 nm. Most of these deviations result from global shifts rather than changes of the local geometry. The single-residue deletion at position 56 affects only the surrounding conformation. Two sites of high electron density and close distances to surrounding oxygen ligands have been found in the Carlsberg enzyme which are probably occupied by calcium ions. Eglin consists of a twisted four-stranded beta-sheet flanked by an alpha-helix and by an exposed proteinase binding loop on opposite sides. Around the reactive site, Leu45I-Asp46I, this loop is mainly stabilized by electrostatic/hydrogen bond interactions with the side chains of two arginine residues which project from the hydrophobic core [Bode, W., Papamokos, E., Musil, D., Seemüller, W. & Fritz, H. (1986) EMBO J. 5, 813-818]. The reactive site loop conformation resembles that found in other 'small' proteinase inhibitors. The scissile peptide bond is not cleaved but its carbonyl group is slightly distorted from planar geometry. Most of the intermolecular contacts are contributed by the nine residues of the reactive-site loop Gly40I-Arg48I.(ABSTRACT TRUNCATED AT 400 WORDS)     

Role of active site binding interactions in 4-chlorobenzoyl-coenzyme A dehalogenase catalysis.

4-Chlorobenzoyl-coenzyme A (4-CBA-CoA) dehalogenase catalyzes the hydrolytic dehalogenation of 4-CBA-CoA to 4-hydroxybenzoyl-CoA (4-HBA-CoA) via a multistep mechanism involving initial attack of Asp145 on C(4) of the substrate benzoyl ring to form a Meisenheimer intermediate (EMc), followed by expulsion of the chloride ion to form an arylated enzyme intermediate (EAr) and then ester hydrolysis in the EAr to form product. This study examines the role of binding interactions in dehalogenase catalysis. The enzyme and substrate groups positioned for favorable binding interaction were identified from the X-ray crystal structure of the enzyme-4-HBA-3'-dephospho-CoA complex. These groups were individually modified (via site-directed mutagenesis or chemical synthesis) for the purpose of disrupting the binding interaction. The changes in the Gibbs free energy of the enzyme-substrate complex (DeltaDeltaG(ES)) and enzyme-transition state complex (DeltaDeltaG) brought about by the modification were measured. Cases where DeltaDeltaG exceeds DeltaDeltaG(ES) are indicative of binding interactions used for catalysis. On the basis of this analysis, we show that the H-bond interactions between the Gly114 and Phe64 backbone amide NHs and the substrate benzoyl C=O group contribute an additional 3.1 kcal/mol of stabilization at the rate-limiting transition state. The binding interactions between the enzyme and the substrate CoA nucleotide moiety also intensify in the rate-limiting transition state, reducing the energy barrier to catalysis by an additional 3.3 kcal/mol. Together, these binding interactions contribute approximately 10(6) to the k(cat)/K(m).     

Increased rigidity of eglin c at acidic pH: evidence from NMR spin relaxation and MD simulations

To gain physical insights into how proteins respond to changes in pH, the picosecond to nanosecond time scale dynamics of the small serine protease inhibitor eglin c have been studied by NMR spin relaxation experiments and MD simulations under two pH solution conditions, pH 7 and 3. Like many proteins, eglin c is destabilized by a lowering of the pH, although it retains enough stability to maintain its native conformation at pH 3. Backbone (15)N relaxation results show comparable global tumbling times (tau(m)) and model-free order parameters (S(2)) under the two pH conditions, indicating that the molecule maintains its overall molecular shape and structure at low pH, although the backbone rigidity is slightly increased (/ = 0.6%). In contrast, the side-chain methyl dynamics, as measured from (2)H relaxation experiments, show a substantial increase in rigidity at lower pH (/ = 14.8%). Molecular dynamics simulations performed at these pH states produce results consistent with NMR measurements, showing that the two methods are in qualitative agreement. Although a full accounting of the physical basis for the concurrent conformational rigidification and destabilization at low pH requires further investigation, the high level of detail in the MD simulations provides a potential molecular mechanism: the breaking of the hydrogen bond between the side chains of Asp46 and Arg53, and changes in electrostatic interactions, appear to allow the binding loop to move closer to the core part of the protein, resulting in a more compact structure at low pH. This more compact structure may be responsible for the increased level of restriction of molecular motion. As these findings show, the stability of a molecular structure is distinct from its conformational rigidity, and the two can even change in opposite directions, against na?ve expectation.