Traian Sulea

Solvated interaction energy (SIE) for scoring protein-ligand binding affinities. 1. Exploring the parameter space.

We present a binding free energy function that consists of force field terms supplemented by solvation terms. We used this function to calibrate the solvation model along with the binding interaction terms in a self-consistent manner. The motivation for this approach was that the solute dielectric-constant dependence of calculated hydration gas-to-water transfer free energies is markedly different from that of binding free energies (J. Comput. Chem. 2003, 24, 954). Hence, we sought to calibrate directly the solvation terms in the context of a binding calculation. The five parameters of the model were systematically scanned to best reproduce the absolute binding free energies for a set of 99 protein-ligand complexes. We obtained a mean unsigned error of 1.29 kcal/mol for the predicted absolute binding affinity in a parameter space that was fairly shallow near the optimum. The lowest errors were obtained with solute dielectric values of Din = 20 or higher and scaling of the intermolecular van der Waals interaction energy by factors ranging from 0.03 to 0.15. The high apparent Din and strong van der Waals scaling may reflect the anticorrelation of the change in solvated potential energy and configurational entropy, that is, enthalpy-entropy compensation in ligand binding (Biophys. J. 2004, 87, 3035-3049). Five variations of preparing the protein-ligand data set were explored in order to examine the effect of energy refinement and the presence of bound water on the calculated results. We find that retaining water in the final protein structure used for calculating the binding free energy is not necessary to obtain good results; that is the continuum solvation model is sufficient. Virtual screening enrichment studies on estrogen receptor and thymidine kinase showed a good ability of the binding free energy function to recover true hits in a collection of decoys.

The Papain-Like Protease from the Severe Acute Respiratory Syndrome Coronavirus Is a Deubiquitinating Enzyme

ABSTRACT The severe acute respiratory syndrome coronavirus papain-like protease (SARS-CoV PLpro) is involved in the processing of the viral polyprotein and, thereby, contributes to the biogenesis of the virus replication complex. Structural bioinformatics has revealed a relationship for the SARS-CoV PLpro to herpesvirus-associated ubiquitin-specific protease (HAUSP), a ubiquitin-specific protease, indicating potential deubiquitinating activity in addition to its function in polyprotein processing (T. Sulea, H. A. Lindner, E. O. Purisima, and R. Menard, J. Virol. 79:4550-4551, 2005). In order to confirm this prediction, we overexpressed and purified SARS-CoV PLpro (amino acids [aa]1507 to 1858) from Escherichia coli. The purified enzyme hydrolyzed ubiquitin-7-amino-4-methylcoumarin (Ub-AMC), a general deubiquitinating enzyme substrate, with a catalytic efficiency of 13,100 M−1s−1, 220-fold more efficiently than the small synthetic peptide substrate Z-LRGG-AMC, which incorporates the C-terminal four residues of ubiquitin. In addition, SARS-CoV PLpro was inhibited by the specific deubiquitinating enzyme inhibitor ubiquitin aldehyde, with an inhibition constant of 210 nM. The purified SARS-CoV PLpro disassembles branched polyubiquitin chains with lengths of two to seven (Ub2-7) or four (Ub4) units, which involves isopeptide bond cleavage. SARS-CoV PLpro processing activity was also detected against a protein fused to the C terminus of the ubiquitin-like modifier ISG15, both in vitro using the purified enzyme and in HeLa cells by coexpression with SARS-CoV PLpro (aa 1198 to 2009). These results clearly establish that SARS-CoV PLpro is a deubiquitinating enzyme, thereby confirming our earlier prediction. This unexpected activity for a coronavirus papain-like protease suggests a novel viral strategy to modulate the host cell ubiquitination machinery to its advantage.

Molecular dynamics-solvated interaction energy studies of protein-protein interactions: the MP1-p14 scaffolding complex.

