Structures of product and inhibitor complexes of Streptomyces griseus protease A at 1.8 A resolution. A model for serine protease catalysis

Docking, thermodynamics and molecular dynamics (MD) studies of a non-canonical protease inhibitor, MP-4, from Mucuna pruriens

Sequence and structural homology suggests that MP-4 protein from Mucuna pruriens belongs to Kunitz-type protease inhibitor family. However, biochemical assays showed that this protein is a poor inhibitor of trypsin. To understand the basis of observed poor inhibition, thermodynamics and molecular dynamics (MD) simulation studies on binding of MP-4 to trypsin were carried out. Molecular dynamics simulations revealed that temperature influences the spectrum of conformations adopted by the loop regions in the MP-4 structure. At an optimal temperature, MP-4 achieves maximal binding while above and below the optimum temperature, its functional activity is hampered due to unfavourable flexibility and relative rigidity, respectively. The low activity at normal temperature is due to the widening of the conformational spectrum of the Reactive Site Loop (RSL) that reduces the probability of formation of stabilizing contacts with trypsin. The unique sequence of the RSL enhances flexibility at ambient temperature and thus reduces its ability to inhibit trypsin. This study shows that temperature influences the function of a protein through modulation in the structure of functional domain of the protein. Modulation of function through appearance of new sequences that are more sensitive to temperature may be a general strategy for evolution of new proteins.

Biochemical and molecular characterization of new keratinoytic protease from Actinomadura viridilutea DZ50

A new extracellular thermostable keratinolytic protease, designated KERDZ, was purified and characterized from a thermophilic actinomycetes Actinomadura viridilutea DZ50 isolated from Algerian fishing port. The isolate exhibited high keratinase production when grown in chicken-feather meal media (18,000U/ml) after 96-h of incubation at 45°C. The enzyme was purified by ammonium sulfate precipitation (35-55%)-dialysis and heat treatment (30min at 75°C) followed by UNO S-1 FPLC cation exchange chromatography and size exclusion HPLC column. The biochemical characterizations carried on include physico-chemical determination and spectroscopic analysis. The MALDI-TOF/MS analysis revealed that the purified enzyme was a monomer with a molecular mass of 19536.10-Da. The sequence of the 25 N-terminal residues of KERDZ showed high homology with those of actinomycetes keratinases. Optimal activity was achieved at pH 11 and 80°C. KERDZ was completely inhibited by PMSF and DFP suggested its belonging to the serine keratinase family. KERDZ displayed higher levels of hydrolysis and catalytic efficiency than bacterial keratinases (KERAK-29, Actinase E, and KERAB) and subtilisins (subtilisin Carlsberg and subtilisin Novo). The kerDZ gene encoding KERDZ was isolated and its DNA sequence was determined. These properties make KERDZ a potential, promising and eco-friendly alternative to the conventional chemicals used for industrial applications.

Quantifying tetrahedral adduct formation and stabilization in the cysteine and the serine proteases

Two new papain inhibitors have been synthesized where the terminal α-carboxyl groups of Z-Phe-Ala-COOH and Ac-Phe-Gly-COOH have been replaced by a proton to give Z-Phe-Ala-H and Ac-Phe-Gly-H. We show that for papain, replacing the terminal carboxylate group of a peptide inhibitor with a hydrogen atom decreases binding 3–4 fold while replacing an aldehyde or glyoxal group with a hydrogen atom decreases binding by 300,000–1,000,000 fold. Thiohemiacetal formation by papain with aldehyde or glyoxal inhibitors is shown to be ~ 10,000 times more effective than hemiacetal or hemiketal formation with chymotrypsin. It is shown using effective molarities, that for papain, thiohemiacetal stabilization is more effective with aldehyde inhibitors than with glyoxal inhibitors. The effective molarity obtained when papain is inhibited by an aldehyde inhibitor is similar to the effective molarity obtained when chymotrypsin is inhibited by glyoxal inhibitors showing that both enzymes can stabilize tetrahedral adducts by similar amounts. Therefore the greater potency of aldehyde and glyoxal inhibitors with papain is not due to greater thiohemiacetal stabilization by papain compared to the hemiketal and hemiacetal stabilization by chymotrypsin, instead it reflects the greater intrinsic reactivity of the catalytic thiol group of papain compared to the catalytic hydroxyl group of chymotrypsin. It is argued that while the hemiacetals and thiohemiacetals formed with the serine and cysteine proteases respectively can mimic the catalytic tetrahedral intermediate they are also analogues of the productive and non-productive acyl intermediates that can be formed with the cysteine and serine proteases.

•

We compare thiohemiacetal and hemiacetal stabilization by papain and chymotrypsin.

•

An aldehyde or glyoxal group increases binding by 300,000–1,000,000 fold.

•

Thiohemiacetal formation is ~ 10,000 fold greater than hemiacetal formation.

•

Thiohemiacetal formation is more effective with aldehyde than glyoxal inhibitors.

•

Both papain and chymotrypsin stabilize tetrahedral adducts by similar amounts.

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.

S46 peptidases are the first exopeptidases to be members of clan PA

The dipeptidyl aminopeptidase BII (DAP BII) belongs to a serine peptidase family, S46. The amino acid sequence of the catalytic unit of DAP BII exhibits significant similarity to those of clan PA endopeptidases, such as chymotrypsin. However, the molecular mechanism of the exopeptidase activity of family S46 peptidase is unknown. Here, we report crystal structures of DAP BII. DAP BII contains a peptidase domain including a typical double β-barrel fold and previously unreported α-helical domain. The structures of peptide complexes revealed that the α-helical domain covers the active-site cleft and the side chain of Asn330 in the domain forms hydrogen bonds with the N-terminus of the bound peptide. These observations indicate that the α-helical domain regulates the exopeptidase activity of DAP BII. Because S46 peptidases are not found in mammals, we expect that our study will be useful for the design of specific inhibitors of S46 peptidases from pathogens.

Poly(2-hydroxyethyl methacrylate) brush surface for specific and oriented adsorption of glycosidases

