Structural changes of active site cleft and different saccharide binding modes in human lysozyme co-crystallized with hexa-N-acetyl-chitohexaose at pH 4.0

Chitosan-coated liposome with lysozyme-responsive properties for on-demand release of levofloxacin

As an anti-infection antibiotic delivery route, a drug-controlled release system based on a specific condition stimulus response can enhance drug stability and bioavailability, reduce antibiotic resistance, achieve on-demand release and improve targeting and utilization efficiency. In this study, chitosan-coated liposomes containing levofloxacin (Lef@Lip@CS) were prepared with lysozyme in body fluids serving as an intelligent "switch" to enable accurate delivery of antibiotics through the catalytic degradation ability of chitosan. Good liposome encapsulation efficacy (64.89 ± 1.86 %) and loading capacity (5.28 ± 0.18 %) were achieved. The controlled-release behavior and morphological characterization before and after enzymatic hydrolysis confirmed that the levofloxacin release rate depended on the lysozyme concentration and the degrees of deacetylation of chitosan. In vitro bacteriostatic experiments showed significant differences in the effects of Lef@Lip@CS before and after enzyme addition, with 6-h inhibition rate of 72.46 % and 100 %, and biofilm removal rates of 51 % and 71 %, respectively. These findings show that chitosan-coated liposomes are a feasible drug delivery system responsive to lysozyme stimulation.

Fluorescence-Polarization-Based Assaying of Lysozyme with Chitooligosaccharide Tracers

Lysozyme is a well-known enzyme found in many biological fluids which plays an important role in the antibacterial protection of humans and animals. Lysozyme assays are used for the diagnosis of a number of diseases and utilized in immunohistochemistry, genetic and cellular engineering studies. The assaying methods are divided into two categories measuring either the concentration of lysozyme as a protein or its activity as an enzyme. While the first category of methods traditionally uses an enzyme-linked immunosorbent assay (ELISA), the methods for the determination of the enzymatic activity of lysozyme use either live bacteria, which is rather inconvenient, or natural peptidoglycans of high heterogeneity and variability, which leads to the low reproducibility of the assay results. In this work, we propose the use of a chemically synthesized substrate of a strictly defined structure to measure in a single experiment both the concentration of lysozyme as a protein and its enzymatic activity by means of the fluorescence polarization (FP) method. Chito-oligosaccharides of different chain lengths were fluorescently labeled and tested leading to the selection of the pentasaccharide as the optimal size tracer and the further optimization of the assay conditions for the accurate (detection limit 0.3 μM) and rapid (<30 min) determination of human lysozyme. The proposed protocol was applied to assay human lysozyme in tear samples and resulted in good correlation with the reference assay. The use of synthetic fluorescently labeled tracer, in contrast to natural peptidoglycan, in FP analysis allows for the development of a reproducible method for the determination of lysozyme activity.

Chitosan-based oral hydrogel formulations of β-galactosidase to improve enzyme supplementation therapy for lactose intolerance

β-Galactosidase supplementation plays an important role in the life of people with lactose intolerance. However, these formulations are rendered ineffective by the low pH and pepsin in the stomach and pancreatic proteases in the intestine. Therefore, it is necessary to develop oral transport systems for carrying this enzyme in the active form up to the intestine, where the lactose digestion occurs. In this research, a new hydrogel was developed that could potentially be used for enzyme supplement therapy. In this regard, the chitosan-based β-Gal formulations described in the manuscript are an alternative long-acting preparation to the so far available preparations that allow for enzyme protection and mucosal targeting. These hydrogels were prepared from chitosan and polyethylene glycol and contained a covalently immobilized β-galactosidase from Aspergillus oryzae. The β-galactosidase in the hydrogel was protected from degradation in a gastric medium at a pH of 2.5 and retained 75 % of its original activity under subsequent intestinal conditions. In the case of a simulated gastric fluid with a pH of 1.5, a copolymer containing methacrylic acid functional groups was sufficient to protect the hybrid hydrogel from the extremely acidic pH. In addition, the surface of the hydrogel was chemically modified with thiol and amidine groups, which increased the binding to intestinal mucin by 20 % compared with the unmodified hydrogel. These results represent a promising approach for oral transport as a reservoir for β-galactosidase in the small intestine to reduce the symptoms of hypolactasia.

Electrostatic interactions of enzymes in non-aqueous conditions: insights from molecular dynamics simulations

Electrostatic interactions of enzymes and their effects on enzyme activity and stability are poorly understood in non-aqueous conditions. Here, we investigate the contribution of the electrostatic interactions on the stability and activity of enzymes in the non-aqueous environment using molecular dynamics simulations. Lipase was selected as active and lysozyme as inactive model enzymes in non-aqueous media. Hexane was used as a common non-aqueous solvent model. In agreement with the previous experiments, simulations show that lysozyme has more structural instabilities than lipase in hexane. The number of hydrogen bonds and salt bridges of both enzymes is dramatically increased in hexane. In contrast to the other opinions, we show that the increase of the electrostatic interactions in non-aqueous media is not so favorable for enzymatic function and stability. In this condition, the newly formed hydrogen bonds and salt bridges can partially denature the local structure of the enzymes. For lysozyme, the changes in electrostatic interactions occur in all domains including the active site cleft, which leads to enzyme inactivation and destabilization. Interestingly, most of the changes in electrostatic interactions of lipase occur far from the active site regions. Therefore, the active site entrance regions remain functional in hexane. The results of this study reveal how the changes in electrostatic interactions can affect enzyme stability and activity in non-aqueous conditions. Moreover, we show for the first time how some enzymes, such as lipase, remain active in a non-aqueous environment.Communicated by Ramaswamy H. Sarma.

A multilayer dynamic perturbation analysis method for predicting ligand-protein interactions

BACKGROUND

Ligand-protein interactions play a key role in defining protein function, and detecting natural ligands for a given protein is thus a very important bioengineering task. In particular, with the rapid development of AI-based structure prediction algorithms, batch structural models with high reliability and accuracy can be obtained at low cost, giving rise to the urgent requirement for the prediction of natural ligands based on protein structures. In recent years, although several structure-based methods have been developed to predict ligand-binding pockets and ligand-binding sites, accurate and rapid methods are still lacking, especially for the prediction of ligand-binding regions and the spatial extension of ligands in the pockets.

RESULTS

In this paper, we proposed a multilayer dynamics perturbation analysis (MDPA) method for predicting ligand-binding regions based solely on protein structure, which is an extended version of our previously developed fast dynamic perturbation analysis (FDPA) method. In MDPA/FDPA, ligand binding tends to occur in regions that cause large changes in protein conformational dynamics. MDPA, examined using a standard validation dataset of ligand-protein complexes, yielded an averaged ligand-binding site prediction Matthews coefficient of 0.40, with a prediction precision of at least 50% for 71% of the cases. In particular, for 80% of the cases, the predicted ligand-binding region overlaps the natural ligand by at least 50%. The method was also compared with other state-of-the-art structure-based methods.

CONCLUSIONS

MDPA is a structure-based method to detect ligand-binding regions on protein surface. Our calculations suggested that a range of spaces inside the protein pockets has subtle interactions with the protein, which can significantly impact on the overall dynamics of the protein. This work provides a valuable tool as a starting point upon which further docking and analysis methods can be used for natural ligand detection in protein functional annotation. The source code of MDPA method is freely available at: https://github.com/mingdengming/mdpa .

Crystal Structure of Human Lysozyme Complexed with N-Acetyl-α-d-glucosamine

