Dissecting structural basis of the unique substrate selectivity of human enteropeptidase catalytic subunit

Enteropeptidase is a key enzyme in the digestion system of higher animals. It initiates enzymatic cascade cleaving trypsinogen activation peptide after a unique sequence DDDDK. Recently, we have found specific activity of human enteropeptidase catalytic subunit (L-HEP) being significantly higher than that of its bovine ortholog (L-BEP). Moreover, we have discovered that L-HEP hydrolyzed several nonspecific peptidic substrates. In this work, we aimed to further characterize species-specific enteropeptidase activities and to reveal their structural basis. First, we compared hydrolysis of peptides and proteins lacking DDDDK sequence by L-HEP and L-BEP. In each case human enzyme was more efficient, with the highest hydrolysis rate observed for substrates with a large hydrophobic residue in P2-position. Computer modeling suggested enzyme exosite residues 96 (Arg in L-HEP, Lys in L-BEP) and 219 (Lys in L-HEP, Gln in L-BEP) to be responsible for these differences in enteropeptidase catalytic activity. Indeed, human-to-bovine mutations Arg96Lys, Lys219Gln shifted catalytic properties of L-HEP toward those of L-BEP. This effect was amplified in case of the double mutation Arg96Lys/Lys219Gln, but still did not cover the full difference in catalytic activities of human and bovine enzymes. To find a missing link, we studied monopeptide benzyl-arginine-β-naphthylamide hydrolysis. L-HEP catalyzed it with an order lower K _m than L-BEP, suggesting the monopeptide-binding S1 site input into catalytic distinction between two enteropeptidase species. Together, our findings suggest structural basis of the unique catalytic properties of human enteropeptidase and instigate further studies of its tentative physiological and pathological roles.     

Conversion of the catalytic specificity of alanine racemase to a -amino acid aminotransferase activity by a double active-site mutation

Alanine racemase depending on pyridoxal 5′-phosphate catalyzes the interconversion between - and -alanine. The enzyme from Bacillus stearothermophilus catalyzes the transamination as a side reaction with both substrates once per 3×10⁷ times of the racemization. In this work, we studied the effects of the mutation of Arg219, and that of Arg219 and Tyr265 on the catalysis of Bacillus alanine racemase. Arg219 interacting with pyridinium nitrogen of the cofactor is conserved in all alanine racemases. The corresponding residue of aminotransferases is an acidic residue, such as glutamate or aspartate. Mutation of Arg219 to a glutamyl residue resulted in a 5.4-fold increase in the forward half transamination activity with -alanine and a 10³-fold decrease in the racemase activity. The double mutation, Arg219→Glu and Tyr265→Ala, completely abolished the racemase activity and increased the forward half transaminase activity 6.6-fold. Arg219 is one of the structural determinants of the catalytic specificity of the alanine racemase.     

Dissecting structural basis of the unique substrate selectivity of human enteropeptidase catalytic subunit

Enteropeptidase is a key enzyme in the digestion system of higher animals. It initiates enzymatic cascade cleaving trypsinogen activation peptide after a unique sequence DDDDK. Recently, we have found specific activity of human enteropeptidase catalytic subunit (L-HEP) being significantly higher than that of its bovine ortholog (L-BEP). Moreover, we have discovered that L-HEP hydrolyzed several nonspecific peptidic substrates. In this work, we aimed to further characterize species-specific enteropeptidase activities and to reveal their structural basis. First, we compared hydrolysis of peptides and proteins lacking DDDDK sequence by L-HEP and L-BEP. In each case human enzyme was more efficient, with the highest hydrolysis rate observed for substrates with a large hydrophobic residue in P2-position. Computer modeling suggested enzyme exosite residues 96 (Arg in L-HEP, Lys in L-BEP) and 219 (Lys in L-HEP, Gln in L-BEP) to be responsible for these differences in enteropeptidase catalytic activity. Indeed, human-to-bovine mutations Arg96Lys, Lys219Gln shifted catalytic properties of L-HEP toward those of L-BEP. This effect was amplified in case of the double mutation Arg96Lys/Lys219Gln, but still did not cover the full difference in catalytic activities of human and bovine enzymes. To find a missing link, we studied monopeptide benzyl-arginine-β-naphthylamide hydrolysis. L-HEP catalyzed it with an order lower K (m) than L-BEP, suggesting the monopeptide-binding S1 site input into catalytic distinction between two enteropeptidase species. Together, our findings suggest structural basis of the unique catalytic properties of human enteropeptidase and instigate further studies of its tentative physiological and pathological roles.     