We present a detailed picture to screen general ligands from simple chemicals for fabricating affinity surface to glycosidase enzymes. The surface was constructed by grafting poly(2-hydroxyethyl methacrylate) (PHEMA) brush on SPR gold chip via surface-initiated atom-transfer radical polymerization, after which poly(methoxyethyl methacrylate) (PMEMA) and poly(oligo(ethylene glycol) methacrylate) (POEGMA) brushes were also prepared for comparison. SPR measurements were adopted to monitor the early-stage adsorption of two glycosidases and three other typical proteins. PHEMA resists the adsorption of lysozyme, bovine serum albumin, and fibrinogen, while it is capable of specifically adsorbing β-glucosidase (GLU) and β-galactosidase (GAL). These are quite different from the nonspecific adsorption of PMEMA and the anti-nonspecific adsorption of POEGMA to the studied proteins, because PHEMA is the acceptor substrate of the glycosidases. About 69.6 and 93.7 ng/cm(2) of GAL and GLU are adsorbed on the PHEMA brush surface, of which more than 49.6 ng/cm(2) is remained after washing with PBS. The specific adsorption process is appropriately described by Freundlich isothermal model rather than Langmuir one, and is also indicated to be spontaneous, endothermic, and entropy driven through thermodynamic studies. Taking into account all stated results above, we propose that molecular recognition takes place between the hydroxyl groups of PHEMA and the active sites of glycosidases, which subsequently enables the oriented adsorption of glycosidases on the brush surface. The adsorbed enzyme can be effectively eluted with 1.0 M aqueous solution of ethanol. Our findings open the door to the further development in the design of novel acceptor substrate-ligand affinity chromatography for enzyme purification.

Crystal structure of human enterovirus 71 3C protease

Human enterovirus 71 (EV71) is the major pathogen that causes hand, foot and mouth disease that particularly affects young children. Growing hand, foot and mouth disease outbreaks were observed worldwide in recent years and caused devastating losses both economically and politically. However, vaccines or effective drugs are unavailable to date. The genome of EV71 consists of a positive sense, single-stranded RNA of ∼7400 bp, encoding a large precursor polyprotein that requires proteolytic processing to generate mature viral proteins. The proteolytic processing mainly depends on EV71 3C protease (3C(pro)) that possesses both proteolysis and RNA binding activities, which enable the protease to perform multiple tasks in viral replication and pathogen-host interactions. The central roles played by EV71 3C(pro) make it an appealing target for antiviral drug development. We determined the first crystal structure of EV71 3C(pro) and analyzed its enzymatic activity. The crystal structure shows that EV71 3C(pro) has a typical chymotrypsin-like fold that is common in picornaviral 3C(pro). Strikingly, we found an important surface loop, also denoted as β-ribbon, which adopts a novel open conformation in EV71 3C(pro). We identified two important residues located at the base of the β-ribbon, Gly123 and His133, which form hinges that govern the intrinsic flexibility of the ribbon. Structure-guided mutagenesis studies revealed that the hinge residues are important to EV71 3C(pro) proteolytic activities. In summary, our work provides the first structural insight into EV71 3C(pro), including a mobile β-ribbon, which is relevant to the proteolytic mechanism. Our data also provides a framework for structure-guided inhibitor design against EV71 3C(pro).

Structure and catalysis of acylaminoacyl peptidase: closed and open subunits of a dimer oligopeptidase

Acylaminoacyl peptidase from Aeropyrum pernix is a homodimer that belongs to the prolyl oligopeptidase family. The monomer subunit is composed of one hydrolase and one propeller domain. Previous crystal structure determinations revealed that the propeller domain obstructed the access of substrate to the active site of both subunits. Here we investigated the structure and the kinetics of two mutant enzymes in which the aspartic acid of the catalytic triad was changed to alanine or asparagine. Using different substrates, we have determined the pH dependence of specificity rate constants, the rate-limiting step of catalysis, and the binding of substrates and inhibitors. The catalysis considerably depended both on the kind of mutation and on the nature of the substrate. The results were interpreted in terms of alterations in the position of the catalytic histidine side chain as demonstrated with crystal structure determination of the native and two mutant structures (D524N and D524A). Unexpectedly, in the homodimeric structures, only one subunit displayed the closed form of the enzyme. The other subunit exhibited an open gate to the catalytic site, thus revealing the structural basis that controls the oligopeptidase activity. The open form of the native enzyme displayed the catalytic triad in a distorted, inactive state. The mutations affected the closed, active form of the enzyme, disrupting its catalytic triad. We concluded that the two forms are at equilibrium and the substrates bind by the conformational selection mechanism.

The effect of the P1 side chain on the binding of optimized carboxylate and activated carbonyl inhibitors of the hepatitis C virus NS3 protease

Peptidyl inhibitors of the hepatitis C virus NS3 protease hold much promise as direct-acting antiviral agents against hepatitis C infection. The optimization of N-terminal cleavage products, found to exhibit activity (product inhibition) against the enzyme, has led to potent tripeptide inhibitors that bear free C-terminal carboxylate groups. An analogous activated carbonyl compound (pentafluoroethyl ketone) bearing a P1 norvaline (Nva) was found to possess comparable activity against hepatitis C virus protease. However, an analogue bearing an aminocyclopropylcarboxylic acid (Acca) P1 residue exhibited very poor activity. 19F-NMR studies indicate that the propensity of the Acca-derived activated carbonyl to form hemiketals is only slightly reduced compared with that of a P1 Nva equivalent. These results, as well as molecular modeling studies, argue against steric hindrance of the nucleophilic attack of Ser-139 accounting for the poor mechanism-based inhibition by the former. We hypothesize that the conformational properties of the respective C-termini in the context of an adaptable active site account for the divergent P1 structure–activity relationships.

Insights into the enzyme-substrate interaction in the norovirus 3C-like protease

