Molecular Dynamics Simulations on Pars Intercerebralis Major Peptide-C (PMP-C) Reveal the Role of Glycosylation and Disulfide Bonds in its Enhanced Structural Stability and Function

Abstract

Fucosylation of Thr 9 in pars intercerebralis major peptide-C (PMP-C) enhances its structural stability and functional ability as a serine protease inhibitor. In order to understand the role of disulfide bonds and glycosylation on the structure and function of PMP-C, we have carried out multiple explicit solvent molecular dynamics (MD) simulations on fucosylated and non-fucosylated forms of PMP-C, both in the presence and absence of the disulfide bonds. Our simulations revealed that there were no significant structural changes in the native disulfide bonded forms of PMP-C due to fucosylation. On the other hand, the non-fucosylated form of PMP-C without disulfide bonds showed larger deviations from the starting structure than the fucosylated form. However, the structural deviations were restricted to the terminal regions while core β-sheet retained its hydrogen bonded structure even in absence of disulfide bonds as well as fucosylation. Interestingly, fucosylation of disulfide bonded native PMP-C led to a decreased thermal flexibility in the residue stretch 29–32 which is known to interact with the active site of the target proteases. Our analysis revealed that disulfide bonds covalently connect the residue stretch 29–32 to the central β-sheet of PMP-C and using a novel network of side chain interactions and disulfide bonds fucosylation at Thr 9 is altering the flexibility of the stretch 29–32 located at a distal site. Thus, our simulations explain for the first time, how presence of disulfide bonds between conserved cysteines and fucosylation enhance the function of PMP-C as a protease inhibitor.     

Direct observations of shifts in the β-sheet register of a protein-peptide complex using explicit solvent simulations

Using explicit solvent molecular dynamics simulations, we were able to obtain direct observations of shifts in the hydrogen-bonding register of an intermolecular β-sheet protein-peptide complex. The β-sheet is formed between the FHA domain of cancer marker protein Ki67 (Ki67FHA) and a peptide fragment of the hNIFK signaling protein. Potential encounter complexes of the Ki67FHA receptor and hNIFK peptide are misregistered states of the β-sheet. Rearrangements of one of these misregistered states to the native state were captured in three independent simulations. All three rearrangements occurred by a common mechanism: an aromatic residue of the peptide (F263) anchors into a transient hydrophobic pocket of the receptor to facilitate the formation of native hydrogen bonds. To our knowledge, these simulations provide the first atomically detailed visualizations of a mechanism by which nature might correct for errors in the alignment of intermolecular β-sheets.     

Conformational properties of the disease-causing Z variant of α1-antitrypsin revealed by theory and experiment

The human serine protease inhibitor (serpin) α-1 antitrypsin (α1-AT) protects tissues from proteases of inflammatory cells. The most common disease-causing mutation in α1-AT is the Z-mutation (E342K) that results in an increased propensity of α1-AT to polymerize in the ER of hepatocytes, leading to a lack of secretion into the circulation. The structural consequences of this mutation, however, remain elusive. We report a comparative molecular dynamics investigation of the native states of wild-type and Z α1-AT, revealing a striking contrast between their structures and dynamics in the breach region at the top of β-sheet A, which is closed in the wild-type simulations but open in the Z form. Our findings are consistent with experimental observations, notably the increased solvent exposure of buried residues in the breach region in Z, as well as polymerization via domain swapping, whereby the reactive center loop is rapidly inserted into an open A-sheet before proper folding of the C-terminal β-strands, allowing C-terminal domain swapping with a neighboring molecule. Taken together, our experimental and simulation data imply that mutations at residue 342 that either stabilize an open form of the top of β-sheet A or increase the local flexibility in this region, may favor polymerization and hence aggregation.     

Structural,mutagenic and in silico studies of xyloglucan fucosylation in Arabidopsis thaliana suggest a water‐mediated mechanism

