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.
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.
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.
ALG-2, a 22-kDa penta-EF-hand protein, is involved in cell death, signal transduction, membrane trafficking, etc., by interacting with various proteins in mammalian cells in a Ca2+-dependent manner. Most known ALG-2-interacting proteins contain proline-rich regions in which either PPYPXnYP (type 1 motif) or PXPGF (type 2 motif) is commonly found. Previous X-ray crystal structural analysis of the complex between ALG-2 and an ALIX peptide revealed that the peptide binds to the two hydrophobic pockets. In the present study, we resolved the crystal structure of the complex between ALG-2 and a peptide of Sec31A (outer shell component of coat complex II, COPII; containing the type 2 motif) and found that the peptide binds to the third hydrophobic pocket (Pocket 3). While amino acid substitution of Phe85, a Pocket 3 residue, with Ala abrogated the interaction with Sec31A, it did not affect the interaction with ALIX. On the other hand, amino acid substitution of Tyr180, a Pocket 1 residue, with Ala caused loss of binding to ALIX, but maintained binding to Sec31A. We conclude that ALG-2 recognizes two types of motifs at different hydrophobic surfaces. Furthermore, based on the results of serial mutational analysis of the ALG-2-binding sites in Sec31A, the type 2 motif was newly defined.
adaptor protein; calcium-binding protein; COPII; crystal structure; EF-hand; motif; protein-protein interaction
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.
Developing antiviral therapies for influenza A virus (IAV) infection is an ongoing process because of the rapid rate of antigenic mutation and the emergence of drug-resistant viruses. The ideal strategy is to develop drugs that target well-conserved, functionally restricted, and unique surface structures without affecting host cell function. We recently identified the antiviral compound, RK424, by screening a library of 50,000 compounds using cell-based infection assays. RK424 showed potent antiviral activity against many different subtypes of IAV in vitro and partially protected mice from a lethal dose of A/WSN/1933 (H1N1) virus in vivo. Here, we show that RK424 inhibits viral ribonucleoprotein complex (vRNP) activity, causing the viral nucleoprotein (NP) to accumulate in the cell nucleus. In silico docking analysis revealed that RK424 bound to a small pocket in the viral NP. This pocket was surrounded by three functionally important domains: the RNA binding groove, the NP dimer interface, and nuclear export signal (NES) 3, indicating that it may be involved in the RNA binding, oligomerization, and nuclear export functions of NP. The accuracy of this binding model was confirmed in a NP-RK424 binding assay incorporating photo-cross-linked RK424 affinity beads and in a plaque assay evaluating the structure-activity relationship of RK424. Surface plasmon resonance (SPR) and pull-down assays showed that RK424 inhibited both the NP-RNA and NP-NP interactions, whereas size exclusion chromatography showed that RK424 disrupted viral RNA-induced NP oligomerization. In addition, in vitro nuclear export assays confirmed that RK424 inhibited nuclear export of NP. The amino acid residues comprising the NP pocket play a crucial role in viral replication and are highly conserved in more than 7,000 NP sequences from avian, human, and swine influenza viruses. Furthermore, we found that the NP pocket has a surface structure different from that of the pocket in host molecules. Taken together, these results describe a promising new approach to developing influenza virus drugs that target a novel pocket structure within NP.
Influenza A virus nucleoprotein (NP) is a highly conserved multifunctional protein that plays an essential role in infection by all subtypes of influenza A virus, making it an attractive target for new antiviral drugs. NP regulates viral polymerase activity and transport of the viral genome into/from the host cell nucleus by forming the viral ribonucleoprotein complex (vRNP). Because NP regulates replication and transcription of the viral genome in addition to its role in nuclear export (all of which are essential for the production of viral progeny), it is a promising drug target. Here, we used the antiviral compound RK424 to identify a novel pocket structure within NP. This structure encompassed three different functional domains that are involved in the above-mentioned replication steps. RK424 inhibits viral genome replication/transcription and nuclear export of NP by destabilizing the NP oligomer and inhibiting the binding of chromosome region maintenance 1 (CRM1) to NP via nuclear export signal (NES) 3, which is located in close proximity to the NP pocket. Taken together, these findings suggest that this small NP pocket is a novel antiviral target.
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.
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
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.
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.
Rabbit hemorrhagic disease virus (RHDV) is a member of the Caliciviridae family (Lagovirus genus). RHDV is highly contagious and attaches to epithelial cells in the digestive or respiratory tract, leading to massive lesions with high mortality rates. A new variant of RHDV (termed RHDVb) recently has emerged, and previously vaccinated rabbits appear to have little protection against this new strain. Similar to human norovirus (Caliciviridae, Norovirus genus), RHDV binds histo-blood group antigens (HBGAs), and this is thought to be important for infection. Here, we report the HBGA binding site on the RHDVb capsid-protruding domain (P domain) using X-ray crystallography. The HBGA binding pocket was located in a negatively charged patch on the side of the P domain and at a dimeric interface. Residues from both monomers contributed to the HBGA binding and involved a network of direct hydrogen bonds and water-mediated interactions. An amino acid sequence alignment of different RHDV strains indicated that the residues directly interacting with the ABH-fucose of the HBGAs (Asp472, Asn474, and Ser479) were highly conserved. This result suggested that different RHDV strains also could bind HBGAs at the equivalent pocket. Moreover, several HBGA binding characteristics between RHDVb and human genogroup II norovirus were similar, which indicated a possible convergent evolution of HBGA binding interactions. Further structural studies with other RHDV strains are needed in order to better understand the HBGA binding mechanisms among the diverse RHDV strains.
IMPORTANCE We identified, for the first time, the HBGA binding site on an RHDVb P domain using X-ray crystallography. Our results showed that RHDVb and human genogroup II noroviruses had similar HBGA binding interactions. Recently, it was discovered that synthetic HBGAs or HBGA-expressing enteric bacteria could enhance human genogroup II norovirus infection in B cells. Considering that RHDVb and genogroup II norovirus similarly interacted with HBGAs, it may be possible that a comparable cell culture system also could work with RHDVb. Taken together, these new findings will extend our understanding of calicivirus HBGA interactions and may help to elucidate the specific roles of HBGAs in calicivirus infections.
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
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.
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.
Schizophrenia (SZ), a chronic mental and heritable disorder characterized by neurophysiological impairment and neuropsychological abnormalities, is strongly associated with D-amino acid oxidase activator (DAOA, G72). Research studies emphasized that overexpression of DAOA may be responsible for improper functioning of neurotransmitters, resulting in neurological disorders like SZ. In the present study, a hybrid approach of comparative modeling and molecular docking followed by inhibitor identification and structure modeling was employed. Screening was performed by two-dimensional similarity search against selected inhibitor, keeping in view the physiochemical properties of the inhibitor. Here, we report an inhibitor compound which showed maximum binding affinity against four selected isoforms of DAOA. Docking studies revealed that Glu-53, Thr-54, Lys-58, Val-85, Ser-86, Tyr-87, Leu-88, Glu-90, Leu-95, Val-98, Ser-100, Glu-112, Tyr-116, Lys-120, Asp-121, and Arg-122 are critical residues for receptor–ligand interaction. The C-terminal of selected isoforms is conserved, and binding was observed on the conserved region of isoforms. We propose that selected inhibitor might be more potent on the basis of binding energy values. Further analysis of this inhibitor through site-directed mutagenesis could be helpful for exploring the details of ligand-binding pockets. Overall, the findings of this study may be helpful in designing novel therapeutic targets to cure SZ.
schizophrenia; bioinformatics; modeling; docking; DAOA; G72; DAO; computer-aided drug designing; phylogenetic analysis; D-amino acid oxidase activator
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.
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.
Methylaspartate ammonia-lyase (MAL; EC 188.8.131.52) 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
The coat protein of bacteriophage MS2 functions as a symmetric dimer to bind an asymmetric RNA hairpin. This implies the existence of two equivalent RNA binding sites related to one another by a 2-fold symmetry axis. In this view the symmetric binding site defined by mutations conferring the repressor-defective phenotype is a composite picture of these two asymmetric sites. In order to determine whether the RNA ligand interacts with amino acid residues on both subunits of the dimer and in the hope of constructing a functional map of the RNA binding site, we performed heterodimer complementation experiments. Taking advantage of the physical proximity of their N- and C-termini, the two subunits of the dimer were genetically fused, producing a duplicated coat protein which folds normally and allows the construction of the functional equivalent of obligatory heterodimers containing all possible pairwise combinations of the repressor-defective mutations. The restoration of repressor function in certain heterodimers shows that a single RNA molecule interacts with both subunits of the dimer and allows the construction of a functional map of the binding site.
The mechanism(s) by which fatty acids are sequestered and transported in muscle have not been fully elucidated. A potential key player in this process is the protein myoglobin (Mb). Indeed, there is a catalogue of empirical evidence supporting direct interaction of globins with fatty acid metabolites; however, the binding pocket and regulation of the interaction remains to be established. In this study, we employed a computational strategy to elucidate the structural determinants of fatty acids (palmitic & oleic acid) binding to Mb. Sequence analysis and docking simulations with a horse (Equus caballus) structural Mb reference reveals a fatty acid-binding site in the hydrophobic cleft near the heme region in Mb. Both palmitic acid and oleic acid attain a “U” shaped structure similar to their conformation in pockets of other fatty acid-binding proteins. Specifically, we found that the carboxyl head group of palmitic acid coordinates with the amino group of Lys45, whereas the carboxyl group of oleic acid coordinates with both the amino groups of Lys45 and Lys63. The alkyl tails of both fatty acids are supported by surrounding hydrophobic residues Leu29, Leu32, Phe33, Phe43, Phe46, Val67, Val68 and Ile107. In the saturated palmitic acid, the hydrophobic tail moves freely and occasionally penetrates deeper inside the hydrophobic cleft, making additional contacts with Val28, Leu69, Leu72 and Ile111. Our simulations reveal a dynamic and stable binding pocket in which the oxygen molecule and heme group in Mb are required for additional hydrophobic interactions. Taken together, these findings support a mechanism in which Mb acts as a muscle transporter for fatty acid when it is in the oxygenated state and releases fatty acid when Mb converts to deoxygenated state.
Porphyromonas gingivalis is a periodontal pathogen whose primary niche is the anaerobic environment of subgingival dental plaque, but initial colonization of the oral cavity is likely to occur on supragingival surfaces that already support robust biofilm communities. Our studies have shown that P. gingivalis adheres to Streptococcus gordonii through interaction of the minor fimbrial antigen Mfa1 with a specific region of the streptococcal SspB polypeptide (residues 1167 to 1193) designated BAR. We show that a synthetic peptide comprising the BAR sequence potently inhibits P. gingivalis adherence to S. gordonii (50% inhibitory concentration = 1.3 μM) and prevents the development of P. gingivalis biofilms. However, a retroinverso peptide that possessed the same side chain topology as that of BAR was inactive, suggesting that interactions of Mfa1 with the peptide backbone of BAR are important for binding. A conformationally constrained analog of BAR inhibited P. gingivalis adherence and biofilm formation but at a lower specific activity than that of BAR. Therefore, to further define the structural features of the Mfa1-BAR interaction, we functionally screened combinatorial libraries of BAR in which active site residues (Asn1182, Thr1184, and Val1185) were replaced with each of the 19 common amino acids. Peptides containing positively charged amino acids at position 1182 or hydrophobic residues at position 1185 bound P. gingivalis more efficiently than did control peptides containing Asn and Val at these positions, suggesting that electrostatic and hydrophobic interactions may contribute to Mfa1-SspB binding. In contrast, replacement of Pro or Gly at these positions was detrimental to adherence, suggesting that perturbation of the BAR secondary structure influences activity. The net effect of substitutions for Thr1184 was less pronounced either positively or negatively than that at the other sites. These results define physicochemical characteristics of the interacting interface of Mfa1 and SspB and suggest that peptides or peptidomimetics with greater specific inhibitory activity than that of BAR can be developed. These compounds may represent potential therapeutics that target some of the first molecular interactions that allow P. gingivalis to colonize the oral cavity.
The shape of the protein surface dictates what interactions are possible with other macromolecules, but defining discrete pockets or possible interaction sites remains difficult. First, there is the problem of defining the extent of the pocket. Second, one has to characterize the shape of each pocket. Third, one needs to make quantitative comparisons between pockets on different proteins. An elegant solution to these problems is to sort all surface and solvent points by Travel Depth, and then collect a hierarchical tree of pockets. The connectivity of the tree is determined via the deepest saddle points between each pair of neighboring pockets. The resulting pocket surfaces tessellate the entire protein surface, producing a complete inventory of pockets. This method of identifying pockets also allows one to easily compute important shape metrics, including the problematic pocket volume, surface area, and mouth size. Pockets are also annotated with their lining residue lists, polarity, and other residue based properties. Using this tree and the various shape metrics pockets can be merged, grouped, or filtered for further analysis. Since this method includes the entire surface it guarantees that any pocket of interest will be found among the output pockets, unlike all previous methods of pocket identification. The resulting hierarchy of pockets is easy to visualize and aids users in higher level analysis. Comparison of pockets is done using the shape metrics, avoiding the complex shape alignment problem. Example applications show that the method facilitates pocket comparison along mutational or time-dependent series. Pockets from families of proteins can be examined using multiple pocket tree alignments to see how ligand binding sites or other pockets have changed with evolution. Our method is called CLIPPERS, for Complete Liberal Inventory of Protein Pockets Elucidating and Reporting on Shape.
Pockets; depth; ligand binding; molecular surface; protein shape
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.
The Toll-like receptor 4/MD-2 receptor complex recognizes endotoxin, a Gram-negative bacterial cell envelope component. Recognition of the most potent hexaacylated form of endotoxin is mediated by the sixth acyl chain that protrudes from the MD-2 hydrophobic pocket and bridges TLR4/MD-2 to the neighboring TLR4 ectodomain, driving receptor dimerization via hydrophobic interactions. In hypoacylated endotoxins all acyl chains could be accommodated within the binding pocket of the human hMD-2. Nevertheless, tetra- and pentaacylated endotoxins activate the TLR4/MD-2 receptor of several species. We observed that amino acid residues 82 and 122, located at the entrance to the endotoxin binding site of MD-2, have major influence on the species-specific endotoxin recognition. We show that substitution of hMD-2 residue V82 with an amino acid residue with a bulkier hydrophobic side chain enables activation of TLR4/MD-2 by pentaacylated and tetraacylated endotoxins. Interaction of the lipid A phosphate group with the amino acid residue 122 of MD-2 facilitates the appropriate positioning of the hypoacylated endotoxin. Moreover, mouse TLR4 contributes to the agonistic effect of pentaacylated msbB endotoxin. We propose a molecular model that explains how the molecular differences between the murine or equine MD-2, which both have sufficiently large hydrophobic pockets to accommodate all five or four acyl chains, influence the positioning of endotoxin so that one of the acyl chains remains outside the pocket and enables hydrophobic interactions with TLR4, leading to receptor activation.