Water-soluble quinoprotein glucose dehydrogenase (PQQGDH-B) from Acinetobacter calcoaceticus has a great potential for application as a glucose sensor constituent. Because this enzyme shows no activity in its monomeric form, correct quaternary structure is essential for the formation of active enzyme. We have previously reported on the increasing of the stability of PQQGDH-B by preventing the subunit dissociation. Previous studies were based on decreasing the entropy of quaternary structure dissociation but not on increasing the interaction between the two subunits. We therefore attempted to introduce a hydrophobic interaction in the dimeric interface to increase the stability of PQQGDH-B.
Amino acid residues Asn340 and Tyr418 face each other at the dimer interface of PQQGDH-B, however no interaction exists between their side chains. We simultaneously substituted Asn340 to Phe and Tyr418 to Phe or Ile, to create the two mutants Asn340Phe/Tyr418Phe and Asn340Phe/Tyr418Ile. Furthermore, residues Leu280, Val282 and Val342 form a hydrophobic region that faces, on the other subunit, residues Thr416 and Thr417, again without any specific interaction. We simultaneously substituted Thr416 and Thr417 to Val, to create the mutant Thr416Val/Thr417Val. The temperatures resulting in lose of half of the initial activity of the constructed mutants were increased by 3–4°C higher over wild type. All mutants showed 2-fold higher thermal stability at 55°C than the wild-type enzyme, without decreasing their catalytic activities. From the 3D models of all the mutant enzymes, the predicted binding energies were found to be significantly greater that in the wild-type enzyme, consistent with the increases in thermal stabilities.
We have achieved via site-directed mutagenesis the improvement of the thermal stability of PQQGDH-B by increasing the dimer interface interaction. Through rational design based on the quaternary structure of the enzyme, we selected residues located at the dimer interface that do not contribute to the intersubunit interaction. By substituting these residues to hydrophobic ones, the thermal stability of PQQGDH-B was increased without decreasing its catalytic activity.
In MS2 assembly of phage particles results from an interaction between a coat protein dimer and a stem-loop of the RNA genome (the operator hairpin). Amino acid residues Thr45, which is universally conserved among the small RNA phages, and Thr59 are part of the specific RNA binding pocket and interact directly with the RNA; the former through a hydrogen bond, the latter through hydrophobic contacts. The crystal structures of MS2 protein capsids formed by mutants Thr45Ala and Thr59Ser, both with and without the 19 nt wild-type operator hairpin bound, are reported here. The RNA hairpin binds to these mutants in a similar way to its binding to wild-type protein. In a companion paper both mutants are shown to be deficient in RNA binding in an in vivo assay, but in vitro the equilibrium dissociation constant is significantly higher than wild-type for the Thr45Ala mutant. The change in binding affinity of the Thr45Ala mutant is probably a direct consequence of removal of direct hydrogen bonds between the protein and the RNA. The properties of the Thr59Ser mutant are more difficult to explain, but are consistent with a loss of non-polar contact.
Peptide binding to class II MHC protein is commonly viewed as a combination of discrete anchor residue preferences for pockets 1, 4, 6/7, and 9. However, previous studies have suggested cooperative effects during the peptide binding process. Investigation of the DRB1*0901 binding motif demonstrated a clear interaction between peptide binding pockets 6 and 9. In agreement with prior studies, pockets 1 and 4 exhibited clear binding preferences. Previously uncharacterized pockets 6 and 7 accommodated a wide variety of residues. However, although it was previously reported that pocket 9 is completely permissive, several substitutions at this position were unable to bind. Structural modeling revealed a probable interaction between pockets 6 and 9 through β9Lys. Additional binding studies with doubly substituted peptides confirmed that the amino acid bound within pocket 6 profoundly influences the binding preferences for pocket 9 of DRB1*0901, causing complete permissiveness of pocket 9 when a small polar residue is anchored in pocket 6 but accepting relatively few residues when a basic residue is anchored in pocket 6. The β9Lys residue is unique to DR9 alleles. However, similar studies with doubly substituted peptides confirmed an analogous interaction effect for DRA1/B1*0301, a β9Glu allele. Accounting for this interaction resulted in improved epitope prediction. These findings provide a structural explanation for observations that an amino acid in one pocket can influence binding elsewhere in the MHC class II peptide binding groove.
HLA-B27 subtypes share many structural features, including their pocket B, which interacts with a conserved Arg residue at the second position of B*2705-bound peptides. Subtypes differ among each other at other locations in the peptide binding site. In this study, metabolic labeling and radiochemical pool sequencing were used to address the following issues: (a) presence of the Arg 2 (R2) motif among peptides bound to the various HLA-B27 subtypes; (b) influence of mutations inside and outside pocket B on this motif; and (c) the degree of similarity among the peptide pools bound to the various B27 subtypes. Sequencing of Arg-labeled peptide pools extracted from B*2701 to B*2706, and from two site-directed mutants of B*2705 with changes outside pocket B, indicated that all of these molecules bind peptides with Arg at position 2. Peptides from several mutants with changes altering the structure of pocket B, and from one mutant at the pocket B rim, also retained the R2 motif. However, this was absent in the peptide pool extracted from the M45 mutant, in which the negative charge of pocket B, conferred to HLA-B27 by Glu45, was canceled. These results indicate that alterations outside pocket B, and even disruption of the network of hydrogen bonds that stabilizes Arg binding in pocket B, do not impair binding of peptides bearing the R2 motif, but a nonconservative substitution at position 45 does. As a substantial fraction of anti-B*2705 cytotoxic T lymphocyte (CTL) clones crossreact with the M45 mutant (Villadangos, J., B. Galocha, D. Lopez, V. Calvo, and J. A. Lopez de Castro. 1992. J. Immunol. 149:505) this result suggest that determinant mimicry by nonidentical peptides may frequently account for unexpected CTL crossreactions. Metabolic labeling with various other amino acids and radiochemical sequencing revealed similarities, but also substantial differences, among the peptide pools from the various HLA-B27 subtypes. This strongly suggests that many peptides bind to multiple subtypes, but significant subsets of peptides bound to a given HLA-B27 subtype do not bind to other subtypes or do so with greatly altered efficiency. These results indicate the importance of polymorphism outside pocket B in modulating peptide binding to HLA-B27.
The MS2 RNA operator capsid offers an unparalleled opportunity to study sequence-specific protein-protein and RNA-protein interactions in molecular detail. RNA molecules encompassing the minimal translational operator recognition elements can be soaked into crystals of RNA-free coat protein shells, allowing the RNA to access the interior of the capsids and make contact with the operator binding sites. Correct interpretation of these structural studies depends critically on functional analysis in solution to confirm that the interactions seen in the crystal are not an artefact of the unusual approach used to generate the RNA-protein complexes. Here we present a series of in vivo and in vitro functional assays, using coat proteins carrying single amino acid substitutions at residues which either interact with the operator RNA or are involved in stabilizing the conformation of the FG loop, the site of the major quasi-equivalent conformational change. Variant operator RNAs have been assayed for coat protein affinity in vitro. The results reveal the robustness of the operator-coat protein interaction and the requirement for both halves of a protein dimer to contact RNA in order to achieve tight binding. They also suggest that there may be a direct link between the conformation of the FG loop and RNA binding.
The α-helical coiled-coil is one of the most common oligomerization motifs found in both native and engineered proteins. To better understand the stability and dynamics of coiled-coil motifs, including those modified by fluorination, several fluorinated and non-fluorinated parallel dimeric coiled-coil protein structures were designed and modeled. We also attempt to investigate how changing the length and geometry of the important stabilizing salt bridges influences the coiled-coil protein structure. Molecular dynamics (MD) and free energy simulations with AMBER employed a particle mesh Ewald treatment of the electrostatics in explicit TIP3P solvent with balanced force field treatments. Preliminary studies with legacy force fields (ff94, ff96, ff99) show a profound instability of the coiled-coil structures in short MD simulation. Significantly better behavior is evident with the more balanced ff99SB and ff03 protein force fields. Overall, the results suggest that the coiled-coil structures can readily accommodate the larger acidic arginine or S-2,7-diaminoheptanedoic acid mutants in the salt bridge, whereas substitution of the smaller L-ornithine residue leads to rapid disruption of the coiled-coil structure on the MD simulation time scale. This structural distortion of the secondary structure allows both the formation of large hydration pockets proximal to the charged groups and within the hydrophobic core. Moreover, the increased structural fluctuations and movement lead to a decrease in the water occupancy lifetimes in the hydration pockets. In contrast, analysis of the hydration in the stable dimeric coiled coils shows high occupancy water sites along the backbone residues with no water occupancy in the hydrophobic core, although transitory water interactions with the salt bridge residues are evident. The simulations of the fluorinated coiled-coils suggest that in some cases fluorination electrostatically stabilizes the intermolecular coiled-coil salt bridges. Structural analyses also reveal different side chain rotamer preferences for leucine compared to 5,5,5,5′,5′,5′-hexafluoroleucine mutants. These observed differences in the side chain rotamer populations suggest differential changes in the side chain conformational entropy upon coiled-coil formation when the protein is fluorinated. The free energy of hydration of the isolated 5,5,5,5′,5′,5′-hexafluoroleucine amino acid is calculated to be 1.1 kcal/mol less stable than leucine; this hydrophobic penalty in the monomer may provide a driving force for coiled-coil dimer formation. Estimation of the ellipticity at 222 nm from a series of snapshots from the MD simulations with DicroCalc show distinct increases in the ellipticity when the coiled-coil is fluorinated which suggests that the helicity in the folded coiled-coils is greater when fluorinated.
computational chemistry; free energy of hydration; 5, 5, 5, 5′, 5′, 5′-hexafluoroleucine; thermodynamic integration; rotamers
The present study reports distinct dynamic consequences for the T- and R-states of human normal adult hemoglobin (Hb A) due to the binding of a heterotropic allosteric effector, inositol hexaphosphate (IHP). A nuclear magnetic resonance (NMR) technique based on modified transverse relaxation optimized spectroscopy (TROSY) has been used to investigate the effect of conformational exchange of Hb A in both deoxy and CO forms, in the absence and presence of IHP, at 14.1 and 21.1 T, and at 37 °C. Our results show that the majority of the polypeptide backbone amino acid residues of deoxy- and carbonmonoxy-forms of Hb A in the absence of IHP is not mobile on the ·s-ms time scale, with the exception of several amino acid residues, that is, · 109Val and · 132Lys in deoxy-Hb A, and R40Lys in HbCO A. The mobility of R40Lys in HbCO A can be explained by the crystallographic data showing that the H-bond between R40Lys and · 146His in deoxy-Hb A is absent in HbCO A. However, the conformational exchange of ·109Val, which is located in the intradimer (R1· 1 or R2· 2) interface, is not consistent with the crystallographic observations that show rigid packing at this site. IHP binding appears to rigidify R40Lys in HbCO A, but does not significantly affect the flexibility of · 109Val in deoxy-Hb A. In the presence of IHP, several amino acid residues, especially those at the interdimer (R1 · 2 or R2 · 1) interface of HbCO A, exhibit significant conformational exchange. The affected residues include the proximal ·92His in the ·-heme pocket, as well as some other residues located in the flexible joint (·C helix-RFG corner) and switch (RC helix-·FG corner) regions that play an important role in the dimer-dimer rotation of Hb during the oxygenation process. These findings suggest that, upon IHP binding, HbCO A undergoes a conformational fluctuation near the R-state but biased toward the T-state, apparently along the trajectory of its allosteric transition, accompanied by structural fluctuations in the heme pocket of the ·-chain. In contrast, no significant perturbation of the dynamic features on the ms-·s time scale has been observed upon IHP binding to deoxy-Hb A. We propose that the allosteric effector-induced quaternary structural fluctuation may contribute to the reduced ligand affinity of ligated hemoglobin. Conformational exchange mapping of the ·-chain of HbCO A observed at 21.1 T shows significantly increased scatter in the chemical exchange contribution to the transverse relaxation rate (Rex) values, relative to those at lower fields, due to the enhanced effect of the local chemical shift anisotropy (CSA) fluctuation. A spring-on-scissors model is proposed to interpret the dynamic phenomena induced by the heterotropic effector, IHP.
Protein kinases are key signaling enzymes that catalyze the transfer of γ-phosphate from an ATP molecule to a phospho-accepting residue in the substrate. Unraveling the molecular features that govern the preference of kinases for particular residues flanking the phosphoacceptor is important for understanding kinase specificities toward their substrates and for designing substrate-like peptidic inhibitors. We applied ANCHORSmap, a new fragment-based computational approach for mapping amino acid side chains on protein surfaces, to predict and characterize the preference of kinases toward Arginine binding. We focus on positions P−2 and P−5, commonly occupied by Arginine (Arg) in substrates of basophilic Ser/Thr kinases. The method accurately identified all the P−2/P−5 Arg binding sites previously determined by X-ray crystallography and produced Arg preferences that corresponded to those experimentally found by peptide arrays. The predicted Arg-binding positions and their associated pockets were analyzed in terms of shape, physicochemical properties, amino acid composition, and in-silico mutagenesis, providing structural rationalization for previously unexplained trends in kinase preferences toward Arg moieties. This methodology sheds light on several kinases that were described in the literature as having non-trivial preferences for Arg, and provides some surprising departures from the prevailing views regarding residues that determine kinase specificity toward Arg. In particular, we found that the preference for a P−5 Arg is not necessarily governed by the 170/230 acidic pair, as was previously assumed, but by several different pairs of acidic residues, selected from positions 133, 169, and 230 (PKA numbering). The acidic residue at position 230 serves as a pivotal element in recognizing Arg from both the P−2 and P−5 positions.
Protein kinases are key signaling enzymes and major drug targets that catalyze the transfer of phosphate group to a phospho-accepting residue in the substrate. Unraveling molecular features that govern the preference of kinases for particular residues flanking the phosphoacceptor (substrate consensus sequence, SCS) is important for understanding kinase-substrates specificities and for designing peptidic inhibitors. Current methods used to predict this set of essential residues usually rely on linking between experimentally determined SCSs to kinase sequences. As such, these methods are less sensitive when specificity is dictated by subtle or kinase-unique sequence/structural features. In this study, we took a different approach for studying kinases specificities, by applying a new fragment-based method for mapping amino acid side chains on protein surfaces. We predicted and characterized the preference of Ser/Thr kinases toward Arginine binding, using the unbound kinase structures. The method produced high quality predictions and was able to provide novel insights and interesting departures from the prevailing views regarding the specificity-determining elements governing specificity toward Arginine. This work paves the way for studying the kinase binding preferences for other amino acids, for predicting protein-peptide structures, for facilitating the design of novel inhibitors, and for re-engineering of kinase specificities.
The A2a adenosine receptor is a member of the G-protein coupled receptor family, and its activation stimulates cyclic AMP production. To determine the residues which are involved in ligand binding, several residues in transmembrane domains 5–7 were individually replaced with alanine and other amino acids. The binding properties of the resultant mutant receptors were determined in transfected COS-7 cells. To study the expression levels in COS-7 cells, mutant receptors were tagged at their amino terminus with a hemagglutinin epitope, which allowed their immunological detection in the plasma membrane by the monoclonal antibody 12CA5. The functional properties of mutant receptors were determined by measuring stimulation of adenylate cyclase. Specific binding of [3H]CGS 21680 (15 nm) and [3H]XAC (4 nm), an A2a agonist and antagonist, respectively, was absent in the following Ala mutants: F182A, H250A, N253A, I274A, H278A, and S281A, although they were well expressed in the plasma membrane. The hydroxy group of Ser-277 is required for high affinity binding of agonists, but not antagonists. An N181S mutant lost affinity for adenosine agonists substituted at N6 or C-2, but not at C-5′. The mutant receptors I274A, S277A, and H278A showed full stimulation of adenylate cyclase at high concentrations of CGS 21680. The functional agonist potencies at mutant receptors that lacked radioligand binding were >30-fold less than those at the wild type receptor. His-250 appears to be a required component of a hydrophobic pocket, and H-bonding to this residue is not essential. On the other hand, replacement of His-278 with other aromatic residues was not tolerated in ligand binding. Thus, some of the residues targeted in this study may be involved in the direct interaction with ligands in the human A2a adenosine receptor. A molecular model based on the structure of rhodopsin, in which the 5’-NH in NECA is hydrogen bonded to Ser-277 and His-278, was developed in order to visualize the environment of the ligand binding site.
Allosteric HIV-1 integrase (IN) inhibitors (ALLINIs) are an important new class of anti-HIV-1 agents. ALLINIs bind at the IN catalytic core domain (CCD) dimer interface occupying the principal binding pocket of its cellular cofactor LEDGF/p75. Consequently, ALLINIs inhibit HIV-1 IN interaction with LEDGF/p75 as well as promote aberrant IN multimerization. Selection of viral strains emerging under the inhibitor pressure has revealed mutations at the IN dimer interface near the inhibitor binding site.
We have investigated the effects of one of the most prevalent substitutions, H171T IN, selected under increasing pressure of ALLINI BI-D. Virus containing the H171T IN substitution exhibited an ~68-fold resistance to BI-D treatment in infected cells. These results correlated with ~84-fold reduced affinity for BI-D binding to recombinant H171T IN CCD protein compared to its wild type (WT) counterpart. However, the H171T IN substitution only modestly affected IN-LEDGF/p75 binding and allowed HIV-1 containing this substitution to replicate at near WT levels. The x-ray crystal structures of BI-D binding to WT and H171T IN CCD dimers coupled with binding free energy calculations revealed the importance of the Nδ- protonated imidazole group of His171 for hydrogen bonding to the BI-D tert-butoxy ether oxygen and establishing electrostatic interactions with the inhibitor carboxylic acid, whereas these interactions were compromised upon substitution to Thr171.
Our findings reveal a distinct mechanism of resistance for the H171T IN mutation to ALLINI BI-D and indicate a previously undescribed role of the His171 side chain for binding the inhibitor.
Electronic supplementary material
The online version of this article (doi:10.1186/s12977-014-0100-1) contains supplementary material, which is available to authorized users.
HIV-1 integrase; Allosteric inhibitors; Aberrant multimerization; Drug resistance
Human DNA topoisomerase I (hTop1p) catalyzes the relaxation of supercoiled DNA and constitutes the cellular target of the antitumor drug camptothecin (CPT). The X-ray crystal structure of the enzyme covalently joined to DNA and bound to the CPT analog Topotecan suggests that there are two classes of mutations that can produce a CPT-resistant enzyme. The first class includes changes in residues that directly interact with the drug, whereas a second class alters interactions with the DNA and thereby destabilizes the drug binding site. The Thr729Ala, that is part of a hydrophobic pocket in the enzyme C-terminal domain, belongs to a third group of mutations that confer CPT resistance, but do not interact directly with the drug or the DNA. To understand the contribution of this residue in drug resistance, we have studied the effect on hTop1p catalysis and CPT sensitivity of four different substitutions in the Thr729 position (Thr729Ala, Thr729Glu, Thr729Lys and Thr729Pro). Tht729Glu and Thr729Lys mutants show severe CPT resistance and furthermore, Thr729Glu shows a remarkable defect in DNA binding. We postulate that the maintenance of the hydrophobic pocket integrity, where Thr729 is positioned, is crucial for drug sensitivity and DNA binding.
C-terminal, aromatic amino acids in the ϕX174 internal scaffolding protein B mediate conformational switches in the viral coat protein. These switches direct the coat protein through early assembly. In addition to the aromatic amino acids, two acidic residues, D111 and E113, form salt bridges with basic, coat protein side chains. Although salt bridge formation did not appear to be critical for assembly, the substitution of an aromatic amino acid for D111 produced a lethal phenotype. This side chain is uniquely oriented toward the center of the coat-scaffolding binding pocket, which is heavily dominated by aromatic ring-ring interactions. Thus, the D111Y substitution may restructure pocket contacts. Previously characterized B− mutants blocked assembly before procapsid formation. However, the D111Y mutant produced an assembled particle, which contained the structural and external scaffolding proteins but lacked protein B and DNA. A suppressor within the external scaffolding protein, which mediates the later stages of particle morphogenesis, restored viability. The unique formation of a postprocapsid particle and the novel suppressor may be indicative of a novel B protein function. However, genetic data suggest that the particle represents the delayed manifestation of an early assembly error. This seemingly late-acting defect was rescued by previously characterized suppressors of early, preprocapsid, B− assembly mutations, which act on the level of coat protein flexibility. Likewise, the newly isolated suppressor in the external scaffolding protein also exhibited a global suppressing phenotype. Thus, the off-pathway product isolated from infected cells may not accurately reflect the temporal nature of the initial defect.
Alfalfa mosaic virus genomic RNAs are infectious only when the viral coat protein binds to the RNA 3´ termini. The crystal structure of an alfalfa mosaic virus RNA-peptide complex reveals that conserved AUGC repeats and Pro-Thr-x-Arg-Ser-x-x-Tyr coat protein amino acids cofold upon interacting. Alternating AUGC residues have opposite orientation, and they base pair in different adjacent duplexes. Localized RNA backbone reversals stabilized by arginine-guanine interactions place the adenosines and guanines in reverse order in the duplex. The results suggest that a uniform, organized 3´ conformation, similar to that found on viral RNAs with transfer RNA-like ends, may be essential for replication.
The cyclic 3', 5' adenosine monophosphate (cAMP) binding pocket of the cAMP receptor protein (CRP) of Escherichia coli was mutagenized to substitute cysteine or glycine for serine 83; cysteine, glycine, isoleucine, or serine for threonine 127; and threonine or alanine for serine 128. Cells that expressed the binding pocket residue-substituted forms of CRP were characterized by measurements of beta-galactosidase activity. Purified wild-type and mutant CRP preparations were characterized by measurement of cAMP binding activity and by their capacity to support lacP activation in vitro. CRP structure was assessed by measurement of sensitivity to protease and DTNB-mediated subunit crosslinking. The results of this study show that cAMP interactions with serine 83, threonine 127 and serine 128 contribute to CRP activation and have little effect on cAMP binding. Amino acid substitutions that introduce hydrophobic amino acid side chain constituents at either position 127 or 128 decrease CRP discrimination of cAMP and cGMP. Finally, cAMP-induced CRP structural change(s) that occur in or near the CRP hinge region result from cAMP interaction with threonine 127; substitution of threonine 127 by cysteine, glycine, isoleucine, or serine produced forms of CRP that contained, independently of cAMP binding, structural changes similar to those of the wild-type CRP:cAMP complex.
ALG-2 (a gene product of PDCD6) belongs to the penta-EF-hand (PEF) protein family and Ca2+-dependently interacts with various intracellular proteins including mammalian Alix, an adaptor protein in the ESCRT system. Our previous X-ray crystal structural analyses revealed that binding of Ca2+ to EF3 enables the side chain of R125 to move enough to make a primary hydrophobic pocket (Pocket 1) accessible to a short fragment of Alix. The side chain of F122, facing a secondary hydrophobic pocket (Pocket 2), interacts with the Alix peptide. An alternatively spliced shorter isoform, designated ALG-2ΔGF122, lacks Gly121Phe122 and does not bind Alix, but the structural basis of the incompetence has remained to be elucidated.
We solved the X-ray crystal structure of the PEF domain of ALG-2ΔGF122 in the Ca2+-bound form and compared it with that of ALG-2. Deletion of the two residues shortened α-helix 5 (α5) and changed the configuration of the R125 side chain so that it partially blocked Pocket 1. A wall created by the main chain of 121-GFG-123 and facing the two pockets was destroyed. Surprisingly, however, substitution of F122 with Ala or Gly, but not with Trp, increased the Alix-binding capacity in binding assays. The F122 substitutions exhibited different effects on binding of ALG-2 to other known interacting proteins, including TSG101 (Tumor susceptibility gene 101) and annexin A11. The X-ray crystal structure of the F122A mutant revealed that removal of the bulky F122 side chain not only created an additional open space in Pocket 2 but also abolished inter-helix interactions with W95 and V98 (present in α4) and that α5 inclined away from α4 to expand Pocket 2, suggesting acquirement of more appropriate positioning of the interacting residues to accept Alix.
We found that the inability of the two-residue shorter ALG-2 isoform to bind Alix is not due to the absence of bulky side chain of F122 but due to deformation of a main-chain wall facing pockets 1 and 2. Moreover, a residue at the position of F122 contributes to target specificity and a smaller side chain is preferable for Alix binding but not favored to bind annexin A11.
TRAP is an 11 subunit RNA binding protein that regulates expression of genes involved in tryptophan biosynthesis and transport in Bacillus subtilis. TRAP is activated to bind RNA by binding up to 11 molecules of l-tryptophan in pockets formed by adjacent subunits. The precise mechanism by which tryptophan binding activates TRAP is not known. Thr30 is in the tryptophan binding pocket. A TRAP mutant in which Thr30 is substituted with Val (T30V) does not bind tryptophan but binds RNA constitutively, suggesting that Thr30 plays a key role in the activation mechanism. We have examined the effects of other substitutions of Thr30. TRAP proteins with small β-branched aliphatic side chains at residue 30 bind RNA constitutively, whereas those with a small polar side chain show tryptophan-dependent RNA binding. Several mutant proteins exhibited constitutive RNA binding that was enhanced by tryptophan. Although the tryptophan and RNA binding sites on TRAP are distinct and are separated by ∼7.5 Å, several substitutions of residues that interact with the bound RNA restored tryptophan binding to T30V TRAP. These observations support the hypothesis that conformational changes in TRAP relay information between the tryptophan and RNA binding sites of the protein.
The seven antigenically distinct serotypes of Clostridium botulinum neurotoxins, the causative agents of botulism, block the neurotransmitter release by specifically cleaving one of the three SNARE proteins and induce flaccid paralysis. The Centers for Disease Control and Prevention (CDC) has declared them as Category A biowarfare agents. The most potent among them, botulinum neurotoxin type A (BoNT/A), cleaves its substrate synaptosome-associated protein of 25 kDa (SNAP-25). An efficient drug for botulism can be developed only with the knowledge of interactions between the substrate and enzyme at the active site. Here, we report the crystal structures of the catalytic domain of BoNT/A with its uncleavable SNAP-25 peptide 197QRATKM202 and its variant 197RRATKM202 to 1.5 Å and 1.6 Å, respectively. This is the first time the structure of an uncleavable substrate bound to an active botulinum neurotoxin is reported and it has helped in unequivocally defining S1 to S5′ sites. These substrate peptides make interactions with the enzyme predominantly by the residues from 160, 200, 250 and 370 loops. Most notably, the amino nitrogen and carbonyl oxygen of P1 residue (Gln197) chelate the zinc ion and replace the nucleophilic water. The P1′-Arg198, occupies the S1′ site formed by Arg363, Thr220, Asp370, Thr215, Ile161, Phe163 and Phe194. The S2′ subsite is formed by Arg363, Asn368 and Asp370, while S3′ subsite is formed by Tyr251, Leu256, Val258, Tyr366, Phe369 and Asn388. P4′-Lys201 makes hydrogen bond with Gln162. P5′-Met202 binds in the hydrophobic pocket formed by the residues from the 250 and 200 loop. Knowledge of interactions between the enzyme and substrate peptide from these complex structures should form the basis for design of potent inhibitors for this neurotoxin.
Botulinum neurotoxins are the most poisonous substance to humans. The ease with which the bacteria can be grown, its potency and persistence have made it a potential bioterrorism agent, and accordingly, botulinum neurotoxin has been declared as Category A agent by the Centers of Disease Control and Prevention. Since it is both a potential bioweapon and a bioterrorism agent, it is imperative to develop counter measures and therapeutics for these neurotoxins, as none are available so far except experimental vaccines and an FDA-approved equine antitoxin. Our work presented here is an important milestone towards achieving this goal. The best antidote can be developed by blocking the active site of any enzyme. The crystal structures of substrate peptide–enzyme complex presented here map the interactions between the two and provide critical information for designing effective drugs against this toxin.
Juvenile hormone (JH) is a key insect developmental hormone that is found at low nanomolar levels in larval insects. The methyl ester of JH is hydrolyzed in many insects by an esterase that shows high specificity for JH. We have previously determined a crystal structure of the JH esterase (JHE) of the tobacco hornworm Manduca sexta (MsJHE) [Wogulis, M., Wheelock, C. E., Kamita, S. G., Hinton, A. C., Whetstone, P. A., Hammock, B. D., and Wilson, D. K. (2006) Biochemistry 45, 4045-4057]. Our molecular modeling indicates that JH fits very tightly within the substrate binding pocket of MsJHE. This tight fit places two non-catalytic amino acid residues, Phe-259 and Thr-314, within the appropriate distance and geometry to potentially interact with the α,β-unsaturated ester and epoxide, respectively, of JH. These residues are highly conserved in numerous biologically active JHEs. Kinetic analyses of mutants of Phe-259 or Thr-314 indicate that these residues contribute to the low KM that MsJHE shows for JH. This low KM, however, comes at the cost of reduced substrate turnover. Neither nucleophilic attack of the resonance stabilized ester by the catalytic serine nor the availability of a water molecule for attack of the acyl-enzyme intermediate appear to be a rate-determining step in the hydrolysis of JH by MsJHE. We hypothesize that the release of the JH acid metabolite from the substrate binding pocket limits the catalytic cycle. Our findings also demonstrate that chemical bond strength does not necessarily correlate with how reactive the bond will be to metabolism.
Agnoprotein is a small multifunctional regulatory protein required for sustaining the productive replication of JC virus (JCV). It is a mostly cytoplasmic protein localizing in the perinuclear area and forms highly stable dimers/oligomers through a Leu/Ile/Phe-rich domain. There have been no three-dimensional structural data available for agnoprotein due to difficulties associated with the dynamic conversion from monomers to oligomers. Here, we report the first nuclear magnetic resonance (NMR) structure of a synthetic agnoprotein peptide spanning amino acids Thr17 to Glu55 where Lys23 to Phe39 encompassing the Leu/Ile/Phe-rich domain forms an amphipathic α-helix. On the basis of these structural data, a number of Ala substitution mutations were made to investigate the role of the α-helix in the structure and function of agnoprotein. Single L29A and L36A mutations exhibited a significant negative effect on both protein stability and viral replication, whereas the L32A mutation did not. In addition, the L29A mutant displayed a highly nuclear localization pattern, in contrast to the pattern for the wild type (WT). Interestingly, a triple mutant, the L29A+L32A+L36A mutant, yielded no detectable agnoprotein expression, and the replication of this JCV mutant was significantly reduced, suggesting that Leu29 and Leu36 are located at the dimer interface, contributing to the structure and stability of agnoprotein. Two other single mutations, L33A and E34A, did not perturb agnoprotein stability as drastically as that observed with the L29A and L36A mutations, but they negatively affected viral replication, suggesting that the role of these residues is functional rather than structural. Thus, the agnoprotein dimerization domain can be targeted for the development of novel drugs active against JCV infection.
IMPORTANCE Agnoprotein is a small regulatory protein of JC virus (JCV) and is required for the successful completion of the viral replication cycle. It forms highly stable dimers and oligomers through its hydrophobic (Leu/Ile/Phe-rich) domain, which has been shown to play essential roles in the stability and function of the protein. In this work, the Leu/Ile/Phe-rich domain has been further characterized by NMR studies using an agnoprotein peptide spanning amino acids T17 to Q54. Those studies revealed that the dimerization domain of the protein forms an amphipathic α-helix. Subsequent NMR structure-based mutational analysis of the region highlighted the critical importance of certain amino acids within the α-helix for the stability and function of agnoprotein. In conclusion, this study provides a solid foundation for developing effective therapeutic approaches against the dimerization domain of the protein to inhibit its critical roles in JCV infection.
17β-Hydroxysteroid dehydrogenase type 1 (17β-HSD1) catalyzes the conversion of estrone to the potent estrogen estradiol. 17β-HSD1 is highly expressed in breast and ovary tissues and represents a prognostic marker for the tumor progression and survival of patients with breast cancer and other estrogen-dependent tumors. Therefore, the enzyme is considered a promising drug target against estrogen-dependent cancers. For the development of novel inhibitors, an improved understanding of the structure-function relationships is essential. In the present study, we examined the role of a cysteine residue, Cys10, in the Rossmann-fold NADPH binding region, for 17β-HSD1 function and tested the sensitivity towards sulfhydryl modifying chemicals. 3D structure modeling revealed important interactions of Cys10 with residues involved in the stabilization of amino acids of the NADPH binding pocket. Analysis of enzyme activity revealed that 17β-HSD1 was irreversibly inhibited by the sulfhydryl modifying agents N-ethylmaleimide (NEM) and dithiocarbamates. Preincubation with increasing concentrations of NADPH protected 17β-HSD1 from inhibition by these chemicals. Cys10Ser mutant 17β-HSD1 was partially protected from inhibition by NEM and dithiocarbamates, emphasizing the importance of Cys10 in the cofactor binding region. Substitution of Cys10 with serine resulted in a decreased protein half-life, without significantly altering kinetic properties. Despite the fact that Cys10 on 17β-HSD1 seems to have limited potential as a target for new enzyme inhibitors, the present study provides new insight into the structure-function relationships of this enzyme.
Methylaspartate ammonia-lyase (MAL; EC 22.214.171.124) catalyzes the reversible addition of ammonia to mesaconate to yield l-threo-(2S,3S)-3-methylaspartate and l-erythro-(2S,3R)-3-methylaspartate as products. In the proposed minimal mechanism for MAL of Clostridium tetanomorphum, Lys-331 acts as the (S)-specific base catalyst and abstracts the 3S-proton from l-threo-3-methylaspartate, resulting in an enolate anion intermediate. This enolic intermediate is stabilized by coordination to the essential active site Mg2+ ion and hydrogen bonding to the Gln-329 residue. Collapse of this intermediate results in the release of ammonia and the formation of mesaconate. His-194 likely acts as the (R)-specific base catalyst and abstracts the 3R-proton from the l-erythro isomer of 3-methylaspartate, yielding the enolic intermediate. In the present study, we have investigated the importance of the residues Gln-73, Phe-170, Gln-172, Tyr-356, Thr-360, Cys-361 and Leu-384 for the catalytic activity of C. tetanomorphum MAL. These residues, which are part of the enzyme surface lining the substrate binding pocket, were subjected to site-directed mutagenesis and the mutant enzymes were characterized for their structural integrity, ability to catalyze the amination of mesaconate, and regio- and diastereoselectivity. Based on the observed properties of the mutant enzymes, combined with previous structural studies and protein engineering work, we propose a detailed catalytic mechanism for the MAL-catalyzed reaction, in which the side chains of Gln-73, Gln-172, Tyr-356, Thr-360, and Leu-384 provide favorable interactions with the substrate, which are important for substrate binding and activation. This detailed knowledge of the catalytic mechanism of MAL can serve as a guide for future protein engineering experiments.
•A detailed catalytic mechanism for methylaspartate ammonia-lyase is proposed•Gln-172, Thr-360 and Cys-361 bind the 1-carboxylate group of 3-methylaspartate•Gln-73 (via a water molecule) and Gln-172 bind the 2-amino group•Tyr-356 and Leu-384 provide stabilizing interactions with the 3-methyl group
Methylaspartate ammonia-lyase; Deamination; Enzyme mechanism; Amino acid; Mutagenesis
Several mutations at Glu81 located on the crystal contact of human acidic fibroblast growth factor were studied in an effort to improve crystal growth. Mutation to Ser and Thr resulted in crystallization of a rather bulky form of the wild type, whereas mutation to Val prohibited crystallization. These results suggest that crystal growth may be controlled by designing a new interface by protein engineering.
An attempt has been made to improve a crystal contact of human acidic fibroblast growth factor (haFGF; 140 amino acids) to control the crystal growth, because haFGF crystallizes only as a thin-plate form, yielding crystals suitable for X-ray but not neutron diffraction. X-ray crystal analysis of haFGF showed that the Glu81 side chain, located at a crystal contact between haFGF molecules, is in close proximity with an identical residue related by crystallographic symmetry, suggesting that charge repulsion may disrupt suitable crystal-packing interactions. To investigate whether the Glu residue affects the crystal-packing interactions, haFGF mutants in which Glu81 was replaced by Ala, Val, Leu, Ser and Thr were constructed. Although crystals of the Ala and Leu mutants were grown as a thin-plate form by the same precipitant (formate) as the wild type, crystals of the Ser and Thr mutants were grown with increased thickness, yielding a larger overall crystal volume. X-ray structural analysis of the Ser mutant determined at 1.35 Å resolution revealed that the hydroxy groups of Ser are linked by hydrogen bonds mediated by the formate used as a precipitant. This approach to engineering crystal contacts may contribute to the development of large protein crystals for neutron crystallography.
acidic fibroblast growth factor; crystal contacts; mutagenesis
RNA modification plays an important role in modulating host-pathogen interaction. Flavivirus NS5 protein encodes N-7 and 2′-O methyltransferase activities that are required for the formation of 5′ type I cap (m7GpppAm) of viral RNA genome. Here we reported, for the first time, that flavivirus NS5 has a novel internal RNA methylation activity. Recombinant NS5 proteins of West Nile virus and Dengue virus (serotype 4; DENV-4) specifically methylates polyA, but not polyG, polyC, or polyU, indicating that the methylation occurs at adenosine residue. RNAs with internal adenosines substituted with 2′-O-methyladenosines are not active substrates for internal methylation, whereas RNAs with adenosines substituted with N6-methyladenosines can be efficiently methylated, suggesting that the internal methylation occurs at the 2′-OH position of adenosine. Mass spectroscopic analysis further demonstrated that the internal methylation product is 2′-O-methyladenosine. Importantly, genomic RNA purified from DENV virion contains 2′-O-methyladenosine. The 2′-O methylation of internal adenosine does not require specific RNA sequence since recombinant methyltransferase of DENV-4 can efficiently methylate RNAs spanning different regions of viral genome, host ribosomal RNAs, and polyA. Structure-based mutagenesis results indicate that K61-D146-K181-E217 tetrad of DENV-4 methyltransferase forms the active site of internal methylation activity; in addition, distinct residues within the methyl donor (S-adenosyl-L-methionine) pocket, GTP pocket, and RNA-binding site are critical for the internal methylation activity. Functional analysis using flavivirus replicon and genome-length RNAs showed that internal methylation attenuated viral RNA translation and replication. Polymerase assay revealed that internal 2′-O-methyladenosine reduces the efficiency of RNA elongation. Collectively, our results demonstrate that flavivirus NS5 performs 2′-O methylation of internal adenosine of viral RNA in vivo and host ribosomal RNAs in vitro.
We report that flavivirus NS5 has a novel internal RNA methylation activity. Recombinant proteins of NS5 and its N-terminal methyltransferase domain of West Nile virus and Dengue virus (DENV) specifically methylates polyA, but not polyG, polyC, or polyU. RNAs with internal adenosines substituted with 2′-O-methyladenosines are not active substrates for internal methylation, suggesting that the internal methylation occurs at the 2′-OH position of adenosine. Mass spectroscopic analysis confirmed that the internal methylation product is 2′-O-methyladenosine. Furthermore, the 2′-O-methyladenosine could also be detected in DENV genomic RNA. The 2′-O methylation of internal adenosine does not require specific RNA sequence context because the DENV methyltransferase can methylate RNAs spanning different regions of viral genome and host ribosomal RNAs at equal efficiencies. Mutagenesis analysis showed that K61-D146-K181-E217 motif of the DENV methyltransferase forms the active site of internal methylation activity; in addition, distinct residues on the surface of the enzyme are critical for the internal methylation activity. Functional analysis showed that internal methylation attenuated viral RNA translation and replication. Overall, our results demonstrate that flavivirus NS5 performs 2′-O methylation of internal adenosine of viral RNA in vivo and host ribosomal RNA in vitro. Such 2′-O-methyladenosine modification may modulate virus-host interaction.
CheY serves as a structural prototype for the response regulator proteins of two-component regulatory systems. Functional roles have previously been defined for four of the five highly conserved residues that form the response regulator active site, the exception being the hydroxy amino acid which corresponds to Thr87 in CheY. To investigate the contribution of Thr87 to signaling, we characterized, genetically and biochemically, several cheY mutants with amino acid substitutions at this position. The hydroxyl group appears to be necessary for effective chemotaxis, as a Thr→Ser substitution was the only one of six tested which retained a Che+ swarm phenotype. Although nonchemotactic, cheY mutants with amino acid substitutions T87A and T87C could generate clockwise flagellar rotation either in the absence of CheZ, a protein that stimulates dephosphorylation of CheY, or when paired with a second site-activating mutation, Asp13→Lys, demonstrating that a hydroxy amino acid at position 87 is not essential for activation of the flagellar switch. All purified mutant proteins examined phosphorylated efficiently from the CheA kinase in vitro but were impaired in autodephosphorylation. Thus, the mutant CheY proteins are phosphorylated to a greater degree than wild-type CheY yet support less clockwise flagellar rotation. The data imply that Thr87 is important for generating and/or stabilizing the phosphorylation-induced conformational change in CheY. Furthermore, the various position 87 substitutions differentially affected several properties of the mutant proteins. The chemotaxis and autodephosphorylation defects were tightly linked, suggesting common structural elements, whereas the effects on self-catalyzed and CheZ-mediated dephosphorylation of CheY were uncorrelated, suggesting different structural requirements for the two dephosphorylation reactions.
Bacteria such as Escherichia coli and Salmonella typhimurium can utilize acetate as the sole source of carbon and energy. Acetate kinase (AckA) and phosphotransacetylase (Pta), key enzymes of acetate utilization pathway, regulate flux of metabolites in glycolysis, gluconeogenesis, TCA cycle, glyoxylate bypass and fatty acid metabolism.
Here we report kinetic characterization of S. typhimurium AckA (StAckA) and structures of its unliganded (Form-I, 2.70 Å resolution) and citrate-bound (Form-II, 1.90 Å resolution) forms. The enzyme showed broad substrate specificity with kcat/Km in the order of acetate > propionate > formate. Further, the Km for acetyl-phosphate was significantly lower than for acetate and the enzyme could catalyze the reverse reaction (i.e. ATP synthesis) more efficiently. ATP and Mg2+ could be substituted by other nucleoside 5′-triphosphates (GTP, UTP and CTP) and divalent cations (Mn2+ and Co2+), respectively. Form-I StAckA represents the first structural report of an unliganded AckA. StAckA protomer consists of two domains with characteristic βββαβαβα topology of ASKHA superfamily of proteins. These domains adopt an intermediate conformation compared to that of open and closed forms of ligand-bound Methanosarcina thermophila AckA (MtAckA). Spectroscopic and structural analyses of StAckA further suggested occurrence of inter-domain motion upon ligand-binding. Unexpectedly, Form-II StAckA structure showed a drastic change in the conformation of residues 230–300 compared to that of Form-I. Further investigation revealed electron density corresponding to a citrate molecule in a pocket located at the dimeric interface of Form-II StAckA. Interestingly, a similar dimeric interface pocket lined with largely conserved residues could be identified in Form-I StAckA as well as in other enzymes homologous to AckA suggesting that ligand binding at this pocket may influence the function of these enzymes.
The biochemical and structural characterization of StAckA reported here provides insights into the biochemical specificity, overall fold, thermal stability, molecular basis of ligand binding and inter-domain motion in AckA family of enzymes. Dramatic conformational differences observed between unliganded and citrate-bound forms of StAckA led to identification of a putative ligand-binding pocket at the dimeric interface of StAckA with implications for enzymatic function.
Acetate metabolism; AckA-Pta pathway; Acetate and Sugar Kinases/Heat shock cognate (Hsc) 70/Actin (ASKHA) superfamily; Conformational rearrangement; Enzyme regulation