The mechanistic underpinnings of the complex process of plant polysaccharide biosynthesis are poorly understood, largely because of the resistance of glycosyltransferase (GT) enzymes to structural characterization. In Arabidopsis thaliana, a glycosyl transferase family 37 (GT37) fucosyltransferase 1 (AtFUT1) catalyzes the regiospecific transfer of terminal 1,2‐fucosyl residues to xyloglucan side chains – a key step in the biosynthesis of fucosylated sidechains of galactoxyloglucan. We unravel the mechanistic basis for fucosylation by AtFUT1 with a multipronged approach involving protein expression, X‐ray crystallography, mutagenesis experiments and molecular simulations. Mammalian cell culture expressions enable the sufficient production of the enzyme for X‐ray crystallography, which reveals the structural architecture of AtFUT1 in complex with bound donor and acceptor substrate analogs. The lack of an appropriately positioned active site residue as a catalytic base leads us to propose an atypical water‐mediated fucosylation mechanism facilitated by an H‐bonded network, which is corroborated by mutagenesis experiments as well as detailed atomistic simulations.     

Atomistic details of effect of disulfide bond reduction on active site of Phytase B from Aspergillus niger: A MD Study

The molecular integrity of the active site of phytases from fungi is critical for maintaining phytase function as efficient catalytic machines. In this study, the molecular dynamics (MD) of two monomers of phytase B from Aspergillus niger, the disulfide intact monomer (NAP) and a monomer with broken disulfide bonds (RAP), were simulated to explore the conformational basis of the loss of catalytic activity when disulfide bonds are broken. The simulations indicated that the overall secondary and tertiary structures of the two monomers were nearly identical but differed in some crucial secondary–structural elements in the vicinity of the disulfide bonds and catalytic site. Disulfide bonds stabilize the β-sheet that contains residue Arg66 of the active site and destabilize the α-helix that contains the catalytic residue Asp319. This stabilization and destabilization lead to changes in the shape of the active–site pocket. Functionally important hydrogen bonds and atomic fluctuations in the catalytic pocket change during the RAP simulation. None of the disulfide bonds are in or near the catalytic pocket but are most likely essential for maintaining the native conformation of the catalytic site.

Abbreviations

PhyB - 2.5 pH acid phophatese from Aspergillus niger, NAP - disulphide intact monomer of Phytase B, RAP - disulphide reduced monomer of Phytase B, Rg - radius of gyration, RMSD - root mean square deviation, MD - molecular dynamics.     

A Molecular Dynamics Study of Human Defensins HBD-1 and HNP-3 in Water

Abstract

Mammalian defensins are crucial components of the innate immune system. They are characterized by three disulfide bridges and exhibit broad spectrum antibacterial activity. The spacing between the cysteines and disulfide connectivities in the two classes of defensins, the α- and β-forms, are different. The structural motif of 3 β-strands appears to be conserved in α and β-defensins despite differences in disulfide connectivities and spacing between cysteines. In this study, Molecular Dynamics Simulations (MDS) have been carried out to study the conformational behavior of α- and β-defensins with and without disulfide bridges. Our results indicate that β-strands in the C-terminal region of HBD-1 and HNP-3 do not unfold during the course of MDS. The segment adopting α-helix in HBD-1 unfolds early during the simulations. The backbone hydrogen bonds in HBD-1 and HNP-3 are broken during MDS. When the disulfide bonds are absent, the N-terminal β-strand unfolds by 20 ns but β-strands are observed in the C-terminal region of HNP-3. HBD-1, without disulfide bridges, unfolds to a greater extent during the course of the MDS. Examination of distances between sulfur atoms of cysteines without disulfide bridges during the simulations indicate that there is no specific preference for native disulfide bridges, which could be the reason for the experimental observation of non-native disulfide bridge formation during chemical synthesis of human α- and β-defensins. Since defensins with non-native disulfide bridges are biologically active, the exact three dimensional structures observed for native HBD-1 and HNP-3 does not appear to be essential for exhibiting antibacterial activity.     

Simulation Studies on the Stabilities of Aggregates Formed by Fibril-Forming Segments of α-Synuclein

Abstract

We performed molecular dynamics simulations for various oligomers with different β-sheet conformations consisting of α-Synuclein 71–82 residues using an all atom force field and explicit water model. Tetramers of antiparallel β-sheet are shown to be stable, whereas parallel sheets are highly unstable due to the repulsive interactions between bulky and polar side chains as well as the weaker backbone hydrogen bonds. We also investigated the stabilities of double antiparallel β-sheets stacked with asymmetric and symmetric geometries. Our results show that this 12 amino acid residue peptide can form stable β-sheet conformers at 320K and higher temperatures. The backbone hydrogen bonds in β-sheet and the steric packing between hydrophobic side chains between β-sheets are shown to give conformational stabilities.     

