The 0.78 A structure of a serine protease: Bacillus lentus subtilisin

Structural Catalytic Core in Subtilisin-like Proteins and Its Comparison to Trypsin-like Serine Proteases and Alpha/Beta-Hydrolases

Subtilisin-like proteins are serine proteases that use two types of catalytic triads: Ser-His-Asp and Ser-Glu-Asp. Here, we investigate the two known families of subtilisin-like proteins, the subtilases (Ser-His-Asp triad) and the serine-carboxyl proteinases (Ser-Glu-Asp triad), and describe the local structural arrangements (cores) that govern the catalytic residues in these proteins. We show the separation of the cores into conserved structural zones, which can be repeatedly found in different structures, and compare the structural cores in subtilisin-like proteins with those in trypsin-like serine proteases and alpha/beta-hydrolases.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Novel Detection Method for Evaluating the Activity of an Alkaline Serine Protease from Bacillus clausii

Until now, the detection methods for serine proteases have been quite time-consuming or cannot indicate the "real" protease activity. Here, a rapid and simple method for determining the "real" activity of serine proteases toward AAPX (a kind of mixed polypeptide substrates, with X representing 20 standard amino acids) was developed. This AAPX method has high reliability, sensitivity, and repeatability and can be used for detecting the serine protease activity spectrophotometrically. Additionally, the site-directed saturation mutagenesis library of alkaline serine protease PRO (BcPRO) from Bacillus clausii was screened with this AAPX method. Three beneficial mutants S99R, S99H, and S99W were identified, and S99W displayed the highest activity. In comparison to wild-type BcPRO, S99W exhibited enhanced catalytic performance toward eight AAPX monomers, and the molecular dynamics simulation revealed the mechanism responsible for its improved activity toward AAPM. Consequently, this work provides an efficient method for detecting, characterizing, mining, and high-throughput screening of serine proteases.

Improving the activity and stability of Bacillus clausii alkaline protease using directed evolution and molecular dynamics simulation

Detergent enzymes have been developed extensively as eco-friendly substitutes for the harmful chemicals in detergent. Among them, alkaline protease accounts for a large share of detergent enzyme sales. Thus, improving the specific activity of alkaline protease could play an important role in reducing the cost of detergent enzymes. In our study, alkaline protease from Bacillus clausii (PRO) was used to construct a mutant library through directed evolution using error-prone PCR, and a variant (G95P) with 9-fold enhancement in specific activity was obtained. After incubation at a pH of 11.0 for 70 h, G95P maintained 67 % of its maximal activity, which was 46 % more than wild-type PRO (WT), indicating that G95P has better alkaline stability than WT. The thermostability of G95P was also enhanced, as G95P achieved 17 % initial activity after incubation for 50 h at 60 °C, while WT lost its activity. The MD simulation results verified that variant G95P possessed improved stability of its Gly95-Gly100 loop region and Arg19-Asp265 salt bridge, leading to enhanced stability and catalytic efficiency. This work enhances the understanding of the structure-function relationship of PRO and provides an academic foundation for improving the enzymatic properties of PRO to satisfy industrial requirements using protein engineering.

On the Case of the Misplaced Hydrogens

Few other elements play a more central role in biology than hydrogen. The interactions, bonding and movement of hydrogen atoms are central to biological catalysis, structure and function. Yet owing to the elusive nature of a single hydrogen atom few experimental and computational techniques can precisely determine its location. This is exemplified in short hydrogen bonds (SHBs) where the location of the hydrogen atom is indicative of the underlying strength of the bonds, which can vary from 1-5 kcal/mol in canonical hydrogen bonds, to an almost covalent nature in single-well hydrogen bonds. Owing to the often-times inferred position of hydrogen, the role of SHBs in biology has remained highly contested and debated. This has also led to discrepancies in computational, biochemical and structural studies of proteins thought to use SHBs in performing chemistry and stabilizing interactions. Herein, we discuss in detail two distinct examples, namely the conserved catalytic triad and the photoreceptor, photoactive yellow protein, where studies of these SHB-containing systems have permitted contextualization of the role these unique hydrogen bonds play in biology.

Low barrier hydrogen bonds in protein structure and function

Low-barrier hydrogen bonds (LBHBs) are a special type of short hydrogen bond (HB) that is characterized by the equal sharing of a hydrogen atom. The existence and catalytic role of LBHBs in proteins has been intensely contested. Advancements in X-ray and neutron diffraction methods has revealed delocalized hydrogen atoms involved in potential LBHBs in a number of proteins, while also demonstrating that short HBs are not necessarily LBHBs. More importantly, a series of experiments on ketosteroid isomerase (KSI) have suggested that LBHBs are significantly stronger than standard HBs in the protein microenvironment in terms of enthalpy, but not free energy. The discrepancy between the enthalpy and free energy of LBHBs offers clues to the challenges, and potential solutions, of the LBHB debate, where the unique strength of LBHBs plays a special role in the kinetic processes of enzyme function and structure, together with other molecular forces in a pre-organized environment.

Combination of coarse-grained molecular dynamics simulations and small-angle X-ray scattering experiments