Human lysozyme is a natural non-specific immune protein that participates in the immune response of infants against bacterial and viral infections. Lysozyme is a well-known hydrolase that cleaves peptidoglycan in bacterial cell walls. Several crystal structures of human lysozyme have been reported, but little is known regarding how it recognizes sugar molecules. In this study, the crystal structures of human lysozyme in its native and two N-acetyl-α-d-glucosamine (α-D-NAG)-bound forms were determined at 1.3 Å and 1.55/1.60 Å resolution, respectively. Human lysozyme formed a typical c-type lysozyme fold and the α-D-NAG molecule was bound to the middle of subsites C and D. The N-acetyl and glucosamine groups of α-D-NAG were stabilized by hydrophobic interactions (Val117, Ala126, and Trp127), hydrogen bonds (Asn64, Asn78, Ala126, and Val128), and water bridges. Conformational changes of Arg80, Tyr81, Val128, and Arg131 of human lysozyme were observed due to the interactions of α-D-NAG with the active-site cleft. The binding configuration of α-D-NAG in human lysozyme was distinct compared with that of other sugar-bound lysozymes. Findings from this structural analysis provide a better understanding of the sugar recognition of human lysozyme during the immune response to microbial pathogens.

Sequence-based prediction of protein binding regions and drug-target interactions

