Free energy perturbation calculations on binding and catalysis after mutating Asn 155 in subtilisin

Site-directed mutagenesis is a very powerful approach to altering the biological functions of proteins, the structural stability of proteins and the interactions of proteins with other molecules. Several experimental studies in recent years have been directed at estimating the changes in catalytic properties, (rates of binding and catalysis) in site-directed mutants of enzymes compared to the native enzymes. Simulation approaches to the study of complex molecules have also become more powerful, in no small measure owing to the increase in computer power. These simulations have often allowed results of experiments to be rationalized and understood mechanistically. A new approach called the free-energy pertubation method, which uses statistical mechanics and molecular dynamics can often be used for quantitative calculation of free energy differences. We have applied such a technique to calculate the differential free energy of binding and free energy of activation for catalysis of a tripeptide substrate by native subtilisin and a subtilisin mutant (Asn 155----Ala 155). Our studies lead to a calculated difference in free energy of binding which is relatively small, but a calculated change in free energy of catalysis which is substantial. These energies are very close to those determined experimentally (J. A. Wells and D. A. Estell, personal communication), which were not known to us until the simulations were completed. This demonstrates the predictive power and utility of theoretical simulation methods in studies of the effects of site-specific mutagenesis on both enzyme binding and catalysis.     

Biosynthesis of selenosubtilisin: A novel way to target selenium into the active site of subtilisin

Glutathione peroxidase （GPx, EC1.11.1.9）, an important anti-oxidative selenoenzyme, can catalyze the reduction of harmful hydroperoxides with concomitant glutathione, thereby protecting cells and other biological issues against oxidative damage. It captures considerable interest in redesign of its function for either the mechanism study or the pharmacological development as an antioxidant. In order to develop a general strategy for specifically targeting and operating selenium in active sites of enzymes, the catalytically essential residue selenocysteine （Sec） was first successfully bioincorporated into the catalytic center of subtilisin by using an auxotrophic expression system. The studies of the catalytic activity and the steady-state kinetics demonstrated that selenosubtilisin is an excellent GPx-like biocatalyst. In comparison with the chemically modified method, biosynthesis exhibits obvious advantages： Sec could be site-directly incorporated into active sites of enzymes to overcome the non-specificity generated by chemical modification. This study provides an important strategy for specifically targeting and operating selenium in the active site of an enzyme.     

A mutant T4 lysozyme displays five different crystal conformations

Phage T4 lysozyme consists of two domains between which is formed the active-site cleft of the enzyme. The crystallographically determined thermal displacement parameters for the protein suggested that the amino terminal of the two domains undergoes 'hinge-bending' motion about an axis passing through the waist of the molecule. Such conformational mobility may be important in allowing access of substrates to the active site of the enzyme. We report here a crystallographic study of a mutant T4 lysozyme which demonstrates further the conformational flexibility of the protein. A mutant form of the enzyme with a methionine residue (Met 6) replaced by isoleucine crystallizes with four independent molecules in the crystal lattice. These four molecules have distinctly different conformations. The mutant protein can also crystallize in standard form with a structure very similar to the wild-type protein. Thus the mutant protein can adopt five different crystal conformations. The isoleucine for methionine substitution at the intersection of the two domains of T4 lysozyme apparently enhances the hinge-bending motion presumed to occur in the wild-type protein, without significantly affecting the catalytic activity or thermal stability of the protein.     

GILT is a critical host factor for Listeria monocytogenes infection