The combination of molecular dynamics (MD) simulations and small-angle X-ray scattering (SAXS), called the MD-SAXS method, is efficient for investigating protein dynamics. To overcome the time-scale limitation of all-atom MD simulations, coarse-grained (CG) representations are often utilized for biomolecular simulations. In this study, we propose a method to combine CG MD simulations with SAXS, termed the CG-MD-SAXS method. In the CG-MD-SAXS method, the scattering factors of CG particles for proteins and nucleic acids are evaluated using high-resolution structural data in the Protein Data Bank, and the excluded volume and the hydration shell are modeled using two adjustable parameters to incorporate solvent effects. To avoid overfitting, only the two parameters are adjusted for an entire structure ensemble. To verify the developed method, theoretical SAXS profiles for various proteins, DNA/RNA, and a protein-RNA complex are compared with both experimental profiles and theoretical profiles obtained by the all-atom representation. In the present study, we applied the CG-MD-SAXS method to the Swi5-Sfr1 complex and three types of nucleosomes to obtain reliable ensemble models consistent with the experimental SAXS data.

Structural complementarity of distance constraints obtained from chemical cross-linking and amino acid coevolution

The analysis of amino acid coevolution has emerged as a practical method for protein structural modeling by providing structural contact information from alignments of amino acid sequences. In parallel, chemical cross-linking/mass spectrometry (XLMS) has gained attention as a universally applicable method for obtaining low-resolution distance constraints to model the quaternary arrangements of proteins, and more recently even protein tertiary structures. Here, we show that the structural information obtained by XLMS and coevolutionary analysis are effectively complementary: the distance constraints obtained by each method are almost exclusively associated with non-coincident pairs of residues, and modeling results obtained by the combination of both sets are improved relative to considering the same total number of constraints of a single type. The structural rationale behind the complementarity of the distance constraints is discussed and illustrated for a representative set of proteins with different sizes and folds.

Precise Probing of Residue Roles by Post-Translational β,γ-C,N Aza-Michael Mutagenesis in Enzyme Active Sites

Biomimicry valuably allows the understanding of the essential chemical components required to recapitulate biological function, yet direct strategies for evaluating the roles of amino acids in proteins can be limited by access to suitable, subtly-altered unnatural variants. Here we describe a strategy for dissecting the role of histidine residues in enzyme active sites using unprecedented, chemical, post-translational side-chain-β,γ C-N bond formation. Installation of dehydroalanine (as a "tag") allowed the testing of nitrogen conjugate nucleophiles in "aza-Michael"-1,4-additions (to "modify"). This allowed the creation of a regioisomer of His (iso-His, His_iso) linked instead through its pros-Nπ atom rather than naturally linked via C4, as well as an aza-altered variant aza-His_iso. The site-selective generation of these unnatural amino acids was successfully applied to probe the contributing roles (e.g., size, H-bonding) of His residues toward activity in the model enzymes subtilisin protease from Bacillus lentus and Mycobacterium tuberculosis pantothenate synthetase.

Neutron Crystallography for the Study of Hydrogen Bonds in Macromolecules

The hydrogen bond (H bond) is one of the most important interactions that form the foundation of secondary and tertiary protein structure. Beyond holding protein structures together, H bonds are also intimately involved in solvent coordination, ligand binding, and enzyme catalysis. The H bond by definition involves the light atom, H, and it is very difficult to study directly, especially with X-ray crystallographic techniques, due to the poor scattering power of H atoms. Neutron protein crystallography provides a powerful, complementary tool that can give unambiguous information to structural biologists on solvent organization and coordination, the electrostatics of ligand binding, the protonation states of amino acid side chains and catalytic water species. The method is complementary to X-ray crystallography and the dynamic data obtainable with NMR spectroscopy. Also, as it gives explicit H atom positions, it can be very valuable to computational chemistry where exact knowledge of protonation and solvent orientation can make a large difference in modeling. This article gives general information about neutron crystallography and shows specific examples of how the method has contributed to structural biology, structure-based drug design; and the understanding of fundamental questions of reaction mechanisms.

Atomic resolution structure of serine protease proteinase K at ambient temperature

Atomic resolution structures (beyond 1.20 Å) at ambient temperature, which is usually hampered by the radiation damage in synchrotron X-ray crystallography (SRX), will add to our understanding of the structure-function relationships of enzymes. Serial femtosecond crystallography (SFX) has attracted surging interest by providing a route to bypass such challenges. Yet the progress on atomic resolution analysis with SFX has been rather slow. In this report, we describe the 1.20 Å resolution structure of proteinase K using 13 keV photon energy. Hydrogen atoms, water molecules, and a number of alternative side-chain conformations have been resolved. The increase in the value of B-factor in SFX suggests that the residues and water molecules adjacent to active sites were flexible and exhibited dynamic motions at specific substrate-recognition sites.

Open and closed states of Candida antarctica lipase B: protonation and the mechanism of interfacial activation

Lipases (EC 3.1.1.3) are ubiquitous hydrolases for the carboxyl ester bond of water-insoluble substrates, such as triacylglycerols, phospholipids, and other insoluble substrates, acting in aqueous as well as in low-water media, thus being of considerable physiological significance with high interest also for their industrial applications. The hydrolysis reaction follows a two-step mechanism, or “interfacial activation,” with adsorption of the enzyme to a heterogeneous interface and subsequent enhancement of the lipolytic activity. Among lipases, Candida antarctica lipase B (CALB) has never shown any significant interfacial activation, and a closed conformation of CALB has never been reported, leading to the conclusion that its behavior was due to the absence of a lid regulating the access to the active site. The lid open and closed conformations and their protonation states are observed in the crystal structure of CALB at 0.91 Å resolution. Having the open and closed states at atomic resolution allows relating protonation to the conformation, indicating the role of Asp145 and Lys290 in the conformation alteration. The findings explain the lack of interfacial activation of CALB and offer new elements to elucidate this mechanism, with the consequent implications for the catalytic properties and classification of lipases.

