Bacteriophage T4 Lysozyme Research Articles

Correlation of Solvent Interaction Analysis Signatures with Thermodynamic Properties and In Silico Calculations of the Structural Effects of Point Mutations in Two Proteins.

The partition behavior of single and double-point mutants of bacteriophage T4 lysozyme (T4 lysozyme) and staphylococcal nuclease A was examined in different aqueous two-phase systems (ATPSs) and studied by Solvent Interaction Analysis (SIA). Additionally, the solvent accessible surface area (SASA) of modeled mutants of both proteins was calculated. The in silico calculations and the in vitro analyses of the staphylococcal nuclease and T4 lysozyme mutants correlate, indicating that the partition analysis in ATPSs provides a valid descriptor (SIA signature) covering various protein features, such as structure, structural dynamics, and conformational stability.

Local Xenon-Protein Interaction Produces Global Conformational Change and Allosteric Inhibition in Lysozyme.

Noble gases have well-established biological effects, yet their molecular mechanisms remain poorly understood. Here, we investigated, both experimentally and computationally, the molecular modes of xenon (Xe) action in bacteriophage T4 lysozyme (T4L). By combining indirect gassing methods with a colorimetric lysozyme activity assay, a reversible, Xe-specific (20 ± 3)% inhibition effect was observed. Accelerated molecular dynamic simulations revealed that Xe exerts allosteric inhibition on the protein by expanding a C-terminal hydrophobic cavity. Xe-induced cavity expansion results in global conformational changes, with long-range transduction distorting the active site where peptidoglycan binds. Interestingly, the peptide substrate binding site that enables lysozyme specificity does not change conformation. Two T4L mutants designed to reshape the C-terminal Xe cavity established a correlation between cavity expansion and enzyme inhibition. This work also highlights the use of Xe flooding simulations to identify new cryptic binding pockets. These results enrich our understanding of Xe-protein interactions at the molecular level and inspire further biochemical investigations with noble gases.

An improved lysozyme substrate preparation method for the biochemistry kinetics laboratory

The mutation and growth of Bacteriophage T4 Lysozyme (T4L) is easily accomplished in an undergraduate biochemistry laboratory and allows students to visualize how mutations affect enzyme stability and kinetics. However, full Michaelis-Menten parameters have been difficult to observe due to lack of a readily available, small molecule substrate. A FITC-tagged peptidoglycan has been used previously, but it is not commercially available and its preparation method was not appropriate for the undergraduate laboratory. Here, we describe a cost-effective, non-toxic, and reproducible means of generating FITC-tagged peptidoglycan fragments. Fluorescence assays monitoring substrate digestion using HEWL and T4L produced signals sufficient for reliable analysis of lysozyme kinetics.

Exploring proten’s conformational space by using encoding layer supervised auto-encoder

Protein function is related to its structure and dynamic change. Molecular dynamics simulation is an important tool for studying protein dynamics by exploring its conformational space, however, conformational sampling is a nontrivial issue, because of the risk of missing key details during sampling. In recent years, deep learning methods, such as auto-encoder, can couple with MD to explore conformational space of protein. After being trained with the MD trajectories, auto-encoder can generate new conformations quickly by inputting random numbers in low dimension space. However, some problems still exist, such as requirements for the quality of the training set, the limitation of explorable area and the undefined sampling direction. In this work, we build a supervised auto-encoder, in which some reaction coordinates are used to guide conformational exploration along certain directions. We also try to expand the explorable area by training through the data generated by the model. Two multi-domain proteins, bacteriophage T4 lysozyme and adenylate kinase, are used to illustrate the method. In the case of the training set consisting of only under-sampled simulated trajectories, the supervised auto-encoder can still explore along the given reaction coordinates. The explored conformational space can cover all the experimental structures of the proteins and be extended to regions far from the training sets. Having been verified by molecular dynamics and secondary structure calculations, most of the conformations explored are found to be plausible. The supervised auto-encoder provides a way to efficiently expand the conformational space of a protein with limited computational resources, although some suitable reaction coordinates are required. By integrating appropriate reaction coordinates or experimental data, the supervised auto-encoder may serve as an efficient tool for exploring conformational space of proteins.