Structure and mechanistic implications of a uroporphyrinogen III synthase-product complex

Uroporphyrinogen III synthase (U3S) catalyzes the asymmetrical cyclization of a linear tetrapyrrole to form the physiologically relevant uroporphyrinogen III (uro'gen III) isomer during heme biosynthesis. Here, we report four apoenzyme and one product complex crystal structures of the Thermus thermophilus (HB27) U3S protein. The overlay of eight crystallographically unique U3S molecules reveals a huge range of conformational flexibility, including a "closed" product complex. The product, uro'gen III, binds between the two domains and is held in place by a network of hydrogen bonds between the product's side chain carboxylates and the protein's main chain amides. Interactions of the product A and B ring carboxylate side chains with both structural domains of U3S appear to dictate the relative orientation of the domains in the closed enzyme conformation and likely remain intact during catalysis. The product C and D rings are less constrained in the structure, consistent with the conformational changes required for the catalytic cyclization with inversion of D ring orientation. A conserved tyrosine residue is potentially positioned to facilitate loss of a hydroxyl from the substrate to initiate the catalytic reaction.     

Crystal structure of human uroporphyrinogen III synthase.

Uroporphyrinogen III synthase, U3S, the fourth enzyme in the porphyrin biosynthetic pathway, catalyzes cyclization of the linear tetrapyrrole, hydroxymethylbilane, to the macrocyclic uroporphyrino gen III, which is used in several different pathways to form heme, siroheme, chlorophyll, F(430) and vitamin B(12). U3S activity is essential in all organisms, and decreased activity in humans leads to the autosomal recessive disorder congenital erythropoetic porphyria. We have determined the crystal structure of recombinant human U3S at 1.85 A resolution. The protein folds into two alpha/beta domains connected by a beta-ladder. The active site appears to be located between the domains, and variations in relative domain positions observed between crystallographically independent molecules indicates the presence of flexibility that may be important in the catalytic cycle. Possible mechanisms of catalysis were probed by mutating each of the four invariant residues in the protein that have titratable side chains. Additionally, six other highly conserved and titratable side chains were also mutated. In no case, however, did one of these mutations abolish enzyme activity, suggesting that the mechanism does not require acid/base catalysis.     

Dissecting the cAMP-inducible allosteric switch in protein kinase A RI��

The regulatory subunits of cAMP‐dependent protein kinase (PKA) are the major receptors of cAMP in most eukaryotic cells. As the cyclic nucleotide binding (CNB) domains release cAMP and bind to the catalytic subunit of PKA, they undergo a major conformational change. The change is mediated by the B/C helix in CNB‐A, which extends into one long helix that now separates the two CNB domains and docks onto the surface of the catalytic subunit. We explore here the role of three key residues on the B/C helix that dock onto the catalytic subunit, Arg226, Leu233, and Met 234. By replacing each residue with Ala, we show that each contributes significantly to creating the R:C interface. By also deleting the second CNB domain (CNB‐B), we show furthermore that CNB‐B is a critical part of the cAMP‐induced conformational switch that dislodges the B/C helix from the surface of the catalytic subunit. Without CNB‐B the K_a for activation by cAMP increases from 80 to 1000 nM. Replacing any of the key interface residues with Ala reduces the K_a to 25–40 nM. Leu233 and M234 contribute to a hydrophobic latch that binds the B/C helix onto the large lobe of the C‐subunit, while Arg226 is part of an electrostatic switch that couples the B/C helix to the phosphate binding cassette where the cAMP docks.     

Orthocaspases are proteolytically active prokaryotic caspase homologues: the case of Microcystis aeruginosa

Caspases are a family of cysteine‐dependent proteases known to be involved in the process of programmed cell death in metazoans. Recently, cyanobacteria were also found to contain caspase‐like proteins, but their existence has only been identified in silico up to now. Here, we present the first experimental characterisation of a prokaryotic caspase homologue. We have expressed the putative caspase‐like gene MaOC1 from the toxic bloom‐forming cyanobacterium Microcystis aeruginosa PCC 7806 in Escherichia coli. Kinetic characterisation showed that MaOC1 is an endopeptidase with a preference for arginine in the P1 position and a pH optimum of 7.5. MaOC1 exhibited high catalytic rates with the k_cat/K_M value for Z‐RR‐AMC substrate of the order 10⁶ M^?1 s^?1. In contrast to plant or metazoan caspase‐like proteins, whose activity is calcium‐dependent or requires dimerisation for activation, MaOC1 was activated by autocatalytic processing after residue Arg219, which separated the catalytic domain and the remaining 55 kDa subunit. The Arg219Ala mutant was resistant to autoprocessing and exhibited no proteolytic activity, confirming that processing of MaOC1 is a prerequisite for its activity. Due to their structural and functional differences to other known caspase‐like proteins, we suggest to name these evolutionary primitive proteins orthocaspases.     

Temperature-sensitive mutants and revertants in the coronavirus nonstructural protein 5 protease (3CLpro) define residues involved in long-distance communication and regulation of protease activity

Positive-strand RNA virus genomes are translated into polyproteins that are processed by viral proteases to yield functional intermediate and mature proteins. Coronaviruses (CoVs) carry genes that encode an nsp5 protease (also known as 3CLpro or Mpro) responsible for 11 maturation cleavages. The nsp5 structure contains two chymotrypsin-like domains (D1 and D2) and a unique domain (D3), and forms functional dimers. However, little is known of interactions or communication across the structure of the protease during nsp5 activity. Using reverse genetic mutagenesis of the CoV murine hepatitis virus (MHV) nsp5, we identified a new temperature-sensitive (ts) mutation in D2 of nsp5 (Ser133Ala) and confirmed a ts residue in D3 (Phe219Leu). Both D2-tsS133A and D3-tsF219L were impaired for viral replication and nsp5-mediated polyprotein processing at the nonpermissive temperature. Passage of tsS133A and tsF219L at the nonpermissive temperature resulted in emergence of multiple second-site suppressor mutations, singly and in combinations. Among the second-site mutations, a D2 His134Tyr change suppressed the ts phenotype of D2-tsS133A and D3-tsF219L, as well as the previously reported D2-tsV148A. Analysis of multiple CoV nsp5 structures, and alignment of nonredundant nsp5 primary sequences, demonstrated that ts and suppressor residues are not conserved across CoVs and are physically distant (>10 Å) from each other, from catalytic and substrate-binding residues, and from the nsp5 dimer interface. These findings demonstrate that long-distance communication pathways between multiple residues and domains of nsp5 play a significant role in nsp5 activity and viral replication, suggesting possible novel targets for non-active site inhibitors of nsp5.     

Molecular dynamics simulations of various coronavirus main proteinases

In this study, two homology models (denoted as MproST and MproSH) of main proteinase (Mpro) from the novel coronavirus associated with severe acute respiratory syndrome (SARS-CoV) were constructed based on the crystal structures of Mpro from transmissible gastroenteritis coronavirus (TGEV) (MproT) and human coronavirus HcoV-229E (MproH), respectively. Both MproST and MproSH exhibit similar folds as their respective template proteins. These homology models reveal three distinct functional domains as well as an intervening loop connecting domains II and III as found in both template proteins. A catalytic cleft containing the substrate binding sites S1 and S2 between domains I and II are also observed. S2 undergoes more significant structural fluctuation than S1 during the 400 ps molecular dynamics simulations because it is located at the open mouth of the catalytic cleft, while S1 is situated in the very bottom of this cleft. The thermal unfolding of these proteins begins at domain III, where the structure is least conserved among these proteins. Mpro may still maintain its proteolytic activity while it is partially unfolded. The electrostatic interaction between Arg40 and Asp186 plays an important role in maintaining the structural integrity of both S1 and S2.     

A Comprehensive Alanine-Scanning Mutagenesis Study Reveals Roles for Salt Bridges in the Structure and Activity of Pseudomonas aeruginosa Elastase

The relationship between salt bridges and stability/enzymatic activity is unclear. We studied this relationship by systematic alanine-scanning mutation analysis using the typical M4 family metalloprotease Pseudomonas aeruginosa elastase (PAE, also known as pseudolysin) as a model. Structural analysis revealed seven salt bridges in the PAE structure. We constructed ten mutants for six salt bridges. Among these mutants, six (Asp189Ala, Arg179Ala, Asp201Ala, Arg205Ala, Arg245Ala and Glu249Ala) were active and four (Asp168Ala, Arg198Ala, Arg253Ala, and Arg279Ala) were inactive. Five mutants were purified, and their catalytic efficiencies (k _cat/K _m), half-lives (t _1/2) and thermal unfolding curves were compared with those of PAE. Mutants Asp189Ala and Arg179Ala both showed decreased thermal stabilities and increased activities, suggesting that the salt bridge Asp189-Arg179 stabilizes the protein at the expense of catalytic efficiency. In contrast, mutants Asp201Ala and Arg205Ala both showed slightly increased thermal stability and slightly decreased activity, suggesting that the salt bridge Asp201-Arg205 destabilizes the protein. Mutant Glu249Ala is related to a C-terminal salt bridge network and showed both decreased thermal stability and decreased activity. Furthermore, Glu249Ala showed a thermal unfolding curve with three discernable states [the native state (N), the partially unfolded state (I) and the unfolded state (U)]. In comparison, there were only two discernable states (N and U) in the thermal unfolding curve of PAE. These results suggest that Glu249 is important for catalytic efficiency, stability and unfolding cooperativity. This study represents a systematic mutational analyses of salt bridges in the model metalloprotease PAE and provides important insights into the structure-function relationship of enzymes.     

Functional analyses of the chitin-binding domains and the catalytic domain of <Emphasis Type="Italic">Brassica juncea</Emphasis> chitinase BjCHI1

We previously isolated a Brassica juncea cDNA encoding BjCHI1, a novel chitinase with two chitin-binding domains. Synthesis of its mRNA is induced by wounding, methyl jasmonate treatment, Aspergillus niger infection and caterpillar Pieris rapae feeding, suggesting that the protein has a role in defense. In that it possesses two chitin-binding domains, BjCHI1 resembles the precursor of Urtica dioica agglutinin but unlike that protein, BjCHI1 retains its chitinase catalytic domain after post-translational processing. To explore the properties of multi-domain BjCHI1, we have expressed recombinant BjCHI1 and two derivatives, which lack one (BjCHI2) or both (BjCHI3) chitin-binding domains, as secreted proteins in Pichia pastoris. Recombinant BjCHI1 and BjCHI2, showed apparent molecular masses on SDS-PAGE larger than calculated, and could be deglycosylated using -mannosidase. Recombinant BjCHI3, without the proline/threonine-rich linker region containing predicted O-glycosylation sites, did not appear to be processed by -mannosidase. BjCHI1s ability to agglutinate rabbit erythrocytes is unique among known chitinases. Both chitin-binding domains are essential for agglutination; this property is absent in recombinant BjCHI2 and BjCHI3. To identify potential catalytic residues, we generated site-directed mutations in recombinant BjCHI3. Mutation E212A showed the largest effect, exhibiting 0 of wild-type specific activity. H211N and R361A resulted in considerable (>91) activity loss, implying these charged residues are also important in catalysis. E234A showed 36 retention of activity and substitution Y269D, 50. The least affected mutants were E349A and D360A, with 73 and 68 retention, respectively. Like Y269, E349 and D360 are possibly involved in substrate binding rather than catalysis.     

Characterization of the catalytic domains of <Emphasis Type="Italic">Trichoderma reesei</Emphasis> endoglucanase I,II, and III,expressed in <Emphasis Type="Italic">Escherichia coli</Emphasis>

The genes encoding the catalytic domains (CD) of the three endoglucanases (EG I; Cel7B, EG II; Cel5A, and EG III; Cel12A) from Trichoderma reesei QM9414 were expressed in Escherichia coli strains Rosetta-gami B (DE3) pLacI or Origami B (DE3) pLacI and were found to produce functional intracellular proteins. Protein production by the three endoglucanase transformants was evaluated as a function of growth temperature. Maximal productivity of EG I-CD at 15°C, EG II-CD at 20°C and EG III at 37°C resulted in yields of 6.9, 72, and 50 mg/l, respectively. The endoglucanases were purified using a simple purification method based on removing E. coli proteins by isoelectric point precipitation. Specific activity toward carboxymethyl cellulose was found to be 65, 49, and 15 U/mg for EG I-CD, EG II-CD, and EG III, respectively. EG II-CD was able to cleave 1,3–1,4-β-d-glucan and soluble cellulose derivatives. EG III was found to be active against cellulose, 1,3–1,4-β-d-glucan and xyloglucan, while EG I-CD was active against cellulose, 1,3–1,4-β-d-glucan, xyloglucan, xylan, and mannan.     

Domain Requirements of the JIL-1 Tandem Kinase for Histone H3 Serine 10 Phosphorylation and Chromatin Remodeling in Vivo

The JIL-1 kinase localizes to Drosophila polytene chromosome interbands and phosphorylates histone H3 at interphase, counteracting histone H3 lysine 9 dimethylation and gene silencing. JIL-1 can be divided into four main domains, including an NH₂-terminal domain, two separate kinase domains, and a COOH-terminal domain. In this study, we characterize the domain requirements of the JIL-1 kinase for histone H3 serine 10 (H3S10) phosphorylation and chromatin remodeling in vivo. We show that a JIL-1 construct without the NH₂-terminal domain is without H3S10 phosphorylation activity despite the fact that it localizes properly to polytene interband regions and that it contains both kinase domains. JIL-1 is a double kinase, and we demonstrate that both kinase domains of JIL-1 are required to be catalytically active for H3S10 phosphorylation to occur. Furthermore, we provide evidence that JIL-1 is phosphorylated at serine 424 and that this phosphorylation is necessary for JIL-1 H3S10 phosphorylation activity. Thus, these data are compatible with a model where the NH₂-terminal domain of JIL-1 is required for chromatin complex interactions that position the kinase domain(s) for catalytic activity in the context of the state of higher order nucleosome packaging and chromatin structure and where catalytic H3S10 phosphorylation activity mediated by the first kinase domain is dependent on autophosphorylation of serine 424 by the second kinase domain. Furthermore, using a lacO repeat tethering system to target mutated JIL-1 constructs with or without catalytic activity, we show that the epigenetic H3S10 phosphorylation mark itself functions as a causative regulator of chromatin structure independently of any structural contributions from the JIL-1 protein.     

Structure of bacteriophage SPN1S endolysin reveals an unusual two‐module fold for the peptidoglycan lytic and binding activity

Bacteriophage SPN1S infects the pathogenic Gram‐negative bacterium Salmonella typhimurium and expresses endolysin for the release of phage progeny by degrading peptidoglycan of the host cell walls. Bacteriophage SPN1S endolysin exhibits high glycosidase activity against peptidoglycans, resulting in antimicrobial activity against a broad range of outer membrane‐permeabilized Gram‐negative bacteria. Here, we report a crystal structure of SPN1S endolysin, indicating that unlike most endolysins from Gram‐negative bacteria background, the α‐helical protein consists of two modular domains, a large and a small domain, with a concave groove between them. Comparison with other structurally homologous glycoside hydrolases indicated a possible peptidoglycan binding site in the groove, and the presence of a catalytic dyad in the vicinity of the groove, one residue in a large domain and the other in a junction between the two domains. The catalytic dyad was further validated by antimicrobial activity assay against outer membrane‐permeabilized Escherichia coli. The three‐helix bundle in the small domain containing a novel class of sequence motif exhibited binding affinity against outer membrane‐permeabilized E. coli and was therefore proposed as the peptidoglycan‐binding domain. These structural and functional features suggest that endolysin from a Gram‐negative bacterial background has peptidoglycan‐binding activity and performs glycoside hydrolase activity through the catalytic dyad.     

Structural insights into the cofactor-assisted substrate recognition of yeast quinone oxidoreductase Zta1

Quinone oxidoreductase (QOR EC1.6.5.5) catalyzes the reduction of quinone to hydroxyquinone using NADPH as a cofactor. Here we present the crystal structure of the ζ-crystallin-like QOR Zta1 from Saccharomycescerevisiae in apo-form at 2.00 Å and complexed with NADPH at 1.59 Å resolution. Zta1 forms a homodimer, with each subunit containing a catalytic and a cofactor-binding domain. Upon NADPH binding to the interdomain cleft, the two domains shift towards each other, producing a better fit for NADPH, and tightening substrate binding. Computational simulation combined with site-directed mutagenesis and enzymatic activity analysis defined a potential quinone-binding site that determines the stringent substrate specificity. Moreover, multiple-sequence alignment and kinetics assays implied that a single-residue change from Arg in lower organisms to Gly in vertebrates possibly resulted in elevation of enzymatic activity of ζ-crystallin-like QORs throughout evolution.     

Development of a biomarker for Geobacter activity and strain composition; Proteogenomic analysis of the citrate synthase protein during bioremediation of U(VI)

Monitoring the activity of target microorganisms during stimulated bioremediation is a key problem for the development of effective remediation strategies. At the US Department of Energy's Integrated Field Research Challenge (IFRC) site in Rifle, CO, the stimulation of Geobacter growth and activity via subsurface acetate addition leads to precipitation of U(VI) from groundwater as U(IV). Citrate synthase (gltA) is a key enzyme in Geobacter central metabolism that controls flux into the TCA cycle. Here, we utilize shotgun proteomic methods to demonstrate that the measurement of gltA peptides can be used to track Geobacter activity and strain evolution during in situ biostimulation. Abundances of conserved gltA peptides tracked Fe(III) reduction and changes in U(VI) concentrations during biostimulation, whereas changing patterns of unique peptide abundances between samples suggested sample‐specific strain shifts within the Geobacter population. Abundances of unique peptides indicated potential differences at the strain level between Fe(III)‐reducing populations stimulated during in situ biostimulation experiments conducted a year apart at the Rifle IFRC. These results offer a novel technique for the rapid screening of large numbers of proteomic samples for Geobacter species and will aid monitoring of subsurface bioremediation efforts that rely on metal reduction for desired outcomes.     

Improving the catalytic activity of a detergent-compatible serine protease by rational design

Serine proteases are among the most important biological additives in various industries such as detergents, leather, animal feed and food. A serine protease gene, Fgapt4, from Fusarium graminearum 2697 was identified, cloned and expressed in Pichia pastoris. The optimal pH and temperature of FgAPT4 were 8.5 and 40°C, respectively. The relative activity was >30% even at 10°C. It had a wide range of pH stability (4.0–12.0) and detergent compatibility. To improve the catalytic activity, a strategy combining molecular docking and evolutionary analysis was adopted. Twelve amino acid residue sites and three loops (A, B and C) were selected as potential hot spots that might play critical roles in the enzyme's functional properties. Twenty-eight mutants targeting changes in individual sites or loops were designed, and mutations with good performance were combined. The best mutant was FgAPT4-M3 (Q70N/D142S/A143S/loop C). The specific activity and catalytic efficiency of FgAPT4-M3 increased by 1.6 (1008.5 vs. 385.9 U/mg) and 2.2-fold (3565.1 vs. 1106.3/s/mM), respectively. Computational analyses showed that the greater flexibility of the substrate pocket may be responsible for the increased catalytic activity. In addition, its application in detergents indicated that FgAPT4-M3 has great potential in washing.     

Design of peptidase‐resistant peptide inhibitors of myosin light chain kinase

Myosin light chain kinase (MLCK) is a key regulator of various forms of cell motility including smooth muscle contraction, cell migration, cytokinesis, receptor capping, secretion, etc. Inhibition of MLCK activity in endothelial and epithelial monolayers using cell‐permeant peptide Arg‐Lys‐Lys‐Tyr‐Lys‐Tyr‐Arg‐Arg‐Lys (PIK, P eptide I nhibitor of K inase) allows protecting the barrier capacity, suggesting a potential medical use of PIK. However, low stability of L ‐PIK in a biological milieu prompts for development of more stable L ‐PIK analogues for use as experimental tools in basic and drug‐oriented biomedical research. Previously, we designed PIK1, H‐(N^αMe)Arg‐Lys‐Lys‐Tyr‐Lys‐Tyr‐Arg‐Arg‐Lys‐NH₂, that was 2.5‐fold more resistant to peptidases in human plasma in vitro than L ‐PIK and equal to it as MLCK inhibitor. In order to further enhance proteolytic stability of PIK inhibitor, we designed the set of six site‐protected peptides based on L ‐PIK and PIK1 degradation patterns in human plasma as revealed by ¹H‐NMR analysis. Implemented modifications increased half‐live of the PIK‐related peptides in plasma about 10‐fold, and these compounds retained 25–100% of L ‐PIK inhibitory activity toward MLCK in vitro. Based on stability and functional activity ranking, PIK2, H‐(N^αMe)Arg‐Lys‐Lys‐Tyr‐Lys‐Tyr‐Arg‐D ‐Arg‐Lys‐NH₂, was identified as the most stable and effective L ‐PIK analogue. PIK2 was able to decrease myosin light chain phosphorylation in endothelial cells stimulated with thrombin, and this effect correlated with the inhibition by PIK2 of thrombin‐induced endothelial hyperpermeability in vitro. Therefore, PIK2 could be used as novel alternative to other cell‐permeant inhibitors of MLCK in cell culture‐based and in vivo studies where MLCK catalytic activity inhibition is required. Copyright © 2016 European Peptide Society and John Wiley & Sons, Ltd.