Ligand-Induced Proton Transfer and Low-Barrier Hydrogen Bond Revealed by X-ray Crystallography

Ligand binding can change the pKa of protein residues and influence enzyme catalysis. Herein, we report three ultrahigh resolution X-ray crystal structures of CTX-M β-lactamase, directly visualizing protonation state changes along the enzymatic pathway: apo protein at 0.79 Å, precovalent complex with nonelectrophilic ligand at 0.89 Å, and acylation transition state (TS) analogue at 0.84 Å. Binding of the noncovalent ligand induces a proton transfer from the catalytic Ser70 to the negatively charged Glu166, and the formation of a low-barrier hydrogen bond (LBHB) between Ser70 and Lys73, with a length of 2.53 Å and the shared hydrogen equidistant from the heteroatoms. QM/MM reaction path calculations determined the proton transfer barrier to be 1.53 kcal/mol. The LBHB is absent in the other two structures although Glu166 remains neutral in the covalent complex. Our data represents the first X-ray crystallographic example of a hydrogen engaged in an enzymatic LBHB, and demonstrates that desolvation of the active site by ligand binding can provide a protein microenvironment conducive to LBHB formation. It also suggests that LBHBs may contribute to stabilization of the TS in general acid/base catalysis together with other preorganized features of enzyme active sites. These structures reconcile previous experimental results suggesting alternatively Glu166 or Lys73 as the general base for acylation, and underline the importance of considering residue protonation state change when modeling protein-ligand interactions. Additionally, the observation of another LBHB (2.47 Å) between two conserved residues, Asp233 and Asp246, suggests that LBHBs may potentially play a special structural role in proteins.

Production and characterization of recombinant perdeuterated cholesterol oxidase

Cholesterol oxidase (CO) is a FAD (flavin adenine dinucleotide) containing enzyme that catalyzes the oxidization and isomerization of cholesterol. Studies directed toward elucidating the catalytic mechanism of CO will provide an important general understanding of Flavin-assisted redox catalysis. Hydrogen atoms play an important role in enzyme catalysis; however, they are not readily visualized in protein X-ray diffraction structures. Neutron crystallography is an ideal method for directly visualizing hydrogen positions at moderate resolutions because hydrogen and deuterium have comparable neutron scattering lengths to other heavy atoms present in proteins. The negative coherent and large incoherent scattering lengths of hydrogen atoms in neutron diffraction experiments can be circumvented by replacing hydrogen atoms with its isotope, deuterium. The perdeuterated form of CO was successfully expressed from minimal medium, purified, and crystallized. X-ray crystallographic structures of the enzyme in the perdeuterated and hydrogenated states confirm that there are no apparent structural differences between the two enzyme forms. Kinetic assays demonstrate that perdeuterated and hydrogenated enzymes are functionally identical. Together, structural and functional studies indicate that the perdeuterated protein is suitable for structural studies by neutron crystallography directed at understanding the role of hydrogen atoms in enzyme catalysis.

Advances in protease engineering for laundry detergents

Proteases are essential ingredients in modern laundry detergents. Over the past 30 years, subtilisin proteases employed in the laundry detergent industry have been engineered by directed evolution and rational design to tailor their properties towards industrial demands. This comprehensive review discusses recent success stories in subtilisin protease engineering. Advances in protease engineering for laundry detergents comprise simultaneous improvement of thermal resistance and activity at low temperatures, a rational strategy to modulate pH profiles, and a general hypothesis for how to increase promiscuous activity towards the production of peroxycarboxylic acids as mild bleaching agents. The three protease engineering campaigns presented provide in-depth analysis of protease properties and have identified principles that can be applied to improve or generate enzyme variants for industrial applications beyond laundry detergents.

Connecting Classical QSAR and LERE Analyses Using Modern Molecular Calculations, LERE-QSAR (VI): Hydrolysis of Substituted Hippuric Acid Phenyl Esters by Trypsin

The reaction mechanism of trypsin was studied by applying DFT and ab initio molecular orbital (MO) calculations to complexes of trypsin with a congeneric series of eight para-substituted hippuric acid phenyl esters, for which a previous quantitative structureactivity relationship (QSAR) study revealed nice linearity of Hammett substitution constant σ(-) with logarithmic values of the MichaelisMenten and catalytic rate constants. Based on the LERE procedure, we performed QSAR analyses on each elementary reaction step during the acylation process. The present calculations showed that the rate-determining step during the acylation process is the transition state (TS) between the enzymesubstrate complex (ES) and tetrahedral intermediate (TET), and that the proton transfer occurs from Ser195 to His57, not between His57 and Asp102. The LERE-QSAR analysis statistically suggested that the variation of overall free-energy changes leading to formation of TS is governed mostly by that of activation energies required to form TS from ES. In spite of a very limited number of congeneric ligands in the current work, it is critically essential to clarify and verify physicochemical meanings of a typical QSAR/Chemoinformatics parameter, Hammett σ(-) based on quantum chemical calculations on the proteinligand kinetics; how Hammett σ(-) behaves in terms of proteinligand interaction energies.

Marvels of enzyme catalysis at true atomic resolution: distortions, bond elongations, hidden flips, protonation states and atom identities