The Glu54 residue of the norovirus 3C-like protease was implicated in proteolysis as a third-member carboxylate of the catalytic triad. The E54L mutant protease cleaved the sequence (133)LSFE/AP between the 3B and 3C regions of norovirus polyprotein, but did not cleave the sequence (198)ATSE/GK between the 3A and 3B. The 3BC junction mutation (3B-L133A or 3B-F135S) hampered the cleavage by the E54L protease, whereas the 3AB junction mutation (3A-A198L, S200F) allowed the E54L protease to digest. These results indicate that the E54L mutant protease is a substrate-specificity mutant and requires large hydrophobic amino acid residues at both P4 and P2 positions of the substrate. It was notable that the 3A-S200F P2 position mutation caused tight interaction between the wild-type protease and the C-terminus of the 3A protein, hence a decreased release rate of the product from the enzyme. This tight binding was dependent on the hydrophobicity of amino acid residues introduced at position 200 of the 3A region and was affected by the mutation in the bII-cII loop of the protease or the mutation of position 198 of 3A corresponding to the P4 position of the substrate. These results suggest that the protease and the substrate sense each other in the process of the proteolysis, being supported by crystal structures.

StoneHinge: hinge prediction by network analysis of individual protein structures

Hinge motions are important for molecular recognition, and knowledge of their location can guide the sampling of protein conformations for docking. Predicting domains and intervening hinges is also important for identifying structurally self-determinate units and anticipating the influence of mutations on protein flexibility and stability. Here we present StoneHinge, a novel approach for predicting hinges between domains using input from two complementary analyses of noncovalent bond networks: StoneHingeP, which identifies domain-hinge-domain signatures in ProFlex constraint counting results, and StoneHingeD, which does the same for DomDecomp Gaussian network analyses. Predictions for the two methods are compared to hinges defined in the literature and by visual inspection of interpolated motions between conformations in a series of proteins. For StoneHingeP, all the predicted hinges agree with hinge sites reported in the literature or observed visually, although some predictions include extra residues. Furthermore, no hinges are predicted in six hinge-free proteins. On the other hand, StoneHingeD tends to overpredict the number of hinges, while accurately pinpointing hinge locations. By determining the consensus of their results, StoneHinge improves the specificity, predicting 11 of 13 hinges found both visually and in the literature for nine different open protein structures, and making no false-positive predictions. By comparison, a popular hinge detection method that requires knowledge of both the open and closed conformations finds 10 of the 13 known hinges, while predicting four additional, false hinges.

Structure of a serine protease poised to resynthesize a peptide bond

The serine proteases are among the most thoroughly studied enzymes, and numerous crystal structures representing the enzyme-substrate complex and intermediates in the hydrolysis reactions have been reported. Some aspects of the catalytic mechanism remain controversial, however, especially the role of conformational changes in the reaction. We describe here a high-resolution (1.46 A) crystal structure of a complex formed between a cleaved form of bovine pancreatic trypsin inhibitor (BPTI) and a catalytically inactive trypsin variant with the BPTI cleavage site ideally positioned in the active site for resynthesis of the peptide bond. This structure defines the positions of the newly generated amino and carboxyl groups following the 2 steps in the hydrolytic reaction. Comparison of this structure with those representing other intermediates in the reaction demonstrates that the residues of the catalytic triad are positioned to promote each step of both the forward and reverse reaction with remarkably little motion and with conservation of hydrogen-bonding interactions. The results also provide insights into the mechanism by which inhibitors like BPTI normally resist hydrolysis when bound to their target proteases.

Iterative-build OMIT maps: map improvement by iterative model building and refinement without model bias

An OMIT procedure is presented that has the benefits of iterative model building density modification and refinement yet is essentially unbiased by the atomic model that is built.

A procedure for carrying out iterative model building, density modification and refinement is presented in which the density in an OMIT region is essentially unbiased by an atomic model. Density from a set of overlapping OMIT regions can be combined to create a composite ‘iterative-build’ OMIT map that is everywhere unbiased by an atomic model but also everywhere benefiting from the model-based information present elsewhere in the unit cell. The procedure may have applications in the validation of specific features in atomic models as well as in overall model validation. The procedure is demonstrated with a molecular-replacement structure and with an experimentally phased structure and a variation on the method is demonstrated by removing model bias from a structure from the Protein Data Bank.

Structural analysis of foot-and-mouth disease virus 3C protease: a viable target for antiviral drugs?

Foot-and-mouth disease virus causes a major global agricultural problem that is difficult to control with existing vaccines. Structural analyses of the viral 3C protease not only have provided fresh insights into the catalytic mechanism of an unusual class of chymotrypsin-like cysteine proteases, but also are generating valuable information to drive the quest for effective antiviral therapies.

Mesophile versus thermophile: insights into the structural mechanisms of kinetic stability

Obtaining detailed knowledge of folding intermediate and transition state (TS) structures is critical for understanding protein folding mechanisms. Comparisons between proteins adapted to survive extreme temperatures with their mesophilic homologs are likely to provide valuable information on the interactions relevant to the unfolding transition. For kinetically stable proteins such as alpha-lytic protease (alphaLP) and its family members, their large free energy barrier to unfolding is central to their biological function. To gain new insights into the mechanisms that underlie kinetic stability, we have determined the structure and high temperature unfolding kinetics of a thermophilic homolog, Thermobifida fusca protease A (TFPA). These studies led to the identification of a specific structural element bridging the N and C-terminal domains of the protease (the "domain bridge") proposed to be associated with the enhanced high temperature kinetic stability in TFPA. Mutagenesis experiments exchanging the TFPA domain bridge into alphaLP validate this hypothesis and illustrate key structural details that contribute to TFPA's increased kinetic thermostability. These results lead to an updated model for the unfolding transition state structure for this important class of proteases in which domain bridge undocking and unfolding occurs at or before the TS. The domain bridge appears to be a structural element that can modulate the degree of kinetic stability of the different members of this class of proteases.

Structural and mechanistic exploration of acid resistance: kinetic stability facilitates evolution of extremophilic behavior