Characterization of an alternative low energy fold for bovine α-lactalbumin formed by disulfide bond shuffling

Bovine α-lactalbumin (αLA) forms a misfolded disulfide bond shuffled isomer, X-αLA. This X-αLA isomer contains two native disulfide bridges (Cys 6-Cys 120 and Cys 28-Cys 111) and two non-native disulfide bridges (Cys 61-Cys 73 and Cys 77-Cys 91). MD simulations have been used to characterize the X-αLA isomer and its formation via disulfide bond shuffling and to compare it with the native fold of αLA. In the simulations of the X-αLA isomer the structure of the α-domain of native αLA is largely retained in agreement with experimental data. However, there are significant rearrangements in the β-domain, including the loss of the native β-sheet and calcium binding site. Interestingly, the energies of X-αLA and native αLA in simulations in the absence of calcium are closely similar. Thus, the X-αLA isomer represents a different low energy fold for the protein. Calcium binding to native αLA is shown to help preserve the structure of the β-domain of the protein limiting possibilities for disulfide bond shuffling. Hence, binding calcium plays an important role in both maintaining the native structure of αLA and providing a mechanism for distinguishing between folded and misfolded species.     

核心岩藻糖基化控制在治疗性抗体研究中的应用

对于人类许多疾病,抗体已成为了一种很重要的治疗制剂。如何进一步提高抗体的效能,是目前研究的主要热点之一。在抗体治疗过程中,抗体依赖性细胞介导的细胞毒作用(antibody dependent cellular cytotoxicity,ADCC)被认为是治疗性抗体临床疗效中一个重要的治疗功能。抗体恒定区Fc片段(crystalline fragment)的糖基化对ADCC起着至关重要的作用,尤其是抗体恒定区的核心岩藻糖基化。近些年来,很多研究报告表明,去除或降低岩藻糖基化的治疗性抗体不论在体内还是在体外都表现出更高的效能。这主要是由于去除或降低岩藻糖基化的治疗性抗体相对于岩藻糖基化的抗体,可以在更低浓度下通过与FcγRIIIa的高亲和力而表现出较强的ADCC作用。因此,去除或降低岩藻糖基化抗体的应用有望成为提高下一代治疗性抗体效能的有效手段。在此综述中,我们主要讨论了控制治疗性抗体恒定区岩藻糖基化的重要性及当前去除或降低岩藻糖基化治疗性抗体的生产控制方法。     

Complexation of two proteic insect inhibitors to the active site of chymotrypsin suggests decoupled roles for binding and selectivity

The crystal structures of two homologous inhibitors (PMP-C and PMP-D2v) from the insect Locusta migratoria have been determined in complex with bovine alpha-chymotrypsin at 2.1- and 3.0-A resolution, respectively. PMP-C is a potent bovine alpha-chymotrypsin inhibitor whereas native PMP-D2 is a weak inhibitor of bovine trypsin. One unique mutation at the P1 position converts PMP-D2 into a potent bovine alpha-chymotrypsin inhibitor. The two peptides have a similar overall conformation, which consists of a triple-stranded antiparallel beta-sheet connected by three disulfide bridges, thus defining a novel family of serine protease inhibitors. They have in common the protease interaction site, which is composed of the classical protease binding loop (position P5 to P'4, corresponding to residues 26-34) and of an internal segment (residues 15-18), held together by two disulfide bridges. Structural divergences between the two inhibitors result in an additional interaction site between PMP-D2v (position P10 to P6, residues 21-25) and the residues 172-175 of alpha-chymotrypsin. This unusual interaction may be responsible for species selectivity. A careful comparison of data on bound and free inhibitors (from this study and previous NMR studies, respectively) suggests that complexation to the protease stabilizes the flexible binding loop (from P5 to P'4).     

Disulfide bridges remain intact while native insulin converts into amyloid fibrils

Amyloid fibrils are β-sheet-rich protein aggregates commonly found in the organs and tissues of patients with various amyloid-associated diseases. Understanding the structural organization of amyloid fibrils can be beneficial for the search of drugs to successfully treat diseases associated with protein misfolding. The structure of insulin fibrils was characterized by deep ultraviolet resonance Raman (DUVRR) and Nuclear Magnetic Resonance (NMR) spectroscopy combined with hydrogen-deuterium exchange. The compositions of the fibril core and unordered parts were determined at single amino acid residue resolution. All three disulfide bonds of native insulin remained intact during the aggregation process, withstanding scrambling. Three out of four tyrosine residues were packed into the fibril core, and another aromatic amino acid, phenylalanine, was located in the unordered parts of insulin fibrils. In addition, using all-atom MD simulations, the disulfide bonds were confirmed to remain intact in the insulin dimer, which mimics the fibrillar form of insulin.     