Automated Path Searching Reveals the Mechanism of Hydrolysis Enhancement by T4 Lysozyme Mutants

Bacteriophage T4 lysozyme (T4L) is a glycosidase that is widely applied as a natural antimicrobial agent in the food industry. Due to its wide applications and small size, T4L has been regarded as a model system for understanding protein dynamics and for large-scale protein engineering. Through structural insights from the single conformation of T4L, a series of mutations (L99A,G113A,R119P) have been introduced, which have successfully raised the fractional population of its only hydrolysis-competent excited state to 96%. However, the actual impact of these substitutions on its dynamics remains unclear, largely due to the lack of highly efficient sampling algorithms. Here, using our recently developed travelling-salesman-based automated path searching (TAPS), we located the minimum-free-energy path (MFEP) for the transition of three T4L mutants from their ground states to their excited states. All three mutants share a three-step transition: the flipping of F114, the rearrangement of α0/α1 helices, and final refinement. Remarkably, the MFEP revealed that the effects of the mutations are drastically beyond the expectations of their original design: (a) the G113A substitution not only enhances helicity but also fills the hydrophobic Cavity I and reduces the free energy barrier for flipping F114; (b) R119P barely changes the stability of the ground state but stabilizes the excited state through rarely reported polar contacts S117OG:N132ND2, E11OE1:R145NH1, and E11OE2:Q105NE2; (c) the residue W138 flips into Cavity I and further stabilizes the excited state for the triple mutant L99A,G113A,R119P. These novel insights that were unexpected in the original mutant design indicated the necessity of incorporating path searching into the workflow of rational protein engineering.

Creation of Cross-Linked Crystals With Intermolecular Disulfide Bonds Connecting Symmetry-Related Molecules Allows Retention of Tertiary Structure in Different Solvent Conditions.

Protein crystals are generally fragile and sensitive to subtle changes such as pH, ionic strength, and/or temperature in their crystallization mother liquor. Here, using T4 phage lysozyme as a model protein, the three-dimensional rigidification of protein crystals was conducted by introducing disulfide cross-links between neighboring molecules in the crystal. The effect of cross-linking on the stability of the crystals was evaluated by microscopic observation and X-ray diffraction. When soaking the obtained cross-linked crystals into a precipitant-free solution, the crystals held their shape without dissolution and diffracted to approximately 1.1 Å resolution, comparable to that of the non-cross-linked crystals. Such cross-linked crystals maintained their diffraction even when immersed in other solutions with pH values from 4 to 10, indicating that the disulfide cross-linking made the packing contacts enforced and resulted in some mechanical strength in response to changes in the preservation conditions. Furthermore, the cross-linked crystals gained stability to permit soaking into solutions containing high concentrations of organic solvents. The results suggest the possibility of obtaining protein crystals for effective drug screening by introducing appropriate cross-linked disulfide bonds.

Choice of force fields and water models for sampling solution conformations of bacteriophage T4 lysozyme

A protein may exist as an ensemble of different conformations in solution, which cannot be represented by a single static structure. Molecular dynamics (MD) simulation has become a useful tool for sampling protein conformations in solution, but force fields and water models are important issues. This work presents a case study of the bacteriophage T4 lysozyme (T4L). We have found that MD simulations using a classic AMBER99SB force field and TIP4P water model cannot well describe hinge-bending domain motion of the wild-type T4L at the timescale of one microsecond. Other combinations, such as a residue-specific force field called RSFF2+ and a dispersion-corrected water model TIP4P-D, are able to sample reasonable solution conformations of T4L, which are in good agreement with experimental data. This primary study may provide candidates of force fields and water models for further investigating conformational transition of T4L.