Kinetically stable proteins are unique in that their stability is determined solely by kinetic barriers rather than by thermodynamic equilibria. To better understand how kinetic stability promotes protein survival under extreme environmental conditions, we analyzed the unfolding behavior and determined the structure of Nocardiopsis alba Protease A (NAPase), an acid-resistant, kinetically stable protease, and compared these results with a neutrophilic homolog, alpha-lytic protease (alphaLP). Although NAPase and alphaLP have the same number of acid-titratable residues, kinetic studies revealed that the height of the unfolding free energy barrier for NAPase is less sensitive to acid than that of alphaLP, thereby accounting for NAPase's improved tolerance of low pH. A comparison of the alphaLP and NAPase structures identified multiple salt-bridges in the domain interface of alphaLP that were relocated to outer regions of NAPase, suggesting a novel mechanism of acid stability in which acid-sensitive electrostatic interactions are rearranged to similarly affect the energetics of both the native state and the unfolding transition state. An acid-stable variant of alphaLP in which a single interdomain salt-bridge is replaced with a corresponding intradomain NAPase salt-bridge shows a dramatic >15-fold increase in acid resistance, providing further evidence for this hypothesis. These observations also led to a general model of the unfolding transition state structure for alphaLP protease family members in which the two domains separate from each other while remaining relatively intact themselves. These results illustrate the remarkable utility of kinetic stability as an evolutionary tool for developing longevity over a broad range of harsh conditions.

The acylaminoacyl peptidase from Aeropyrum pernix K1 thought to be an exopeptidase displays endopeptidase activity

Mammalian acylaminoacyl peptidase, a member of the prolyl oligopeptidase family of serine peptidases, is an exopeptidase, which removes acylated amino acid residues from the N terminus of oligopeptides. We have investigated the kinetics and inhibitor binding of the orthologous acylaminoacyl peptidase from the thermophile Aeropyrum pernix K1 (ApAAP). Complex pH-rate profiles were found with charged substrates, indicating a strong electrostatic effect in the surroundings of the active site. Unexpectedly, we have found that oligopeptides can be hydrolysed beyond the N-terminal peptide bond, demonstrating that ApAAP exhibits endopeptidase activity. It was thought that the enzyme is specific for hydrophobic amino acids, in particular phenylalanine, in accord with the non-polar S1 subsite of ApAAP. However, cleavage after an Ala residue contradicted this notion and demonstrated that P1 residues of different nature may bind to the S1 subsite depending on the remaining peptide residues. The crystal structures of the complexes formed between the enzyme and product-like inhibitors identified the oxyanion-binding site unambiguously and demonstrated that the phenylalanine ring of the P1 peptide residue assumes a position different from that established in a previous study, using 4-nitrophenylphosphate. We have found that the substrate-binding site extends beyond the S2 subsite, being capable of binding peptides with a longer N terminus. The S2 subsite displays a non-polar character, which is unique among the enzymes of this family. The S3 site was identified as a hydrophobic region that does not form hydrogen bonds with the inhibitor P3 residue. The enzyme-inhibitor complexes revealed that, upon ligand-binding, the S1 subsite undergoes significant conformational changes, demonstrating the plasticity of the specificity site.

Structural and mutagenic analysis of foot-and-mouth disease virus 3C protease reveals the role of the beta-ribbon in proteolysis

The 3C protease (3C(pro)) from foot-and-mouth disease virus (FMDV), the causative agent of a widespread and economically devastating disease of domestic livestock, is a potential target for antiviral drug design. We have determined the structure of a new crystal form of FMDV 3C(pro), a chymotrypsin-like cysteine protease, which reveals features that are important for catalytic activity. In particular, we show that a surface loop which was disordered in previous structures adopts a beta-ribbon structure that is conformationally similar to equivalent regions on other picornaviral 3C proteases and some serine proteases. This beta-ribbon folds over the peptide binding cleft and clearly contributes to substrate recognition. Replacement of Cys142 at the tip of the beta-ribbon with different amino acids has a significant impact on enzyme activity and shows that higher activity is obtained with more hydrophobic side chains. Comparison of the structure of FMDV 3C(pro) with homologous enzyme-peptide complexes suggests that this correlation arises because the side chain of Cys142 contacts the hydrophobic portions of the P2 and P4 residues in the peptide substrate. Collectively, these findings provide compelling evidence for the role of the beta-ribbon in catalytic activity and provide valuable insights for the design of FMDV 3C(pro) inhibitors.

Simulation of diffusion time of small molecules in protein crystals

A simple model for evaluation of diffusion times of small molecule into protein crystals has been developed, which takes into account the physical and chemical properties both of protein crystal and the diffusing molecules. The model also includes consideration of binding and the binding affinity of a ligand to the protein. The model has been validated by simulation of experimental set-ups of several examples found in the literature. These experiments cover a wide range of situations: from small to relatively large diffusing molecules, crystals having low, medium, or high protein density, and different size. The reproduced experiments include ligand exchange in protein crystals by soaking techniques. Despite the simplifying assumptions of the model, theoretical and experimental data are in agreement with available data, with experimental diffusion times ranging from a few seconds to several hours. The method has been used successfully for planning intermediate cryotrapping experiments in maltodextrin phosphorylase crystals.

Hierarchical PLS Modeling for Predicting the Binding of a Comprehensive Set of Structurally Diverse Protein−Ligand Complexes

A new approach is presented for predicting ligand binding to proteins using hierarchical partial-least-squares regression to latent structures (Hi-PLS). Models were based on information from the 2002 release of the PDBbind database containing (after in-house refinement) high-resolution X-ray crystallography and binding affinity (Kd or Ki) data for 612 protein-ligand complexes. The complexes were characterized by four different descriptor blocks: three-dimensional (3D) structural descriptors of the proteins, protein-ligand interactions according to the Validate scoring function, binding site surface areas, and ligand 2D and 3D descriptors. These descriptor blocks were used in Hi-PLS models, generated using both linear and nonlinear terms, to relate the characterizations to pKd/i. The results show that each of the four descriptor blocks contributed to the model, and the predictions of pKd/i of the internal test set gave a root-mean-square error of prediction (RMSEP) of 1.65. The data were further divided according to the structural classification of the proteins, and Hi-PLS models were constructed for the resulting subclasses. The models for the four subclasses differed considerably in terms of both their ability to predict pKd/i (with RMSEPs ranging from 0.8 to 1.56) and the descriptor block that had the strongest influence. The models were validated with an external test set of 174 complexes from the 2003 release of the PDBbind database. The overall results show that the presented Hi-PLS methodology could facilitate the difficult task of predicting binding affinity.

