Structure and dynamics of the active site of hen egg-white lysozyme from atomic resolution neutron crystallography

Hen egg-white lysozyme (HEWL) is a widely used model protein in crystallographic studies and its enzymatic mechanism has been extensively investigated for decades. Despite this, the interaction between the reaction intermediate and the catalytic Asp52, as well as the orientation of Asn44 and Asn46 side chains, remain ambiguous. Here, we report the crystal structures of perdeuterated HEWL and D2O buffer-exchanged HEWL from 0.91 and 1.1 Å resolution neutron diffraction data, respectively. These structures were obtained at room temperature and acidic pH, representing the active state of the enzyme. The unambiguous assignment of hydrogen positions based on the neutron scattering length density maps elucidates the roles of Asn44, Asn46, Asn59, and nearby water molecules in the stabilization of Asp52. Additionally, the identification of hydrogen positions reveals unique details of lysozyme's folding, hydrogen (H)/deuterium (D) exchange, and side chain disorder.

Hen egg white lysozyme encapsulated in ZIF-8 for performing promiscuous enzymatic Mannich reaction

Hen egg white lysozyme (HEWL) was exploited for the synthesis of β-amino carbonyl compounds through a direct and three-component Mannich reaction in aqueous, confirming high chemoselectivity toward imine. In order to further extend the applications of the enzyme, HEWL was encapsulated using a metal-organic framework (MOF). The reactivity, stereoselectivity, and reusability of the encapsulated enzyme were investigated. The reaction was significantly enhanced as compared to the non-encapsulated enzyme. A mutated version of the enzyme, containing Asp52Ala (D52A), lacking important catalytical residue, has lost the bacterial site activity against Micrococcus luteus (M. luteus) while the D52A variant displayed an increased rate of the Mannich reaction, indicating a different catalytical residue involved in the promiscuous reaction. Based on site-directed mutagenesis, molecular docking, and molecular dynamic studies, it was proposed that π-stacking, H-bond interactions, and the presence of water in the active site may play crucial roles in the mechanism of the reaction.

Identifying pseudoenzymes using functional annotation: pitfalls of common practice

Pseudoenzymes are proteins that are evolutionary related to enzymes but lack relevant catalytic activity. They are usually evolved from enzymatic ancestors that have lost their catalytic activities. The loss of catalytic function is one extreme amongst the other evolutionary changes that can occur to enzymes, like the changing of substrate specificity or the reaction catalysed. However, the loss of catalytic function events remain poorly characterised, except for some notable examples, like the pseudokinases. In this review, we aim to analyse current knowledge related to pseudoenzymes across a large number of enzymes families. This aims to be a review of the data available in biological databases, rather than a more traditional literature review. In particular, we use UniProtKB as the source for functional annotation and M-CSA (Mechanism and Catalytic Site Atlas) for information on the catalytic residues of enzymes. We show that explicit annotation of lack of activity is not exhaustive in UniProtKB and that a protocol using lack of catalytic annotation as an indication for lack of function can be an adequate alternative, after some corrections. After identifying pseudoenzymes related to enzymes in M-CSA, we were able to comment on their prevalence across enzyme families, and on the correlation between lack of catalytic function and the mutation of catalytic residues. These analyses challenge two common ideas in the emerging literature: that pseudoenzymes are ubiquitous across enzyme families and that mutations in the catalytic residues of enzyme homologues are always a good indication of lack of activity.

The binding of precipitant ions in the tetragonal crystals of hen egg white lysozyme

The bonds between lysozyme molecules and precipitant ions in single crystals grown with chlorides of several metals are analysed on the basis of crystal structure data. Crystals of tetragonal hen egg lysozyme (HEWL) were grown with chlorides of several alkali and transition metals (LiCl, NaCl, KCl, NiCl₂ and CuCl₂) as precipitants and the three-dimensional structures were determined at 1.35 Å resolution by X-ray diffraction method. The positions of metal and chloride ions attached to the protein were located, divided into three groups and analysed. Some of them, in accordance with the recently proposed and experimentally confirmed crystal growth model, provide connections in protein dimers and octamers that are precursor clusters in the crystallization lysozyme solution. The first group, including Cu⁺², Ni⁺² and Na⁺¹ cations, binds specifically to the protein molecule. The second group consists of metal and chloride ions bound inside the dimers and octamers. The third group of ions can participate in connections between the octamers that are suggested as building units during the crystal growth. The arrangement of chloride and metal ions associated with lysozyme molecule at all stages of the crystallization solution formation and crystal growth is discussed.Communicated by Ramaswamy H. Sarma.

X-ray crystallographic analysis of time-dependent binding of guanidine hydrochloride to HEWL: First steps during protein unfolding

Time-dependent binding of guanidine hydrochloride (GuHCl) to hen egg-white lysozyme (HEWL), and effects of this binding on the protein structure have been investigated by solving X-ray structures of crystalline complexes. The complexes have been prepared by soaking, for different periods of time, native lysozyme crystals in solutions containing 2.5M GuHCl. In the refined structures, the number of water molecules in the protein's first solvent shell has progressively decreased from 152 to 115, showing protein's preference for guanidinium over water. Guanidinium ions preferentially hydrogen bond with the backbone carbonyl oxygen atoms. In their van der Waals interactions, they do not show any preference for apolar residues. Guanidinium ions have replaced water molecules that form cages around exposed hydrophobic residues. Guanidinium binding has decreased the average length of water-water hydrogen bond by 0.1Å. The hydrogen bonds between main chain atoms have been weakened by GuHCl, and this may be the reason for increased potency of GuHCl compared to urea. Guanidinium binding destabilizes the β-domain by causing loss of hydrogen bonds involving Asn 59 side chain. Interestingly, this loss is almost identical to that observed in structures of amyloidogenic variants of human lysozyme. Compounds preventing this loss could be anti-amyloidogenic.

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.

