Protein cofactor-dependent acquisition of novel catalytic activity by the RNase P ribonucleoprotein of E. coli

Escherichia coli RNase P derivatives were evolved in vitro for DNA cleavage activity. Ribonucleoproteins sampled after ten generations of selection show a >400-fold increase in the first-order rate constant (k(cat)) on a DNA substrate, reflecting a significant improvement in the chemical cleavage step. This increase is offset by a reduction in substrate binding, as measured by K(M). We trace the catalytic enhancement to two ubiquitous A-->U sequence changes at positions 136 and 333 in the M1 RNA component, positions that are phylogenetically conserved in the Eubacteria. Furthermore, although the mutations are located in different folding domains of the catalytic RNA, the first in the substrate binding domain, the second near the catalytic core, their effect on catalytic activity is significantly influenced by the presence of the C5 protein. The activity of the evolved ribonucleoproteins on both pre-4.5 S RNA and on an RNA oligo substrate remain at wild-type levels. In contrast, improved DNA cleavage activity is accompanied by a 500-fold decrease in pre-tRNA cleavage efficiency (k(cat)/K(M)). The presence of the C5 component does not buffer this tradeoff in catalytic activities, despite the in vivo role played by the C5 protein in enhancing the substrate versatility of RNase P. The change at position 136, located in the J11/12 single-stranded region, likely alters the geometry of the pre-tRNA-binding cleft and may provide a functional explanation for the observed tradeoff. These results thus shed light both on structure/function relations in E. coli RNase P and on the crucial role of proteins in enhancing the catalytic repertoire of RNA.     

Substitution of the 2'-hydroxyl group at position 2.1 by an amino group interferes with Mg(2+) binding and efficient cleavage by hammerhead ribozyme.

Recently we have demonstrated that hammerhead ribozymes can be fully substituted with 2'-amino pyrimidines without detriment to the catalytic activity, provided that positions 2.2 and/or 2.1 are not modified. We now report on the potential molecular mechanisms by which 2'-amino groups at these positions inhibit the ribozyme cleavage activity. In the presence of Mg(2+), the 2'-amino modification at positions 2.2 and/or 2.1 had no significant effect on substrate binding. Detailed analysis of the ribozyme initial cleavage rates in the presence of various Mg(2+) concentrations indicated that Mg(2+) binding is inhibited by the 2'-amino group at position 2.1. Furthermore, preannealed substrate molecules to the modified ribozyme are not effectively cleaved upon Mg(2+) addition, indicating an alteration of the ribozyme cleavage step. Surprisingly, the cleavage activity of the modified ribozymes was substantially increased when Mg(2+) ions were replaced by the thiophilic Mn(2+) ions, whereas only a moderate cleavage enhancement occurred with its unmodified version. Taken together, our findings indicate that changes in the sugar at position 2.1 alter Mg(2+)-promoting ribozyme cleavage.     

Site-directed cross-linking of mRNA analogues to 16S ribosomal RNA; a complete scan of cross-links from all positions between '+1' and '+16' on the mRNA, downstream from the decoding site.