Effects of soman inhibition and of structural differences on cholinesterase molecular dynamics: a neutron scattering study

Incoherent elastic neutron scattering experiments on members of the cholinesterase family were carried out to investigate how molecular dynamics is affected by covalent inhibitor binding and by differences in primary and quaternary structure. Tetrameric native and soman-inhibited human butyrylcholinesterase (HuBChE) as well as native dimeric Drosophila melanogaster acetylcholinesterase (DmAChE) hydrated protein powders were examined. Atomic mean-square displacements (MSDs) were found to be identical for native HuBChE and for DmAChE in the whole temperature range examined, leading to the conclusion that differences in activity and substrate specificity are not reflected by a global modification of subnanosecond molecular dynamics. MSDs of native and soman-inhibited HuBChE were identical below the thermal denaturation temperature of the native enzyme, indicating a common mean free-energy surface. Denaturation of the native enzyme is reflected by a relative increase of MSDs consistent with entropic stabilization of the unfolded state. The results suggest that the stabilization of HuBChE phosphorylated by soman is due to an increase in free energy of the unfolded state due to a decrease in entropy.

Potent Inhibitors of the Hepatitis C Virus NS3 Protease: Design and Synthesis of Macrocyclic Substrate-Based β-Strand Mimics

The virally encoded NS3 protease is essential to the life cycle of the hepatitis C virus (HCV), an important human pathogen causing chronic hepatitis, cirrhosis of the liver, and hepatocellular carcinoma. The design and synthesis of 15-membered ring beta-strand mimics which are capable of inhibiting the interactions between the HCV NS3 protease enzyme and its polyprotein substrate will be described. The binding interactions between a macrocyclic ligand and the enzyme were explored by NMR and molecular dynamics, and a model of the ligand/enzyme complex was developed.

Extraordinarily slow binding of guanosine to the Tetrahymena group I ribozyme: implications for RNA preorganization and function

The Tetrahymena ribozyme derived from the self-splicing group I intron binds a 5'-splice site analog (S) and guanosine (G), catalyzing their conversion to a 5'-exon analog (P) and GA. Herein, we show that binding of guanosine is exceptionally slow, limiting the reaction at near neutral pH. Our results implicate a conformational rearrangement on guanosine binding, likely because the binding site is not prearranged in the absence of ligand. The fast accommodation of guanosine (10(2) to 10(3) x s(-1)) and prior structural data suggest local rather than global rearrangements, raising the possibility that folding of this and perhaps other large RNAs is not fully cooperative. Guanosine binding is accelerated by addition of residues that form helices, referred to as P9.0 and P10, immediately 5' and 3' to the guanosine. These rate enhancements provide evidence for binding intermediates that have the adjacent helices formed before accommodation of guanosine into its binding site. Because the ability to form the P9.0 and P10 helices distinguishes the guanosine at the correct 3'-splice site from other guanosine residues, the faster binding of the correct guanosine can enhance specificity of 3'-splice site selection. Thus, paradoxically, the absence of a preformed binding site and the resulting slow guanosine binding can contribute to splicing specificity by providing an opportunity for the adjacent helices to increase the rate of binding of the guanosine specifying the 3'-splice site.

High activity of human butyrylcholinesterase at low pH in the presence of excess butyrylthiocholine

Butyrylcholinesterase is a serine esterase, closely related to acetylcholinesterase. Both enzymes employ a catalytic triad mechanism for catalysis, similar to that used by serine proteases such as alpha-chymotrypsin. Enzymes of this type are generally considered to be inactive at pH values below 5, because the histidine member of the catalytic triad becomes protonated. We have found that butyrylcholinesterase retains activity at pH <or= 5, under conditions of excess substrate activation. This low-pH activity appears with wild-type butyrylcholinesterase as well as with all mutants we examined: A328G, A328I, A328F, A328Y, A328W, E197Q, L286W, V288W and Y332A (residue A328 is at the bottom of the active-site gorge, near the pi-cation-binding site; E197 is next to the active-site serine S198; L286 and V288 form the acyl-binding pocket; and Y332 is a component of the peripheral anionic site). For example, the kcat value at pH 5.0 for activity in the presence of excess substrate was 32900 +/- 4400 min(-1) for wild-type, 55200 +/- 1600 min(-1) for A328F, and 28 700 +/- 700 min(-1) for A328W. This activity is titratable, with pKa values of 6.0-6.6, suggesting that the catalytic histidine is protonated at pH 5. The existence of activity when the catalytic histidine is protonated indicates that the catalytic-triad mechanism of butyrylcholinesterase does not operate for catalysis at low pH. The mechanism explaining the catalytic behaviour of butyrylcholinesterase at low pH in the presence of excess substrate remains to be elucidated.

Inhibiting viral proteases: challenges and opportunities

Inhibitor design against viral targets must take into account the peculiar characteristics of viral biology-in particular, the plasticity of their replicative machinery. This includes maturational cleavage of the polyprotein, which is mediated by virally encoded proteases. Designing against a movable target is particularly challenging, but at the same time it offers new opportunities. Here we describe our experience with the NS3/4A (NS: nonstructural) serine protease of human hepatitis C virus (HCV). By extensive use of combinatorial peptide libraries, various inhibitor types were generated, including product inhibitors, serine traps, P-P' inhibitors, and prime side inhibitors. The latter represent a first case for a serine protease. A key finding, derived from structural studies utilizing these inhibitors, was that NS3 is an induced-fit protease, requiring both the NS4A cofactor protein and the substrate to fully activate its catalytic machinery. In the absence of cofactor and/or substrate, NS3 exists in solution as a large conformational ensemble, which can be matched by a correspondingly large set of peptide inhibitors, each one stabilizing a given conformer. In the perspective of inhibiting viral proteases in general, we suggest that combinatorial ligand ensembles may be a powerful tool, to contrast the adaptive potential of the viral quasispecies.

Water structure theory and some implications for drug design