Although general principles of enzyme catalysis are fairly well understood nowadays, many important details of how exactly the substrate is bound and processed in an enzyme remain often invisible and as such elusive. In fortunate cases, structural analysis of enzymes can be accomplished at true atomic resolution thus making possible to shed light on otherwise concealed fine-structural traits of bound substrates, intermediates, cofactors and protein groups. We highlight recent structural studies of enzymes using ultrahigh-resolution X-ray protein crystallography showcasing its enormous potential as a tool in the elucidation of enzymatic mechanisms and in unveiling fundamental principles of enzyme catalysis. We discuss the observation of seemingly hyper-reactive, physically distorted cofactors and intermediates with elongated scissile substrate bonds, the detection of 'hidden' conformational and chemical equilibria and the analysis of protonation states with surprising findings. In delicate cases, atomic resolution is required to unambiguously disclose the identity of atoms as demonstrated for the metal cluster in nitrogenase. In addition to the pivotal structural findings and the implications for our understanding of enzyme catalysis, we further provide a practical framework for resolution enhancement through optimized data acquisition and processing.

Atomic resolution structure of a lysine-specific endoproteinase fromLysobacter enzymogenessuggests a hydroxyl group bound to the oxyanion hole

Lysobacter enzymogeneslysyl endoproteinase (LysC) is a trypsin-type serine protease with a high pH optimum that hydrolyses all Lys-Xaa peptide bonds. The high specificity of LysC renders it useful for biotechnological purposes. The K30R variant of a related lysyl endoproteinase fromAchromobacter lyticushas favourable enzymatic properties that might be transferrable to LysC. To visualize structural differences in the substrate-binding sites, the crystal structures of wild-type and the K30R variant of LysC were determined. The mutation is located at a distance of 12 Å from the catalytic triad and subtly changes the surface properties of the substrate-binding site. The high pH optimum of LysC can be attributed to electrostatic effects of an aromatic Tyr/His stack on the catalytic aspartate and is a general feature of this enzyme subfamily. LysC crystals in complex with the covalent inhibitorNα-p-tosyl-lysyl chloromethylketone yielded data to 1.1 and 0.9 Å resolution, resulting in unprecedented precision of the active and substrate-binding sites for this enzyme subfamily. Error estimates on bond lengths and difference electron density indicate that instead of the expected oxyanion a hydroxyl group binds to the partially solvent-exposed oxyanion hole. Protonation of the alkoxide catalytic intermediate might be a recurring feature during serine protease catalysis.

Distant homology modeling of LCAT and its validation through in silico targeting and in vitro and in vivo assays

LCAT (lecithin:cholesterol acyltransferase) catalyzes the transacylation of a fatty acid of lecithin to cholesterol, generating a cholesteryl ester and lysolecithin. The knowledge of LCAT atomic structure and the identification of the amino acids relevant in controlling its structure and function are expected to be very helpful to understand the enzyme catalytic mechanism, as involved in HDL cholesterol metabolism. However - after an early report in the late ‘90 s - no recent advance has been made about LCAT three-dimensional structure. In this paper, we propose an LCAT atomistic model, built following the most up-to-date molecular modeling approaches, and exploiting newly solved crystallographic structures. LCAT shows the typical folding of the α/β hydrolase superfamily, and its topology is characterized by a combination of α-helices covering a central 7-strand β-sheet. LCAT presents a Ser/Asp/His catalytic triad with a peculiar geometry, which is shared with such other enzyme classes as lipases, proteases and esterases. Our proposed model was validated through different approaches. We evaluated the impact on LCAT structure of some point mutations close to the enzyme active site (Lys218Asn, Thr274Ala, Thr274Ile) and explained, at a molecular level, their phenotypic effects. Furthermore, we devised some LCAT modulators either designed through a de novo strategy or identified through a virtual high-throughput screening pipeline. The tested compounds were proven to be potent inhibitors of the enzyme activity.

Structure and protective efficacy of the Staphylococcus aureus autocleaving protease EpiP

Despite the global medical needs associated with Staphylococcus aureus infections, no licensed vaccines are currently available. We identified and characterized a protein annotated as an epidermin leader peptide processing serine protease (EpiP), as a novel S. aureus vaccine candidate. In addition, we determined the structure of the recombinant protein (rEpiP) by X-ray crystallography. The crystal structure revealed that rEpiP was cleaved somewhere between residues 95 and 100, and we found that the cleavage occurs through an autocatalytic intramolecular mechanism. The protein expressed by S. aureus cells also appeared to undergo a similar processing event. To determine whether the protein acts as a serine protease, we mutated the hypothesized catalytic serine 393 residue to alanine, generating rEpiP-S393A. The crystal structure of this mutant protein showed that the polypeptide chain was not cleaved and was not interacting stably with the active site. Indeed, rEpiP-S393A was shown to be impaired in its protease activity. Mice vaccinated with rEpiP were protected from S. aureus infection (34% survival, P=0.0054). Moreover, the protective efficacy generated by rEpiP and rEpiP-S393A was comparable, implying that the noncleaving mutant could be used for vaccination purposes.

How round is a protein? Exploring protein structures for globularity using conformal mapping

We present a new algorithm that automatically computes a measure of the geometric difference between the surface of a protein and a round sphere. The algorithm takes as input two triangulated genus zero surfaces representing the protein and the round sphere, respectively, and constructs a discrete conformal map f between these surfaces. The conformal map is chosen to minimize a symmetric elastic energy E_S(f) that measures the distance of f from an isometry. We illustrate our approach on a set of basic sample problems and then on a dataset of diverse protein structures. We show first that E_S(f) is able to quantify the roundness of the Platonic solids and that for these surfaces it replicates well traditional measures of roundness such as the sphericity. We then demonstrate that the symmetric elastic energy E_S(f) captures both global and local differences between two surfaces, showing that our method identifies the presence of protruding regions in protein structures and quantifies how these regions make the shape of a protein deviate from globularity. Based on these results, we show that E_S(f) serves as a probe of the limits of the application of conformal mapping to parametrize protein shapes. We identify limitations of the method and discuss its extension to achieving automatic registration of protein structures based on their surface geometry.

