NMR assignment of backbone and side chain resonances for a dockerin-containing C-terminal fragment of the putative μ-toxin from Clostridium perfringens

The μ-toxin of Clostridium perfringens, termed CpGH84A, is modular hydrolytic enzyme that contributes to the pathogenicity of this organism. Backbone and side chain ¹H, ¹³C, and ¹⁵N resonance assignments have been determined for the C-terminal 15.5 kDa FIVAR-Doc modular pair of CpGH84A.     

Structure and functional investigation of ligand binding by a family 35 carbohydrate binding module (CtCBM35) of β-mannanase of family 26 glycoside hydrolase from Clostridium thermocellum

Functional attributes of recombinant CtCBM35 (family 35 carbohydrate binding module) of β-mannanase of family 26 Glycoside Hydrolase from Clostridium thermocellum were deduced by biochemical and in silico approaches. Ligand-binding analysis of expressed CtCBM35 analyzed by affinity-gel electrophoresis and fluorescence spectroscopy exhibited association constants K _a ～ 1.2·10⁵ and 3.0·10⁵ M^?1 with locust bean galactomannan and mannotriose, respectively. However, CtCBM35 showed low ligand-binding affinity with insoluble ivory nut mannan with K _a of 5.0·10^?5 M^?1. Unfolding transition analysis by fluorescence spectroscopy explained the conformational changes of CtCBM35 in the presence of guanidine hydrochloride (5 M) and urea (6.25 M). This explained that CtCBM35 has good conformational stability and requires higher free energy of denaturation to invoke unfolding. The three-dimensional (3-D) model of CtCBM35 from C. thermocellum generated by Modeller9v8 displayed predominance of β-sheets arranged as β-jelly-roll fold. The secondary structure of CtCBM35 by PredictProtein showed the presence of two α-helices (3%), 12 β-sheets (45%), and 15 random coils (52%). Secondary structural element analysis of cloned, expressed, and purified recombinant CtCBM35 by circular dichroism also corroborated the in silico predicted secondary structure. Multiple sequence alignment of CtCBM35 showed conserved residues (Tyr123, Gly124, and Phe125), which are commonly observed in mannan specific CBMs. Docking analysis of CtCBM35 with manno-oligosaccharide displayed the involvement of Tyr26, Gln29, Asn43, Trp66, Tyr68, Leu69, Arg76, and Leu127 residues, making polar contact with the ligand molecules. Ligand docking analysis of CtCBM35 exhibiting higher binding affinity with mannotriose and galactomannan (Man-Gal-Man moiety) substantiated the affinity binding and fluorescence results, displaying similar values of K _a.     

Characterization of a glycoside hydrolase family 42 ��-galactosidase from Deinococcus geothermalis

A putative recombinant β-galactosidase from Deinococcus geothermalis was purified as a single 79 kDa band of 42 U activity/mg using His-Trap affinity chromatography. The molecular mass of the native enzyme was a 158 kDa dimer. The catalytic residues E151 and E325 of β-galactosidase from D. geothermalis were conserved in all aligned GH family 42 β-galactosidases, indicating that this enzyme is also a GH family 42 β-galactosidase. Maximal activity of the enzyme was at pH 6.5 and 60°C. It has a unique hydrolytic activity for p-nitrophenyl(pNP)-β-D-galactopyranoside (k (cat)/K (m) = 69 s(-1) mM(-1)), pNP-β-D-fucopyranoside (13), oNP-β-D-galactopyranoside (9.5), oNP-β-D-fucopyranoside (2.6), lactose (0.97), and pNP-α-L-arabinopyranoside (0.78), whereas no activity, or less than 2% of the pNP-β-D-galactopyranoside activity, for other pNP- and oNP-glycosides.     

Complete backbone and DENQ side chain NMR assignments in proteins from a single experiment: implications to structure–function studies

