1 Ludwig Institute for Cancer Research, PO Box 2008, Royal Melbourne Hospital, Parkville, VIC 3050, 2 Cytopia Pty Ltd, 7th Floor, Daly Wing,St. Vincent's Hospital, 41 Victoria Parade, Fitzroy, VIC 3065, Australia and 4 Cooperative Research Centre for Cellular Growth Factors, Royal Melbourne Hospital, Parkville, VIC 3050, Australia
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Keywords: computational combinatorial chemistry/inhibitor design/ras protein
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Computational combinatorial inhibitor design is one of the newly developed structure-based de novo design methods (Caflisch et al., 1993; Caflisch and Karplus, 1995
, 1996
; Caflisch, 1996
; 1996; Joseph-McCarthy et al., 1997). It consists of three steps: first, the positions of functional chemical groups are identified on the key regions of a target protein using an exhaustive `multiple copy simultaneous search' (MCSS) approach (Caflisch et al., 1993
); second, the MCSS minima of a specific functional group are used to construct links to the other functional groups; finally, potential inhibitors are selected according to their calculated or estimated binding affinities to the target protein. Previously, we had developed a novel scheme that allows the application of computational combinatorial inhibitor design to identify peptide inhibitors that target the protein surface where the details of the peptide-protein binding site are not clearly defined (Zeng and Treutlein, 1999
; Zeng, 2000
). In our scheme, the third step involves the sequence alignment of designed peptides and evaluation of the probability of a specific amino acid type at each position of the peptides. Peptide sequences of potential inhibitors were based on the amino acid preferences at each position, rather than on approximately calculated binding affinities (Caflisch and Karplus, 1995
, 1996
; Caflisch, 1996
). Our method is therefore also suitable to design combinatorial libraries for specific drug targets and to identify and analyze amino acid sequence patterns of potential inhibitor peptides.
Here we describe the successful application of our approach to the optimization of peptides that inhibit the association between Ras and its downstream target protein Raf. Ras acts as a molecular switch between its active state where guananine triphosphate (GTP) is bound and its inactive form where a guananine diphosphate (GDP) is bound (Barbacid, 1987). Only the GTP-bound form of Ras binds to the downstream targets, one of which is the threonine/serine kinase Raf (Barbacid, 1987
). Once Ras binds to Raf, it stimulates the activation of the MEKMAPKK signaling pathway, transmitting signaling events into the nucleus to control cell growth and differentiation. An oncogenic form of Ras exists only in the active form, resulting in a constitutive activation of Ras-mediating signaling events and promoting aberrant growth in more than 30% of all human tumors (Barbacid, 1987
). Antagonists of the RasRaf interactions that are likely to inhibit the Ras-stimulated signal transduction pathway are thus of great potential value to anti-cancer therapy (Maruta and Burgess, 1994
).
RasRaf interactions have been studied extensively by site-directed mutagenesis [see Campbell et al. (1998) and references therein]. The important regions of Raf that bind to Ras protein have been identified as a 78 amino acid N-terminal region, the so-called `Ras binding domain' (RBD) (Chuang et al., 1994; Scheffler et al., 1994
) and a zinc-finger domain, the so-called `cysteine-rich domain' (CRD) (Drugan et al., 1996
). Whereas the binding mechanism between the CRD and Ras remains unknown, details of the interaction between RBD and Ras have been revealed from the crystal structure of the complex between the RBD and Rap1A (Nassar et al., 1995
), a protein with a >80% sequence homology to Ras. Previously, we have carried out extensive computational and experimental studies on the structural and binding properties of the RasRaf interactions (Zeng et al., 1998
, 1999a
,Zeng et al., b
). Based on the available structural information (Zeng et al., 1999b
), potential peptide inhibitors have been designed using the protocol we have described previously (Zeng and Treutlein, 1999
). Here we describe the application of our previous results to optimize a previously published inhibitor and its experimental characterization by in vitro radioassay, enzyme-linked immunosorbent assay (ELISA) and BIAcore techniques. Our optimized peptides should act as lead compounds for the development of novel pharmaceuticals (Bures and Martin, 1998
; Schneider et al., 1999
).
![]() |
Materials and methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The details of our computational combinatorial inhibitor design approach have been described elsewhere (Zeng and Treutlein, 1999; Zeng, 2000
). Our procedure consists of three steps. In the first step, the MCSS method is used to determine energetically favorable positions and orientations of functional groups on the surface of the complex between Ras and a Ras-binding helix (RBH). The RBH helix (residues 8489) contains two critical residues (Lys84 and Arg89) for RasRaf interaction, as mutants at these positions have been shown to abolish the RasRaf binding in vivo completely (Nassar et al., 1996
; Zeng et al., 1999a
). The structure of the RasRBD complex was taken from our 2 ns molecular dynamics simulation (Zeng et al., 1999b
) and shown in Figure 1
. In the complex, residues Lys84 and Arg89 of RBH form direct salt bridges to Asp33 and Asp38 of the Ras effector loop (residues 3040), respectively, consistent with mutagenesis studies (Nassar et al., 1996
; Zeng et al., 1999a
). Our basic design strategy was to extend the RBH in order to allow interactions with all the important regions such as the effector loop, switch II (residues 6076) (Wittinghofer and Nassar, 1996
) and an `interaction surface' (residues 58, 5256) that was identified previously (Zeng and Treutlein, 1999
). Inclusion of the RBH into our design scheme significantly reduces the number of peptide conformations that need to be considered on the surface of Ras. The functional groups chosen for the MCSS calculations were N-methylacetamide (NMA), benzene, propane, phenol, methanol, acetate ion, methylammonium, methylguanidinium and water, representing functional moieties of the peptide backbone, hydrophobic residues, polar uncharged residues, charged residues and solvent. All parameters for the functional groups were taken from the CHARMM22 all-hydrogen atom force field (MacKerell et al., 1998
). Amongst the MCSS minima on the surface of Ras, only those within a 4.0 Å distance from Ras were selected for further analysis.
|
In the third step, subsets of functional groups were selected for each C position of the polyglycine peptide according to the distance from each replica of functional groups to the C
atom of the polyglycine representative. The amino acid type was randomly chosen from the available subsets and the missing side chain atoms were grown to connect to the C
atoms by inserting CH2 groups. All completely built side chains were energy minimized to insure feasible side chain geometry. The sequences of functional groups were thus transformed into a peptide. All the procedures were performed using the program XPLOR (Brunger, 1992
).
The peptides were then connected to the RBH and those with significant geometric distortion (i.e. large bond and angle energies) were ignored. Our procedure resulted in 100 peptides on the surface of Ras. We aligned the sequences of these peptides using the program CLUSTALX (Thompson et al., 1994). The multiple sequence alignment was used to determine the probability of an amino acid type at each position of the peptide. The optimal sequences of the Ras binding peptides were constructed by placing amino acid types with high probabilities at each position of the peptide.
Inclusion of RBH for construction of peptide inhibitors
The designed peptides that are expected to inhibit RasRaf interaction contains both RBH and the peptides designed using our computational combinatorial inhibitor design approach. However, the RBH itself hardly affects RasRaf binding (Barnard et al., 1995, 1998
) as it might not adopt its wild-type helical conformation. Therefore, our design strategy for RasRaf inhibitor is to identify the lead peptide using first our design method and then to include specific residues of the RBH which directly form salt bridges or hydrogen bonds in the crystal structure of the RasRBD complex (Nassar et al., 1995
; Zeng et al., 1999b
). The resulting inhibitors can therefore maintain the important binding determinants between the RBH and the effector loop.
Reversed sequences of designed peptides
Our method described above primarily identifies favorite interaction sites of amino acids. Those sites are then linked together to form a peptide. The locations of the N- and C-termini are initially determined by our design strategy to extend the RBH towards the `interaction surface': the N-terminus would be located on an area that coincides with the RBH binding site, the C-terminus would be located on the `interaction surface'. However, the linkage algorithm which links the amino acids is not precise enough to exclude the possibility that a peptide with its amino acid sequence reversed would bind to the surface patch with significantly higher affinity, although it would bind in a reversed way on to the surface, i.e. with its C-terminus at the position of the RBH and its N-terminus at the `interaction surface'. This imprecision of our algorithm makes it necessary also to take into account reversed peptide sequences and sequence patterns when searching for an optimized inhibitor peptide.
Expression and preparation of Ras, GSTRBDRaf and designed peptides
V-Ha-Ras was purified from Escherichia coli according to the procedure of Gibbs et al. (1984). The concentration of Ras was determined with the Bio-Rad protein assay. The RBD was produced in E.coli as a GST fusion protein, which was affinity-purified on GSH beads (Fridman et al., 1994). The concentrations of purified GSTRBD constructs were determined by scanning Coomassie Brilliant Blue-stained SDSPAGE gels using densitometry and measuring the band densities relative to BSA standards with ImageQuant 3.3. The peptides were purchased from Chiron Technology (Melbourne, Australia) purified to apparent homogeneity by micro- preparation RP-HPLC using a Brownlee RP3000 (300x 201 mm i.d.) column eluted at a flow-rate of 100 µl/min with a linear 60 min gradient between a primary solvent of 0.15% (v/v) TFA and a secondary solvent of 60% aqueous CH3CN0.125% (v/v) TFA. The column temperature was 45°C and detection was effected by measuring the absorbance at 215 nm.
Enzyme-linked immunosorbent assay (ELISA)
Purified RasGTP complex (100 µl, 10 µg/ml) was plated overnight on to 96-well microtiter plates (Nunc) coated with a buffer solution containing 10 mM Na2CO3 and 30 mM NaHCO3 (pH 9.6). The plates were washed twice with a `blocking' buffer solution of 0.05% Tween, 20 mM Tris, 2 mM MgCl2, 150 mM NaCl, 1 mM DTT and 0.5% BSA. The wells were then blocked in the blocking buffer at 37°C for 1 h. GSTRBDRaf (10 µg/ml) was diluted into PBS containing 0.05% Tween, 20 mM Tris, 2 mM MgCl2, 150 mM NaCl and 1 mM DTT. The peptides (100 µM) were incubated with GSTRBDRaf (10 µg/mol) at 37°C for 15 min prior to addition to the wells. The plates were washed three times with PBS, then incubated at 37°C for 2 h with a goat anti-GST antibody (Amersham Pharmacia Biotech) diluted 2000-fold in PBS. The washing step was repeated and the plates were incubated with a 1:1000 dilution of rabbit anti-goat antibody conjugated to horseradish peroxidase (HRP) (Amersham Pharmacia Biotech) for 1 h. The wells were washed three times, followed by development with ABST substrate (ZYMED, San Francisco, CA) diluted 100-fold in a buffer solution containing 0.1 M sodium citrate buffer and 0.03% hydrogen peroxide. Optical densities (416 nm) were read after 15 min in a microtiter plate reader. The amount bound was quantified by comparison with an appropriate calibration curve (3010 µg/ml).
Radioassay for RasGTP/RBD binding
Binding between Ras and RBDGST was carried out as described previously (Fridman et al., 1994; Zeng et al., 1999a
). Briefly, v-Ha-Ras was labelled with [
32P]GTP (Du Pont). A 2.5 µg amount of the radioactively labelled Ras (100 µl) was incubated with GSTRBD bound to GSH beads in presence of 100 µM of peptides for 15 min at room temperature on a rotator. The mixture was washed four times with 500 µl of the binding buffer [50 mM TrisHCl (pH 7.5), 1 mM DTT, 10 mM MgCl2, 0.5 mg/ml BSA and 1 mM ATP]. The GSH beads were recovered by centrifugation at 104g for 5min. The radioactivity bound to the beads in the pellet was measured with a scintillation counter. Binding was expressed as the percentage of the radioactivity associated with RBD. Non-specific binding of Ras to the beads was determined by incubation with GST alone.
Biosensor experiments
The binding of GSTRas to the peptides was also analyzed on the BIAcore optical biosensor (Pharmacia Biosensor, Sweden) using peptide-derivatized CM 5 sensor chips. The immobilization of peptides on the sensor chip surface was performed essentially as described previously (Catimel et al., 1999). The peptides at a concentration of 150 µg/ml in 20 mM sodium acetate (pH 4.5) was coupled at 1 µl/min on to sensor chips that activated by a pulse of 0.2 M N-ethyl-N'-(dimethylaminopropyl)carbodiimide and 0.05 M N-hydroxysuccinimide (45 µl, 1 µl/min), in order to yield an increase in the response level of 8001500 response units.
Known Ras/Raf inhibitors
To date, only three inhibitors of RasRaf interaction have been identified. However, even these bind to Ras with only low affinity and are unsuitable for practical use. While sulindac sulfide has been shown to inhibit RasRaf interaction with an IC50 of 400 µM (Hermann et al., 1998), two peptide inhibitors identified from the Ras effectors (i.e. GAP, c-raf-1, etc.) indicate only ~20% inhibition at a concentration of 100 µM (Clark et al., 1996). More recently, a peptide sequence CCAVFRL derived from the Ras-binding domain of c-raf-1 has been shown to interfere with RasRaf association (Barnard et al., 1998
). This peptide was determined by a random screening approach without any consideration of the RasRaf complex structure. In fact, this short peptide is located ~15 Å away from the binding interface in the crystal structure of the RasRBD complex (Nassar et al., 1995
; Zeng et al., 1999b
), as shown in Figure 1
. As we will demonstrate later, our approach can identify this peptide and its derivatives, providing a way towards de novo design of inhibitor that cannot be revealed from the complex structure.
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
In our approach, the RBH was intended to interact with the effector loop of Ras protein as observed in the X-ray structure of the RBDRap1A complex (Nassar et al., 1995). The extended part of the inhibitor was designed to bind to the `interaction surface' and the switch II region. Based on the minimum energy positions of functional groups obtained from the MCSS calculations (Zeng and Treutlein, 1999
), 300 sequences were generated using the protocol described in Materials and methods. After growing the side chains and merging to the Ras-binding helix (RBH), 196 peptides with significant bond distortion (i.e. bond stretch energy >30 kcal/mol) were excluded. The remaining 104 peptides were kept for further investigation in the study. Sequence alignment using the CLUSTALX program shows distinct regions for RBH and extended regions, consistent with our design strategy. The probabilities of different amino acids at each position of the 104 peptides were derived from the alignment of the peptide sequences. Two additional linker residues, namely glycine and phenylalanine, were inserted between the RBH and extended region such that the bulky side chain of phenylalanine could prevent cross-interaction between the regions. Figure 2
shows the probabilities of amino acids at each position of the extended part. Only the amino acids with probabilities
0.30 were considered for the design of the inhibitor and labeled in the figure. While position P1 is clearly defined by Tyr by its hydrophobic interaction to the aliphatic part of RasArg41, position P2 is dominated by Asp/Glu owing to its internal electrostatic interaction with Arg89 and Lys84 of RBH. The residues at positions P3, P4 and P5 interact with Glu37, Asp38, Ser39 and Tyr40 within the Ras effector loop and consequently are controlled by amino acids Arg/Lys or Phe, Arg/Lys and Ser/Thr, respectively. Positions P6P8 are located on the `interaction surface', demanding hydrophobic Phe or hydrophobic/polar residue Tyr to interact residues Leu56, Lys5 and Val7 of Ras, respectively. The last residue (position P9) is located between Lys5 of the interaction surface and Thr74 of the switch II, resulting in uniform distribution amongst Phe, Tyr and Asp/Glu. This position is therefore not clearly defined and can be occupied by any amino acid, labeled X in Figure 2
. By joining the RBH and extended parts with two linker residues (glycine and phenylalanine), the motif of inhibitor sequences can thus be defined as KALKVRGFYDRRTFFFX.
|
|
Table I lists the peptide inhibitors derived from our designed peptides. We first identify the lead peptide based on the extended part of the designed inhibitor sequence. Three modifications are introduced: first, we replace the first two residues (Y and D) by a Leu as these residues mainly interact with the RBH; second, we attach two cysteines to the C-terminus, considering their importance demonstrated previously (Barnard et al., 1998
); third, we substitute the hydrophobic `cluster' XFFF into AVFL in order to dilute the peptide in water. The resulting lead peptide has a sequence of CCAVFLTRRL (peptide 1), altering RasRaf interaction by 23.4% from ELISA and 12.5% from radioassay.
|
The inhibitors are further optimized by extending the lead peptides (i.e. peptides 1 and 4) to inclusion of specific residues of RBH. These residues are defined in Materials and methods and highlighted in bold in Table I. In addition, a fragment of DYFG between the RBH and extended parts is ignored because it mainly interacts with the RBH and precludes the cross-interaction between the RBH and extended regions. Extending peptide 4, peptide 5 increases inhibition up to 39.0% from ELISA, but decreases it to 8.3% from radioassay. Although peptides 14 have similar magnitudes of inhibition from radioassay, peptide 1 inhibits RasRaf binding as twice as much of peptide 4 from the ELISA assay. By elongating peptide 1 to the RBH, however, peptide 6 decreases the inhibition to 10.4% from ELISA and to 4.2% from radioassay.
Binding of GSTRas protein to sensor chip-immobilized peptides
As both ELISA and radioassay experiments showed that peptides 1 and 3 inhibit RasRaf interaction, BIAcore analyses were also carried out to detect the direct binding of GSTRas protein to the peptides. The peptides were sensor chip-immobilized at levels of ~800 and ~1500 response units, respectively. Figure 4 shows the binding curves for these interactions. The relative response units clearly indicate that both peptides bind to Ras, in agreement with ELISA and radioassay results (see Table I
). However, the biosensor response for peptide 1 is shown to be significantly greater than the response of peptide 4, consistent with the ELISA experiment where peptide 1 inhibits RasRaf binding twice as much as peptide 4.
|
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
In total, six peptides were evaluated using ELISA and radioassay, as well as in BIAcore assay. The sequences correspond to the RBH and extended components in the designed inhibitors. The peptides from the extended part were shown to inhibit RasRaf interaction up to 12.015.0 and 8.923.4% from radioassay and ELISA, respectively. BIAcore analysis suggested that the inhibition mechanism is mainly due to the direct binding of peptides to Ras. Relative biosensor responses can be correlated with data from competition experiments (ELISA). Extension of these peptides to include the RBH improved the inhibition up to 39.0% in ELISA, but reduced the inhibition to 8% in radioassay. This difference might be due to different assay protocols and conditions used in ELISA and radioassay (see Materials and methods). As indicated from the results for peptides 14 in Table I, the magnitude of inhibition obtained from the ELISA experiment is more sensitive to the condition than those from radioassay. Regardless of the uncertainty in the assays, our results demonstrate that relatively small molecules can disrupt the formation of the RasRaf signalling complex in vitro.
It is also interesting that extension of peptides by charged residues, i.e. peptide 1 (sequence CCAVFLTRL) vs peptide 6 (sequence CCAVFLTRRLRA), indeed decreases the inhibition in both ELISA and radioassay. In principle, the difference of binding affinities between two peptides is the difference of the free energy change of two peptides in the complexes with the target protein and the desolvation free energies of the two peptides in solution (Zeng et al., 1999a). When charged residues are properly introduced into a peptide to interact strongly with the target protein, it could result in a larger desolvation free energy in solution. If the desolvation free energy decreases more than the binding free energy, the charged residues could indeed decrease the binding affinity of the initial peptide. It seems likely that a proper increase in hydrophobicity should be considered when charged residues are inserted/substituted within the peptides, in order to maintain the desolvation free energy of the peptides in solution.
Previous studies have demonstrated that the multiple conformations in the effector loop and switch II regions of Ras are important for target recognition (i.e. RafRBD) (Sydor et al., 1998; Terada et al., 1999
), so that the dynamics of Ras and peptides should be an important aspect in the design of Ras inhibitors. Therefore, molecular dynamics simulation on the complex between Ras and each designed peptide solvated in the explicit solvent molecules and considering the protein flexibility in the design algorithm (Stultz and Karplus, 1999
) are expected to improve the design results.
From the results of the bioassays we conclude that our computational combinatorial design approach combined with experimental measurements (ELISA or radioassay) as a guide in sequence pattern identification and optimization of Ras inhibitors shows promising results. With the rapid development of computer technology, our method could be used for the economical design of combinatorial libraries, which could be screened for more optimal inhibitor peptides.
![]() |
Notes |
---|
5 Present address: Uppsala School of Engineering, Uppsala University, Uppsala, Sweden
![]() |
Acknowledgments |
---|
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Barbacid,M. (1987) Annu. Rev. Biochem., 56, 779827.[ISI][Medline]
Barnard,D., Diaz,B., Hettich,L., Chuang,E., Zhang,X.F., Avruch,K. and Marshall,M. (1995) Oncogene, 10, 1283[ISI][Medline]
Barnard,D., Sun,H., Baker,L. and Marshall,M. (1998) Biochem. Biophys. Res. Commun., 247, 176180.[ISI][Medline]
Bohacek,R. and McMartin,C. (1997) Curr. Opin. Chem. Biol., 1, 157161.[ISI][Medline]
Brunger,A.T. (1992) X-PLOR Version 3.1, A System for X-Ray Crystallography and NMR. Yale University Press, New Haven, CT.
Bures,M. and Martin,Y. (1998) Curr. Opin. Chem. Biol., 2, 376380.[ISI][Medline]
Caflisch,A. (1996) J. Comput.-Aided Mol. Des. 10, 372396.
Caflisch,A. and Karplus,M. (1995) Perspect. Drug. Discov. Des., 3, 51.[ISI]
Caflisch,A. and Karplus,M. (1996) J. Comput.-Aided Mol. Des., 10, 372396.
Caflisch,A., Miranker,A. and Karplus,M. (1993) J. Med. Chem., 36, 2142.[ISI][Medline]
Campbell,S., Khosravi-Far,R., Rossman,K., Clark,G. and Der,C. (1998) Oncogene, 17, 13951413.[ISI][Medline]
Catimel,B., Domagala,T., Nerrie,M., Weinstock,J., White,S., Abud,H., Heath,J. and Nice,E. (1999) Protein Pept. Lett., 6, 319340.[ISI]
Chuang,E., Barnard,D., Hettich,L., Zhang,X.-F., Avruch,J. and Marshall,M.S. (1994) Mol. Cell. Biol., 14, 53185325.[Abstract]
Clark,G., Drugan,J., Terrell,R., Bradham,C., Der,C., Bell,R. and Campbell,S. (1996) Proc. Natl Acad. Sci. USA, 93, 15771581.
Drugan J., Khosravi-Far,R., White,M., Der,C., Sung,Y., Hwang,Y. and Campbell,S. (1996) J. Biol. Chem., 271, 233.
Fridman,M., Tikoo,A., Varga,M. Murphy,A., Nur-El-Kamal,M. and Maruta,H. (1994) J. Biol. Chem., 269, 3010530108.
Gibbs,J., Sigal,I., Poe,M. and Scolnic,E. (1984) Proc. Natl Acad. Sci. USA, 81, 57045708.[Abstract]
Greer,J., Erickson,J., Baldwin,J. and Varney,M. (1994) J. Med. Chem., 37, 1035.[ISI][Medline]
Herrmann,C., Block,C., Geisen,C., Haas,K., Weber,C., Winde,G., Moroy,T. and Muller,O. (1998) Oncogene, 17, 17691776.[ISI][Medline]
Hilpert,K. et al. (1994) J. Med. Chem., 37, 3889.[ISI][Medline]
Josephy-McCarthy,D., Hogle,J. and Karplus,M. (1997) Proteins: Struct. Funct. Genet., 29, 3258.[ISI][Medline]
MacKerell,A. et al. (1998) J. Phys. Chem. B, 102, 35863616.[ISI]
Maruta,H. and Burgess,A. (1994) Bioessays, 16, 489496.[ISI][Medline]
Nassar,N., Horn,G., Herrmann,C., Scherer,A., McCormick,F. and Wittinghofer,A. (1995) Nature, 375, 554560.[ISI][Medline]
Nassar,N., Horn,G., Herrmann,C., Block,C., Janknecht,R. and Wittinghofer,A. (1996) Nature Struct. Biol., 3, 723729.[ISI][Medline]
Scheffler,J. et al. (1994) J. Biol. Chem., 269, 2234022346.
Schneider,G., Schrodl,W., Gerd,W., Muller,J., Nissen,E., Ronspeck,W., Wrede,P. and Kunze,R. (1999) Proc. Natl Acad. Sci. USA, 95, 1217912184.
Stultz,C. and Karplus,M. (1999) Proteins, 37, 512529.[ISI][Medline]
Sydor,J., Engelhard,M., Wittinghofer,A., Goody,R. and Herrmann C. (1998) Biochemistry, 37, 1429214299.[ISI][Medline]
Terada,T. et al. (1999) J. Mol. Biol., 286, 219232.[ISI][Medline]
Thompson,J., Higgins,D. and Gibson,T. (1994) Nucleic Acids Res., 22, 46734680.[Abstract]
Varghese,J., Smith,P., Sollis,S., Blick,T., Sahasrabudhe,A., McKimm-Breschkin,J. and Colman,P. (1998) Structure, 6, 735746.[ISI][Medline]
Von Itzstein,M., Dyson,J., Oliver,S., White,H., Wu,W., Kok,G. and Pegg,M. (1996) J. Med. Chem., 39, 388391.[ISI][Medline]
Wittinghofer,A. and Nassar,N. (1996) Trends Biochem. Sci., 21, 488491.[ISI][Medline]
Zeng,J. (2000) Comb. Chem. High Throughput Screen., 3, 355363.[ISI][Medline]
Zeng,J. and Treutlein,H.R. (1999) Protein Eng., 12, 457468.
Zeng,J., Treutlein,H.R. and Simonson,T. (1998) Proteins: Struct. Funct. Genet., 31, 186200.[Medline]
Zeng,J., Fridman,M., Maruta,H., Treutlein,H.R. and Simonson,T. (1999a) Protein Sci., 8, 5061.[Abstract]
Zeng,J., Treutlein,H.R. and Simonson,T. (1999b) Proteins: Struct. Funct. Genet., 35, 89100.[ISI][Medline]
Received October 3, 2000; revised October 30, 2000; accepted October 31, 2000.