A novel approach to the site-selective dual labelling of a protein via chemoselective cysteine modification

Local protein microenvironment is used to control the outcome of reaction between cysteine residues and 2,5-dibromohexanediamide. The differential reactivity is exploited to introduce two orthogonal reactive handles onto the surface of a double cysteine mutant of superfolder green fluorescent protein in a regioselective manner. Subsequent elaboration with commonly used thiol and alkyne containing reagents affects site-selective protein dual labelling.

In vitro Ca(2+)-dependent maturation of milk-clotting recombinant Epr: minor extracellular protease: from Bacillus licheniformis

The minor extracellular protease (Epr) is secreted into the culture medium during Bacillus licheniformis, strain USC13, stationary phase of growth. Whereas, B. subtilis Epr has been reported to be involved in swarming; the B. licheniformis protease is also involved in milk-clotting as shown by the curd forming ability of culture broths expressing this protein. The objectives of this study are the characterization of recombinant B. licheniformis Epr (minor extracellular protease) and the determination of its calcium-dependent activation process. In this work, we have cloned and expressed B. licheniformis Epr in Escherichia coli. We were also able to construct a tridimensional model for Epr based on its homology to Thermococcus kodakarensis pro-tk-subtilisin 2e1p, fervidolysin from Fervidobacterium pennivorans 1rv6, and B. lentus 1GCI subtilisin. Recombinant Epr was accumulated into inclusion bodies; after protein renaturation, Epr undergoes an in vitro calcium-dependent activation, similar to that described for tk protease. The recombinant Epr is capable of producing milk curds with the same clotting activity previously described for the native B. licheniformis Epr enzyme although further rheological and industrial studies should be carried out to confirm its real applicability. This work represents for the first time that Epr may be successfully expressed in a non-bacilli microorganism.

Enumerating pathways of proton abstraction based on a spatial and electrostatic analysis of residues in the catalytic site

The pathways of proton abstraction (PA), a key aspect of most catalytic reactions, is often controversial and highly debated. Ultrahigh-resolution diffraction studies, molecular dynamics, quantum mechanics and molecular mechanic simulations are often adopted to gain insights in the PA mechanisms in enzymes. These methods require expertise and effort to setup and can be computationally intensive. We present a push button methodology--Proton abstraction Simulation (PRISM)--to enumerate the possible pathways of PA in a protein with known 3D structure based on the spatial and electrostatic properties of residues in the proximity of a given nucleophilic residue. Proton movements are evaluated in the vicinity of this nucleophilic residue based on distances, potential differences, spatial channels and characteristics of the individual residues (polarity, acidic, basic, etc). Modulating these parameters eliminates their empirical nature and also might reveal pathways that originate from conformational changes. We have validated our method using serine proteases and concurred with the dichotomy in PA in Class A β-lactamases, both of which are hydrolases. The PA mechanism in a transferase has also been corroborated. The source code is made available at www.sanchak.com/prism.

Protonation-state determination in proteins using high-resolution X-ray crystallography: effects of resolution and completeness

A bond-distance analysis has been undertaken to determine the protonation states of ionizable amino acids in trypsin, subtilisin and lysozyme. The diffraction resolutions were 1.2 Å for trypsin (97% complete, 12% H-atom visibility at 2.5σ), 1.26 Å for subtilisin (100% complete, 11% H-atom visibility at 2.5σ) and 0.65 Å for lysozyme (PDB entry 2vb1; 98% complete, 30% H-atom visibility at 3σ). These studies provide a wide diffraction resolution range for assessment. The bond-length e.s.d.s obtained are as small as 0.008 Å and thus provide an exceptional opportunity for bond-length analyses. The results indicate that useful information can be obtained from diffraction data at around 1.2-1.3 Å resolution and that minor increases in resolution can have significant effects on reducing the associated bond-length standard deviations. The protonation states in histidine residues were also considered; however, owing to the smaller differences between the protonated and deprotonated forms it is much more difficult to infer the protonation states of these residues. Not even the 0.65 Å resolution lysozyme structure provided the necessary accuracy to determine the protonation states of histidine.

Gly or Ala substitutions for Pro(210)Thr(211)Asn(212) at the β8-β9 turn of subtilisin Carlsberg increase the catalytic rate and decrease thermostability