Steric zipper formed by hydrophobic peptide fragment of Syrian hamster prion protein

Steric zippers, where the residues of two neighboring β-sheet layers are tightly interdigitated, have been proposed as fundamental structural units of amyloid fibrils by Eisenberg and co-workers. The steric zipper formed by polypeptides containing the palindromic sequence AGAAAAGA has a distinctive feature that the distance between two interdigitated β-sheet layers is comparable to the interstrand distance of the individual β-sheet. This structural motif is of great interest in the study of prion disease because the AGAAAAGA sequence is highly conserved in prion proteins of different species. In this work, the amyloid fibrils formed by the polypeptides of PrP(113-127), viz. Ac-AGAAAAGAVVGGLGG-NH(2), are taken as the model compound to investigate the biophysical principles governing the steric zipper formation. The target fibrils adopt the structural motif of class 7 steric zipper, which is formed by stacking of antiparallel β-sheet layers with residue 117 + k forming backbone hydrogen bonds to residue 120 - k. Implication of our results in the infectivity of scrapie prion is briefly discussed.     

Comparative simulation studies of native and single-site mutant human beta-defensin-1 peptides

Human defensins play important roles in a broad range of biological functions, such as microbial defense and immunity. Yet, little is known about their molecular properties, i.e. secondary structure stability, structural variability, important side chain interactions, surface charge distribution, and resistance to thermal fluctuations, and how these properties are related to their functions. To assess these factors, we studied the native human β-defensin-1 monomer and dimer as well as several single-site mutants using molecular dynamics simulations. The results showed that disulfide bonds are important determinants in maintaining the defensins’ structural integrity, as no structural transitions were observed at 300?K and only minor structural unfolding was detected upon heating to 500?K. The α-helix was less thermally stable than the core β-sheet structure held together by hydrogen bonds and hydrophobic interactions. The monomer α-helix stability was directly correlated, whereas the end-to-end distance was inversely correlated to the experimentally measured β-defensin-1 chemotactic activity, in the order: mutant 2 (Gln24Glu)?>?mutant 3 (Lys31Ala)?=?wild type?>?mutant 1 (Asn4Ala). The structural stability of the β-defensin-1 dimer species exhibited an inverse correlation to their chemotactic activity. In dimers formed by mutants 2 and 3, we observed sliding of one monomer upon the surface of the other in the absence of unbinding. This dynamic sliding feature may enhance the molecular oligomerization of β-defensin-1 peptides contributing to their antibacterial activity. It could also help these peptides orient correctly in the CC chemokine receptor 6 binding site, thereby initiating their chemotactic activity. In agreement with this notion, the remarkable sliding behavior was observed only for the mutants with the highest chemotactic activity.     

Engineering the thermostability of a xylanase from Aspergillus oryzae by an enhancement of the interactions between the N-terminus extension and the β-sheet A2 of the enzyme

A mesophilic Aspergillus oryzae xylanase (AoXyn11A) belongs to glycoside hydrolase family 11. Hydrogen bonds and a disulfide bridge were introduced between the N-terminus extension and the β-sheet A2 of AoXyn11A, which were located in the corresponding region of a hyperthermostable xylanase. The mutants were designated as AoXyn11A^C5 and AoXyn11A^C5–C32, respectively. The thermostabilities of AoXyn11A and the mutants were assessed by the molecular dynamics simulations. After being incubated at 55 °C for 30 min, AoXyn11A^C5–C32 retained 49 % of its original activity, AoXyn11A^C5 retained 12 % and AoXyn11A retained 3 %. The interactions between the N-terminus extension and the β-sheet A2 were analyzed in depth: there was enhancement of the interactions between the N-terminus extension and the β-sheet A2 of AoXyn11A that improved its thermostability.     

Unusually stable β-sheet formation in an ionic self-complementary oligopeptide