Using the MP1-p14 scaffolding complex from the mitogen-activated protein kinase signaling pathway as model system, we explored a structure-based computational protocol to probe and characterize binding affinity hot spots at protein-protein interfaces. Hot spots are located by virtual alanine-scanning consensus predictions over three different energy functions and two different single-structure representations of the complex. Refined binding affinity predictions for select hot-spot mutations are carried out by applying first-principle methods such as the molecular mechanics generalized Born surface area (MM-GBSA) and solvated interaction energy (SIE) to the molecular dynamics (MD) trajectories for mutated and wild-type complexes. Here, predicted hot-spot residues were actually mutated to alanine, and crystal structures of the mutated complexes were determined. Two mutated MP1-p14 complexes were investigated, the p14(Y56A)-mutated complex and the MP1(L63A,L65A)-mutated complex. Alternative ways to generate MD ensembles for mutant complexes, not relying on crystal structures for mutated complexes, were also investigated. The SIE function, fitted on protein-ligand binding affinities, gave absolute binding affinity predictions in excellent agreement with experiment and outperformed standard MM-GBSA predictions when tested on the MD ensembles of Ras-Raf and Ras-RalGDS protein-protein complexes. For wild-type and mutant MP1-p14 complexes, SIE predictions of relative binding affinities were supported by a yeast two-hybrid assay that provided semiquantitative relative interaction strengths. Results on the MP1-mutated complex suggested that SIE predictions deteriorate if mutant MD ensembles are approximated by just mutating the wild-type MD trajectory. The SIE data on the p14-mutated complex indicated feasibility for generating mutant MD ensembles from mutated wild-type crystal structure, despite local structural differences observed upon mutation. For energetic considerations, this would circumvent costly needs to produce and crystallize mutated complexes. The sensitized protein-protein interface afforded by the p14(Y56A) mutation identified here has practical applications in screening-based discovery of first-generation small-molecule hits for further development into specific modulators of the mitogen-activated protein kinase signaling pathway.

Deubiquitination, a New Function of the Severe Acute Respiratory Syndrome Coronavirus Papain-Like Protease?