A comparison of the primary structures among psychrophilic, mesophilic, and thermophilic subtilases revealed that the turn between the β8 and β9 strands (β8-β9 turn, BPN' numbering) of psychrophilic subtilases are more flexible than those of their mesophilic and thermophilic counterparts. To investigate the relationship between structure of this turn and enzyme activity as well as thermostability of mesophilic subtilisin Carlsberg (sC), we analyzed 6 mutants of sC with a single, double, or triple Gly or Ala substitutions for Pro(210)Thr(211)Asn(212) at the β8-β9 turn. Among the single Gly substitutions, the P210G substitution most significantly (1.5-fold) increased the specific activity on N-succinyl-Ala-Ala-Pro-Phe-p-nitroanilide (AAPF) substrate and 12-fold decreased the thermostability. All mutants tested showed the increased k(cat) for the AAPF substrate and reduced thermostability compared with the wild-type sC. The k(cat) values of the P210G, P210G/T211G, and P210G/T211G/N212G mutants were 1.5-, 1.7-, and 1.8-fold higher than that of the wild-type sC. There were significant positive correlations between k(cat) and thermal inactivation rates as well as k(cat) and K(m) of the wild-type and mutants. These results demonstrate that the structure of β8-β9 turn, despite its distance from the active site, has significant effects on the catalytic rate and thermostability of sC through a global network of intramolecular interactions and suggest that the lack of flexibility of this turn stabilizes the wild-type sC against thermal inactivation in compensation for some loss of catalytic activity.

Structural insight into serine protease Rv3671c that Protects M. tuberculosis from oxidative and acidic stress

Rv3671c, a putative serine protease, is crucial for persistence of Mycobacterium tuberculosis in the hostile environment of the phagosome. We show that Rv3671c is required for M. tuberculosis resistance to oxidative stress in addition to its role in protection from acidification. Structural and biochemical analyses demonstrate that the periplasmic domain of Rv3671c is a functional serine protease of the chymotrypsin family and, remarkably, that its activity increases on oxidation. High-resolution crystal structures of this protease in an active strained state and in an inactive relaxed state reveal that a solvent-exposed disulfide bond controls the protease activity by constraining two distant regions of Rv3671c and stabilizing it in the catalytically active conformation. In vitro biochemical studies confirm that activation of the protease in an oxidative environment is dependent on this reversible disulfide bond. These results suggest that the disulfide bond modulates activity of Rv3671c depending on the oxidative environment in vivo.

The catalytic aspartate is protonated in the Michaelis complex formed between trypsin and an in vitro evolved substrate-like inhibitor: a refined mechanism of serine protease action

The mechanism of serine proteases prominently illustrates how charged amino acid residues and proton transfer events facilitate enzyme catalysis. Here we present an ultrahigh resolution (0.93 Å) x-ray structure of a complex formed between trypsin and a canonical inhibitor acting through a substrate-like mechanism. The electron density indicates the protonation state of all catalytic residues where the catalytic histidine is, as expected, in its neutral state prior to the acylation step by the catalytic serine. The carboxyl group of the catalytic aspartate displays an asymmetric electron density so that the O_δ2–C_γ bond appears to be a double bond, with O_δ2 involved in a hydrogen bond to His-57 and Ser-214. Only when Asp-102 is protonated on O_δ1 atom could a density functional theory simulation reproduce the observed electron density. The presence of a putative hydrogen atom is also confirmed by a residual mF_obs − DF_calc density above 2.5 σ next to O_δ1. As a possible functional role for the neutral aspartate in the active site, we propose that in the substrate-bound form, the neutral aspartate residue helps to keep the pK_a of the histidine sufficiently low, in the active neutral form. When the histidine receives a proton during the catalytic cycle, the aspartate becomes simultaneously negatively charged, providing additional stabilization for the protonated histidine and indirectly to the tetrahedral intermediate. This novel proposal unifies the seemingly conflicting experimental observations, which were previously seen as either supporting the charge relay mechanism or the neutral pK_a histidine theory.

Atomic resolution crystallography of a complex of triosephosphate isomerase with a reaction-intermediate analog: new insight in the proton transfer reaction mechanism

Enzymes achieve their catalytic proficiency by precisely positioning the substrate and catalytic residues with respect to each other. Atomic resolution crystallography is an excellent tool to study the important details of these geometric active-site features. Here, we have investigated the reaction mechanism of triosephosphate isomerase (TIM) using atomic resolution crystallographic studies at 0.82-A resolution of leishmanial TIM complexed with the well-studied reaction-intermediate analog phosphoglycolohydroxamate (PGH). Remaining unresolved aspects of the reaction mechanism of TIM such as the protonation state of the first reaction intermediate and the properties of the hydrogen-bonding interactions in the active site are being addressed. The hydroxamate moiety of PGH interacts via unusually short hydrogen bonds of its N1-O1 moiety with the carboxylate group of the catalytic glutamate (Glu167), for example, the distance of N1(PGH)-OE2(Glu167) is 2.69 +/- 0.01 A and the distance of O1(PGH)-OE1(Glu167) is 2.60 +/- 0.01 A. Structural comparisons show that the side chain of the catalytic base (Glu167) can move during the reaction cycle in a small cavity, located above the hydroxamate plane. The structure analysis suggests that the hydroxamate moiety of PGH is negatively charged. Therefore, the bound PGH mimics the negatively charged enediolate intermediate, which is formed immediately after the initial proton abstraction from DHAP by the catalytic glutamate. The new findings are discussed in the context of the current knowledge of the TIM reaction mechanism.

Crystal structure of glyceraldehyde-3-phosphate dehydrogenase 1 from methicillin-resistant Staphylococcus aureus MRSA252 provides novel insights into substrate binding and catalytic mechanism

The dreaded pathogen Staphylococcus aureus is one of the causes of morbidity and mortality worldwide. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), one of the key glycolytic enzymes, is irreversibly oxidized under oxidative stress and is responsible for sustenance of the pathogen inside the host. With an aim to elucidate the catalytic mechanism and identification of intermediates involved, we describe in this study different crystal structures of GAPDH1 from methicillin-resistant S. aureus MRSA252 (SaGAPDH1) in apo and holo forms of wild type, thioacyl intermediate, and ternary complexes of active-site mutants with physiological substrate d-glyceraldehyde-3-phosphate (G3P) and coenzyme NAD(+). A new phosphate recognition site, "new P(i)" site, similar to that observed in GAPDH from Thermotoga maritima, is reported here, which is 3.40 A away from the "classical P(i)" site. Ternary complexes discussed are representatives of noncovalent Michaelis complexes in the ground state. d-G3P is bound to all the four subunits of C151S.NAD and C151G.NAD in more reactive hydrate (gem-di-ol) form. However, in C151S+H178N.NAD, the substrate is bound to two chains in aldehyde form and in gem-di-ol form to the other two. This work reports binding of d-G3P to the C151G mutant in an inverted manner for the very first time. The structure of the thiaocyl complex presented here is formed after the hydride transfer. The C3 phosphate of d-G3P is positioned at the "P(s)" site in the ternary complexes but at the "new P(i)" site in the thioacyl complex and C1-O1 bond points opposite to His178 disrupting the alignment between itself and NE2 of His178. A new conformation (Conformation I) of the 209-215 loop has also been identified, where the interaction between phosphate ion at the "new P(i)" site and conserved Gly212 is lost. Altogether, inferences drawn from the kinetic analyses and crystal structures suggest the "flip-flop" model proposed for the enzyme mechanism.