Genetically Encoded Biosensor-Based Screening for Directed Bacteriophage T4 Lysozyme Evolution.

Lysozyme is widely used as a model protein in studies of structure–function relationships. Recently, lysozyme has gained attention for use in accelerating the degradation of secondary sludge, which mainly consists of bacteria. However, a high-throughput screening system for lysozyme engineering has not been reported. Here, we present a lysozyme screening system using a genetically encoded biosensor. We first cloned bacteriophage T4 lysozyme (T4L) into a plasmid under control of the araBAD promoter. The plasmid was expressed in Escherichia coli with no toxic effects on growth. Next, we observed that increased soluble T4L expression decreased the fluorescence produced by the genetic enzyme screening system. To investigate T4L evolution based on this finding, we generated a T4L random mutation library, which was screened using the genetic enzyme screening system. Finally, we identified two T4L variants showing 1.4-fold enhanced lytic activity compared to native T4L. To our knowledge, this is the first report describing the use of a genetically encoded biosensor to investigate bacteriophage T4L evolution. Our approach can be used to investigate the evolution of other lysozymes, which will expand the applications of lysozyme.

Mapping out the Degree of Freedom of Hosted Enzymes in Confined Spatial Environments

Summary The integration of enzymes with solid materials is crucial for promoting their industrialization. Understanding the enzyme behavior upon association within a confined space, though of fundamental importance for biocomposite development, remains a persistent, unresolved challenge. Here, we present a comprehensive elucidation of the spatial environment’s impact on hosted enzymes’ degree of freedom and, consequently, their accompanying reactivity. Site-directed spin labeling in combination with electron paramagnetic resonance spectroscopy allows the direct detection of host-guest interactions at atomic resolution, while the tailorable synthesis of covalent organic frameworks (COFs) enables an evaluation of factors affecting such interactions. Specifically, lysozyme is found to be more constrained and less active along with increasing hydrophilicity of the COFs. These results support the establishment of a connection between the hydrophilicity of the spatial environment and the resulting biocomposites’ reactivity, enabling the prediction of the performance of unknown biocomposites. This study provides a unique insight into the mechanistic pathways underpinning biocatalysis.

Evaluation of Feasibility of Using the Bacteriophage T4 Lysozyme to Improve the Hydrolysis and Biochemical Methane Potential of Secondary Sludge

Anaerobic digestion (AD) of secondary sludge is a rate-limiting step due to the bacterial cell wall. In this study, experiments were performed to characterize secondary sludges from three wastewater treatment plants (WWTPs), and to investigate the feasibility of using bacteriophage lysozymes to speed up AD by accelerating the degradation of bacterial cell walls. Protein was the main organic material (67.7% of volatile solids in the sludge). The bacteriophage T4 lysozyme (T4L) was tested for hydrolysis and biochemical methane potential. Variations in the volatile suspended solid (VSS) concentration and biogas production were monitored. The VSS reduction efficiencies by hydrolysis using T4L for 72 h increased and ranged from 17.8% to 26.4%. Biogas production using T4L treated sludges increased and biogas production was increased by as much as 82.4%. Biogas production rate also increased, and the average reaction rate coefficient of first-order kinetics was 0.56 ± 0.02/d, which was up to 47.5% higher compared to the untreated samples at the maximum. Alphaproteobacteria, Betaproteobacteria, Flavobacteriia, Gammaproteobacteria, and Sphingobacteriia were major microbial classes in all sludges. The interpretation of the microbial community structure indicated that T4L treatment is likely to increase the rate of cell wall digestion.

The pKa values of the catalytic residues in the retaining glycoside hydrolase T26H mutant of T4 lysozyme.