Asp48 function in the hydrogen-bonding network involving Asp52 of hen egg-white lysozyme

To characterize the hydrogen-bonding network in lysozyme, we focused on the residue of Asp48 located at the active site in hen egg-white lysozyme. We constructed a mutant lysozyme (D48A) and analyzed using (GlcNAc)3 and chitin-affinity chromatography. The substrate binding of subsites D–F in D48A and the activity against (GlcNAc)5 were decreased. The parameters of substrate binding and rate constants obtained from computer simulations confirmed these changes. In the crystal structure, (GlcNAc)4 was located at the same position as wildtype. However, the side chains of Arg45 and Thr47 at subsites E–F were moved by the replacement. Further, the loss of the hydrogen bond between Asp48 and Ser50 changed the hydrogen-bonding network, and this resulted in an alteration of the side chain of Asn59. This result suggests that the hydrogen-bonding network plays a crucial in the function of Asp52 and of transglycosylation at subsites E–F.

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.

Evolution of the genetic code by incorporation of amino acids that improved or changed protein function

Fifty years have passed since the genetic code was deciphered, but how the genetic code came into being has not been satisfactorily addressed. It is now widely accepted that the earliest genetic code did not encode all 20 amino acids found in the universal genetic code as some amino acids have complex biosynthetic pathways and likely were not available from the environment. Therefore, the genetic code evolved as pathways for synthesis of new amino acids became available. One hypothesis proposes that early in the evolution of the genetic code four amino acids-valine, alanine, aspartic acid, and glycine-were coded by GNC codons (N = any base) with the remaining codons being nonsense codons. The other sixteen amino acids were subsequently added to the genetic code by changing nonsense codons into sense codons for these amino acids. Improvement in protein function is presumed to be the driving force behind the evolution of the code, but how improved function was achieved by adding amino acids has not been examined. Based on an analysis of amino acid function in proteins, an evolutionary mechanism for expansion of the genetic code is described in which individual coded amino acids were replaced by new amino acids that used nonsense codons differing by one base change from the sense codons previously used. The improved or altered protein function afforded by the changes in amino acid function provided the selective advantage underlying the expansion of the genetic code. Analysis of amino acid properties and functions explains why amino acids are found in their respective positions in the genetic code.

Binding structures of tri-N-acetyl-β-glucosamine in hen egg white lysozyme using molecular dynamics with a polarizable force field

Lysozyme is a well-studied enzyme that hydrolyzes the β-(1,4)-glycosidic linkage of N-acetyl-β-glucosamine (NAG)(n) oligomers. The active site of hen egg-white lysozyme (HEWL) is believed to consist of six subsites, A-F that can accommodate six sugar residues. We present studies exploring the use of polarizable force fields in conjunction with all-atom molecular dynamics (MD) simulations to analyze binding structures of complexes of lysozyme and NAG trisaccharide, (NAG)(3). MD trajectories are applied to analyze structures and conformation of the complex as well as protein-ligand interactions, including the hydrogen-bonding network in the binding pocket. Two binding modes (ABC and BCD) of (NAG)(3) are investigated independently based on a fixed-charge model and a polarizable model. We also apply molecular mechanics with generalized born and surface area (MM-GBSA) methods based on MD using both nonpolarizable and polarizable force fields to compute binding free energies. We also study the correlation between root-mean-squared deviation and binding free energies of the wildtype and W62Y mutant; we find that for this prototypical system, approaches using the MD trajectories coupled with implicit solvent models are equivalent for polarizable and fixed-charge models.

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.

Solvation properties of N-acetyl-β-glucosamine: molecular dynamics study incorporating electrostatic polarization

N-Acetyl-β-glucosamine (NAG) is an important moiety of glycoproteins and is involved in many biological functions. However, conformational and dynamical properties of NAG molecules in aqueous solution, the most common biological environment, remain ambiguous due to limitations of experimental methods. Increasing efforts are made to probe structural properties of NAG and NAG-containing macromolecules, like peptidoglycans and polymeric chitin, at the atomic level using molecular dynamics simulations. In this work, we develop a polarizable carbohydrate force field for NAG and contrast simulation results of various properties using this novel force field and an analogous nonpolarizable (fixed charge) model. Aqueous solutions of NAG and its oligomers are investigated; we explore conformational properties (rotatable bond geometry), electrostatic properties (dipole moment distribution), dynamical properties (self-diffusion coefficient), hydrogen bonding (water bridge structure and dynamics), and free energy of hydration. The fixed-charge carbohydrate force field exhibits deviations from the gas phase relative rotation energy of exocyclic hydroxymethyl side chain and of chair/boat ring distortion. The polarizable force field predicts conformational properties in agreement with corresponding first-principles results. NAG-water hydrogen bonding pattern is studied through radial distribution functions (RDFs) and correlation functions. Intermolecular hydrogen bonding between solute and solvent is found to stabilize NAG solution structures while intramolecular hydrogen bonds define glycosidic linkage geometry of NAG oligomers. The electrostatic component of hydration free energy is highly dependent on force field atomic partial charges, influencing a more favorable free energy of hydration in the fixed-charge model compared to the polarizable model.