The development of theories of water structure has been hindered in the past by the difficulty of experimental measurement. Both measurement and computer modelling studies have now reached the stage where theoretical treatments of water structure are converging to a broadly acceptable model. In current understanding, water is a mixture of randomly hydrogen-bonded molecules and larger structures comprised of tetrahedral oxygen centres which, when hydrogen-bonded to each other, lead to five-membered and other rings which can aggregate to form three-dimensional structures. Evidence is taken from studies of the ices, from clathrates and other solid solutions, as well as from liquid solutions, that certain motifs occur very frequently and have relatively high stability, such as the (H2O)20 cavity-forming structure known from studies on clathrates. The implications of recent models of water structure for an understanding of biological events, including the interactions of drugs with receptors, are profound. It is becoming clear that modelling of aqueous solutions of any molecule must consider the explicit interactions with water molecules, which should not be regarded as a continuum: water itself is not a continuum. Solute molecules which possess hydrogen-bonding groups will provoke the formation of further hydrogen-bonding chains of water molecules: if these can form rings, such ringswilltend to persist longerthan chains, givingthesolute a secondary identity of associated water which may play a role in molecular recognition. Solutes that do not have hydrogen-bonding capability, or regions of solutes which are non-polar, may also produce partial cage-like water structures that are characteristic of the solute. The classification of many solutes as structure makers or structure breakers has relevance to the interactions between ligands and large biomolecules such as proteins. While it is generally accepted that sulfate and urea, respectively structure maker and breaker, may alter protein conformation through effects on water, it has not been recognised that bioactive ligands, which also change the conformation of proteins, may do so by a related, but more selective, mechanism. Very early studies of cell contents suggested that the associated water might be different from bulk water, a concept that lost support in the mid-20th century. Current theories of water structure may invite a reappraisal of this position, given the observation that structuring may extend for many molecular diameters from an ordered surface.

Efficient docking of peptides to proteins without prior knowledge of the binding site

Reliability in docking of ligand molecules to proteins or other targets is an important challenge for molecular modeling. Applications of the docking technique include not only prediction of the binding mode of novel drugs, but also other problems like the study of protein-protein interactions. Here we present a study on the reliability of the results obtained with the popular AutoDock program. We have performed systematical studies to test the ability of AutoDock to reproduce eight different protein/ligand complexes for which the structure was known, without prior knowledge of the binding site. More specifically, we look at factors influencing the accuracy of the final structure, such as the number of torsional degrees of freedom in the ligand. We conclude that the Autodock program package is able to select the correct complexes based on the energy without prior knowledge of the binding site. We named this application blind docking, as the docking algorithm is not able to "see" the binding site but can still find it. The success of blind docking represents an important finding in the era of structural genomics.

Structure of coronavirus main proteinase reveals combination of a chymotrypsin fold with an extra alpha-helical domain

The key enzyme in coronavirus polyprotein processing is the viral main proteinase, M(pro), a protein with extremely low sequence similarity to other viral and cellular proteinases. Here, the crystal structure of the 33.1 kDa transmissible gastroenteritis (corona)virus M(pro) is reported. The structure was refined to 1.96 A resolution and revealed three dimers in the asymmetric unit. The mutual arrangement of the protomers in each of the dimers suggests that M(pro) self-processing occurs in trans. The active site, comprised of Cys144 and His41, is part of a chymotrypsin-like fold that is connected by a 16 residue loop to an extra domain featuring a novel alpha-helical fold. Molecular modelling and mutagenesis data implicate the loop in substrate binding and elucidate S1 and S2 subsites suitable to accommodate the side chains of the P1 glutamine and P2 leucine residues of M(pro) substrates. Interactions involving the N-terminus and the alpha-helical domain stabilize the loop in the orientation required for trans-cleavage activity. The study illustrates that RNA viruses have evolved unprecedented variations of the classical chymotrypsin fold.

Biomimetic self-assembly of a functional asymmetrical electronic device

This paper introduces a biomimetic strategy for the fabrication of asymmetrical, three-dimensional electronic devices modeled on the folding of a chain of polypeptide structural motifs into a globular protein. Millimeter-size polyhedra-patterned with logic devices, wires, and solder dots-were connected in a linear string by using flexible wire. On self-assembly, the string folded spontaneously into two domains: one functioned as a ring oscillator, and the other one as a shift register. This example demonstrates that biomimetic principles of design and self-organization can be applied to generate multifunctional electronic systems of complex, three-dimensional architecture.

The role of hydrophobic active-site residues in substrate specificity and acyl transfer activity of penicillin acylase

Penicillin acylase of Escherichia coli catalyses the hydrolysis and synthesis of beta-lactam antibiotics. To study the role of hydrophobic residues in these reactions, we have mutated three active-site phenylalanines. Mutation of alphaF146, betaF24 and betaF57 to Tyr, Trp, Ala or Leu yielded mutants that were still capable of hydrolysing the chromogenic substrate 2-nitro-5-[(phenylacetyl)amino]-benzoic acid. Mutations on positions alphaF146 and betaF24 influenced both the hydrolytic and acyl transfer activity. This caused changes in the transferase/hydrolase ratios, ranging from a 40-fold decrease for alphaF146Y and alphaF146W to a threefold increase for alphaF146L and betaF24A, using 6-aminopenicillanic acid as the nucleophile. Further analysis of the betaF24A mutant showed that it had specificity constants (kcat/Km) for p-hydroxyphenylglycine methyl ester and phenylglycine methyl ester that were similar to the wild-type values, whereas the specificity constants for p-hydroxyphenylglycine amide and phenylglycine amide had decreased 10-fold, due to a decreased kcat value. A low amidase activity was also observed for the semisynthetic penicillins amoxicillin and ampicillin and the cephalosporins cefadroxil and cephalexin, for which the kcat values were fivefold to 10-fold lower than the wild-type values. The reduced specificity for the product and the high initial transferase/hydrolase ratio of betaF24A resulted in high yields in acyl transfer reactions.

Primary metabolism and its control in streptomycetes: a most unusual group of bacteria