T4 phage lysozyme (T4L) is an enzyme that cleaves bacterial cell wall peptidoglycan. Remarkably, the single substitution of the active site Thr26 to a His (T26H) converts T4L from an inverting to a retaining glycoside hydrolase with transglycosylase activity. It has been proposed that T26H-T4L follows a double displacement mechanism with His26 serving as a nucleophile to form a covalent glycosyl-enzyme intermediate (Kuroki et al., PNAS 1999; 96:8949-8954). To gain further insights into this or alternative mechanisms, we used NMR spectroscopy to measure the acid dissociation constants (pKa values) and/or define the ionization states of the Asp, Glu, His, and Arg residues in the T4L mutant. Most notably, the pKa value of the putative nucleophile His26 is 6.8 ± 0.1, whereas that of the general acid Glu11 is 4.7 ± 0.1. If the proposed mechanism holds true, then T26H-T4L follows a reverse protonation pathway in which only a minor population of the free enzyme is in its catalytically competent ionization state with His26 deprotonated and Glu11 protonated. Our studies also confirm that all arginines in T26H-T4L, including the active site Arg145, are positively charged under neutral pH conditions. BRIEF STATEMENT: The replacement of a single amino acid changes T4 lysozyme from an inverting to a retaining glycoside hydrolase. Using NMR spectroscopy, we measured the pKa values of the ionizable residues in the active site of this mutant enzyme. Along with previously reported data, these results provide important constraints for understanding the catalytic mechanisms by which the wild-type and mutant form of T4 lysozyme cleave bacterial peptidoglycan.

Biomolecular Simulations of Halogen Bonds with a GROMOS Force Field.

Halogen bonds (XBs) are non-covalent interactions in which halogens (X), acting as electrophiles, interact with Lewis bases. XBs are able to mediate protein-ligand recognition and therefore play an important role in rational drug design. In this context, the development of molecular modeling tools that can tackle XBs is paramount. XBs are predominantly explained by the existence of a positive region on the electrostatic potential of X named the σ-hole. Typically, with molecular mechanics force fields, this region is modeled using a charged extra point (EP) linked to X along the R-X covalent bond axis. In this work, we developed the first EP-based strategy for GROMOS force fields (specifically GROMOS 54A7) using bacteriophage T4 lysozyme in complex with both iodobenzene and iodopentafluorobenzene as a prototype system. Several EP parametrization schemes were tested by adding a virtual interaction site to ligand topologies retrieved from the Automated Topology Builder (ATB) and Repository. Contrary to previous approaches using other force fields, our analysis is based on the capability of each parametrization scheme to sample XBs during MD simulations. Our results indicate that the implementation of an EP at a distance from iodine corresponding to Rmin provides a good qualitative description of XBs in MD simulations, supporting the compatibility of our approach with the GROMOS 54A7 force field.

Amino-Acid Network Clique Analysis of Protein Mutation Non-Additive Effects: A Case Study of Lysozme.

Optimizing amino-acid mutations in enzyme design has been a very challenging task in modern bio-industrial applications. It is well known that many successful designs often hinge on extensive correlations among mutations at different sites within the enzyme, however, the underpinning mechanism for these correlations is far from clear. Here, we present a topology-based model to quantitively characterize non-additive effects between mutations. The method is based on the molecular dynamic simulations and the amino-acid network clique analysis. It examines if the two mutation sites of a double-site mutation fall into to a 3-clique structure, and associates such topological property of mutational site spatial distribution with mutation additivity features. We analyzed 13 dual mutations of T4 phage lysozyme and found that the clique-based model successfully distinguishes highly correlated or non-additive double-site mutations from those additive ones whose component mutations have less correlation. We also applied the model to protein Eglin c whose structural topology is significantly different from that of T4 phage lysozyme, and found that the model can, to some extension, still identify non-additive mutations from additive ones. Our calculations showed that mutation non-additive effects may heavily depend on a structural topology relationship between mutation sites, which can be quantitatively determined using amino-acid network k-cliques. We also showed that double-site mutation correlations can be significantly altered by exerting a third mutation, indicating that more detailed physicochemical interactions should be considered along with the network clique-based model for better understanding of this elusive mutation-correlation principle.