Listeria monocytogenes is a gram-positive, intracellular, food-borne pathogen that can cause severe illness in humans and animals. On infection, it is actively phagocytosed by macrophages; it then escapes from the phagosome, replicates in the cytosol, and subsequently spreads from cell to cell by a non-lytic mechanism driven by actin polymerization. Penetration of the phagosomal membrane is initiated by the secreted haemolysin listeriolysin O (LLO), which is essential for vacuolar escape in vitro and for virulence in animal models of infection. Reduction is required to activate the lytic activity of LLO in vitro, and we show here that reduction by the enzyme gamma-interferon-inducible lysosomal thiol reductase (GILT, also called Ifi30) is responsible for the activation of LLO in vivo. GILT is a soluble thiol reductase expressed constitutively within the lysosomes of antigen-presenting cells, and it accumulates in macrophage phagosomes as they mature into phagolysosomes. The enzyme is delivered by a mannose-6-phosphate receptor-dependent mechanism to the endocytic pathway, where amino- and carboxy-terminal pro-peptides are cleaved to generate a 30-kDa mature enzyme. The active site of GILT contains two cysteine residues in a CXXC motif that catalyses the reduction of disulphide bonds. Mice lacking GILT are deficient in generating major histocompatibility complex class-II-restricted CD4(+) T-cell responses to protein antigens that contain disulphide bonds. Here we show that these mice are resistant to L. monocytogenes infection. Replication of the organism in GILT-negative macrophages, or macrophages expressing an enzymatically inactive GILT mutant, is impaired because of delayed escape from the phagosome. GILT activates LLO within the phagosome by the thiol reductase mechanism shared by members of the thioredoxin family. In addition, purified GILT activates recombinant LLO, facilitating membrane permeabilization and red blood cell lysis. The data show that GILT is a critical host factor that facilitates L. monocytogenes infection.     

The alpha-lytic protease pro-region does not require a physical linkage to activate the protease domain in vivo

alpha-Lytic protease, an extracellular serine protease of Lysobacter enzymogenes 495, is synthesized as a pre-pro-protein. Previously it has been shown that when expressed in Escherichia coli, the protein is autocatalytically processed in the periplasmic space, and that the functional protease domain accumulates extracellularly. Engineered proteins lacking the 166 amino-acid pro-region were enzymatically inactive and remained cell-associated. By independently expressing the pro- and protease domains in vivo, evidence is provided here that direct covalent linkage is not required for production of active protease. We postulate that the pro-region acts as a template to promote the folding of the protease domain into an active configuration. Our results, combined with recent experiments on the evolutionarily unrelated subtilisin E (ref. 3), suggest that the ability of the pro-region of these bacterial proteases to facilitate folding of their protease domains is not a curiosity of a single system, but may reflect a general property of extracellular bacterial serine proteases.     

Suppression of c-ras transformation by GTPase-activating protein

The ras genes are required for normal cell growth and mediate transformation by oncogenes encoding protein tyrosine kinases. Normal ras can transform cells in vitro and in vivo, but mutationally activated ras does so much more efficiently, and highly transforming mutant versions of ras have been isolated from a variety of human and animal tumours. The ras genes encode membrane-associated, guanine nucleotide-binding proteins that are active when GTP is bound and inactive when GDP is bound. The slow intrinsic GTPase activity of normal mammalian Ras proteins can be greatly accelerated by the GTPase-activating protein (GAP), which is predominantly cytoplasmic. This activity of GAP, which can increase with cell density in contact-inhibited cells, suggests that it functions as a negative, upstream regulator of ras. Other studies, however, show that GAP interacts with a region of ras-encoded protein implicated in ras effector function, which raises the possibility that GAP might also be a downstream target of ras. Mutationally activated ras-encoded proteins also interact with GAP, although they are resistant to its catalytic activity. In an attempt to define the role of GAP in ras-mediated transformation, we examined the effects on transformation of normal or mutant ras when cells overexpress GAP. We found that GAP suppresses transformation of NIH 3T3 cells by normal Ha-ras (c-ras) but does not inhibit transformation by activated Ha-ras (v-ras). These results support the hypothesis that GAP functions as a negative regulator of normal ras and make it unlikely that GAP alone is the ras target.     

Prediction of electrostatic effects of engineering of protein charges

Accurate prediction of electrostatic effects on catalytic activity is an essential component of protein design. Site-directed mutagenesis of charged groups in subtilisin of Bacillus amyloliquefaciens has provided experimental measurements of electrostatic interactions which may be used to test such theoretical methods. The pKa of the histidine of the active site has been perturbed by +0.08 to -1.0 units by modifying one or two residues. Electrostatic effects in proteins can be modelled by the algorithm of Warwicker and Watson, which uses classical electrostatics and considers both the charge position and the shape of the molecule. Here we report that the algorithm can model several pKa shifts in subtilisin to fair accuracy.     