Streptomycetes are Gram-positive bacteria with a unique capacity for the production of a multitude of varied and complex secondary metabolites. They also have a complex life cycle including differentiation into at least three distinct cell types. Whilst much attention has been paid to the pathways and regulation of secondary metabolism, less has been paid to the pathways and the regulation of primary metabolism, which supplies the precursors. With the imminent completion of the total genome sequence of Streptomyces coelicolor A3(2), we need to understand the pathways of primary metabolism if we are to understand the role of newly discovered genes. This review is written as a contribution to supplying these wants. Streptomycetes inhabit soil, which, because of the high numbers of microbial competitors, is an oligotrophic environment. Soil nutrient levels reflect the fact that plant-derived material is the main nutrient input; i.e. it is carbon-rich and nitrogen- and phosphate-poor. Control of streptomycete primary metabolism reflects the nutrient availability. The variety and multiplicity of carbohydrate catabolic pathways reflects the variety and multiplicity of carbohydrates in the soil. This multiplicity of pathways has led to investment by streptomycetes in pathway-specific and global regulatory networks such as glucose repression. The mechanism of glucose repression is clearly different from that in other bacteria. Streptomycetes feed by secreting complexes of extracellular enzymes that break down plant cell walls to release nutrients. The induction of these enzyme complexes is often coordinated by inducers that bear no structural relation to the substrate or product of any particular enzyme in the complex; e.g. a product of xylan breakdown may induce cellulase production. Control of amino acid catabolism reflects the relative absence of nitrogen catabolites in soil. The cognate amino acid induces about half of the catabolic pathways and half are constitutive. There are reduced instances of global carbon and nitrogen catabolite control of amino acid catabolism, which again presumably reflects the relative rarity of the catabolites. There are few examples of feedback repression of amino acid biosynthesis. Again this is taken as a reflection of the oligotrophic nature of the streptomycete ecological niche. As amino acids are not present in the environment, streptomycetes have rarely invested in feedback repression. Exceptions to this generalization are the arginine and branched-chain amino acid pathways and some parts of the aromatic amino acid pathways which have regulatory systems similar to Escherichia coli and Bacillus subtilis and other copiotrophic bacteria.

Structures of prolyl oligopeptidase substrate/inhibitor complexes. Use of inhibitor binding for titration of the catalytic histidine residue

Structure determination of the inactive S554A variant of prolyl oligopeptidase complexed with an octapeptide has shown that substrate binding is restricted to the P4-P2' region. In addition, it has revealed a hydrogen bond network of potential catalytic importance not detected in other serine peptidases. This involves a unique intramolecular hydrogen bond between the P1' amide and P2 carbonyl groups and another between the P2' amide and Nepsilon2 of the catalytic histidine 680 residue. It is argued that both hydrogen bonds promote proton transfer from the imidazolium ion to the leaving group. Another complex formed with the product-like inhibitor benzyloxycarbonyl-glycyl-proline, indicating that the carboxyl group of the inhibitor forms a hydrogen bond with the Nepsilon2 of His(680). Because a protonated histidine makes a stronger interaction with the carboxyl group, it offers a possibility of the determination of the real pK(a) of the catalytic histidine residue. This was found to be 6.25, lower than that of the well studied serine proteases. The new titration method gave a single pK(a) for prolyl oligopeptidase, whose reaction exhibited a complex pH dependence for k(cat)/K(m), and indicated that the observed pK(a) values are apparent. The procedure presented may be applicable for other serine peptidases.

A personal account of the role of peptide research in drug discovery: the case of hepatitis C

Although peptides themselves are not usually the end products of a drug discovery effort, peptide research often plays a key role in many aspects of this process. This will be illustrated by reviewing the experience of peptide research carried out at IRBM in the course of our study of hepatitis C virus (HCV). The target of our work is the NS3/4A protease, which is essential for maturation of the viral polyprotein. After a thorough examination of its substrate specificity we fine-tuned several substrate-derived peptides for enzymology studies, high-throughput screening and as fluorescent probes for secondary binding assays. In the course of these studies we made the key observation: that the protease is inhibited by its own cleavage products. Single analog and combinatorial optimization then derived potent peptide inhibitors. The crucial role of the NS4A cofactor was also addressed. NS4A is a small transmembrane protein, whose central domain is the minimal region sufficient for enzyme activation. Structural studies were performed with a peptide corresponding to the minimal activation domain, with a series of product inhibitors and with both. We found that NS3/4A is an induced fit enzyme, requiring both the cofactor and the substrate to acquire its bioactive conformation; this explained some puzzling results of 'serine-trap' type inhibitors. A more complete study on NS3 activation, however, requires the availability of the full-length NS4A protein. This was prepared by native chemical ligation, after sequence engineering to enhance its solubility; structural studies are in progress. Current work is focused on the P' region of the substrate, which, at variance with the P region, is not used for ground state binding to the enzyme and might give rise to inhibitors showing novel interactions with the enzyme.

New bivalent thrombin inhibitors with N(alpha)(methyl)arginine at the P1-position

A series of bivalent thrombin inhibitors was synthesized, consisting of a d-phenylalanyl-prolyl-N(alpha)(methyl)arginyl active site blocking segment, a fibrinogen recognition exosite inhibitor part, and a peptidic linker connecting these fragments. The methylation of the P1 amino acid led to a moderate decrease in affinity compared with the unmethylated analog. In addition, it prevented the thrombin catalyzed proteolysis, independent of the P1' amino acid used. This is a significant advantage compared to the original hirulogs, which strictly require a proline as P1' amino acid to reduce the cleavage C-terminal to the arginyl residue. Several analogs were prepared by incorporation of different P1' amino acids found in natural thrombin substrates. The most potent inhibitor was I-11 [dCha-Pro-N(Me)Arg-Thr-(Gly)5-DYEPIPEEA-Cha-dGlu] with a Ki of 37 pM. I-11 is highly selective and no inhibition of the related serine proteases trypsin, factor Xa and plasmin was observed. The stability of I-11 in human plasma in vitro was strongly improved compared to hirulog-1. In addition, a significantly reduced plasma clearance of I-11 was observed after intravenous injection in rats. Results from molecular modeling suggest that a strong reorganization of the hydrogen bonds in the active site of thrombin may result in the proteolytic stability found in this inhibitor series.