Using disruptive insertional mutagenesis to identify the in situ structure-function landscape of the Shigella translocator protein IpaB.

Bacterial type III secretion systems (T3SS) are used to inject proteins into mammalian cells to subvert cellular functions. The Shigella T3SS apparatus (T3SA) is comprised of a basal body, cytoplasmic sorting platform and exposed needle with needle "tip complex" (TC). TC maturation occurs when the translocator protein IpaB is recruited to the needle tip where both IpaD and IpaB control secretion induction. IpaB insertion into the host membrane is the first step of translocon pore formation and secretion induction. We employed disruptive insertional mutagenesis, using bacteriophage T4 lysozyme (T4L), within predicted IpaB loops to show how topological features affect TC functions (secretion control, translocon formation and effector secretion). Insertions within the N-terminal half of IpaB were most likely to result in a loss of steady-state secretion control, however, all but the two that were not recognized by the T3SA retained nearly wild-type hemolysis (translocon formation) and invasiveness levels (effector secretion). In contrast, all but one insertion in the C-terminal half of IpaB maintained secretion control but were impaired for hemolysis and invasion. These nature of the data suggest the latter mutants are defective in a post-secretion event, most likely due to impaired interactions with the second translocator protein IpaC. Intriguingly, only two insertion mutants displayed readily detectable T4L on the bacterial surface. The data create a picture in which the makeup and structure of a functional T3SA TC is highly amenable to physical perturbation, indicating that the tertiary structure of IpaB within the TC is more plastic than previously realized.

Neutron crystallographic studies of T4 lysozyme at cryogenic temperature.

Bacteriophage T4 lysozyme (T4L) has been used as a paradigm for seminal biophysical studies on protein structure, dynamics, and stability. Approximately 700 mutants of this protein and their respective complexes have been characterized by X-ray crystallography; however, despite the high resolution diffraction limits attained in several studies, no hydrogen atoms were reported being visualized in the electron density maps. To address this, a 2.2 Å-resolution neutron data set was collected at 80 K from a crystal of perdeuterated T4L pseudo-wild type. We describe a near complete atomic structure of T4L, which includes the positions of 1737 hydrogen atoms determined by neutron crystallography. The cryogenic neutron model reveals explicit detail of the hydrogen bonding interactions in the protein, in addition to the protonation states of several important residues.

Neutron structure of the T26H mutant of T4 phage lysozyme provides insight into the catalytic activity of the mutant enzyme and how it differs from that of wild type.

T4 phage lysozyme is an inverting glycoside hydrolase that degrades the murein of bacterial cell walls by cleaving the β‐1,4‐glycosidic bond. The substitution of the catalytic Thr26 residue to a histidine converts the wild type from an inverting to a retaining enzyme, which implies that the original general acid Glu11 can also act as an acid/base catalyst in the hydrolysis. Here, we have determined the neutron structure of the perdeuterated T26H mutant to clarify the protonation states of Glu11 and the substituted His26, which are key in the retaining reaction. The 2.09‐Å resolution structure shows that the imidazole group of His26 is in its singly protonated form in the active site, suggesting that the deprotonated Nɛ2 atom of His26 can attack the anomeric carbon of bound substrate as a nucleophile. The carboxyl group of Glu11 is partially protonated and interacts with the unusual neutral state of the guanidine moiety of Arg145, as well as two heavy water molecules. Considering that one of the water‐binding sites has the potential to be occupied by a hydronium ion, the bulk solvent could be the source for the protonation of Glu11. The respective protonation states of Glu11 and His26 are consistent with the bond lengths determined by an unrestrained refinement of the high‐resolution X‐ray structure of T26H at 1.04‐Å resolution. The detail structural information, including the coordinates of the deuterium atoms in the active site, provides insight into the distinctively different catalytic activities of the mutant and wild type enzymes.