Resonance assignment is the first and the most crucial step in all nuclear magnetic resonance (NMR) investigations on structure–function relationships in biological macromolecules. Often, the assignment exercise has to be repeated several times when specific interactions with ligands, substrates etc., have to be elucidated for understanding the functional mechanisms. While the protein backbone serves to provide a scaffold, the side chains interact directly with the ligands. Such investigations will be greatly facilitated, if there are rapid methods for obtaining exhaustive information with minimum of NMR experimentation. In this context, we present here a pulse sequence which exploits the recently introduced technique of parallel detection of multiple nuclei, e.g. ¹H and ¹³C, and results in two 3D-data sets simultaneously. These yield complete backbone resonance assignment (¹H^N, ¹⁵N, ¹³CO, ¹Hα/¹³Cα, and ¹Hβ/¹³Cβ chemical shifts) and side chain assignment of D, E, N and Q residues. Such an exhaustive assignment has the potential of yielding accurate 3D structures using one or more of several algorithms which calculate structures of the molecules very reliably on the basis of NMR chemical shifts alone. The side chain assignments of D, E, N, and Q will be extremely valuable for interaction studies with different ligands; D and E side chains are known to be involved in majority of catalytic activities. Utility of this experiment has been demonstrated with Ca²⁺ bound M-crystallin, which contains largely D, E, N and Q residues at the metal binding sites.     

Structural and biochemical characterization of glycoside hydrolase family 79 β-glucuronidase from Acidobacterium capsulatum

We present the first structure of a glycoside hydrolase family 79 β-glucuronidase from Acidobacterium capsulatum, both as a product complex with β-D-glucuronic acid (GlcA) and as its trapped covalent 2-fluoroglucuronyl intermediate. This enzyme consists of a catalytic (β/α)(8)-barrel domain and a β-domain with irregular Greek key motifs that is of unknown function. The enzyme showed β-glucuronidase activity and trace levels of β-glucosidase and β-xylosidase activities. In conjunction with mutagenesis studies, these structures identify the catalytic residues as Glu(173) (acid base) and Glu(287) (nucleophile), consistent with the retaining mechanism demonstrated by (1)H NMR analysis. Glu(45), Tyr(243), Tyr(292)-Gly(294), and Tyr(334) form the catalytic pocket and provide substrate discrimination. Consistent with this, the Y292A mutation, which affects the interaction between the main chains of Gln(293) and Gly(294) and the GlcA carboxyl group, resulted in significant loss of β-glucuronidase activity while retaining the side activities at wild-type levels. Likewise, although the β-glucuronidase activity of the Y334F mutant is ~200-fold lower (k(cat)/K(m)) than that of the wild-type enzyme, the β-glucosidase activity is actually 3 times higher and the β-xylosidase activity is only 2.5-fold lower than the equivalent parameters for wild type, consistent with a role for Tyr(334) in recognition of the C6 position of GlcA. The involvement of Glu(45) in discriminating against binding of the O-methyl group at the C4 position of GlcA is revealed in the fact that the E45D mutant hydrolyzes PNP-β-GlcA approximately 300-fold slower (k(cat)/K(m)) than does the wild-type enzyme, whereas 4-O-methyl-GlcA-containing oligosaccharides are hydrolyzed only 7-fold slower.     

Utilization of lysine 13C-methylation NMR for protein–protein interaction studies

Chemical modification is an easy way for stable isotope labeling of non-labeled proteins. The reductive ¹³C-methylation of the amino group of the lysine side-chain by ¹³C-formaldehyde is a post-modification and is applicable to most proteins since this chemical modification specifically and quickly proceeds under mild conditions such as 4 °C, pH 6.8, overnight. ¹³C-methylation has been used for NMR to study the interactions between the methylated proteins and various molecules, such as small ligands, nucleic acids and peptides. Here we applied lysine ¹³C-methylation NMR to monitor protein–protein interactions. The affinity and the intermolecular interaction sites of methylated ubiquitin with three ubiquitin-interacting proteins were successfully determined using chemical-shift perturbation experiments via the ¹H–¹³C HSQC spectra of the ¹³C-methylated-lysine methyl groups. The lysine ¹³C-methylation NMR results also emphasized the importance of the usage of side-chain signals to monitor the intermolecular interaction sites, and was applicable to studying samples with concentrations in the low sub-micromolar range.     

Cloning, expression, and characterization of a glycoside hydrolase family 50 β-agarase from a marine Agarivorans isolate