Identifying drug–target interactions (DTIs) is important for drug discovery. However, searching all drug–target spaces poses a major bottleneck. Therefore, recently many deep learning models have been proposed to address this problem. However, the developers of these deep learning models have neglected interpretability in model construction, which is closely related to a model’s performance. We hypothesized that training a model to predict important regions on a protein sequence would increase DTI prediction performance and provide a more interpretable model. Consequently, we constructed a deep learning model, named Highlights on Target Sequences (HoTS), which predicts binding regions (BRs) between a protein sequence and a drug ligand, as well as DTIs between them. To train the model, we collected complexes of protein–ligand interactions and protein sequences of binding sites and pretrained the model to predict BRs for a given protein sequence–ligand pair via object detection employing transformers. After pretraining the BR prediction, we trained the model to predict DTIs from a compound token designed to assign attention to BRs. We confirmed that training the BRs prediction model indeed improved the DTI prediction performance. The proposed HoTS model showed good performance in BR prediction on independent test datasets even though it does not use 3D structure information in its prediction. Furthermore, the HoTS model achieved the best performance in DTI prediction on test datasets. Additional analysis confirmed the appropriate attention for BRs and the importance of transformers in BR and DTI prediction. The source code is available on GitHub (https://github.com/GIST-CSBL/HoTS).

Insights on peptide topology in the computational design of protein ligands: the example of lysozyme binding peptides

Herein, we compared the ability of linear and cyclic peptides generated in silico to target different protein sites: internal pockets and solvent-exposed sites. We selected human lysozyme (HuL) as a model target protein combined with the computational evolution of linear and cyclic peptides. The sequence evolution of these peptides was based on the PARCE algorithm. The generated peptides were screened based on their aqueous solubility and HuL binding affinity. The latter was evaluated by means of scoring functions and atomistic molecular dynamics (MD) trajectories in water, which allowed prediction of the structural features of the protein-peptide complexes. The computational results demonstrated that cyclic peptides constitute the optimal choice for solvent exposed sites, while both linear and cyclic peptides are capable of targeting the HuL pocket effectively. The most promising binders found in silico were investigated experimentally by surface plasmon resonance (SPR), nuclear magnetic resonance (NMR), and electrospray ionization mass spectrometry (ESI-MS) techniques. All tested peptides displayed dissociation constants in the micromolar range, as assessed by SPR; however, both NMR and ESI-MS suggested multiple binding modes, at least for the pocket binding peptides. A detailed NMR analysis confirmed that both linear and cyclic pocket peptides correctly target the binding site they were designed for.

Generic nature of the condensed states of proteins

Proteins undergoing liquid-liquid phase separation are being discovered at an increasing rate. Since at the high concentrations present in the cell most proteins would be expected to form a liquid condensed state, this state should be considered to be a fundamental state of proteins along with the native state and the amyloid state. Here we discuss the generic nature of the liquid-like and solid-like condensed states, and describe a wide variety of biological functions conferred by these condensed states.

A lysozyme with altered substrate specificity facilitates prey cell exit by the periplasmic predator Bdellovibrio bacteriovorus

Lysozymes are among the best-characterized enzymes, acting upon the cell wall substrate peptidoglycan. Here, examining the invasive bacterial periplasmic predator Bdellovibrio bacteriovorus, we report a diversified lysozyme, DslA, which acts, unusually, upon (GlcNAc-) deacetylated peptidoglycan. B. bacteriovorus are known to deacetylate the peptidoglycan of the prey bacterium, generating an important chemical difference between prey and self walls and implying usage of a putative deacetyl-specific “exit enzyme”. DslA performs this role, and ΔDslA strains exhibit a delay in leaving from prey. The structure of DslA reveals a modified lysozyme superfamily fold, with several adaptations. Biochemical assays confirm DslA specificity for deacetylated cell wall, and usage of two glutamate residues for catalysis. Exogenous DslA, added ex vivo, is able to prematurely liberate B. bacteriovorus from prey, part-way through the predatory lifecycle. We define a mechanism for specificity that invokes steric selection, and use the resultant motif to identify wider DslA homologues.

The bacterial periplasmic predator Bdellovibrio bacteriovorus deacetylates the peptidoglycan of the prey bacterium early upon invasion. Here, the authors identify and characterize a Bdellovibrio lysozyme that acts specifically on deacetylated peptidoglycan and is important for periplasmic exit.

Development of a High Throughput Screen for the Identification of Inhibitors of Peptidoglycan O-Acetyltransferases, New Potential Antibacterial Targets

The O-acetylation of peptidoglycan occurs in many Gram-negative and most Gram-positive pathogens and this modification to the essential wall polymer controls the lytic activity of the autolysins, particularly the lytic transglycosylases, and inhibits that of the lysozymes of innate immunity systems. As such, the peptidoglycan O-acetyltransferases PatA/B and OatA are recognized as virulence factors. In this study, we present the high throughput screening of small compound libraries to identify the first known inhibitors of these enzymes. The fluorometric screening assay developed involved monitoring the respective O-acetyltransferases as esterases using 4-methylumbelliferylacetate as substrate. Pilot screens of 3921 compounds validated the usefulness of the HTS protocol. A number of potential inhibitors were identified amongst a total of 145,000 low molecular-weight compounds, some of which were common to both enzymes, while others were unique to each. After eliminating a number of false positives in secondary screens, dose response curves confirmed the apparent specificity of a benzothiazolyl-pyrazolo-pyridine as an inhibitor of Neisseria gonorrhoeae PatB, and several coumarin-based compounds as inhibitors of both this PatB and OatA from Staphylococcus aureus. The benzothiazolyl-pyrazolo-pyridine was determined to be a non-competitive inhibitor of PatB with a Ki of 126 µM. At 177 µg/mL and close to its solubility limit, this compound caused a 90% reduction in growth of N. gonorrhoeae, while growth of Escherichia coli, a bacterium that lacks PatB and, hence, does not produce O-acetylated peptidoglycan, was unaffected. These data provide preliminary proof of concept that peptidoglycan O-acetyltransferases would serve as useful antibacterial targets.

Nanomedical Relevance of the Intermolecular Interaction Dynamics-Examples from Lysozymes and Insulins

Insulin and lysozyme share the common features of being prone to aggregate and having biomedical importance. Encapsulating lysozyme and insulin in micellar nanoparticles probably would prevent aggregation and facilitate oral drug delivery. Despite the vivid structural knowledge of lysozyme and insulin, the environment-dependent oligomerization (dimer, trimer, and multimer) and associated structural dynamics remain elusive. The knowledge of the intra- and intermolecular interaction profiles has cardinal importance for the design of encapsulation protocols. We have employed various biophysical methods such as NMR spectroscopy, X-ray crystallography, Thioflavin T fluorescence, and atomic force microscopy in conjugation with molecular modeling to improve the understanding of interaction dynamics during homo-oligomerization of lysozyme (human and hen egg) and insulin (porcine, human, and glargine). The results obtained depict the atomistic intra- and intermolecular interaction details of the homo-oligomerization and confirm the propensity to form fibrils. Taken together, the data accumulated and knowledge gained will further facilitate nanoparticle design and production with insulin or lysozyme-related protein encapsulation.

Lysozyme's lectin-like characteristics facilitates its immune defense function

Interactions between human lysozyme (HL) and the lipopolysaccharide (LPS) of Klebsiella pneumoniae O1, a causative agent of lung infection, were identified by surface plasmon resonance. To characterize the molecular mechanism of this interaction, HL binding to synthetic disaccharides and tetrasaccharides representing one and two repeating units, respectively, of the O-chain of this LPS were studied. pH-dependent structural rearrangements of HL after interaction with the disaccharide were observed through nuclear magnetic resonance. The crystal structure of the HL-tetrasaccharide complex revealed carbohydrate chain packing into the A, B, C, and D binding sites of HL, which primarily occurred through residue-specific, direct or water-mediated hydrogen bonds and hydrophobic contacts. Overall, these results support a crucial role of the Glu35/Asp53/Trp63/Asp102 residues in HL binding to the tetrasaccharide. These observations suggest an unknown glycan-guided mechanism that underlies recognition of the bacterial cell wall by lysozyme and may complement the HL immune defense function.

Structural, Functional and Phylogenetic Analysis of Sperm Lysozyme-Like Proteins

Sperm lysozyme-like proteins belonging to c-type lysozyme family evolved in multiple forms. Lysozyme-like proteins, viz., LYZL2, LYZL3 or SLLP1, LYZL4, LYZL5 and LYZL6 are expressed in the testis of mammals. Not all members of LYZL family have been uniformly and unambiguously identified in the genome and proteome of mammals. Some studies suggested a role of SLLP1 and LYZL4 in fertilization; however, the function of other LYZL proteins is unknown. We identified all known forms of LYZL proteins in buffalo sperm by LC-MS/MS. Cloning and sequence analysis of the Lyzl cDNA showed 38-50% identity at amino acid level among the buffalo LYZL paralogs, complete conservation of eight cysteines and other signature sequences of c-type lysozyme family. Catalytic residues in SLLP1, LYZL4 and LYZL5 have undergone replacement. The substrate binding residues showed significant variation in LYZL proteins. Residues at sites 62, 101, 114 in LYZL4; 101 in SLLP1; 37, 62, and 101 in LYZL6 were more variable among diverse species. Sites 63 and 108 occupied by tryptophan were least tolerant to variation. Site 37 also showed lower tolerance to substitution in SLLP1, LYZL4 and LYZL5, but more variable in non-testicular lysozymes. Models of LYZL proteins were created by homology modeling and the substrate binding pockets were analyzed in term of binding energies and contacting residues of LYZL proteins with tri-N-acetylglucosamine (NAG)3 in the A-B-C and B-C-D binding mode. Except LYZL6, LYZL proteins did not show significant difference in binding energies in comparison to hen egg white lysozyme in the A-B-C mode. (NAG)3 binding energy in the B-C-D mode was higher by 1.3-2.2 kcal/mol than in A-B-C mode. Structural analysis indicated that (NAG)3 was involved in making more extensive interactions including hydrogen bonding with LYZL proteins in B-C-D mode than in A-B-C mode. Despite large sequence divergence among themselves and with respect to c-type lysozymes, substrate binding residues as well as hydrogen bonding network between (NAG)3 and proteins were mostly conserved. LYZL5 in buffalo and other mammalian species contained additional 10-12 amino acid sequence at c-terminal that matched with ankyrin repeat domain-containing protein 27. Phylogenetic analysis indicated LYZL2 to be most ancient among all the LYZL proteins and that the evolution of LYZL proteins occurred through several gene duplications preceding the speciation of mammals from other vertebrates as distant as reptiles and amphibians.

Molecular characterization of a c-type lysozyme from the desert locust, Schistocerca gregaria (Orthoptera: Acrididae)

Lysozymes are bacteriolytic peptides that are implicated in the insect nonspecific innate immune responses. In this study, a full-length cDNA encoding a c-type lysozyme from Schistocerca gregaria (SgLys) has been cloned and characterized from the fat body of immune-challenged 5(th) instar. The deduced mature lysozyme is 119 amino acid residues in length, has a calculated molecular mass of 13.4 kDa and an isoelectric point (Ip) of 9.2. SgLys showed high identities with other insect lysozymes, ranging from 41.5% to 93.3% by BLASTp search in NCBI. Eukaryotic in vitro expression of the SgLys ORF (rSgLys) with an apparent molecular mass of ∼16 kDa under SDS-PAGE is close to the calculated molecular weight of the full-length protein. rSgLys displayed growth inhibitory activity against Gram-negative and Gram-positive bacteria. 3D structure modeling of SgLys, based on comparison with that of silkworm lysozyme, and sequence comparison with the helix-loop-helix (α-hairpin) structure of hen egg white lysozyme (HEWL) were employed to interpret the antibacterial potencies. Phylogenetic alignments indicate that SgLys aligns well with insect c-type lysozymes that expressed principally in fat body and hemocytes and whose role has been defined as immune-related. Western blot analysis showed that SgLys expression was highest at 6-12 h post-bacterial challenge and subsequently decreased with time. Transcriptional profiles of SgLys were determined by semi-quantitative RT-PCR analysis. SgLys transcript was upregulated at the highest level in fat body, hemocytes, salivary gland, thoracic muscles, and epidermal tissue. It was expressed in all developmental stages from egg to adult. These data indicate that SgLys is a predominant acute-phase protein that is expressed and upregulated upon immune challenge.

The RpfC (Rv1884) atomic structure shows high structural conservation within the resuscitation-promoting factor catalytic domain

The first structure of the catalytic domain of RpfC (Rv1884), one of the resuscitation-promoting factors (RPFs) from Mycobacterium tuberculosis, is reported. The structure was solved using molecular replacement once the space group had been correctly identified as twinned P21 rather than the apparent C2221 by searching for anomalous scattering sites in P1. The structure displays a very high degree of structural conservation with the previously published structures of the catalytic domains of RpfB (Rv1009) and RpfE (Rv2450). This structural conservation highlights the importance of the versatile domain composition of the RPF family.

Hot spots in a network of functional sites

It is of significant interest to understand how proteins interact, which holds the key phenomenon in biological functions. Using dynamic fluctuations in high frequency modes, we show that the Gaussian Network Model (GNM) predicts hot spot residues with success rates ranging between S 8–58%, C 84–95%, P 5–19% and A 81–92% on unbound structures and S 8–51%, C 97–99%, P 14–50%, A 94–97% on complex structures for sensitivity, specificity, precision and accuracy, respectively. High specificity and accuracy rates with a single property on unbound protein structures suggest that hot spots are predefined in the dynamics of unbound structures and forming the binding core of interfaces, whereas the prediction of other functional residues with similar dynamic behavior explains the lower precision values. The latter is demonstrated with the case studies; ubiquitin, hen egg-white lysozyme and M2 proton channel. The dynamic fluctuations suggest a pseudo network of residues with high frequency fluctuations, which could be plausible for the mechanism of biological interactions and allosteric regulation.

TraG encoded by the pIP501 type IV secretion system is a two-domain peptidoglycan-degrading enzyme essential for conjugative transfer

pIP501 is a conjugative broad-host-range plasmid frequently present in nosocomial Enterococcus faecalis and Enterococcus faecium isolates. We focus here on the functional analysis of the type IV secretion gene traG, which was found to be essential for pIP501 conjugative transfer between Gram-positive bacteria. The TraG protein, which localizes to the cell envelope of E. faecalis harboring pIP501, was expressed and purified without its N-terminal transmembrane helix (TraGΔTMH) and shown to possess peptidoglycan-degrading activity. TraGΔTMH was inhibited by specific lytic transglycosylase inhibitors hexa-N-acetylchitohexaose and bulgecin A. Analysis of the TraG sequence suggested the presence of two domains which both could contribute to the observed cell wall-degrading activity: an N-terminal soluble lytic transglycosylase domain (SLT) and a C-terminal cysteine-, histidine-dependent amidohydrolases/peptidases (CHAP) domain. The protein domains were expressed separately, and both degraded peptidoglycan. A change of the conserved glutamate residue in the putative catalytic center of the SLT domain (E87) to glycine resulted in almost complete inactivity, which is consistent with this part of TraG being a predicted lytic transglycosylase. Based on our findings, we propose that TraG locally opens the peptidoglycan to facilitate insertion of the Gram-positive bacterial type IV secretion machinery into the cell envelope.

Difucosylation of chitooligosaccharides by eukaryote and prokaryote α1,6-fucosyltransferases

BACKGROUND

The synthesis of eukaryotic N-glycans and the rhizobia Nod factor both involve α1,6-fucosylation. These fucosylations are catalyzed by eukaryotic α1,6-fucosyltransferase, FUT8, and rhizobial enzyme, NodZ. The two enzymes have similar enzymatic properties and structures but display different acceptor specificities: FUT8 and NodZ prefer N-glycan and chitooligosaccharide, respectively. This study was conducted to examine the fucosylation of chitooligosaccharides by FUT8 and NodZ and to characterize the resulting difucosylated chitooligosaccharides in terms of their resistance to hydrolysis by glycosidases.

METHODS

The issue of whether FUT8 or NodZ catalyzes the further fucosylation of chitooligosaccharides that had first been monofucosylated by the other. The oligosaccharide products from the successive reactions were analyzed by normal-phase high performance liquid chromatography, mass spectrometry and nuclear magnetic resonance. The effect of difucosylation on sensitivity to glycosidase digestion was also investigated.

RESULTS

Both FUT8 and NodZ are able to further fucosylate the monofucosylated chitooligosaccharides. Structural analyses of the resulting oligosaccharides showed that the reducing terminal GlcNAc residue and the third GlcNAc residue from the non-reducing end are fucosylated via α1,6-linkages. The difucosylation protected the oligosaccharides from extensive degradation to GlcNAc by hexosamidase and lysozyme, and also even from defucosylation by fucosidase.

CONCLUSIONS

The sequential actions of FUT8 and NodZ on common substrates effectively produce site-specific-difucosylated chitooligosaccharides. This modification confers protection to the oligosaccharides against various glycosidases.

GENERAL SIGNIFICANCE

The action of a combination of eukaryotic and bacterial α1,6-fucosyltransferases on chitooligosaccharides results in the formation of difucosylated products, which serves to stabilize chitooligosaccharides against the action of glycosidases.

Changes in Lysozyme Flexibility upon Mutation Are Frequent, Large and Long-Ranged

We investigate changes in human c-type lysozyme flexibility upon mutation via a Distance Constraint Model, which gives a statistical mechanical treatment of network rigidity. Specifically, two dynamical metrics are tracked. Changes in flexibility index quantify differences within backbone flexibility, whereas changes in the cooperativity correlation quantify differences within pairwise mechanical couplings. Regardless of metric, the same general conclusions are drawn. That is, small structural perturbations introduced by single point mutations have a frequent and pronounced affect on lysozyme flexibility that can extend over long distances. Specifically, an appreciable change occurs in backbone flexibility for 48% of the residues, and a change in cooperativity occurs in 42% of residue pairs. The average distance from mutation to a site with a change in flexibility is 17–20 Å. Interestingly, the frequency and scale of the changes within single point mutant structures are generally larger than those observed in the hen egg white lysozyme (HEWL) ortholog, which shares 61% sequence identity with human lysozyme. For example, point mutations often lead to substantial flexibility increases within the β-subdomain, which is consistent with experimental results indicating that it is the nucleation site for amyloid formation. However, β-subdomain flexibility within the human and HEWL orthologs is more similar despite the lowered sequence identity. These results suggest compensating mutations in HEWL reestablish desired properties.

The functional importance of protein dynamics is universally accepted, making the study of dynamical similarities and differences among proteins of the same function an intriguing problem. While some metrics are likely to be conserved across family, differences are also very common. In previous works we have used a Distance Constraint Model to quantify flexibility differences across sets of orthologous proteins, which reproduce this diversity. In the same manner, this work investigates changes occurring upon individual point mutations. Somewhat surprisingly, the small structural perturbations caused by mutation lead to changes throughout the protein. These changes can be quite large, actually surpassing the scale for differences between ortholog pairs. Moreover, changes in flexibility frequently occur at sites far from the mutation site. These results underscore the sensitivity of protein dynamics in connection with allostery, and help explain why differences across protein families are so common.

Robust identification of binding hot spots using continuum electrostatics: application to hen egg-white lysozyme

Binding hot spots, protein regions with high binding affinity, can be identified by using X-ray crystallography or NMR spectroscopy to screen libraries of small organic molecules that tend to cluster at such hot spots. FTMap, a direct computational analogue of the experimental screening approaches, uses 16 different probe molecules for global sampling of the surface of a target protein on a dense grid and evaluates the energy of interaction using an empirical energy function that includes a continuum electrostatic term. Energy evaluation is based on the fast Fourier transform correlation approach, which allows for the sampling of billions of probe positions. The grid sampling is followed by off-grid minimization that uses a more detailed energy expression with a continuum electrostatics term. FTMap identifies the hot spots as consensus clusters formed by overlapping clusters of several probes. The hot spots are ranked on the basis of the number of probe clusters, which predicts their binding propensity. We applied FTMap to nine structures of hen egg-white lysozyme (HEWL), whose hot spots have been extensively studied by both experimental and computational methods. FTMap found the primary hot spot in site C of all nine structures, in spite of conformational differences. In addition, secondary hot spots in sites B and D that are known to be important for the binding of polysaccharide substrates were found. The predicted probe–protein interactions agree well with those seen in the complexes of HEWL with various ligands and also agree with an NMR-based study of HEWL in aqueous solutions of eight organic solvents. We argue that FTMap provides more complete information on the HEWL binding site than previous computational methods and yields fewer false-positive binding locations than the X-ray structures of HEWL from crystals soaked in organic solvents.

Crystal structure of a charge engineered human lysozyme having enhanced bactericidal activity

Human lysozyme is a key component of the innate immune system, and recombinant forms of the enzyme represent promising leads in the search for therapeutic agents able to treat drug-resistant infections. The wild type protein, however, fails to participate effectively in clearance of certain infections due to inherent functional limitations. For example, wild type lysozymes are subject to electrostatic sequestration and inactivation by anionic biopolymers in the infected airway. A charge engineered variant of human lysozyme has recently been shown to possess improved antibacterial activity in the presence of disease associated inhibitory molecules. Here, the 2.04 Å crystal structure of this variant is presented along with an analysis that provides molecular level insights into the origins of the protein's enhanced performance. The charge engineered variant's two mutated amino acids exhibit stabilizing interactions with adjacent native residues, and from a global perspective, the mutations cause no gross structural perturbations or loss of stability. Importantly, the two substitutions dramatically expand the negative electrostatic potential that, in the wild type enzyme, is restricted to a small region near the catalytic residues. The net result is a reduction in the overall strength of the engineered enzyme's electrostatic potential field, and it appears that the specific nature of this remodeled field underlies the variant's reduced susceptibility to inhibition by anionic biopolymers.

Crystal structures of g-type lysozyme from Atlantic cod shed new light on substrate binding and the catalytic mechanism

Crystal structures of Atlantic cod lysozyme have been solved with and without ligand bound in the active site to 1.7 and 1.9 A resolution, respectively. The structures reveal the presence of NAG in the substrate binding sites at both sides of the catalytic Glu73, hence allowing the first crystallographic description of the goose-type (g-type) lysozyme E-G binding sites. In addition, two aspartic acid residues suggested to participate in catalysis (Asp101 and Asp90) were mutated to alanine. Muramidase activity data for two single mutants and one double mutant demonstrates that both residues are involved in catalysis, but Asp101 is the more critical of the two. The structures and activity data suggest that a water molecule is the nucleophile completing the catalytic reaction, and the roles of the aspartic acids are to ensure proper positioning of the catalytic water.

The Chitopentaose Complex of a Mutant Hen Egg-White Lysozyme Displays No Distortion of the –1 Sugar Away from a 4C1 Chair Conformation

Glycosidase inhibitors frequently reflect either the charge or the ‘flattened’ shape of the oxocarbenium-ion like transition state. Much of the impetus for such inhibitory strategies derives from historical studies on ligand binding to hen egg white lysozyme (HEWL); not least those suggesting that product complexes of the enzyme showed distortion of the pyranosides in the –1 subsite. Ironically, while distortion is undoubtedly a defining feature of glycosidases, product complexes themselves are rarely distorted. Here we show that the chitopentaose product complex of a mutant E35Q HEWL, solved at 1.8 Å resolution, is bound with all sugars in 4C1 conformation.

Peptidoglycan degradation by specialized lytic transglycosylases associated with type III and type IV secretion systems

Specialized lytic transglycosylases are muramidases capable of locally degrading the peptidoglycan meshwork of Gram-negative bacteria. Specialized lytic transglycosylase genes are present in clusters encoding diverse macromolecular transport systems. This paper reports the analysis of selected members of the specialized lytic transglycosylase family from type III and type IV secretion systems. These proteins were analysedin vivoby assaying their ability to complement the DNA transfer defect of the conjugative F-like plasmid R1-16 lacking a functional P19 protein, the specialized lytic transglycosylase of this type IV secretion system. Heterologous complementation was accomplished using IpgF from the plasmid-encoded type III secretion system ofShigella sonneiand TrbN from the type IV secretion system of the conjugative plasmid RP4. In contrast, neither VirB1 proteins (Agrobacterium tumefaciens,Brucella suis) nor IagB (Salmonella enterica) could functionally replace P19.In vitro, IpgF, IagB, both VirB1 proteins, HP0523 (Helicobacter pylori) and P19 displayed peptidoglycanase activity in zymogram analyses. Using an established test system and a newly developed assay it was shown that IpgF degraded peptidoglycan in solution. IpgF was active only after removal of the chaperonin GroEL, which co-purified with IpgF and inhibited its enzymic activity. A mutant IpgF protein in which the predicted catalytic amino acid, Glu42, was replaced by Gln, was completely inactive. IpgF-catalysed peptidoglycan degradation was optimal at pH 6 and was inhibited by the lytic transglycosylase inhibitors hexa-N-acetylchitohexaose and bulgecin A.

Structural and functional modeling of human lysozyme reveals a unique nonapeptide, HL9, with anti-HIV activity

We previously reported that lysozyme accounts for anti-HIV activity associated with the beta-core fraction of human chorionic gonadotropin [Lee-Huang, S., Huang, P. L., Sun, Y., Kung, H. F., Blithe, D. L. & Chen, H. C. (1999) Proc Natl Acad Sci U S A 96, 2678-81]. To define the structural and sequence requirements for anti-HIV activity, we carried out peptide fragmentation and activity mapping of human lysozyme. We identified two peptides that consist of 18 and 9 amino acids of human lysozyme (HL18 and HL9), corresponding to residues 98-115 and 107-115. HL18 and HL9 are potent inhibitors of HIV-1 infection and replication with EC(50)s of 50 to 55 nM, comparable to intact lysozyme. Scrambling the sequence or substitution of key arginine or tryptophan residues results in loss of antiviral activity. HL9, with the sequence RAWVAWRNR, is the smallest peptide we identified with full anti-HIV activity. It forms a pocket with its basic residues on the surface of the molecule. HL9 exists as an alpha-helix in native human lysozyme, in a region of the protein distinct from the muramidase catalytic site. Monte Carlo peptide folding energy minimizing simulation modeling and CD studies indicate that helical propensity does not correlate with antiviral activity. HL9 blocks HIV-1 viral entrance and replication, and modulates gene expression of HIV-infected cells, affecting pathways involved in survival, stress, TGFbeta, p53, NFkappaB, protein kinase C and hedgehog signaling.

Pattern similarity study of functional sites in protein sequences: lysozymes and cystatins

BACKGROUND

Although it is generally agreed that topography is more conserved than sequences, proteins sharing the same fold can have different functions, while there are protein families with low sequence similarity. An alternative method for profile analysis of characteristic conserved positions of the motifs within the 3D structures may be needed for functional annotation of protein sequences. Using the approach of quantitative structure-activity relationships (QSAR), we have proposed a new algorithm for postulating functional mechanisms on the basis of pattern similarity and average of property values of side-chains in segments within sequences. This approach was used to search for functional sites of proteins belonging to the lysozyme and cystatin families.

RESULTS

Hydrophobicity and beta-turn propensity of reference segments with 3-7 residues were used for the homology similarity search (HSS) for active sites. Hydrogen bonding was used as the side-chain property for searching the binding sites of lysozymes. The profiles of similarity constants and average values of these parameters as functions of their positions in the sequences could identify both active and substrate binding sites of the lysozyme of Streptomyces coelicolor, which has been reported as a new fold enzyme (Cellosyl). The same approach was successfully applied to cystatins, especially for postulating the mechanisms of amyloidosis of human cystatin C as well as human lysozyme.

CONCLUSION

Pattern similarity and average index values of structure-related properties of side chains in short segments of three residues or longer were, for the first time, successfully applied for predicting functional sites in sequences. This new approach may be applicable to studying functional sites in un-annotated proteins, for which complete 3D structures are not yet available.

The structure of a resuscitation-promoting factor domain from Mycobacterium tuberculosis shows homology to lysozymes

Resuscitation-promoting factor (RPF) proteins reactivate stationary-phase cultures of (G+C)-rich Gram-positive bacteria including the causative agent of tuberculosis, Mycobacterium tuberculosis. We report the solution structure of the RPF domain from M. tuberculosis Rv1009 (RpfB) solved by heteronuclear multidimensional NMR. Structural homology with various glycoside hydrolases suggested that RpfB cleaved oligosaccharides. Biochemical studies indicate that a conserved active site glutamate is important for resuscitation activity. These data, as well as the presence of a clear binding pocket for a large molecule, indicate that oligosaccharide cleavage is probably the signal for revival from dormancy.

NMR studies of protein hydration and TEMPOL accessibility

Understanding the mechanisms of the interaction between a protein surface and its outer molecular environment is of primary relevance for the rational design of new drugs and engineered proteins. Protein surface accessibility is emerging as a new dimension of Structural Biology, since NMR methods have been developed to follow how molecules, even those different from physiological ligands, preferentially approach specific regions of the protein surface. Hen egg-white lysozyme, a paradigmatic example of the state of the art of protein structure and dynamics, has been selected as a model system to study protein surface accessibility. Bound water and soluble spin-labels have been used to investigate the interaction of this enzyme, both free and bound to the inhibitor (NAG)(3), with its molecular environment. No tightly bound water molecules were found inside the enzyme active site, which, conversely, appeared as the most exposed to visits from the soluble paramagnetic probe TEMPOL. From the presented set of data, an integrated view of lysozyme surface accessibility towards water and TEMPOL molecules is obtained.

Low-temperature-induced structural changes in human lysozyme elucidated by three-dimensional NMR spectroscopy

The three-dimensional solution structures of human lysozyme were determined at 35 and 4 degrees C using the heteronuclear multidimensional NMR spectroscopy, which were compared with each other to clarify the structural response of this enzyme to lowering of the temperature. Together with the data of the temperature dependence experiments of the lytic activity against Micrococcus luteus, we consider the implication of the observed structural change for the low-temperature-induced reduction of the activity of human lysozyme. The structures of human lysozyme determined at the two temperatures are found to be similar, both of which comprise four alpha-helices (A- to D-helices) and three antiparallel beta-strands (beta(1)-beta(3)), leading to the constructions of the alpha- and beta-domains as previously identified in the X-ray crystal structure. A significant structural change was observed for the "active site lobe" comprising the loop region connecting C- and D-helices and the following D-helix, which moves toward the active site cleft located between the alpha- and beta-domains so as to obstruct the cleft according to the temperature lowering. It further appeared that the total volume as well as the accessible surface area of human lysozyme decreases with lowering of the temperature, suggesting that the internal cavity of this enzyme shrinks under low temperature environment. Because in human lysozyme the region comprising the active site lobe is responsible for turnover of the enzymatic reaction against the substrate, the low-temperature-induced structural change of the active site lobe presumably controls the efficiency of the lytic activity under low temperatures.

A helix-loop-helix peptide at the upper lip of the active site cleft of lysozyme confers potent antimicrobial activity with membrane permeabilization action

Recently, we have found that partially unfolded lysozyme exerts broad spectrum antimicrobial action in vitro against Gram-negative and Gram-positive bacteria independent of its catalytic activity. In parallel, an internal peptide (residues 98-112) of hen egg white lysozyme, obtained after digestion with clostripain, possessed broad spectrum antimicrobial action in vitro. This internal peptide is part of a helix-loop-helix domain (87-114 sequence of hen lysozyme) located at the upper lip of the active site cleft of lysozyme. The helix-loop-helix (HLH) structures are known motifs commonly found in membrane-active and DNA-binding proteins. To evaluate the contribution of the HLH peptide to the antimicrobial properties of lysozyme, the HLH sequence and its secondary structure derivatives of chicken and human lysozyme were synthesized and tested for antimicrobial activity against several bacterial strains. We found that the full HLH peptide of both chicken and human lysozymes was potently microbicidal against both Gram-positive and Gram-negative bacteria and the fungus Candida albicans. The N-terminal helix of HLH was specifically bactericidal to Gram-positive bacteria, whereas the C-terminal helix was bactericidal to all tested strains. Outer and inner membrane permeabilization studies, as well as measurements of transmembrane electrochemical potentials, provided evidence that HLH peptide and its C-terminal helix domain kill Gram-negative bacteria by crossing the outer membrane via self-promoted uptake and causing damage to the inner membrane through channel formation. The results are discussed in terms of proposed mechanisms for the catalytically independent antimicrobial activity of lysozyme that offer a new strategy for the design of potential antimicrobial drugs in the treatment of infectious diseases.

Crystal structure of the lytic transglycosylase from bacteriophage lambda in complex with hexa-N-acetylchitohexaose

The three-dimensional structure of the lytic transglycosylase from bacteriophage lambda, also known as bacteriophage lambda lysozyme, complexed to the hexasaccharide inhibitor, hexa-N-acetylchitohexaose, has been determined by X-ray crystallography at 2.6 A resolution. The unit cell contains two molecules of the lytic transglycosylase with two hexasaccharides bound. Each enzyme molecule is found to interact with four N-acetylglucosamine units from one hexasaccharide (subsites A-D) and two N-acetylglucosamine units from the second hexasaccharide (subsites E and F), resulting in all six subsites of the active site of this enzyme being filled. This crystallographic structure, therefore, represents the first example of a lysozyme in which all subsites are occupied, and detailed protein-oligosaccharide interactions are now available for this bacteriophage lytic transglycosylase. Examination of the active site furthermore reveals that of the two residues that have been implicated in the reaction mechanism of most other c-type lysozymes (Glu35 and Asp52 in hen egg white lysozyme), only a homologous Glu residue is present. The lambda lytic transglycosylase is therefore functionally closely related to the Escherichia coli Slt70 and Slt35 lytic transglycosylases and goose egg white lysozyme which also lack the catalytic aspartic acid.

Hydrogen bonds with pi-acceptors in proteins: frequencies and role in stabilizing local 3D structures

A comprehensive structural analysis of X--H...pi hydrogen bonding in proteins is performed based on 592 published high-resolution crystal structures (< or = 1.6 A). All potential donors and acceptors are considered, including acidic C--H groups. The sample contains 1311 putative X--H...pi hydrogen bonds with N--H, O--H or S--H donors, that is about one per 10.8 aromatic residues. By far the most efficient pi-acceptor is the side-chain of Trp, which accepts one X--H...pi hydrogen bond per 5.7 residues. The focus of the analysis is on recurrent structural patterns involving regular secondary structure elements. Numerous examples are found where peptide X--H...pi interactions are functional in stabilization of helix termini, strand ends, strand edges, beta-bulges and regular turns. Side-chain X--H...pi hydrogen bonds are formed in considerable numbers in alpha-helices and beta-sheets. Geometrical data on various types of X--H...pi hydrogen bonds are given.

Crystal structures of apo- and holo-bovine alpha-lactalbumin at 2. 2-A resolution reveal an effect of calcium on inter-lobe interactions

High affinity binding of Ca(2+) to alpha-lactalbumin (LA) stabilizes the native structure and is required for the efficient generation of native protein with correct disulfide bonds from the reduced denatured state. A progressive increase in affinity of LA conformers for Ca(2+) as they develop increasingly native structures can account for the tendency of the apo form to assume a molten globule state and the large acceleration of folding by Ca(2+). To investigate the effect of calcium on structure of bovine LA, x-ray structures have been determined for crystals of the apo and holo forms at 2.2-A resolution. In both crystal forms, which were grown at high ionic strength, the protein is in a similar global native conformation consisting of alpha-helical and beta-subdomains separated by a cleft. Even though alternative cations and Ca(2+) liganding solvent molecules are absent, removal of Ca(2+) has only minor effects on the structure of the metal-binding site and a structural change was observed in the cleft on the opposite face of the molecule adjoining Tyr(103) of the helical lobe and Gln(54) of the beta-lobe. Changes include increased separation of the lobes, loss of a buried solvent molecule near the Ca(2+)-binding site, and the replacement of inter- and intra-lobe H-bonds of Tyr(103) by interactions with new immobilized water molecules. The more open cleft structure in the apo protein appears to be an effect of calcium binding transmitted via a change in orientation of helix H3 relative to the beta-lobe to the inter-lobe interface. Calcium is well known to promote the folding of LA. The results from the comparison of apo and holo structures of LA provide high resolution structural evidence that the acceleration of folding by Ca(2+) is mediated by an effect on interactions between the two subdomains.

Rapid evolution in plant chitinases: molecular targets of selection in plant-pathogen coevolution

Many pathogen recognition genes, such as plant R-genes, undergo rapid adaptive evolution, providing evidence that these genes play a critical role in plant-pathogen coevolution. Surprisingly, whether rapid adaptive evolution also occurs in genes encoding other kinds of plant defense proteins is unknown. Unlike recognition proteins, plant chitinases attack pathogens directly, conferring disease resistance by degrading chitin, a component of fungal cell walls. Here, we show that nonsynonymous substitution rates in plant class I chitinase often exceed synonymous rates in the plant genus Arabis (Cruciferae) and in other dicots, indicating a succession of adaptively driven amino acid replacements. We identify individual residues that are likely subject to positive selection by using codon substitution models and determine the location of these residues on the three-dimensional structure of class I chitinase. In contrast to primate lysozymes and plant class III chitinases, structural and functional relatives of class I chitinase, the adaptive replacements of class I chitinase occur disproportionately in the active site cleft. This highly unusual pattern of replacements suggests that fungi directly defend against chitinolytic activity through enzymatic inhibition or other forms of chemical resistance and identifies target residues for manipulating chitinolytic activity. These data also provide empirical evidence that plant defense proteins not involved in pathogen recognition also evolve in a manner consistent with rapid coevolutionary interactions.

Protein-carbohydrate interactions in human lysozyme probed by combining site-directed mutagenesis and affinity labeling

The synergism between apolar and polar interactions in the carbohydrate recognition by human lysozyme (HL) was probed by site-directed mutagenesis and affinity labeling. The three-dimensional structures of the Tyr63-->Leu mutant HL labeled with 2',3'-epoxypropyl beta-glycoside of N,N'-diacetylchitobiose (L63-HL/NAG-NAG-EPO complex) and the Asp102-->Glu mutant HL labeled with the 2',3'-epoxypropyl beta-glycoside of N-acetyllactosamine were revealed by X-ray diffraction at 2.23 and 1.96 A resolution, respectively. Compared to the wild-type HL labeled with the 2', 3'-epoxypropyl beta-glycoside of N,N'-diacetylchitobiose, the N-acetylglucosamine residue at subsite B of the L63-HL/NAG-NAG-EPO complex markedly moved away from the 63rd residue, with substantial loss of hydrogen-bonding interactions. Evidently, the stacking interaction with the aromatic side chain of Tyr63 is essential in positioning the N-acetylglucosamine residue in the productive binding mode. On the other hand, the position of the galactose residue in subsite B of HL is almost unchanged by the mutation of Asp102 to Glu. Most hydrogen bonds, including the one between the carboxylate group of Glu102 and the axial 4-OH group of the galactose residue, were maintained by local movement of the backbone from residues 102-104. In both structures, the conformation of the disaccharide was conserved, reflecting an intrinsic conformational rigidity of the disaccharides. The structural analysis suggested that CH-pi interactions played an important role in the recognition of the carbohydrate residue at subsite B of HL.

Crystal structure of chitosanase from Bacillus circulans MH-K1 at 1.6-A resolution and its substrate recognition mechanism

Chitosanase from Bacillus circulans MH-K1 is a 29-kDa extracellular protein composed of 259 amino acids. The crystal structure of chitosanase from B. circulans MH-K1 has been determined by multiwavelength anomalous diffraction method and refined to crystallographic R = 19.2% (R(free) = 23.5%) for the diffraction data at 1.6-A resolution collected by synchrotron radiation. The enzyme has two globular upper and lower domains, which generate the active site cleft for the substrate binding. The overall molecular folding is similar to chitosanase from Streptomyces sp. N174, although there is only 20% identity at the amino acid sequence level between both chitosanases. However, there are three regions in which the topology is remarkably different. In addition, the disulfide bridge between Cys(50) and Cys(124) joins the beta1 strand and the alpha7 helix, which is not conserved among other chitosanases. The orientation of two backbone helices, which connect the two domains, is also different and is responsible for the differences in size and shape of the active site cleft in these two chitosanases. This structural difference in the active site cleft is the reason why the enzymes specifically recognize different substrates and catalyze different types of chitosan degradation.

Molecular modeling of an active loop structure in lysozyme. Sequence effects or crystal packing?

The Val99-Gly 104 variable region in egg white lysozyme is part of the active site cleft and of the epitope recognized by some monoclonal antibodies. In general, this loop is found in a conformation inflected towards the active site (proximal conformational) such as in free hen lysozyme (HEL). But in a lysozyme such as Japanese quail's (JEL), the loop turns away from the active site cleft (distal conformation). In order to differentiate sequence effects from crystal packing, we generated and refined loop conformations for the 99-104 variable region in lysozyme, then estimated their relative conformational free energies. Some of the results indicate that (i) the flexibility of the 99-104 segment is much greater for HEL than for JEL sequences when unconstrained by the crystal lattice, (ii) for JEL, only distal structures are favored, while for HEL the states span the zone between proximal and distal regions, and (iii) epitopes elucidated from crystal structures may not always be conserved in solution. For the JEL loop, model building shows that an energy-costly distal to proximal transition appears necessary. Finally, analysis of available structural data indicates that changes of humidity, temperature and pressure on loop conformation are negligible.

Molecular dynamics simulation of alpha-lactalbumin and calcium binding c-type lysozyme

Alpha-lactalbumins (LAs) and c-type lysozymes (LYZs) are two classes of proteins which have a 35-40% sequence homology and share a common three dimensional fold but perform different functions. Lysozymes bind and cleave the glycosidic bond linkage in sugars, where as, alpha-lactalbumin does not bind sugar but participates in the synthesis of lactose. Alpha-lactalbumin is a metallo-protein and binds calcium, where as, only a few of the LYZs bind calcium. These proteins consist of two domains, an alpha-helical and a beta-strand domain, separated by a cleft. Calcium is bound at a loop situated at the bottom of the cleft and is important for the structural integrity of the protein. Calcium is an ubiquitous intracellular signal in higher eukaryotes and structural changes induced on calcium binding have been observed in a number of proteins. In the present study, molecular dynamics simulations of equine LYZ and human LA, with and without calcium, were carried out. We detail the differences in the dynamics of equine LYZ and human LA, and discuss it in the light of experimental data already available and relate it to the behavior of the functionally important regions of both the proteins. These simulations bring out the role of calcium in the conformation and dynamics of these metallo-proteins. In the calcium bound LA, the region of the protein around the calcium binding site is not only frozen but the atomic fluctuations are found to increase away from the binding site and peak at the exposed sites of the protein. This channeling of fluctuations away from the metal binding site could serve as a general mechanism by which the effect of metal binding at a site is transduced to other parts of the protein and could play a key role in protein-ligand and/or protein-protein interaction.

Dual affinity labeling of the active site of human lysozyme with an N-acetyllactosamine derivative: first ligand assisted recognition of the second ligand

Among the three kinds of the 2',3'-epoxypropyl beta-glycoside of disaccharides (GlcNAc-beta1,4-GlcNAc, Gal-beta1,4-GlcNAc, and Man-beta1,4-GlcNAc), the derivative of N-acetyllactosamine (Gal-beta1,4-GlcNAc-Epo) caused the dual labeling of human lysozyme (HL) most efficiently. The labeled HL was crystallized and analyzed by X-ray diffraction methodology. The X-ray analysis located the two Gal-beta1,4-GlcNAc-Epo moieties inside the catalytic cleft of HL. The attachment sites were the side-chain carboxylate groups of the catalytic residues Glu35 and Asp53 in HL. The first Gal-beta1, 4-GlcNAc-Epo moiety occupied virtually the same position as observed in the HL labeled with single Gal-beta1,4-GlcNAc-Epo molecule. The second Gal-beta1,4-GlcNAc-Epo moiety was recognized via the carbohydrate-carbohydrate interaction with the first Gal-beta1, 4-GlcNAc-Epo moiety in addition to the protein-carbohydrate interaction with the "right-side" catalytic cleft of HL through a number of hydrogen bonds including water-mediated ones as well as many van der Waals contacts. The two N-acetylglucosamine residues stacked with each other, while the two rings of galactose residues approximately shared the same plane. The dual labeling with two Gal-beta1,4-GlcNAc-Epo molecules was supposed to have occurred sequentially, which was accompanied with the alteration to the pKa of Glu35 derived from the esterification of Asp53 in the first labeling. Both asymmetric carbons in the connection parts between HL and N-acetyllactosamine moieties showed the same stereoconfiguration derived from the reaction with (2'R) stereoisomer concerning the epoxide group in the labeling reagent. The results demonstrated that the HL labeled with single Gal-beta1,4-GlcNAc-Epo was functional as a novel N-acetyllactosamine-binding protein, and the second labeling was performed by way of the first-ligand assisted recognition of the second ligand.

Structural and functional effect of Trp-62-->Gly and Asp-101-->Gly substitutions on substrate-binding modes of mutant hen egg-white lysozymes

In order to clarify the structural role of subsite B of hen egg-white lysozyme in hydrolytic activity towards a carbohydrate substrate, we analysed the structures of Trp-62-->Gly and Asp-101-->Gly mutant hen lysozymes, which have no side chain at positions 62 or 101, complexed with a substrate analogue, (N-acetyl-d-glucosamine)3 [(GlcNAc)3], using X-ray crystallography. The overall protein structures in the mutant lysozyme complexes were almost identical to those in the wild type. In the crystals of all the mutant complexes, the (GlcNAc)3 molecule, which is an inhibitor of wild-type lysozyme, had no inhibitory effect, but was hydrolysed as a substrate. One of the products, (GlcNAc)2, the reducing end of which is an alpha-anomer, was bound in an unproductive binding mode, protruding from the active-site cleft, and was able to act as an inhibitor. Hydrolysis of the synthetic substrate by the mutants occurred in a beta-anomer-retaining manner, and so the alpha-anomer product was converted from the beta-anomer product. Thus the interactions of Asp-101 and Trp-62 in subsite B are not essential for the catalytic mechanism, but co-operatively enhance the affinity of the substrate in the productive binding mode, other than the inhibitor in the unproductive mode.

Structural analysis of mutant hen egg-white lysozyme preferring a minor binding mode

Trp62 in hen egg-white lysozyme has general features observed in protein-carbohydrate interactions, a stacking interaction toward nonpolar surface of the substrate sugar residue B and a hydrogen bonding network with the residue C. Our previous report (I. Kumagai, K. Maenaka, F. Sunada, S. Takeda, K. Miura, Eur. J. Biochem. 212 (1993) 151-156.) showed that the substitution of Trp62 changed the substrate binding modes; especially, the Trp62His mutant exhibited the drastic change of the binding mode and preferred to a minor binding mode of the wild-type enzyme. In order to clarify the relationship between functional and structural changes of the Trp62His mutant, we analyzed the structure of the Trp62His mutant hen lysozyme complexed with the substrate analogue, (GlcNAc)3, by X-ray crystallography. The overall protein structure in the mutant lysozyme complex was almost identical to that in the wild-type. His62 shared almost the same plane as the indole ring of Trp62 of the wild-type. Although the (GlcNAc)3 molecule which is an inhibitor against the wild-type lysozyme was cocrystallized, the Trp62His mutant did not put it in the sites A-B-C but hydrolyzed it as a substrate. One of the products, (GlcNAc)2, whose reducing end is alpha-anomer, was bound in another binding mode sticking out from the active-site cleft. The hydrolytic activity against the synthetic substrate showed that the mutant was a beta-anomer retaining enzyme, so the alpha-anomer product was converted from the beta-anomer product. Therefore, the Trp62His mutant showed the remarkable change of the substrate binding modes not by alteration of the catalytic system but possibly by subtle rearrangement of general features of protein-carbohydrate interactions between His62 and the sugar residues B and C.

Molecular divergence of lysozymes and alpha-lactalbumin

The vast number of proteins that sustain the currently living organisms have been generated from a relatively small number of ancestral genes that has involved a variety of processes. Lysozyme is an ancient protein whose origin goes back an estimated 400 to 600 million years. This protein was originally a bacteriolytic defensive agent and has been adapted to serve a digestive function on at least two occasions, separated by nearly 40 million years. The origins of the related goose type and T4 phage lysozyme that are distinct from the more common C type are obscure. They share no discernable amino acid sequence identity and yet they possess common secondary and tertiary structures. Lysozyme C gene also gave rise, after gene duplication 300 to 400 million years ago, to a gene that currently codes for alpha-lactalbumin, a protein expressed only in the lactating mammary gland of all but a few species of mammals. It is required for the synthesis of lactose, the sugar secreted in milk. alpha-Lactalbumin shares only 40% identity in amino acid sequence with lysozyme C, but it has a closer spatial structure and gene organization. Although structurally similar, functionally they are quite distinct. Specific amino acid substitutions in alpha-lactalbumin account for the loss of the enzyme activity of lysozyme and the acquisition of the features necessary for its role in lactose synthesis. Evolutionary implications are as yet unclear but are being unraveled in many laboratories.

Molecular modelling studies of lysozyme catalysed hydrolysis of synthetic substrates

Kinetic data for the lysozyme catalysed hydrolysis of aryl chitooligosides were surveyed. Both electron-donating and electron-withdrawing substituents on the departing aryl aglycones enhance the rate of hydrolyses. The parallel pH-rate profiles implicate that identical catalytic residues are involved in the hydrolytic fission of the glycosyl-aryloxy bond of these two groups of synthetic substrates. Molecular modelling studies of lysozyme complexes with aryl diN-acetyl chitobiosides and their intermediates were performed. The two synthetic substrates bearing aryl aglycones with opposite electronic effects bind to the active site of lysozyme in different conformations. Based on the energetic and geometric considerations, the oxocarbonium ion whose pyranose ring D in a sofa conformation is the most plausible reaction intermediate for the lysozyme catalysed hydrolysis of the synthetic substrates. The modelling study also suggests that considerable conformational changes of both the lysozyme binding site and the chitobiosyl group accompany the formation of the glycosyl enzyme intermediate. In particular, the chitobiosyl group undergoes a dislocation of the pyranose ring C from the subsite C and a constraint of the pyranose ring D to form a boat conformer.

Origin of carbohydrate recognition specificity of human lysozyme revealed by affinity labeling

In order to reveal the origin of carbohydrate recognition specificity of human lysozyme by clarifying the difference in the binding mode of ligands in the active site, the inactivation of human lysozyme by 2',3'-epoxypropyl beta-glycoside derivatives of the disaccharides, N,N'-diacetylchitobiose [GlcNAc-beta-(1-->4)-GlcNAc] and N-acetyllactosamine [Gal-beta-(1-->4)-GlcNAc], was investigated and the three-dimensional structures of the affinity-labeled enzymes were determined by X-ray crystallography at 1.7 A resolution. Under the conditions comprising 2.0 x 10(-3) M labeling reagent and 1.0 x 10(-5) M human lysozyme at pH 5.4, 37 degrees C, the reaction time required to reduce the lytic activity against Micrococcus luteus cells to 50% of its initial activity was lengthened by 3.7 times through the substitution of the nonreducing end sugar residue, GlcNAc to Gal. The refined structure of human lysozyme labeled by 2',3'-epoxypropyl beta-glycoside derivatives of N,N'-diacetylchitobiose (HL/NAG-NAG-EPO complex) indicated that the interaction mode of the N,N'-diacetylchitobiose moiety in substites B and C in this study was essentially the same as in the case of the complex of human lysozyme with the free ligand. On the other hand, the hydrogen-bonding pattern and the stacking interaction at subsite B were remarkably different between the HL/NAG-NAG-EPO complex and human lysozyme labeled by the 2',3'-epoxypropyl beta-glycoside of N-acetyllactosamine (HL/GAL-NAG-EPO complex). The reduced number of possible hydrogen bonds as well as the less favorable stacking between the side chain of Tyr63 in human lysozyme and the galactose residue in the HL/GAL-NAG-EPO complex reasonably explained the less efficient ability of the 2',3'-epoxypropyl beta-glycoside of N-acetyllactosamine as compared to that of N,N'-diacetylchitobiose as an affinity labeling reagent toward human lysozyme.

Functional site in alpha-lactalbumin encompasses a region corresponding to a subsite in lysozyme and parts of two adjacent flexible substructures

Aromatic cluster 1 of alpha-lactalbumin (LA), a substructure adjacent to the cleft, is important for its interaction with galactosyltransferase (GT) and effects on glucose binding in the lactose synthase complex [Grobler, J. A., Wang, M., Pike, A. K., & Brew, K. (1994) J. Biol. Chem. 269, 5106-5114]. The full extent of the functional region in LA has been probed by mutagenesis of residues that are near aromatic cluster 1 or within the cleft that corresponds to the active site in the homologous type c lysozymes. The conserved residues Val42, Gln54, and Ile59, which correspond to residues of lysozyme that act in substrate binding in subsites C to E, together with residues adjacent to aromatic cluster 1, were found to be not required for activity. In contrast, replacing Leu110, a component of the region corresponding to lysozyme subsite F, with His or Glu greatly reduces the affinity of LA for GT while the introduction of Arg lowers the synergism of LA and glucose binding to GT and also reduces the affinity of LA for GT. Substitutions for Ala106, which is adjacent to Leu110 in the structure, also perturb activity. The region of the cleft corresponding to subsite F is important for function in LA as well as in lysozyme since other components of this subsite, His32 and Phe31, are also crucial for LA activity. The qualitatively different effects of various substitutions for Leu110 may be mediated by their influence on His32 or by changes in the structure of the lactose synthase complex.