Thermostability of subtilisin nattokinase obtained by site-directed mutagenesis

To study the thermostability of Nattokinase(subtilisin NAT,NK),three double mutant plasmids(pET-28a-NKG61C/S98C,pET-28a-NKT22C/S87C,pET-28a-NKS24C/S87C)were constructed by site-directed mutagenesis.Target enzymes were detected using SDS-PAGE and disulfide bond formation was detected using Western blotting analysis.Thermostability was tested by rates of inactivation at certain temperature.The results showed that disulfide bond was not formed within two cysteines and the thermostability of three double mutants was not increased compared with the wild-type NK.The thermostability of NK performed in Ca2+was stronger than in ethylenediaminetetraacetic acid(EDTA).But when the temperature reached 62℃,the enzymes rapidly denatured and inactivated even in the presence of Ca2+.Although the thermostability of mutants was not increased,this study shows a tendency of improving thermostability of NK in protein engineering.     

Calculation of electrostatic potentials in an enzyme active site

To be able to calculate the contributions of individual amino acids to the electrostatic field of a protein would be of considerable value in designing proteins of enhanced or altered function and stability. Recent studies on the serine protease subtilisin provide direct measurements of the electrostatic potential in the active site of the enzyme produced by two charged amino acids. We have used these results to test a recently developed method for the calculation of electrostatic interactions between two specific sites on a protein. The extent of agreement between the theoretical and experimental results suggests that the continuum solvent model used in the calculations reproduces the essential features of the interaction.     

Involvement of the 'leucine zipper' region in the oligomerization and transforming activity of human c-myc protein

c-Myc plays a part in the regulation of important cellular processes such as growth, differentiation and neoplastic transformation. Although c-myc gene structure and expression are well characterized, the function and biochemical properties of the protein are less well understood. Human c-myc is a 439-amino acid phosphoprotein which binds DNA in vitro and belongs to a discrete subset of nuclear proteins. Using the human c-myc mutants generated by linker-insertion and deletion mutagenesis, we have defined regions of the protein that are important for its transforming activities and its nuclear localization. Here, we show that human c-myc exists as an oligomer in vitro and use mutant proteins to localize the oligomerization domain to a carboxyl-terminal peptide containing the 'leucine zipper' motif. The 'leucine zipper' describes a structure found in a number of DNA-binding proteins that contains leucines occurring at intervals of every seventh amino acid in a region predicted to be alpha-helical. The 'leucine zipper' might mediate dimerization by intermolecular interdigitation of the leucine side-chains. We show that a c-myc mutant, which is inactive but can oligomerize, dominantly inhibits the cotransforming activity with wild-type c-myc of rat embryo cells, whereas inactive mutants which cannot oligomerize properly because of deletions in the oligomerization domain are recessive.     

Receptor binding redefined by a structural switch in a mutant human insulin

Crystal structures of insulin have been determined in various distinct forms, the relevance of which to receptor recognition has long been the subject of speculation. Recently the crystal structure of an inactive insulin analogue has been determined and, surprisingly, found to have a conformation identical to native insulin. On this basis Dodson and colleagues have suggested that the known insulin crystal structures reflect an inactive conformation, and that a change in conformation is required for activity--specifically, the carboxy terminal residues of the B-chain are proposed to separate from the amino terminal residues of the A-chain. Here we report the solution structure of an active insulin mutant, determined by two-dimensional NMR, which supports this hypothesis. In the mutant, the carboxy terminal beta-turn and beta-strand of the B-chain are destabilized and do not pack across the rest of the molecule. We suggest that analogous detachment of the carboxy terminal region of the B-chain occurs in native insulin on binding to its receptor. Our finding that partial unfolding of the B-chain exposes an alternative protein surface rationalizes the receptor-binding properties of a series of anomalous insulin analogues, including a mutant insulin associated with diabetes mellitus in man.     

The structural basis of agonist-induced activation in constitutively active rhodopsin