The gene for a thermostable β-agarase from Agarivorans sp. JA-1 was cloned and sequenced. It comprised an open reading frame of 2,988 base pairs, which encode a protein of 109,450 daltons consisting of 995 amino acid residues. A comparison of the entire sequence showed that the enzyme has 98.8% sequence similarities to β-agarase from Vibrio sp. JT1070, indicating that it belongs to the family glycoside hydrolase (GH)-50. The gene corresponding to a mature protein of 976 amino acids was inserted and expressed in Escherichia coli. The recombinant β-agarase was purified to homogeneity. It had maximal activity at 40°C and pH 8.0 in the presence of 1 mM NaCl and 1 mM CaCl₂. The enzyme hydrolyzed agarose as well as neoagarohexaose and neoagarotetraose to yield neoagarobiose as the main product. Thus, the enzyme would be useful for the industrial production of neoagarobiose.     

Expression of a Clostridium perfringens genome-encoded putative N-acetylmuramoyl–l-alanine amidase as a potential antimicrobial to control the bacterium

Clostridium perfringens is a gram-positive, spore-forming anaerobic bacterium that plays a substantial role in non-foodborne human, animal, and avian diseases as well as human foodborne disease. Previously discovered C. perfringens bacteriophage lytic enzyme amino acid sequences were utilized to identify putative prophage lysins or autolysins by BLAST analyses encoded by the genomes of C. perfringens isolates. A predicted N-acetylmuramoyl–l-alanine amidase or MurNAc–LAA (also known as peptidoglycan aminohydrolase, NAMLA amidase, NAMLAA, amidase 3, and peptidoglycan amidase; EC 3.5.1.28) was identified that would hydrolyze the amide bond between N-acetylmuramoyl and l-amino acids in certain cell wall glycopeptides. The gene encoding this protein was subsequently cloned from genomic DNA of a C. perfringens isolate by polymerase chain reaction, and the gene product (PlyCpAmi) was expressed to determine if it could be utilized as an antimicrobial to control the bacterium. By spot assay, lytic zones were observed for the purified amidase and the E. coli expression host cellular lysate containing the amidase gene. Turbidity reduction and plate counts of C. perfringens cultures were significantly reduced by the expressed protein and observed morphologies for cells treated with the amidase appeared vacuolated, non-intact, and injured compared to the untreated cells. Among a variety of C. perfringens strains, there was little gene sequence heterogeneity that varied from 1 to 21 nucleotide differences. The results further demonstrate that it is possible to discover lytic proteins encoded in the genomes of bacteria that could be utilized to control bacterial pathogens.     

Structural and enzymatic characterization of a glycoside hydrolase family 31 α-xylosidase from Cellvibrio japonicus involved in xyloglucan saccharification

The desire for improved methods of biomass conversion into fuels and feedstocks has re-awakened interest in the enzymology of plant cell wall degradation. The complex polysaccharide xyloglucan is abundant in plant matter, where it may account for up to 20% of the total primary cell wall carbohydrates. Despite this, few studies have focused on xyloglucan saccharification, which requires a consortium of enzymes including endo-xyloglucanases, α-xylosidases, β-galactosidases and α-L-fucosidases, among others. In the present paper, we show the characterization of Xyl31A, a key α-xylosidase in xyloglucan utilization by the model Gram-negative soil saprophyte Cellvibrio japonicus. CjXyl31A exhibits high regiospecificity for the hydrolysis of XGOs (xylogluco-oligosaccharides), with a particular preference for longer substrates. Crystallographic structures of both the apo enzyme and the trapped covalent 5-fluoro-β-xylosyl-enzyme intermediate, together with docking studies with the XXXG heptasaccharide, revealed, for the first time in GH31 (glycoside hydrolase family 31), the importance of a PA14 domain insert in the recognition of longer oligosaccharides by extension of the active-site pocket. The observation that CjXyl31A was localized to the outer membrane provided support for a biological model of xyloglucan utilization by C. japonicus, in which XGOs generated by the action of a secreted endo-xyloglucanase are ultimately degraded in close proximity to the cell surface. Moreover, the present study diversifies the toolbox of glycosidases for the specific modification and saccharification of cell wall polymers for biotechnological applications.     

Biotechnological potential of novel glycoside hydrolase family 70 enzymes synthesizing α-glucans from starch and sucrose