Ultrahigh (0.93A) resolution structure of manganese peroxidase from Phanerochaete chrysosporium: implications for the catalytic mechanism

Manganese peroxidase (MnP) is an extracellular heme enzyme produced by the lignin-degrading white-rot fungus Phanerochaete chrysosporium. MnP catalyzes the peroxide-dependent oxidation of Mn(II) to Mn(III). The Mn(III) is released from the enzyme in complex with oxalate, enabling the oxalate-Mn(III) complex to serve as a diffusible redox mediator capable of oxidizing lignin, especially under the mediation of unsaturated fatty acids. One heme propionate and the side chains of Glu35, Glu39 and Asp179 have been identified as Mn(II) ligands in our previous crystal structures of native MnP. In our current work, new 0.93A and 1.05A crystal structures of MnP with and without bound Mn(II), respectively, have been solved. This represents only the sixth structure of a protein of this size at 0.93A resolution. In addition, this is the first structure of a heme peroxidase from a eukaryotic organism at sub-Angstrom resolution. These new structures reveal an ordering/disordering of the C-terminal loop, which is likely required for Mn binding and release. In addition, the catalytic Arg42 residue at the active site, normally thought to function only in the peroxide activation process, also undergoes ordering/disordering that is coupled to a transient H-bond with the Mn ligand, Glu39. Finally, these high-resolution structures also reveal the exact H atoms in several parts of the structure that are relevant to the catalytic mechanism.

N-terminal deletion of Tk1689, a subtilisin-like serine protease from Thermococcus kodakaraensis, copes with its cytotoxicity in Escherichia coli

Tk1689, a subtilisin-like serine protease from Thermococcus kodakaraensis, was found to be toxic to the host cells when produced in the pro-protein form (Pro-Tk1689) in Eschericia coli. Cytotoxic effect of Pro-Tk1689 was reduced when signal and pro-peptide both were removed and the protein was produced in the mature form (MP-Tk1689). The mature protein was produced in E. coli in the soluble form. Recombinant MP-Tk1689 was catalytically active and exhibited optimum activity at 55°C and pH 7. Specific activity of the enzyme was 700 U/mg. The enzyme displayed a half life of 80 min at 60°C.

Challenging a paradigm: theoretical calculations of the protonation state of the Cys25-His159 catalytic diad in free papain

A central mechanistic paradigm of cysteine proteases is that the His-Cys catalytic diad forms an ion-pair NH(+)/S(-) already in the catalytically active free enzyme. Most molecular modeling studies of cysteine proteases refer to this paradigm as their starting point. Nevertheless, several recent kinetics and X-ray crystallography studies of viral and bacterial cysteine proteases depart from the ion-pair mechanism, suggesting general base catalysis. We challenge the postulate of the ion-pair formation in free papain. Applying our QM/SCRF(VS) molecular modeling approach, we analyzed all protonation states of the catalytic diad in free papain and its SMe derivative, comparing the predicted and experimental pK(a) data. We conclude that the His-Cys catalytic diad in free papain is fully protonated, NH(+)/SH. The experimental pK(a) = 8.62 of His159 imidazole in free papain, obtained by NMR-controlled titration and originally interpreted as the NH(+)/S(-) <==> N/S(-) NH(+)/S(-) <==> N/S(-) equilibrium, is now assigned to the NH(+)/SH <==> N/SH NH(+)/SH <==> N/SH equilibrium.

Combined high-resolution neutron and X-ray analysis of inhibited elastase confirms the active-site oxyanion hole but rules against a low-barrier hydrogen bond

To help resolve long-standing questions regarding the catalytic activity of the serine proteases, the structure of porcine pancreatic elastase has been analyzed by high-resolution neutron and X-ray crystallography. To mimic the tetrahedral transition intermediate, a peptidic inhibitor was used. A single large crystal was used to collect room-temperature neutron data to 1.65 A resolution and X-ray data to 1.20 A resolution. Another crystal provided a low-temperature X-ray data set to 0.94 A resolution. The neutron data are to higher resolution than previously reported for a serine protease and the X-ray data are comparable with other studies. The neutron and X-ray data show that the hydrogen bond between His57 and Asp102 (chymotrypsin numbering) is 2.60 A in length and that the hydrogen-bonding hydrogen is 0.80-0.96 A from the histidine nitrogen. This is not consistent with a low-barrier hydrogen which is predicted to have the hydrogen midway between the donor and acceptor atom. The observed interaction between His57 and Asp102 is essentially a short but conventional hydrogen bond, sometimes described as a short ionic hydrogen bond. The neutron analysis also shows that the oxygen of the oxopropyl group of the inhibitor is present as an oxygen anion rather than a hydroxyl group, supporting the role of the "oxyanion hole" in stabilizing the tetrahedral intermediate in catalysis.