The noncatalytic beta-propeller domain of prolyl oligopeptidase enhances the catalytic capability of the peptidase domain

Prolyl oligopeptidase, which is involved in memory disorders, is a member of a new family of serine peptidases. In addition to the peptidase domain, the enzyme contains a beta-propeller, which excludes large peptides from the active site. The enzyme is inhibited with thiol reagents, possibly by reacting with Cys-255 located close to the substrate binding site. This assumption was tested with the Cys-255 --> Thr, Cys-255 --> Ala, and Cys-255 --> Ser variants of prolyl oligopeptidase. In contrast to the wild type enzyme, the Cys-255 --> Thr variant was not inhibited with N-ethylmaleimide, indicating that Cys-255, of the 16 free cysteine residues, exclusively accounts for the enzyme inhibition. Unlike the wild type enzyme that showed a doubly bell-shaped pH rate profile, the modified enzyme displayed a single bell-shaped pH dependence with benzyloxycarbonyl-Gly-Pro-naphthylamide. It was the high pH form of the enzyme that virtually disappeared with all three enzyme variants. A substantial reduction was also observed in k(cat)/K(m) for the aminobenzoyl-Ser-Pro-Phe(NO(2))-Ala-OH substrate. The high pK(a) (9.77) of Cys-255 determined by titration with N-ethylmaleimide excluded the possibility that ionization of the thiol group was responsible for generation of the two active enzyme forms. The impaired activity of the enzyme variants could be rationalized in terms of weaker binding, which manifests itself in high K(m) for substrates and high K(i) for inhibitors, like benzyloxycarbonyl-Gly-Pro-OH and aminobenzoyl-Ser-d-Pro-Phe(NO(2))-Ala-OH. It was concluded that, besides selecting substrates by size, the beta-propeller domain containing Cys-255 remarkably contributed to catalysis of the peptidase domain.

Structure of porcine pancreatic elastase complexed with FR901277, a novel macrocyclic inhibitor of elastases, at 1.6 A resolution

Human leukocyte elastase (HLE) is a serine protease that contributes to tissue destruction in various disease states-for example, in emphysema. FR901277 is a natural product isolated from the culture filtrate of Streptomyces resistomicificus and is a potent inhibitor of both HLE and porcine pancreatic elastase (PPE). FR901277 consists of four normal amino acids and three unusual amino acids, and is a unique bicyclic peptide compound. The crystal structure of PPE complexed with FR901277 has been determined at 1.6 A resolution. The Ogamma atom of Ser-195 in PPE did not form a covalent bond with FR901277, but formed a hydrogen bond with the Nvarepsilon atom of His-57. On the other hand, the portion from L-Orn(1) through dehydroxyThr(3) in FR901277 formed an antiparallel beta-sheet structure with the backbone of the active site in PPE. The S4 through S2' binding subsites in PPE were all occupied by the hydrophobic side chains of the inhibitor molecule. Especially, the ethylidene moiety of FR901277 occupied the S1 specific pocket, indicating a CH/pi interaction. In addition, the isopropyl side chain of L-Val(7) was located at the enzyme surface between the S2 and S1' pockets with several van der Waals contacts. However, the amino acid (4) residue was not involved in a significant interaction with PPE. Comparison of inhibitor structures in different environments showed that FR901277 has a highly rigid bicyclic framework; however, it can slightly change its conformation according to the circumstances. The binding mode of FR901277 at the active site of PPE was directly applicable to that in HLE, after consideration of induced fit. The structure of the PPE-FR901277 complex provided much information regarding potential sites for modification of the physicochemical properties of FR901277.

Backbone dipoles generate positive potentials in all proteins: origins and implications of the effect

Asymmetry in packing the peptide amide dipole results in larger positive than negative regions in proteins of all folding motifs. The average side chain potential in 305 proteins is 109 +/- 30 mV (2. 5 +/- 0.7 kcal/mol/e). Because the backbone has zero net charge, the non-zero potential is unexpected. The larger oxygen at the negative and smaller proton at the positive end of the amide dipole yield positive potentials because: 1) at allowed phi and psi angles residues come off the backbone into the positive end of their own amide dipole, avoiding the large oxygen; and 2) amide dipoles with their carbonyl oxygen surface exposed and amine proton buried make the protein interior more positive. Twice as many amides have their oxygens exposed than their amine protons. The distribution of acidic and basic residues shows the importance of the bias toward positive backbone potentials. Thirty percent of the Asp, Glu, Lys, and Arg are buried. Sixty percent of buried residues are acids, only 40% bases. The positive backbone potential stabilizes ionization of 20% of the acids by >3 pH units (-4.1 kcal/mol). Only 6.5% of the bases are equivalently stabilized by negative regions. The backbone stabilizes bound anions such as phosphates and rarely stabilizes bound cations.

Structural characterization of the interactions of optimized product inhibitors with the N-terminal proteinase domain of the hepatitis C virus (HCV) NS3 protein by NMR and modelling studies

The interactions of peptide inhibitors, obtained by the optimization of N-terminal cleavage products of natural substrates, with the protease of human hepatitis C virus (HCV) are characterized by NMR and modelling studies. The S-binding region of the enzyme and the bound conformation of the ligands are experimentally determined. The NMR data are then used as the experimental basis for modelling studies of the structure of the complex. The S-binding region involves the loop connecting strands E2 and F2, and appears shallow and solvent-exposed. The ligand binds in an extended conformation, forming an antiparallel beta-sheet with strand E2 of the protein, with the P1 carboxylate group in the oxyanion hole.

Structural and energetic determinants of the S1-site specificity in serine proteases

In recent years the number of determined three-dimensional structures of serine proteases that are accompanied by detailed mutational studies has grown rapidly. In particular, spatial structures have been described for enzymes involved in processes of critical medical significance, often related to severe pathophysiological diseases. There has also been significant progress in the understanding of the structural grounds for the substrate specificity of serine proteases. This review is concerned mainly with primary structural determinants of the S1 specificity, the crucial component of substrate selectivity, often in relation to more distant specificity elements, which cooperatively influence the S1 site.