A 16-residue amphiphilic oligopeptide (EAK16) with every other residue alanine and also containing glutamic acid and lysine (Ac-NH-AEAEAKAKAEAEAKAK-CONH₂) is able to form an unusually stable β-sheet structure. The β-sheet structure is stable at very low concentrations in water and at high temperatures. Various pH changes at 1.5, 3, 7, and 11 had little effect on the stability of the β-sheet structure. The β-sheet structure was not altered significantly even in the presence of 0.1% SDS, 7 molar guanidine hydrochloride, or 8 molar urea. One of the structural characteristics of the EAK16 is its ionic self-complementarity in that ionic bonds and hydrogen bonds between Glu and Lys can form readily between two oligopeptide β-sheet structures. This structural feature is probably one of the factors that promotes its extreme stability. This is the first example of such an extended ionic self-complementarity in a protein structure. EAK16 and its related peptides may have applications as useful biomaterials. It also offers a good model for studying the mechanism of β-sheet formation. Because the oligopeptide can self-assemble to form a membranous structure, it may have relevance to origin of life research. © 1994 John Wiley & Sons, Inc.     

The effect of dissolved oxygen on the production and the glycosylation profile of recombinant human erythropoietin produced from CHO cells

Human recombinant erythropoietin (rHuEPO) was produced from Chinese hamster ovary (CHO) cells transfected with the human EPO gene. The cells were grown in batch cultures in controlled bioreactors in which the set-points for dissolved oxygen varied between 3% and 200%. The cell-specific growth rate and final cell yield was significantly lower under hyperoxic conditions (200% DO). However, there was no significant difference in growth rates at other oxygen levels compared to control cultures run under a normoxic condition (50% DO). The specific productivity of EPO was significantly lower at a DO set-point of 3% and 200% but maintained a consistently high value between 10% to 100% DO. The EPO produced under all conditions as analyzed by two-dimensional electrophoresis showed a molecular weight range of 33 to 37 kDa and a low isoelectric point range of 3.5 to 5.0. This corresponds to a highly glycosylated and sialylated protein with a profile showing at least seven distinct isoforms. The glycan pattern of isolated samples of EPO was analyzed by weak anion exchange (WAX) HPLC and by normal-phase HPLC incorporating sequential digestion with exoglycosidase arrays. Assigned structures were confirmed by mass spectrometry (MALDI-MS). The most prominent glycan structures were core fucosylated tetranntenary with variable sialylation. However, significant biantennary, triantennary, and non-fucosylated glycans were also identified. Detailed analysis of these glycan structures produced under variable dissolved oxygen levels did not show consistently significant variations except for the ratio of fucosylated to non-fucosylated isoforms. Maximum core fucosylation (80%) was observed at 50% and 100% DO, whereas higher or lower DO levels resulted in reduced fucosylation. This observation of lower fucosylation at high or low DO levels is consistent with previous data reported for glycoprotein production in insect cells.     

Role of invariant Thr80 in human immunodeficiency virus type 1 protease structure, function, and viral infectivity

Sequence variability associated with human immunodeficiency virus type 1 (HIV-1) is useful for inferring structural and/or functional constraints at specific residues within the viral protease. Positions that are invariant even in the presence of drug selection define critically important residues for protease function. While the importance of conserved active-site residues is easily understood, the role of other invariant residues is not. This work focuses on invariant Thr80 at the apex of the P1 loop of HIV-1, HIV-2, and simian immunodeficiency virus protease. In a previous study, we postulated, on the basis of a molecular dynamics simulation of the unliganded protease, that Thr80 may play a role in the mobility of the flaps of protease. In the present study, both experimental and computational methods were used to study the role of Thr80 in HIV protease. Three protease variants (T80V, T80N, and T80S) were examined for changes in structure, dynamics, enzymatic activity, affinity for protease inhibitors, and viral infectivity. While all three variants were structurally similar to the wild type, only T80S was functionally similar. Both T80V and T80N had decreased the affinity for saquinavir. T80V significantly decreased the ability of the enzyme to cleave a peptide substrate but maintained infectivity, while T80N abolished both activity and viral infectivity. Additionally, T80N decreased the conformational flexibility of the flap region, as observed by simulations of molecular dynamics. Taken together, these data indicate that HIV-1 protease functions best when residue 80 is a small polar residue and that mutations to other amino acids significantly impair enzyme function, possibly by affecting the flexibility of the flap domain.     