G-protein-coupled receptors (GPCRs) comprise the largest family of membrane proteins in the human genome and mediate cellular responses to an extensive array of hormones, neurotransmitters and sensory stimuli. Although some crystal structures have been determined for GPCRs, most are for modified forms, showing little basal activity, and are bound to inverse agonists or antagonists. Consequently, these structures correspond to receptors in their inactive states. The visual pigment rhodopsin is the only GPCR for which structures exist that are thought to be in the active state. However, these structures are for the apoprotein, or opsin, form that does not contain the agonist all-trans retinal. Here we present a crystal structure at a resolution of 3 ? for the constitutively active rhodopsin mutant Glu 113 Gln in complex with a peptide derived from the carboxy terminus of the α-subunit of the G protein transducin. The protein is in an active conformation that retains retinal in the binding pocket after photoactivation. Comparison with the structure of ground-state rhodopsin suggests how translocation of the retinal β-ionone ring leads to a rotation of transmembrane helix 6, which is the critical conformational change on activation. A key feature of this conformational change is a reorganization of water-mediated hydrogen-bond networks between the retinal-binding pocket and three of the most conserved GPCR sequence motifs. We thus show how an agonist ligand can activate its GPCR.     

Enhancing protein C interaction with thrombin results in a clot-activated anticoagulant.

Human protein C is a vitamin K-dependent plasma glycoprotein that circulates as an inactive zymogen. At the endothelial cell surface, thrombin in complex with the integral membrane protein thrombomodulin converts protein C to its active form by specific cleavage of an activation peptide. The activated form of protein C has potent anticoagulant activity as a feedback regulator of thrombin generation (reviewed in refs 4-6), and also has profibrinolytic, anti-ischaemic and anti-inflammatory properties. Protein C is effective in the treatment of model and human thrombotic diseases but, except when it has been used to treat genetic or acquired deficiencies and microvascular thrombosis, it is administered as the activated enzyme, which has a short biological half-life. We have altered two putative inhibitory acidic residues near the thrombin cleavage site, which results in a 30-fold increase in substrate utilization by alpha-thrombin. We combined these changes with a genetically altered glycoform to generate a zymogen protein C with a 60-fold increased cleavage rate by free alpha-thrombin, independent of its cofactor thrombomodulin. We show that this 'proform' of protein C, unlike the natural circulating zymogen, can be activated by thrombin generated in clotting human plasma, resulting in an inhibition of further clot formation. Our data therefore show that we have engineered a site-activated agent, which only has anticoagulant activity when significant amounts of thrombin are being generated.     

产纤溶酶菌种的鉴定及纤溶酶的分离纯化

为了获得纯纤溶酶并探讨其酶学性质,从广泛收集的豆豉中筛选到一株具有高产纤溶酶能力的菌株,经鉴定为枯草芽孢杆菌.将该菌突变株DC-12N8的发酵液经离心除菌、硫酸铵分级沉淀、DEAE-Sepharose Fast Flow和CM-Sepharose Fast Flow离子交换层析、Sephadex G-75凝胶过滤,获得电泳纯的豆豉纤溶酶,每升发酵液可得到4.5mg活性酶蛋白,每克酶蛋白活力达1.18mol/s,纯度为发酵原液的53.6倍,回收率为11.7%.SDS-聚丙烯酰胺凝胶电泳表明,该酶是单链蛋白质,其分子质量约为28ku.     

Idealization of the hydrophobic segment of the alkaline phosphatase signal peptide

Proteins secreted by prokaryotic cells are synthesized as precursors containing an amino-terminal extension sequence or signal peptide. Although these signal peptides share little primary sequence homology, recent studies suggest that they function via common pathways during the transport process and that a common element may reside in their secondary structural characteristics. We are investigating the role of an idealized hydrophobic sequence with high potential for alpha-helix formation in the Escherichia coli alkaline phosphatase signal peptide. Here, amino-acid substitutions were made using site-directed mutagenesis to produce a mutant signal sequence containing nine consecutive leucine residues in the hydrophobic core segment. Transport studies with this mutant precursor indicate that mature alkaline phosphatase is correctly targeted to the E. coli periplasm and that processing of the precursor to the mature form of the enzyme is extremely rapid. In contrast, processing is slowed when the mutant signal sequence is lengthened by the insertion of five additional leucine residues and one serine.     