Transglucosidases belonging to the glycoside hydrolase (GH) family 70 are promising enzymatic tools for the synthesis of α-glucans with defined structures from renewable sucrose and starch substrates. Depending on the GH70 enzyme specificity, α-glucans with different structures and physicochemical properties are produced, which have found diverse (potential) commercial applications, e.g. in food, health and as biomaterials. Originally, the GH70 family was established only for glucansucrase enzymes of lactic acid bacteria that catalyze the synthesis of α-glucan polymers from sucrose. In recent years, we have identified 3 novel subfamilies of GH70 enzymes (designated GtfB, GtfC and GtfD), inactive on sucrose but converting starch/maltodextrin substrates into novel α-glucans. These novel starch-acting enzymes considerably enlarge the panel of α-glucans that can be produced. They also represent very interesting evolutionary intermediates between sucrose-acting GH70 glucansucrases and starch-acting GH13 α-amylases. Here we provide an overview of the repertoire of GH70 enzymes currently available with focus on these novel starch-acting GH70 enzymes and their biotechnological potential. Moreover, we discuss key developments in the understanding of structure-function relationships of GH70 enzymes in the light of available three-dimensional structures, and the protein engineering strategies that were recently applied to expand their natural product specificities.     

Decoupled side chain and backbone dynamics for proton translocation – M2 of influenza A

The 97 amino acid bitopic membrane protein M2 of influenza A forms a tetrameric bundle in which two of the monomers are covalently linked via a cysteine bridge. In its tetrameric assembly the protein conducts protons across the viral envelope and within intracellular compartments during the infectivity cycle of the virus. A key residue in the translocation of the protons is His-37 which forms a planar tetrad in the configuration of the bundle accepting and translocating the incoming protons from the N terminal side, exterior of the virus, to the C terminal side, inside the virus. With experimentally available data from NMR spectroscopy of the transmembrane domains of the tetrameric M2 bundle classical MD simulations are conducted with the protein bundle in different protonation stages in respect to His-37. A full correlation analysis (FCA) of the data sets with the His-37 tetrad either in a fully four times unprotonated or protonated state, assumed to mimic high and low pH in vivo, respectively, in both cases reveal asymmetric backbone dynamics. His-37 side chain rotation dynamics is increased at full protonation of the tetrad compared to the dynamics in the fully unprotonated state. The data suggest that proton translocation can be achieved by decoupled side chain or backbone dynamics.

Open image in new window

Graphical abstract Visualization of the tetrameric bundle of the transmembrane domains of M2 of influenza A after 200 ns of MD simulations (upper left). The four histidine residues 37 are either not protonated as in M2⁰ or fully protonated is in M2⁴⁺. The asymmetric dynamics of the backbones are shown after a full correlation analysis (FCA) in blue (lower left). The rotational dynamics of the χ2 dihedral angles of the histidines in M2⁰ (upper right) are less than those in M2⁴⁺ (lower right)

     

NMR assignment of the apo-form of a Desulfovibrio gigas protein containing a novel Mo–Cu cluster

We report the 98% assignment of the apo-form of an orange protein, containing a novel Mo–Cu cluster isolated from Desulfovibrio gigas. This protein presents a region where backbone amide protons exchange fast with bulk solvent becoming undetectable. These residues were assigned using ¹³C-detection experiments.     

Structural and functional analyses of a glycoside hydrolase family 5 enzyme with an unexpected β-fucosidase activity

We present characterization of PbFucA, a family 5 glycoside hydrolase (GH5) from Prevotella bryantii B(1)4. While GH5 members typically are xylanases, PbFucA shows no activity toward xylan polysaccharides. A screen against a panel of p-nitrophenol coupled sugars identifies PbFucA as a β-D-fucosidase. We also present the 2.2 ? resolution structure of PbFucA and use structure-based mutational analysis to confirm the role of catalytically essential residues. A comparison of the active sites of PbFucA with those of family 5 and 51 glycosidases reveals that while the essential catalytic framework is identical between these enzymes, the steric contours of the respective active site clefts are distinct and likely account for substrate discrimination. Our results show that members of this cluster of orthologous group (COG) 5520 have β-D-fucosidase activities, despite showing an overall sequence and structural similarity to GH-5 xylanases.     

Structural and biochemical insights into the substrate-binding mechanism of a novel glycoside hydrolase family 134 β-mannanase