A new coronavirus has been identified as the infectious agent of severe acute respiratory syndrome (SARS). Although SARS was successfully contained by quarantine measures, the reconstructed itinerary of the virus through 30 countries and its high mortality rate illustrate the global threat that this newly emerging disease represents (17). During the expression of the SARS coronavirus (SCoV) genome, two viral cysteine proteases, a papain-like protease (PLpro) and a chymotrypsin-like protease (3CLpro), process the encoded polyprotein precursor to release most of the proteins required for virus replication. PLpro refers to a domain of nonstructural protein 3, whose boundaries are defined by homology to the papain-like fold (7). The PLpro domain can be regarded as the catalytic core behind PLpro-mediated cleavages, even though processing by PLpros has been reported to be modulated by additional amino acid residues outside of these boundaries (15, 18). The SCoV utilizes a single PLpro, whereas most coronaviruses contain two paralogous enzymes, termed PL1pro and PL2pro (14). PLpros in general are not as well characterized as 3CLpros and have not generated as much interest as pharmaceutical targets. However, further structure-to-function annotations might refocus the attention on PLpros. For this purpose, we mined the current Protein Data Bank (PDB) content using the Structure Prediction Meta Server (http://www.bioinfo.pl/meta), which assembles state-of-the-art fold recognition methods and provides a consensus sequence-to-structure hyperscore with the 3D-Jury method (4). The only highly reliable prediction (3D-Jury score > 100) obtained for the SCoV PLpro sequence from K1632 to E1847 was the structure of the catalytic core domain of the herpesvirus-associated ubiquitin-specific protease (HAUSP), also known as USP7 (8). There is compelling evidence for the relevance of a structural relationship between SCoV PLpro and HAUSP (Fig. ​(Fig.1).1). First, HAUSP is a cysteine protease with a finger domain inserted between the two subdomains of a papain-like fold (8). This structure mirrors the general architecture proposed for the members of the coronaviral PLpro family, which now includes SCoV PLpro, with a Zn ribbon domain inserted in the middle of a papain-like protease domain (7). Our structural bioinformatics data support a circularly permuted rather than a classical Zn ribbon domain for SCoV PLpro, as was recently reported for HAUSP and related enzymes (12). Second, as a deubiquitinating enzyme, HAUSP recognizes the C-terminal ubiquitin sequence, LRGG, which matches the narrow specificity profile of SCoV PLpro (LXGG) derived from the three PLpro-processing sites of the polyprotein (6, 16). The modeled SCoV PLpro binding site is highly complementary to the LXGG sequence and establishes extensive hydrogen bonding with the substrate main chain. In particular, significant occlusions of the S1 subsite (due to N1649 and L1702) and S2 subsite (due to Y1804 and Y1813) account for the strict specificity for diglycine at substrate positions P1 and P2 (Fig. ​(Fig.1C).1C). The structural signatures for strict specificity are present in human coronavirus (HCoV 229E) PL1pro and PL2pro as well as in mouse hepatitis virus (MHV) PL2pro but not in MHV PL1pro (Fig. ​(Fig.1D).1D). This fact correlates with the available specificity data showing a preference for the large arginine residue at the P2 position in the case of MHV PL1pro (2, 3, 9) and with the observation that the irreversible inhibitor E-64d, which contains a bulky P2 residue (leucine), inhibits MHV PL1pro (5, 11) but not MHV PL2pro (10). FIG. 1. Structure-to-function relationships between SCoV PLpro and HAUSP. (A) Three-dimensional model of SCoV PLpro. The protease domain is rendered in cyan, and the circularly permuted Zn ribbon domain is shown in red. The zinc ion is shown as a magenta sphere, ... Structural similarities to HAUSP suggest that, in addition to polyprotein processing activity, SCoV PLpro might possess deubiquitinating activity (including deconjugation of other ubiquitin-like modifiers), as was observed for an adenoviral protease (1). This unexpected activity prediction raises provocative hypotheses regarding the ability of the SARS virus to evade cellular defense mechanisms. For example, it is tempting to speculate that ISG15 deconjugation by PLpro allows the SARS virus to counteract protein ISGylation, an interferon-induced process that may contribute to innate immunity to viral infection (13). The deubiquitination function would greatly impact the value of PLpro as a therapeutic target and provide a framework for the development of antivirals to treat SARS. Strategies for the design of inhibitors of SCoV PLpro must also take into consideration the potentially overlapping specificity of this protease with those of cellular deubiquitinating enzymes. Our finding suggests the performance of follow-up experiments that will increase the understanding of the functional roles of papain-like proteases in the viral life cycles of SCoV and related viruses.

Characterization of a Pseudomonad 2-Nitrobenzoate Nitroreductase and Its Catabolic Pathway-Associated 2-Hydroxylaminobenzoate Mutase and a Chemoreceptor Involved in 2-Nitrobenzoate Chemotaxis

Pseudomonas fluorescens strain KU-7 is a prototype microorganism that metabolizes 2-nitrobenzoate (2-NBA) via the formation of 3-hydroxyanthranilate (3-HAA), a known antioxidant and reductant. The initial two steps leading to the sequential formation of 2-hydroxy/aminobenzoate and 3-HAA are catalyzed by a NADPH-dependent 2-NBA nitroreductase (NbaA) and 2-hydroxylaminobenzoate mutase (NbaB), respectively. The 216-amino-acid protein NbaA is 78% identical to a plasmid-encoded hypothetical conserved protein of Polaromonas strain JS666; structurally, it belongs to the homodimeric NADH:flavin mononucleotide (FMN) oxidoreductase-like fold family. Structural modeling of complexes with the flavin, coenzyme, and substrate suggested specific residues contributing to the NbaA catalytic activity, assuming a ping-pong reaction mechanism. Mutational analysis supports the roles of Asn40, Asp76, and Glu113, which are predicted to form the binding site for a divalent metal ion implicated in FMN binding, and a role in NADPH binding for the 10-residue insertion in the beta5-alpha2 loop. The 181-amino-acid sequence of NbaB is 35% identical to the 4-hydroxylaminobenzoate lyases (PnbBs) of various 4-nitrobenzoate-assimilating bacteria, e.g., Pseudomonas putida strain TW3. Coexpression of nbaB with nbaA in Escherichia coli produced a small amount of 3-HAA from 2-NBA, supporting the functionality of the nbaB gene. We also showed by gene knockout and chemotaxis assays that nbaY, a chemoreceptor NahY homolog located downstream of the nbaA gene, is responsible for strain KU-7 being attracted to 2-NBA. NbaY is the first chemoreceptor in nitroaromatic metabolism to be identified, and this study completes the gene elucidation of 2-NBA metabolism that is localized within a 24-kb chromosomal locus of strain KU-7.

High incidence of ubiquitin‐like domains in human ubiquitin‐specific proteases

Ubiquitin‐specific proteases (USPs) emerge as key regulators of numerous cellular processes and account for the bulk of human deubiquitinating enzymes (DUBs). Their modular structure, mostly annotated by sequence homology, is believed to determine substrate recognition and subcellular localization. Currently, a large proportion of known human USP sequences are not annotated either structurally or functionally, including regions both within and flanking their catalytic cores. To extend the current understanding of human USPs, we applied consensus fold recognition to the unannotated content of the human USP family. The most interesting discovery was the marked presence of reliably predicted ubiquitin‐like (UBL) domains in this family of enzymes. The UBL domain thus appears to be the most frequently occurring domain in the human USP family, after the characteristic catalytic domain. The presence of multiple UBL domains per USP protein, as well as of UBL domains embedded in the USP catalytic core, add to the structural complexity currently recognized for many DUBs. Possible functional roles of the newly uncovered UBL domains of human USPs, including proteasome binding, and substrate and protein target specificities, are discussed. Proteins 2007.

Improvement of the Thermostability and Activity of a Pectate Lyase by Single Amino Acid Substitutions, Using a Strategy Based on Melting-Temperature-Guided Sequence Alignment

ABSTRACT In the vast number of random mutagenesis experiments that have targeted protein thermostability, single amino acid substitutions that increase the apparent melting temperature (Tm) of the enzyme more than 1 to 2°C are rare and often require the creation of a large library of mutated genes. Here we present a case where a single beneficial mutation (R236F) of a hemp fiber-processing pectate lyase of Xanthomonas campestris origin (PLXc) produced a 6°C increase in Tm and a 23-fold increase in the half-life at 45°C without compromising the enzymes catalytic efficiency. This success was based on a variation of sequence alignment strategy where a mesophilic amino acid sequence is matched with the sequences of its thermophilic counterparts that have established Tm values. Altogether, two-thirds of the nine targeted single amino acid substitutions were found to have effects either on the thermostability or on the catalytic activity of the enzyme, evidence of a high success rate of mutation without the creation of a large gene library and subsequent screening of clones. Combination of R236F with another beneficial mutation (A31G) resulted in at least a twofold increase in specific activity while preserving the improved Tm value. To understand the structural basis for the increased thermal stability or activity, the variant R236F and A31G R236F proteins and wild-type PLXc were purified and crystallized. By structure analysis and computational methods, hydrophobic desolvation was found to be the driving force for the increased stability with R236F.

The crystal structure of the formiminotransferase domain of formiminotransferase-cyclodeaminase: implications for substrate channeling in a bifunctional enzyme.

BACKGROUND The bifunctional enzyme formiminotransferase-cyclodeaminase (FTCD) contains two active sites at different positions on the protein structure. The enzyme binds a gamma-linked polyglutamylated form of the tetrahydrofolate substrate and channels the product of the transferase reaction from the transferase active site to the cyclodeaminase active site. Structural studies of this bifunctional enzyme and its monofunctional domains will provide insight into the mechanism of substrate channeling and the two catalytic reactions. RESULTS The crystal structure of the formiminotransferase (FT) domain of FTCD has been determined in the presence of a product analog, folinic acid. The overall structure shows that the FT domain comprises two subdomains that adopt a novel alpha/beta fold. Inspection of the folinic acid binding site reveals an electrostatic tunnel traversing the width of the molecule. The distribution of charged residues in the tunnel provides insight into the possible mode of substrate binding and channeling. The electron density reveals that the non-natural stereoisomer, (6R)-folinic acid, binds to the protein; this observation suggests a mechanism for product release. In addition, a single molecule of glycerol is bound to the enzyme and indicates a putative binding site for formiminoglutamate. CONCLUSIONS The structure of the FT domain in the presence of folinic acid reveals a possible novel mechanism for substrate channeling. The position of the folinic acid and a bound glycerol molecule near to the sidechain of His82 suggests that this residue may act as the catalytic base required for the formiminotransferase mechanism.

Catalytic mechanism of heparinase II investigated by site-directed mutagenesis and the crystal structure with its substrate.

Heparinase II (HepII) is an 85-kDa dimeric enzyme that depolymerizes both heparin and heparan sulfate glycosaminoglycans through a β-elimination mechanism. Recently, we determined the crystal structure of HepII from Pedobacter heparinus (previously known as Flavobacterium heparinum) in complex with a heparin disaccharide product, and identified the location of its active site. Here we present the structure of HepII complexed with a heparan sulfate disaccharide product, proving that the same binding/active site is responsible for the degradation of both uronic acid epimers containing substrates. The key enzymatic step involves removal of a proton from the C5 carbon (a chiral center) of the uronic acid, posing a topological challenge to abstract the proton from either side of the ring in a single active site. We have identified three potential active site residues equidistant from C5 and located on both sides of the uronate product and determined their role in catalysis using a set of defined tetrasaccharide substrates. HepII H202A/Y257A mutant lost activity for both substrates and we determined its crystal structure complexed with a heparan sulfate-derived tetrasaccharide. Based on kinetic characterization of various mutants and the structure of the enzyme-substrate complex we propose residues participating in catalysis and their specific roles.

Evolutionary Reshaping of Fungal Mating Pathway Scaffold Proteins

ABSTRACT Scaffold proteins play central roles in the function of many signaling pathways. Among the best-studied examples are the Ste5 and Far1 proteins of the yeast Saccharomyces cerevisiae. These proteins contain three conserved modules, the RING and PH domains, characteristic of some ubiquitin-ligating enzymes, and a vWA domain implicated in protein-protein interactions. In yeast, Ste5p regulates the mating pathway kinases while Far1p coordinates the cellular polarity machinery. Within the fungal lineage, the Basidiomycetes and the Pezizomycetes contain a single Far1-like protein, while several Saccharomycotina species, belonging to the CTG (Candida) clade, contain both a classic Far1-like protein and a Ste5-like protein that lacks the vWA domain. We analyzed the function of C. albicans Ste5p (Cst5p), a member of this class of structurally distinct Ste5 proteins. CST5 is essential for mating and still coordinates the mitogen-activated protein (MAP) kinase (MAPK) cascade elements in the absence of the vWA domain; Cst5p interacts with the MEK kinase (MEKK) C. albicans Ste11p (CaSte11p) and the MAPK Cek1 as well as with the MEK Hst7 in a vWA domain-independent manner. Cst5p can homodimerize, similar to Ste5p, but can also heterodimerize with Far1p, potentially forming heteromeric signaling scaffolds. We found direct binding between the MEKK CaSte11p and the MEK Hst7p that depends on a mobile acidic loop absent from S. cerevisiae Ste11p but related to the Ste7-binding region within the vWA domain of Ste5p. Thus, the fungal lineage has restructured specific scaffolding modules to coordinate the proteins required to direct the gene expression, polarity, and cell cycle regulation essential for mating. IMPORTANCE The mitogen-activated protein (MAP) kinase cascade is an extensively used signaling module in eukaryotic cells, and the ability to regulate these modules is critical for ensuring proper responses to a wide variety of stimuli. One way that cells regulate this signaling module is through scaffold proteins that insulate related pathways against cross talk, improve signaling efficiency, and ensure that signals are connected to the correct response. The Ste5 scaffold of the S. cerevisiae mating response is a well-studied representative of this class of proteins. Using bioinformatics, structural modeling, and molecular genetic approaches, we have investigated the equivalent scaffold in the pathogenic yeast Candida albicans. We show that the C. albicans protein is structurally distinct from that of Saccharomyces cerevisiae but still provides similar functions. Increases in pathway complexity have been associated with changes in scaffold connectivity, and overall, the tethering capacity of the scaffolds has been more conserved than their structural organization. The mitogen-activated protein (MAP) kinase cascade is an extensively used signaling module in eukaryotic cells, and the ability to regulate these modules is critical for ensuring proper responses to a wide variety of stimuli. One way that cells regulate this signaling module is through scaffold proteins that insulate related pathways against cross talk, improve signaling efficiency, and ensure that signals are connected to the correct response. The Ste5 scaffold of the S. cerevisiae mating response is a well-studied representative of this class of proteins. Using bioinformatics, structural modeling, and molecular genetic approaches, we have investigated the equivalent scaffold in the pathogenic yeast Candida albicans. We show that the C. albicans protein is structurally distinct from that of Saccharomyces cerevisiae but still provides similar functions. Increases in pathway complexity have been associated with changes in scaffold connectivity, and overall, the tethering capacity of the scaffolds has been more conserved than their structural organization.