Serpin-resistant mutants of human tissue-type plasminogen activator

Tissue-type plasminogen activator (t-PA) converts the inactive zymogen, plasminogen, into the powerful protease, plasmin, which then degrades the fibrin meshwork of thrombi. To prevent systemic activation of plasminogen, plasma contains several inhibitors of t-PA, the most important of which is plasminogen activator inhibitor-1 (PAI-1), a member of the serpin superfamily. As the ability to produce serpin-resistant variants of t-PA could increase the potential of this enzyme as a thrombolytic agent, we have used the known three-dimensional structure of the complex between trypsin and bovine pancreatic trypsin inhibitor (BPTI) to model the interactions between the active site of human t-PA and PAI-1. On the basis of this model we then altered by site-directed mutagenesis those amino acids of t-PA predicted to make contact with PAI-1 but not with the substrate plasminogen. We report here that although the resulting mutants have enzymatic properties similar to those of wild-type t-PA, they display significant resistance to inhibition by PAI-1. For example, following incubation with an amount of the serpin that completely inhibits the wild-type enzyme, one variant retains 95% of its initial activity. This mutant is also resistant to inhibition by the complex mixture of serpins present in human plasma.     

Visualizing transient events in amino-terminal autoprocessing of HIV-1 protease

HIV-1 protease processes the Gag and Gag-Pol polyproteins into mature structural and functional proteins, including itself, and is therefore indispensable for viral maturation. The mature protease is active only as a dimer with each subunit contributing catalytic residues. The full-length transframe region protease precursor appears to be monomeric yet undergoes maturation via intramolecular cleavage of a putative precursor dimer, concomitant with the appearance of mature-like catalytic activity. How such intramolecular cleavage can occur when the amino and carboxy termini of the mature protease are part of an intersubunit beta-sheet located distal from the active site is unclear. Here we visualize the early events in N-terminal autoprocessing using an inactive mini-precursor with a four-residue N-terminal extension that mimics the transframe region protease precursor. Using paramagnetic relaxation enhancement, a technique that is exquisitely sensitive to the presence of minor species, we show that the mini-precursor forms highly transient, lowly populated (3-5%) dimeric encounter complexes that involve the mature dimer interface but occupy a wide range of subunit orientations relative to the mature dimer. Furthermore, the occupancy of the mature dimer configuration constitutes a very small fraction of the self-associated species (accounting for the very low enzymatic activity of the protease precursor), and the N-terminal extension makes transient intra- and intersubunit contacts with the substrate binding site and is therefore available for autocleavage when the correct dimer orientation is sampled within the encounter complex ensemble.     

Direct control of shoot meristem activity by a cytokinin-activating enzyme

The growth of plants depends on continuous function of the meristems. Shoot meristems are responsible for all the post-embryonic aerial organs, such as leaves, stems and flowers. It has been assumed that the phytohormone cytokinin has a positive role in shoot meristem function. A severe reduction in the size of meristems in a mutant that is defective in all of its cytokinin receptors has provided compelling evidence that cytokinin is required for meristem activity. Here, we report a novel regulation of meristem activity, which is executed by the meristem-specific activation of cytokinins. The LONELY GUY (LOG) gene of rice is required to maintain meristem activity and its loss of function causes premature termination of the shoot meristem. LOG encodes a novel cytokinin-activating enzyme that works in the final step of bioactive cytokinin synthesis. Revising the long-held idea of multistep reactions, LOG directly converts inactive cytokinin nucleotides to the free-base forms, which are biologically active, by its cytokinin-specific phosphoribohydrolase activity. LOG messenger RNA is specifically localized in shoot meristem tips, indicating the activation of cytokinins in a specific developmental domain. We propose the fine-tuning of concentrations and the spatial distribution of bioactive cytokinins by a cytokinin-activating enzyme as a mechanism that regulates meristem activity.