Polarizable atomic multipole X-ray refinement: application to peptide crystals

Recent advances in computational chemistry have produced force fields based on a polarizable atomic multipole description of biomolecular electrostatics. In this work, the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field is applied to restrained refinement of molecular models against X-ray diffraction data from peptide crystals. A new formalism is also developed to compute anisotropic and aspherical structure factors using fast Fourier transformation (FFT) of Cartesian Gaussian multipoles. Relative to direct summation, the FFT approach can give a speedup of more than an order of magnitude for aspherical refinement of ultrahigh-resolution data sets. Use of a sublattice formalism makes the method highly parallelizable. Application of the Cartesian Gaussian multipole scattering model to a series of four peptide crystals using multipole coefficients from the AMOEBA force field demonstrates that AMOEBA systematically underestimates electron density at bond centers. For the trigonal and tetrahedral bonding geometries common in organic chemistry, an atomic multipole expansion through hexadecapole order is required to explain bond electron density. Alternatively, the addition of interatomic scattering (IAS) sites to the AMOEBA-based density captured bonding effects with fewer parameters. For a series of four peptide crystals, the AMOEBA-IAS model lowered R(free) by 20-40% relative to the original spherically symmetric scattering model.

Protonate3D: assignment of ionization states and hydrogen coordinates to macromolecular structures

A new method, called Protonate3D, is presented for the automated prediction of hydrogen coordinates given the 3D coordinates of the heavy atoms of a macromolecular structure. Protonate3D considers side-chain "flip," rotamer, tautomer, and ionization states of all chemical groups, ligands, and solvent, provided suitable templates are available in a parameter file. The energy model includes van der Waals, Coulomb, solvation, rotamer, tautomer, and titration effects. The results of computational validation experiments suggest that Protonate3D can accurately predict the location of hydrogen atoms in macromolecular structures.

The crystal structures of the psychrophilic subtilisin S41 and the mesophilic subtilisin Sph reveal the same calcium-loaded state

We determine and compare the crystal structure of two proteases belonging to the subtilisin superfamily: S41, a cold-adapted serine protease produced by Antarctic bacilli, at 1.4 A resolution and Sph, a mesophilic serine protease produced by Bacillus sphaericus, at 0.8 A resolution. The purpose of this comparison was to find out whether multiple calcium ion binding is a molecular factor responsible for the adaptation of S41 to extreme low temperatures. We find that these two subtilisins have the same subtilisin fold with a root mean square between the two structures of 0.54 A. The final models for S41 and Sph include a calcium-loaded state of five ions bound to each of these two subtilisin molecules. None of these calcium-binding sites correlate with the high affinity known binding site (site A) found for other subtilisins. Structural analysis of the five calcium-binding sites found in these two crystal structures indicate that three of the binding sites have two side chains of an acidic residue coordinating the calcium ion, whereas the other two binding sites have either a main-chain carbonyl, or only one acidic residue side chain coordinating the calcium ion. Thus, we conclude that three of the sites are of high affinity toward calcium ions, whereas the other two are of low affinity. Because Sph is a mesophilic subtilisin and S41 is a psychrophilic subtilisin, but both crystal structures were found to bind five calcium ions, we suggest that multiple calcium ion binding is not responsible for the adaptation of S41 to low temperatures.

Enzymatic proteolysis of a surface-bound alpha-helical polypeptide

In this work, we studied the interactions of enzymes with model substrate surfaces using label-free techniques. Our model system was based on serine proteases (a class of enzymes that digests proteins) and surface-bound polypeptide substrates. While previous studies have focused on bulk media factors such as pH, ionic strength, and surfactants, this study focuses on the role of the surface-bound substrate itself. In particular, we assess how the substrate density of a polypeptide with an alpha-helical secondary structure influences surface reactivity. An alpha-helical secondary structure was chosen based on literature indicating that stable alpha-helices can resist enzymatic digestion. To investigate the protease resistance of a surface-bound a-helix, we designed an a-helical polypeptide (SS-polypeptide, where SS = disulfide), used it to form films of varying surface coverage and then measured responses of the films to enzymatic exposure. Using quartz-crystal microbalance with dissipation (QCM-D), angle-resolved X-ray photoelectron spectroscopy (AR-XPS), grazing-angle infrared spectroscopy (GAIRS), and other techniques, we characterized the degradation of films to determine how the lateral packing density of the surface-bound SS-polypeptide substrate affected surface proteolysis. Characterization of pure SS-polypeptide films indicated dense packing of helices that maintained their helical structure and were generally oriented normal to the surface. We found that films of pure SS-polypeptide significantly resisted enzymatic digestion, while incorporation of very minor amounts of a diluent in such films resulted in rapid digestion. In part, this may be due to the need for the enzyme to bind several peptides along the peptide substrate within the cleft for digestion to occur. Only SS-polypeptide films that were densely packed and did not permit catalytic access to multiple peptides (e.g., terminal peptides only) were resistant to enzymatic proteolysis.