mRNA analogues containing 4-thiouridine residues at selected sites were used to extend our analysis of photo-induced cross-links between mRNA and 16S RNA to cover the entire downstream range between positions +1 and +16 on the mRNA (position +1 is the 5'-base of the P-site codon). No tRNA-dependent cross-links were observed from positions +1, +2, +3 or +5. Position +4 on the mRNA was cross-linked in a tRNA-dependent manner to 16S RNA at a site between nucleotides ca 1402-1415 (most probably to the modified residue 1402), and this was absolutely specific for the +4 position. Similarly, the previously observed cross-link to nucleotide 1052 was absolutely specific for the +6 position. The previously observed cross-links from +7 to nucleotide 1395 and from +11 to 532 were however seen to a lesser extent with certain types of mRNA sequence from neighbouring positions (+6 to +10, and +10 to +13, respectively); no tRNA-dependent cross-links to other sites on 16S RNA were found from these positions, and no cross-linking was seen from positions +14 to +16. In each case the effect of a second cognate tRNA (at the ribosomal A-site) on the level of cross-linking was studied, and the specificity of each cross-link was confirmed by translocation experiments with elongation factor G, using appropriate mRNA analogues.     

NAIM and site-specific functional group modification analysis of RNase P RNA: magnesium dependent structure within the conserved P1-P4 multihelix junction contributes to catalysis

The tRNA processing endonuclease ribonuclease P contains an essential and highly conserved RNA molecule (RNase P RNA) that is the catalytic subunit of the enzyme. To identify and characterize functional groups involved in RNase P RNA catalysis, we applied self-cleaving ribozyme-substrate conjugates, on the basis of the RNase P RNA from Escherichia coli, in nucleotide analogue interference mapping (NAIM) and site-specific modification experiments. At high monovalent ion concentrations (3 M) that facilitate protein-independent substrate binding, we find that the ribozyme is largely insensitive to analogue substitution and that concentrations of Mg2+ (1.25 mM) well below that necessary for optimal catalytic rate (>100 mM) are required to produce interference effects because of modification of nucleotide bases. An examination of the pH dependence of the reaction rate at 1.25 mM Mg2+ indicates that the increased sensitivity to analogue interference is not due to a change in the rate-limiting step. The nucleotide positions detected by NAIM under these conditions are located exclusively in the catalytic domain, consistent with the proposed global structure of the ribozyme, and predominantly occur within the highly conserved P1-P4 multihelix junction. Several sensitive positions in J3/4 and J2/4 are proximal to a previously identified site of divalent metal ion binding in the P1-P4 element. Kinetic analysis of ribozymes with site-specific N7-deazaadenosine and deazaguanosine modifications in J3/4 was, in general, consistent with the interference results and also permitted the analysis of sites not accessible by NAIM. These results show that, in this region only, modification of the N7 positions of A62, A65, and A66 resulted in measurable effects on reaction rate and modification at each position displayed distinct sensitivities to Mg2+ concentration. These results reveal a restricted subset of individual functional groups within the catalytic domain that are particularly important for substrate cleavage and demonstrate a close association between catalytic function and metal ion-dependent structure in the highly conserved P1-P4 multihelix junction.     

Improved thrombin binding aptamer by incorporation of a single unlocked nucleic acid monomer

A 15-mer DNA aptamer (named TBA) adopts a G-quadruplex structure that strongly inhibits fibrin-clot formation by binding to thrombin. We have performed thermodynamic analysis, binding affinity and biological activity studies of TBA variants modified by unlocked nucleic acid (UNA) monomers. UNA-U placed in position U3, U7 or U12 increases the thermodynamic stability of TBA by 0.15–0.50 kcal/mol. In contrast, modification of any position within the two G-quartet structural elements is unfavorable for quadruplex formation. The intramolecular folding of the quadruplexes is confirmed by T_m versus ln c analysis. Moreover, circular dichroism and thermal difference spectra of the modified TBAs displaying high thermodynamic stability show bands that are characteristic for antiparallel quadruplex formation. Surface plasmon resonance studies of the binding of the UNA-modified TBAs to thrombin show that a UNA monomer is allowed in many positions of the aptamer without significantly changing the thrombin-binding properties. The biological effect of a selection of the modified aptamers was tested by a thrombin time assay and showed that most of the UNA-modified TBAs possess anticoagulant properties, and that the construct with a UNA-U monomer in position 7 is a highly potent inhibitor of fibrin-clot formation.     

Thermodynamic and biological evaluation of a thrombin binding aptamer modified with several unlocked nucleic acid (UNA) monomers and a 2'-C-piperazino-UNA monomer

Thrombin binding aptamer is a DNA 15-mer which forms a G-quadruplex structure and possess promising anticoagulant properties due to specific interactions with thrombin. Herein we present the influence of a single 2'-C-piperazino-UNA residue and UNA residues incorporated in several positions on thermodynamics, kinetics and biological properties of the aptamer. 2'-C-Piperazino-UNA is characterized by more efficient stabilization of quadruplex structure in comparison to regular UNA and increases thermodynamic stability of TBA by 0.28-0.44 kcal/mol in a position depending manner with retained quadruplex topology and molecularity. The presence of UNA-U in positions U3, U7, and U12 results in the highest stabilization of G-quadruplex structure (ΔΔG(37)(°)=-1.03kcal/mol). On the contrary, the largest destabilization mounting to 1.79 kcal/mol was observed when UNA residues were placed in positions U7, G8, and U9. Kinetic studies indicate no strict correlation between thermodynamic stability of modified variants and their binding affinity to thrombin. Most of the studied variants bind thrombin, albeit with decreased affinity in reference to unmodified TBA. Thrombin time assay studies indicate three variants as being as potent as TBA in fibrin clotting inhibition.     

Intracellular ribozyme-catalyzed trans-cleavage of RNA monitored by fluorescence resonance energy transfer

Small catalytic RNAs like the hairpin ribozyme are proving to be useful intracellular tools; however, most attempts to demonstrate trans-cleavage of RNA by ribozymes in cells have been frustrated by rapid cellular degradation of the cleavage products. Here, we describe a fluorescence resonance energy transfer (FRET) assay that directly monitors cleavage of target RNA in tissue-culture cells. An oligoribonucleotide substrate was modified to inhibit cellular ribonuclease degradation without interfering with ribozyme cleavage, and donor (fluorescein) and acceptor (tetramethylrhodamine) fluorophores were introduced at positions flanking the cleavage site. In simple buffers, the intact substrate produces a strong FRET signal that is lost upon cleavage, resulting in a red-to-green shift in dominant fluorescence emission. Hairpin ribozyme and fluorescent substrate were microinjected into murine fibroblasts under conditions in which substrate cleavage can occur only inside the cell. A strong FRET signal was observed by fluorescence microscopy when substrate was injected, but rapid decay of the FRET signal occurred when an active, cognate ribozyme was introduced with the substrate. No acceleration in cleavage rates was observed in control experiments utilizing a noncleavable substrate, inactive ribozyme, or an active ribozyme with altered substrate specificity. Subsequently, the fluorescent substrates were injected into clonal cell lines that expressed cognate or noncognate ribozymes. A decrease in FRET signal was observed only when substrate was microinjected into cells expressing its cognate ribozyme. These results demonstrate trans-cleavage of RNA within mammalian cells, and provide an experimental basis for quantitative analysis of ribozyme activity and specificity within the cell.     

Binary hammerhead ribozymes with improved catalytic activity

A new design of binary hammerhead ribozymes displaying high catalytic activity and nucleolytic stability is described. These catalytic structures consist of two partially complementary oligoribonucleotides, capable of assembling into the hammerhead-like structure without tetraloop II on binding to the RNA target. A series of these binary ribozymes targeting the translation initiation region of multiple drug resistance gene mdr1 mRNA was synthesized and assessed in terms of catalytic activity under single and multiple reaction turnover conditions. Enhanced nuclease resistance of the binary ribozymes was achieved by incorporation of 2'-modified nucleotides at selected positions, along with addition of a 3'-3'-linked thymidine cap. The new binary ribozymes exhibit higher RNA cleavage activity than their full-length analogs because of faster dissociation of cleavage products. Furthermore, an excess of one of the ribozyme strands provides the possibility to unfold structured regions of the target RNA and facilitate productive complex formation.     

Structural and functional analysis of methylation and 5'-RNA sequence requirements of short capped RNAs by the methyltransferase domain of dengue virus NS5

The N-terminal 33 kDa domain of non-structural protein 5 (NS5) of dengue virus (DV), named NS5MTase(DV), is involved in two of four steps required for the formation of the viral mRNA cap (7Me)GpppA(2'OMe), the guanine-N7 and the adenosine-2'O methylation. Its S-adenosyl-l-methionine (AdoMet) dependent 2'O-methyltransferase (MTase) activity has been shown on capped (7Me+/-)GpppAC(n) RNAs. Here we report structural and binding studies using cap analogues and capped RNAs. We have solved five crystal structures at 1.8 A to 2.8 A resolution of NS5MTase(DV) in complex with cap analogues and the co-product of methylation S-adenosyl-l-homocysteine (AdoHcy). The cap analogues can adopt several conformations. The guanosine moiety of all cap analogues occupies a GTP-binding site identified earlier, indicating that GTP and cap share the same binding site. Accordingly, we show that binding of (7Me)GpppAC(4) and (7Me)GpppAC(5) RNAs is inhibited in the presence of GTP, (7Me)GTP and (7Me)GpppA but not by ATP. This particular position of the cap is in accordance with the 2'O-methylation step. A model was generated of a ternary 2'O-methylation complex of NS5MTase(DV), (7Me)GpppA and AdoMet. RNA-binding increased when (7Me+/-)GpppAGC(n-1) starting with the consensus sequence GpppAG, was used instead of (7Me+/-)GpppAC(n). In the NS5MTase(DV)-GpppA complex the cap analogue adopts a folded, stacked conformation uniquely possible when adenine is the first transcribed nucleotide at the 5' end of nascent RNA, as it is the case in all flaviviruses. This conformation cannot be a functional intermediate of methylation, since both the guanine-N7 and adenosine-2'O positions are too far away from AdoMet. We hypothesize that this conformation mimics the reaction product of a yet-to-be-demonstrated guanylyltransferase activity. A putative Flavivirus RNA capping pathway is proposed combining the different steps where the NS5MTase domain is involved.     

3'-5' Exonucleolytic activity of DNA polymerases: structural features that allow kinetic discrimination between ribo- and deoxyribonucleotide residues

We have determined rates for the excision of nucleotides from the 3' termini of chimeric DNA-RNA oligonucleotides using the Klenow fragment (KF) and two other DNA polymerases, from phages T4 and T7. For these studies, we synthesized DNA-RNA chimeric oligonucleotides with RNA residues in defined positions. When a ribonucleotide residue was placed at the 3' terminus, all three DNA polymerases removed it at the same rate as they did for substrates composed solely of deoxynucleotide residues. There was a decrease in the excision rate, however, when a ribonucleotide residue was located at the second or third position from the 3' terminus. When both the second and third positions were occupied by ribonucleotide residues, the excision rate for the 3' terminal nucleotide was reduced even further and was almost identical to the rate observed when the DNA polymerases encountered single-stranded RNA. The magnitude of the effect of ribonucleotide residues on the excision rate was lower when Mn(2+) replaced Mg(2+) as the essential divalent cation. Two KF mutations, Y423A and N420A, selectively affected the excision rates for the chimeric substrates. Specifically, Y423A totally abolished the rate reduction when there was a single ribonucleotide residue immediately preceding the 3' terminus, whereas N420A diminished, but did not eliminate, the rate reduction relative to that of wild-type KF when the single ribonucleotide residue occupied either the second or third position from the 3' terminus. These results are consistent with the structure of a KF-ss DNA complex from which it can be deduced, by modeling, that a 2' OH group on the second sugar from the 3' terminus would sterically clash with the Tyr 423 side chain, and a 2' OH group on the third sugar would clash with the side chain of Asn 420. The corresponding mutations in T4 DNA polymerase did not affect the rate of hydrolysis of the chimeric oligonucleotides. Thus, there appears to be a major difference in the kinetic behavior of KF and T4 DNA polymerase with respect to the exonuclease reaction. These results are discussed with respect to their possible biological relevance to DNA replication.     

Template nucleotide moieties required for de novo initiation of RNA synthesis by a recombinant viral RNA-dependent RNA polymerase

The recombinant RNA-dependent RNA polymerase of the bovine viral diarrhea virus specifically requires a cytidylate at the 3' end for the de novo initiation of RNA synthesis (C. C. Kao, A. M. Del Vecchio, and W. Zhong, Virology 253:1-7, 1999). Using RNAs containing nucleotide analogs, we found that the N3 and C4-amino group at the initiation cytidine were required for RNA synthesis. However, the ribose C2'-hydroxyl of the initiating cytidylate can accept several modifications and retain the ability to direct synthesis. The only unacceptable modification is a protonated C2'-amino group. Quite strikingly, the recognition of the functional groups for the initiation cytidylate and other template nucleotides are different. For example, a C5-methyl group in cytidine can direct RNA synthesis at all template positions except at the initiation cytidylate and C2'-amino modifications are tolerated better after the +11 position. When a 4-thiouracil (4sU) base analog that allows only imperfect base pairing with the nascent RNA is placed at different positions in the template, the efficiency of synthesis is correlated with the calculated stability of the template-nascent RNA duplex adjacent to the position of the 4sU. These results define the requirements for the specific interactions required for the initiation of RNA synthesis and will be compared to the mechanisms of initiation by other RNA-dependent and DNA-dependent RNA polymerases.     

Mapping of the functional phosphate groups in the catalytic core of deoxyribozyme 10-23

The RNA phosphodiester bond cleavage activity of a series of 16 thio-deoxyribozymes 10-23, containing a P-stereorandom single phosphorothioate linkage in predetermined positions of the catalytic core from P1 to P16, was evaluated under single-turnover conditions in the presence of either 3 mM Mg(2+) or 3 mM Mn(2+). A metal-specificity switch approach permitted the identification of nonbridging phosphate oxygens (proR(P) or proS(P)) located at seven positions of the core (P2, P4 and P9-13) involved in direct coordination with a divalent metal ion(s). By contrast, phosphorothioates at positions P3, P6, P7 and P14-16 displayed no functional relevance in the deoxyribozyme-mediated catalysis. Interestingly, phosphorothioate modifications at positions P1 or P8 enhanced the catalytic efficiency of the enzyme. Among the tested deoxyribozymes, thio-substitution at position P5 had the largest deleterious effect on the catalytic rate in the presence of Mg(2+), and this was reversed in the presence of Mn(2+). Further experiments with thio-deoxyribozymes of stereodefined P-chirality suggested direct involvement of both oxygens of the P5 phosphate and the proR(P) oxygen at P9 in the metal ion coordination. In addition, it was found that the oxygen atom at C6 of G(6) contributes to metal ion binding and that this interaction is essential for 10-23 deoxyribozyme catalytic activity.     

Poliovirus RNA-dependent RNA polymerase (3Dpol): structural, biochemical, and biological analysis of conserved structural motifs A and B

We have constructed a structural model for poliovirus RNA-dependent RNA polymerase (3D(pol)) in complex with a primer-template (sym/sub) and ATP. Residues found in conserved structural motifs A (Asp-238) and B (Asn-297) are involved in nucleotide selection. Asp-238 appears to couple binding of nucleotides with the correct sugar configuration to catalytic efficiency at the active site of the enzyme. Asn-297 is involved in selection of ribonucleoside triphosphates over 2'-dNTPs, a role mediated most likely via a hydrogen bond between the side chain of this residue and the 2'-OH of the ribonucleoside triphosphate. Substitutions at position 238 or 297 of 3D(pol) produced derivatives exhibiting a range of catalytic efficiencies when assayed in vitro for poly(rU) polymerase activity or sym/sub elongation activity. A direct correlation existed between activity on sym/sub and biological phenotypes; a 2.5-fold reduction in polymerase elongation rate produced virus with a temperature-sensitive growth phenotype. These data permit us to propose a detailed, structural model for nucleotide selection by 3D(pol), confirm the biological relevance of the sym/sub system, and provide additional evidence for kinetic coupling between RNA synthesis and subsequent steps in the virus life cycle.     

The mRNA codon environment at the P and E sites of human ribosomes deduced from photo crosslinking with pUUUGUU

Three mRNA analogs--derivatives of hexaribonucleotide pUUUGUU comprising phenylalanine and valine codons with a perfluoroarylazido group attached to the C5 atom of the uridine residue at the first, second, or third position--were used for photocrosslinking with 80S ribosomes from human placenta. The mRNA analogs were positioned on the ribosome with tRNA recognizing these codons: UUU was at the P site if tRNA(Phe) was used, while tRNA(Val) was used to put there the GUU codon (UUU at the E site). Thus, the crosslinking group of mRNA analog might occupy positions -3 to +3 with respect to the first nucleotide of the codon at the P site. Irradiation of the complexes with soft UV light (lambda > 280 nm) resulted in the crosslinking of pUUUGUU derivatives with 18S RNA and proteins in the ribosome small subunit. The crosslinking with rRNA was observed only in the presence of tRNA. The photoactivatable group in positions -1 to +3 binds to G1207, while that in positions -2 or -3 binds to G961 of 18S RNA. In all cases, we observed crosslinking with S2 and S3 proteins irrespective of the presence of tRNA in the complex. Crosslinking with S23 and S26 proteins was observed mainly in the presence of tRNA when modified nucleotide occupied the +1 position (for both proteins) or the -3 position (for S26 protein). The crosslinking with S5/S7 proteins was substantial when modified nucleotide was in the -3 position, this crosslinking was not observed in the absence of tRNA.     

Model of the regulation of activity of immobilized enzymes (amylases) in soil

The preservation of activity of extracellular enzymes in soil is presently associated with their immobilization on organic or inorganic carriers. Enzyme immobilization results, however, in a significant decrease in enzymatic activity. In the present work, the mechanism responsible for promotion of the catalytic activity was revealed, as well as the favorable effect of low-molecular alkylhydrozybenzenes of the class of alkylresorcinols, which are common in soil organic matter, on stability of immobilized enzymes (exemplified by amylases) by their post-translational modification. Optimal conditions (enzyme to sorbent ratio, pH optimum, CaCl₂ concentration, and sorption time) for amylase sorption on a biological sorbent (yeast cell walls) were determined and decreased activity of the immobilized enzyme compared to its dissolved state was confirmed. Alkylresorcinols (C₇AHB) at concentrations of 1.6 to 80 mM were found to cause an increase of amylase activity both in the case of already sorbed enzymes (by 30%) and in the case of a free dissolved enzyme with its subsequent immobilization (by 50–60%). In both cases, the optimal C₇AHB concentration was 16 mM. Amylase stability was determined for C7AHB-modified and unmodified enzymes immobilized on the biological sorbent after two cycles of freezing (–20°C) and thawing (4°C). Inverse dependence was revealed between increasing stability of C₇AHB-modified enzymes and an increase in their activity, as well as higher stability of immobilized modified amylases than of the dissolved modified enzyme. Investigation of the effect of C₇HOB-modification in the preservation of activity in immobilized amylases after four freeze–thaw cycles revealed: (1) better preservation of activity by the modified immobilized enzymes compared to immobilized ones; (2) differences in the dynamics of activity loss within compared pairs, with activity of immobilized amylases decreasing after the second cycle to a lower level (42%) than activity of the modified immobilized enzymes after the fourth cycle (48%). These results demonstrate that in the preservation of activity of extracellular enzymes in soil both stabilization mechanisms are of importance: immobilization on organic carriers and modification of the enzyme conformation by low-molecular compounds with the functions of chemical chaperones.     

Mutations in a nonconserved sequence of the Tetrahymena ribozyme increase activity and specificity

The RNA substrate-binding site of the Tetrahymena ribozyme is connected to the catalytic core by the joining region J1/2. Although J1/2 is not conserved among group I introns, small insertions or deletions in this sequence have dramatic effects, enhancing the turnover number and sequence specificity of ribozyme-catalyzed RNA cleavage. Measurements of rate constants for individual steps in the reaction have revealed the basis of these improvements. Ironically, the higher turnover and specificity both result from decreased affinity for RNA, rather than better cleavage. These results provide evidence that the nonconserved J1/2 sequence positions the RNA substrate to optimize tertiary interactions and ensure cleavage at the position corresponding to the 5' splice site. The wild-type RNA is well adapted to its biological function, and its limitations in multiple turnover can be corrected by mutation.     

Three conserved guanosines approach the reaction site in native and minimal hammerhead ribozymes

Native hammerhead ribozymes contain RNA domains that enable high catalytic activity under physiological conditions, where minimal hammerheads show little activity. However, little is known about potential differences in native versus minimal ribozyme folding. Here, we present results of photocross-linking analysis of native and minimal hammerheads containing photoreactive nucleobases 6-thioguanosine, 2,6-diaminopurine, 4-thiouridine, and pyrrolocytidine, introduced at specific sites within the catalytic core. Under conditions where catalytic activity is observed, the two substrate nucleobases spanning the cleavage site approach and stack upon G8 and G12 of the native hammerhead, two conserved nucleobases that show similar behavior in minimal constructs, have been implicated in general acid-base catalysis, and are >15 A from the cleavage site in the crystal structures. Pyrrolocytidine at cleavage site position 17 forms an efficient crosslink to G12, and the crosslinked RNA retains catalytic activity. Multiple cross-linked species point to a structural rearrangement within the U-turn, positioning residue G5 in the vicinity of cleavage site position 1.1. Intriguing crosslinks were triggered by nucleotide analogues at positions distal to the crosslinked residues; for example, 6-thioguanosine at position 5 induced a crosslink between G12 and C17, suggesting an intimate functional communication among these three nucleobases. Together, these results support a model in which the native hammerhead folds to an active structure similar to that of the minimal ribozyme, and significantly different from the crystallographic structures.     

Stabilizing basic fibroblast growth factor using protein engineering

Using site directed mutagenesis, each of the four cysteines present at amino acid residues 26, 70, 88, and 93 of the mature protein of human basic fibroblast growth factor (bFGF) was individually changed to serine. The biological activity and heparin binding ability was retained when the serine was substituted for the cysteine residue at either 70 or 88 of the bFGF protein. This finding indicates that the cysteines at these positions are not essential for expressing biological activity. The substitution of the residues at these positions, especially at position 88, reduced the heterogeneity recognized as several peaks of bFGF eluted from a heparin affinity column, even after oxidation with hydrogen peroxide, suggesting that the cysteines at these positions are exposed to the surface of the molecule to form disulfide bonds that induce heterologous conformations. Furthermore, under acidic conditions, these modified bFGFs are revealed to be more stable in maintaining their activity. These facts suggest that this protein has been successfully modified by protein engineering.