Mannan is one of the major constituent groups of hemicellulose, which is a renewable resource from higher plants. β-Mannanases are enzymes capable of degrading lignocellulosic biomass. Here, an endo-β-mannanase from Rhizopus microsporus (RmMan134A) was cloned and expressed. The recombinant RmMan134A showed maximal activity at pH?5.0 and 50?°C, and exhibited high specific activity towards locust bean gum (2337?U/mg). To gain insight into the substrate-binding mechanism of RmMan134A, four complex structures (RmMan134A–M3, RmMan134A-M4, RmMan134A-M5 and RmMan134A-M6) were further solved. These structures showed that there were at least seven subsites (?3 to +4) in the catalytic groove of RmMan134A. Mannose in the ?1 subsite hydrogen bonded with His113 and Tyr131, revealing a unique conformation. Lys48 and Val159 formed steric hindrance, which impedes to bond with galactose branches. In addition, the various binding modes of RmMan134A–M5 indicated that subsites ?2 to +2 are indispensable during the hydrolytic process. The structure of RmMan134A–M4 showed that mannotetrose only binds at subsites +1 to +4, and RmMan134A could therefore not hydrolyze mannan oligosaccharides with degree of polymerization ≤4. Through rational design, the specific activity and optimal conditions of RmMan134A were significantly improved. The purpose of this paper is to investigate the structure and function of fungal GH family 134 β-1,4-mannanases, and substrate-binding mechanism of GH family 134 members.     

Site-specific ¹⁹F NMR chemical shift and side chain relaxation analysis of a membrane protein labeled with an unnatural amino acid

Site‐specific ¹⁹F chemical shift and side chain relaxation analysis can be applied on large size proteins. Here, one‐dimensional ¹⁹F spectra and T₁, T₂ relaxation data were acquired on a SH3 domain in aqueous buffer containing 60% glycerol, and a nine‐transmembrane helices membrane protein diacyl‐glycerol kinase (DAGK) in dodecyl phosphochoine (DPC) micelles. The high quality of the data indicates that this method can be applied to site‐specifically analyze side chain internal mobility of membrane proteins or large size proteins.     

Generation and characterization of an inter-generic bivalent alpha domain fusion protein αCS from Clostridium perfringens and Staphylococcus aureus for concurrent diagnosis and therapeutic applications

Aim: To evaluate an inter‐generic recombinant alpha domain fusion protein for simultaneous detection and neutralization of Clostridium perfringens and Staphylococcus aureus alpha toxins. Methods and Results: Truncated portions of clostridial and staphylococcal alpha haemolysin genes were PCR amplified and linked to each other through a hydrophilic flexible Glycine linker sequence using overlap‐extension PCR to form a chimeric gene αCS. The recombinant αCS fusion protein was expressed and characterized for its toxicity, cell binding capacity and haemolysis inhibition properties. The fusion protein was nontoxic and effectively retarded staphylococcal alpha haemolysis, probably by competitively interacting with putative staphylococcal alpha haemolysin receptors on erythrocytes. Murine hyperimmune polysera raised against r‐αCS specifically detected 42‐kDa and 33‐kDa proteins when culture supernatants of Cl. perfringens (clostridial alpha toxin) and Staph. aureus (staphylococcal alpha toxin), respectively, were analysed in Western blot. The polyclonal antisera effectively diminished the haemolytic action of both the wild‐type toxins in vitro. Conclusions: The r‐αCS fusion protein was nontoxic competitive inhibitor of staphylococcal alpha haemolysin. The protein elicited specific immune response against Cl. perfringens and Staph. aureus alpha toxins. The antisera also neutralized the toxicities of both the native wild‐type toxins in vitro. Significance of the Study: The bivalent recombinant αCS protein could be a novel intervention in the field of diagnostics and therapeutics against Cl. perfringens and Staph. aureus infections, particularly, in case of co‐infections like gangrenous ischaemia, gangrenous mastitis, etc.     

IRAA: A statistical tool for investigating a protein–protein interaction interface from multiple structures