14-3-3zeta C-terminal stretch changes its conformation upon ligand binding and phosphorylation at Thr232

14-3-3 proteins are abundant binding proteins involved in many biologically important processes. 14-3-3 proteins bind to other proteins in a phosphorylation-dependent manner and function as scaffold molecules modulating the activity of their binding partners. In this work, we studied the conformational changes of 14-3-3 C-terminal stretch, a region implicated in playing a role in the regulation of 14-3-3. Time-resolved fluorescence and molecular dynamics were used to investigate structural changes of the C-terminal stretch induced by phosphopeptide binding and phosphorylation at Thr232, a casein kinase I phosphorylation site located within this region. A tryptophan residue placed at position 242 was exploited as an intrinsic fluorescence probe of the C-terminal stretch dynamics. Other tryptophan residues were mutated to phenylalanine. Time-resolved fluorescence measurements revealed that phosphopeptide binding changes the conformation and increases the flexibility of 14-3-3zeta C-terminal stretch, demonstrating that this region is directly involved in ligand binding. Phosphorylation of 14-3-3zeta at Thr232 resulted in inhibition of phosphopeptide binding and suppression of 14-3-3-mediated enhancement of serotonin N-acetyltransferase activity. Time-resolved fluorescence of Trp242 also revealed that phosphorylation at Thr232 induces significant changes of the C-terminal stretch conformation. In addition, molecular dynamic simulations suggest that phosphorylation at Thr232 induces a more extended conformation of 14-3-3zeta C-terminal stretch and changes its interaction with the rest of the 14-3-3 molecule. These results indicate that the conformation of the C-terminal stretch plays an important role in the regulation of 14-3-3 binding properties.     

Oxidative folding and structural analyses of a Kunitz-related inhibitor and its disulfide intermediates: functional implications

Tick-derived protease inhibitor (TdPI) is a tight-binding Kunitz-related inhibitor of human tryptase β with a unique structure and disulfide-bond pattern. Here we analyzed its oxidative folding and reductive unfolding by chromatographic and disulfide analyses of acid-trapped intermediates. TdPI folds through a stepwise generation of heterogeneous populations of one-disulfide, two-disulfide, and three-disulfide intermediates, with a major accumulation of the nonnative three-disulfide species IIIa. The rate-limiting step of the process is disulfide reshuffling within the three-disulfide population towards a productive intermediate that oxidizes directly into the native four-disulfide protein. TdPI unfolds through a major accumulation of the native three-disulfide species IIIb and the subsequent formation of two-disulfide and one-disulfide intermediates. NMR characterization of the acid-trapped and further isolated IIIa intermediate revealed a highly disordered conformation that is maintained by the presence of the disulfide bonds. Conversely, the NMR structure of IIIb showed a native-like conformation, with three native disulfide bonds and increased flexibility only around the two free cysteines, thus providing a molecular basis for its role as a productive intermediate. Comparison of TdPI with a shortened variant lacking the flexible prehead and posthead segments revealed that these regions do not contribute to the protein conformational stability or the inhibition of trypsin but are important for both the initial steps of the folding reaction and the inhibition of tryptase β. Taken together, the results provide insights into the mechanism of oxidative folding of Kunitz inhibitors and pave the way for the design of TdPI variants with improved properties for biomedical applications.     

Conversion of trypsin to a functional threonine protease

The hydroxyl group of a serine residue at position 195 acts as a nucleophile in the catalytic mechanism of the serine proteases. However, the chemically similar residue, threonine, is rarely used in similar functional context. Our structural modeling suggests that the Ser 195 --> Thr trypsin variant is inactive due to negative steric interaction between the methyl group on the beta-carbon of Thr 195 and the disulfide bridge formed by cysteines 42 and 58. By simultaneously truncating residues 42 and 58 and substituting Ser 195 with threonine, we have successfully converted the classic serine protease trypsin to a functional threonine protease. Substitution of residue 42 with alanine and residue 58 with alanine or valine in the presence of threonine 195 results in trypsin variants that are 10(2) -10(4) -fold less active than wild type in kcat/KM but >10(6)-fold more active than the Ser 195 --> Thr single variant. The substitutions do not alter the substrate specificity of the enzyme in the P1'- P4' positions. Removal of the disulfide bridge decreases the overall thermostability of the enzyme, but it is partially rescued by the presence of threonine at position 195.