Metal ions-binding T4 lysozyme as an intramolecular protein purification tag compatible with X-ray crystallography.

Phage T4 lysozyme is a well folded and highly soluble protein that is widely used as an insertion tag to improve solubility and crystallization properties of poorly behaved recombinant proteins. It has been used in the fusion protein strategy to facilitate crystallization of various proteins including multiple G protein-coupled receptors, lipid kinases, or sterol binding proteins. Here, we present a structural and biochemical characterization of its novel, metal ions-binding mutant (mbT4L). We demonstrate that mbT4L can be used as a purification tag in the immobilized-metal affinity chromatography and that, in many respects, it is superior to the conventional hexahistidine tag. In addition, structural characterization of mbT4L suggests that mbT4L can be used as a purification tag compatible with X-ray crystallography.

PDB2Graph: A toolbox for identifying critical amino acids map in proteins based on graph theory

The integrative and cooperative nature of protein structure involves the assessment of topological and global features of constituent parts. Network concept takes complete advantage of both of these properties in the analysis concomitantly. High compatibility to structural concepts or physicochemical properties in addition to exploiting a remarkable simplification in the system has made network an ideal tool to explore biological systems. There are numerous examples in which different protein structural and functional characteristics have been clarified by the network approach. Here, we present an interactive and user-friendly Matlab-based toolbox, PDB2Graph, devoted to protein structure network construction, visualization, and analysis. Moreover, PDB2Graph is an appropriate tool for identifying critical nodes involved in protein structural robustness and function based on centrality indices. It maps critical amino acids in protein networks and can greatly aid structural biologists in selecting proper amino acid candidates for manipulating protein structures in a more reasonable and rational manner. To introduce the capability and efficiency of PDB2Graph in detail, the structural modification of Calmodulin through allosteric binding of Ca(2+) is considered. In addition, a mutational analysis for three well-identified model proteins including Phage T4 lysozyme, Barnase and Ribonuclease HI, was performed to inspect the influence of mutating important central residues on protein activity.

Analysis of the solution conformations of T4 lysozyme by paramagnetic NMR spectroscopy.

A large number of crystal structures of bacteriophage T4 lysozyme (T4-L) have shown that it contains two subdomains, which can arrange in a compact conformation (closed state) or, in mutants of T4-L, more extended structures (open state). In solution, wild-type T4-L displays only a single set of nuclear magnetic resonance (NMR) signals, masking any conformational heterogeneity. To probe the conformational space of T4-L, we generated a site-specific lanthanide binding site by attaching 4-mercaptomethyl dipicolinic acid via a disulfide bond to Cys44 in the triple-mutant C54T/C97A/S44C of T4-L and measured pseudocontact shifts (PCS) and magnetically induced residual dipolar couplings (RDC). The data indicate that, in solution and in the absence of substrate, the structure of T4-L is on average more open than suggested by the closed conformation of the crystal structure of wild-type T4-L. A slightly improved fit was obtained by assuming a population-weighted two-state model involving an even more open conformation and the closed state, but paramagnetic relaxation enhancements measured with Gd(3+) argue against such a conformational equilibrium. The fit could not be improved by including a third conformation picked from the hundreds of crystal structures available for T4-L mutants.

Functional heterogeneity as reflected by topological parameters in a classical protein molecular model: t4 phage lysozyme.

A systematic comparison with the Wild-Type (WT) of one-point mutants of bacteriophage T4 lysozyme was carried out using as difference markers the topological parameters of the protein contact networks corresponding to each crystallographic structure. The investigation concerned changes at the resolution level of single residue along the protein sequence. The results were correlated with (reported) changes in functional properties and (observed) changes in the information provided by the energy dissipation algorithm of the "Turbine" software simulation tool. The critical factor leading to significant difference among mutants and WT is in most cases associated to the sensitivity towards mutation of relatively short windows in the amino acidic sequence not necessarily contiguous to the active site.