Understanding protein–protein interactions (PPIs) is fundamental to infer how different molecular systems work. A major component to model molecular recognition is the buried surface area (BSA), that is, the area that becomes inaccessible to solvent upon complex formation. To date, many attempts tried to connect BSA to molecular recognition principles, and in particular, to the underlying binding affinity. However, the most popular approach to calculate BSA is to use a single (or in some cases few) bound structures, consequently neglecting a wealth of structural information of the interacting proteins derived from ensembles corresponding to their unbound and bound states. Moreover, the most popular method inherently assumes the component proteins to bind as rigid entities. To address the above shortcomings, we developed a Monte Carlo method-based Interface Residue Assessment Algorithm (IRAA), to calculate a combined distribution of BSA for a given complex. Further, we apply our algorithm to human ACE2 and SARS-CoV-2 Spike protein complex, a system of prime importance. Results show a much broader distribution of BSA compared to that obtained from only the bound structure or structures and extended residue members of the interface with implications to the underlying biomolecular recognition. We derive that specific interface residues of ACE2 and of S-protein are consistently highly flexible, whereas other residues systematically show minor conformational variations. In effect, IRAA facilitates the use of all available structural data for any biomolecular complex of interest, extracting quantitative parameters with statistical significance, thereby providing a deeper biophysical understanding of the molecular system under investigation.     

155 3D QSAR and protein–protein interaction studies on neuraminidase against Clostridium perfringens: An approach toward target identification using structure-based drug designing

The rapid onset of resistance to new drugs and emergence of antibiotic resistant bacteria has led to resurgence in life-threatening bacterial infections. These problems have revitalized interest in antibiotics and lead to new research. To gain further insight between structural and biological activity of Clostridium perfringens, a gram-positive anaerobe responsible for food poisoning, mynecrosis in wound infections, and enterotoxemia in humans. We have considered various in silico approaches for developing new drug leads, based on small ligand structure. The importance of neuraminidases in the virulence of C. perfringens makes it a potent target for the studies of drug designing against this microbe. Natural products or their direct derivatives play crucial roles in many diseases. In the present study, 3D QSAR analysis using kNN-MFA method was performed on a series of pterocarpan derivatives as Clostridial neuraminidase inhibitors. Twenty-five compounds using random selection and sphere exclusion method for the division of dataset into training and test set were chosen. kNN-MFA methodology with stepwise, simulated annealing and genetic algorithm was used for model building and four predictive models have been generated. The most significant model has a high internal predictivity of 64.80% (q ²?=?0.6480) and an external predictivity of 95.46% (r ²?=?0.9546). Model showed that electrostatic and steric interactions play important role in determining neuraminidase inhibitory activity. The kNN-MFA plots provide further understanding of the relationship between structural features of substituted pterocarpan derivatives and their activities which were applied for designing new compounds as inhibitors. The drug likeliness of these compounds was checked through ADME property prediction and their interaction with the neuraminidase was checked by molecular docking studies. Moreover, analysis of protein–protein interaction network of NanI in C. perfringens was done.     

Estimation of binding parameters for the protein–protein interaction using a site-directed spin labeling and EPR spectroscopy

Sensitivity of the electron paramagnetic resonance (CW EPR) to molecular tumbling provides potential means for studying processes of molecular association. It uses spin-labeled macromolecules, whose CW EPR spectra may change upon binding to other macromolecules. When a spin-labeled molecule is mixed with its liganding partner, the EPR spectrum constitutes a linear combination of spectra of the bound and unbound ligand (as seen in our example of spin-labeled cytochrome c ₂ interacting with cytochrome bc ₁ complex). In principle, the fraction of each state can be extracted by the numerical decomposition of the spectrum; however, the accuracy of such decomposition may often be compromised by the lack of the spectrum of the fully bound ligand, imposed by the equilibrium nature of molecular association. To understand how this may affect the final estimation of the binding parameters, such as stoichiometry and affinity of the binding, a series of virtual titration experiments was conducted. Our non-linear regression analysis considered a case in which only a single class of binding sites exists, and a case in which classes of both specific and non-specific binding sites co-exist. The results indicate that in both models, the error due to the unknown admixture of the unbound ligand component in the EPR spectrum causes an overestimation of the bound fraction leading to the bias in the dissociation constant. At the same time, the stoichiometry of the binding remains relatively unaffected, which overall makes the decomposition of the EPR spectrum an attractive method for studying protein–protein interactions in equilibrium. Our theoretical treatment appears to be valid for any spectroscopic techniques dealing with overlapping spectra of free and bound component.