Mutant Lysozymes 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.

Mechanistic and biophysical insight into the inhibitory and disaggregase role of antibiotic moxifloxacin on human lysozyme amyloid formation

Lysozyme amyloidosis is a systemic non-neuropathic disease caused by the accumulation of amyloids of mutant lysozyme. Presently, therapeutic interventions targeting lysozyme amyloidosis, remain elusive with only therapy available for lysozyme amyloidosis being supportive management. In this work, we examined the effects of moxifloxacin, a synthetic fluoroquinolone antibiotic on the amyloid formation of human lysozyme. The ability of moxifloxacin to interfere with lysozyme amyloid aggregation was examined using various biophysical methods like Rayleigh light scattering, Thioflavin T fluorescence assay, transmission electron microscopy and docking method. The reduction in scattering and ThT fluorescence along with extended lag phase in presence of moxifloxacin, suggest that the antibiotic inhibits and impedes the lysozyme fibrillation in concentration dependent manner. From ANS experiment, we deduce that moxifloxacin is able to decrease the hydrophobicity of the protein molecule thereby preventing aggregation. Our CD and DLS results show that moxifloxacin stabilizes the protein in its native monomeric structure, thus also showing retention of lytic activity upto 69% and inhibition of cytotoxicity at highest concentration of moxifloxacin. The molecular docking showed that moxifloxacin forms a stable complex of −7.6 kcal/mol binding energy and binds to the aggregation prone region of lysozyme thereby stabilising it and preventing aggregation. Moxifloxacin also showed disaggregase potential by disrupting fibrils and decreasing the β-sheet content of the fibrils. Our current study, thus highlight the anti-amyloid and disaggregase property of an antibiotic moxifloxacin and hence sheds light on the future of antibiotics against protein aggregation, a hallmark event in many neurodegenerative diseases.

Ligand Gaussian Accelerated Molecular Dynamics 2 (LiGaMD2): Improved Calculations of Ligand Binding Thermodynamics and Kinetics with Closed Protein Pocket.

Ligand binding thermodynamics and kinetics are critical parameters for drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics from molecular simulations due to limited simulation timescales. Protein dynamics, especially in the ligand binding pocket, often plays an important role in ligand binding. Based on our previously developed Ligand Gaussian accelerated molecular dynamics (LiGaMD), here we present LiGaMD2 in which a selective boost potential was applied to both the ligand and protein residues in the binding pocket to improve sampling of ligand binding and dissociation. To validate the performance of LiGaMD2, the T4 lysozyme (T4L) mutants with open and closed pockets bound by different ligands were chosen as model systems. LiGaMD2 could efficiently capture repetitive ligand dissociation and binding within microsecond simulations of all T4L systems. The obtained ligand binding kinetic rates and free energies agreed well with available experimental values and previous modeling results. Therefore, LiGaMD2 provides an improved approach to sample opening of closed protein pockets for ligand dissociation and binding, thereby allowing for efficient calculations of ligand binding thermodynamics and kinetics.

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.

Symmetry-Adapted Restraints for Binding Free Energy Calculations.

Binding free energy calculations rely critically on a precise definition of the bound state and well-designed ligand restraints to ensure that binding free energy calculations converge rapidly and yield estimates of well-defined thermodynamic quantities. The distance-to-bound-configuration (DBC) is a single variable that can precisely delineate the bound state of a ligand including translational, rotational and conformational degrees of freedom and has been successfully used to capture binding modes with complex geometries. DBC is defined as the root-mean-square deviation (RMSD) of ligand coordinates in the frame of reference of the binding site. In the special case where the ligand features symmetry-equivalent atoms, a standard RMSD arbitrarily distinguishes equivalent poses, mixing equivalent and nonequivalent degrees of freedom, and preventing the precise delineation of the bound state ensemble, which negates the benefits of defining a flat-bottom binding restraint. To remedy this, we introduce a symmetry-adapted DBC coordinate where the RMSD is minimized over permutations of equivalent ligand atoms. This coordinate is implemented in a portable software library, the Collective Variables Module. We tested the approach by computing the absolute binding free energy of benzene to the engineered site of a mutant lysozyme (L99A/M102H) using alchemical free energy perturbation. We found that the symmetry-adapted restraint leads to well-behaved convergence of both the decoupling free energy in the binding site and the restrained free energy in the gas phase, recovering the affinity computed using a classic center-of-mass restraint. Thus, symmetry-adapted DBC seamlessly generalizes the benefits of DBC restraints to the case of symmetric ligands. The underlying symmetric RMSD coordinate can also be used for analyzing or biasing simulations in other contexts than affinity predictions.

Rationally designed protein cross-linked hydrogel for bone regeneration via synergistic release of magnesium and zinc ions

The development of recombinant protein cross-linked injectable hydrogels with good mechanical strength and effective drug loading capacity for bone regeneration is extremely attractive and rarely reported. Here, we report the fabrication of a smart hydrogel delivery system by incorporating a rationally designed T4 lysozyme mutant (T4M) to mediate the localized delivery and synergistic release of Mg2+ and Zn2+ for bone repair. Apart from its intrinsic antibacterial properties, T4M bears abundant free amine groups on its surface to function as effective covalent crosslinkers to strengthen the hydrogel network as well as exhibits specific binding affinity to multivalent cations such as Zn2+. Moreover, the integrin receptor-binding Arg-Gly-Asp (RGD) sequence was introduced onto the C-terminus of T4 lysozyme to improve its cellular affinity and further facilitate rapid tissue regeneration. The final composite hydrogel displays excellent injectability, improved mechanical properties, antibacterial activity, and unique bioactivities. The effective loading of Mg2+/Zn2+ in the hydrogels could mediate the sequential and sustained release of Mg2+ and Zn2+, thereby resulting in synergistic enhancement on bone regeneration through modulation of the MAPK signaling pathway. We believe that the strategy proposed in this paper opens up a new route for developing protein cross-linked smart delivery systems for tissue regeneration.

Comprehensive Characterization of Ligand Unbinding Mechanisms and Kinetics for T4 Lysozyme Mutants

Comprehensive Characterization of Ligand Unbinding Mechanisms and Kinetics for T4 Lysozyme Mutants

Ligand unbinding mechanisms and kinetics for T4 lysozyme mutants from τRAMD simulations.

The protein-ligand residence time, τ, influences molecular function in biological networks and has been recognized as an important determinant of drug efficacy. To predict τ, computational methods must overcome the problem that τ often exceeds the timescales accessible to conventional molecular dynamics (MD) simulation. Here, we apply the τ-Random Acceleration Molecular Dynamics (τRAMD) method to a set of kinetically characterized complexes of T4 lysozyme mutants with small, engineered binding cavities. τRAMD yields relative ligand dissociation rates in good accordance with experiments across this diverse set of complexes that differ with regard to measurement temperature, ligand identity, protein mutation and binding cavity. τRAMD thereby allows a comprehensive characterization of the ligand egress routes and determinants of τ. Although ligand dissociation by multiple egress routes is observed, we find that egress via the predominant route determines the value of τ. We also find that the presence of a greater number of metastable states along egress pathways leads to slower protein-ligand dissociation. These physical insights could be exploited in the rational optimization of the kinetic properties of drug candidates.

Effect of glyoxal and 1-methylisatin on stress-induced fibrillation of Hen Egg White Lysozyme: Insight into the anti-amyloidogenic property of the compounds with possible therapeutic implications

Systemic amyloidosis is a hereditary disorder that mostly arises as a result of specific point mutations to the wild type gene of lysozyme, forming mutant lysozyme variants leading to aggregation of the protein. The small monomeric protein Hen Egg White Lysozyme (HEWL) is a structural homolog of Human Lysozyme and is widely used as a model protein to investigate protein aggregation. In the present study, we have investigated the effect of 1-methylisatin, an indole derivative and glyoxal, a reactive dicarbonyl compound, on stress-induced aggregation of HEWL. Interaction of the compounds with HEWL induced changes in structure and surface hydrophobicity of the protein as evident from CD spectroscopy, tryptophan fluorescence and ANS binding studies. Additional experiments (Thioflavin T fluorescence, AFM imaging and DLS studies) demonstrate that stress induces amyloid-like fibrillation of HEWL, however, prior modification of the protein with glyoxal or 1-methylisatin significantly reduces its susceptibility to aggregation. High resolution mass spectrometric analysis indicated that 1-methylisatin primarily complexes with the protein in the form of a dimer. On the other hand, glyoxal-mediated modification of the protein induces formation of glycated adducts (carboxymethyllysine, hydroimidazolone). The results highlight possible clinical implications of the compounds in treatment of systemic amyloidosis and protein conformational disorder.

Discovering Protein Conformational Flexibility through Artificial-Intelligence-Aided Molecular Dynamics.

Proteins sample a variety of conformations distinct from their crystal structure. These structures, their propensities, and the pathways for moving between them contain an enormous amount of information about protein function that is hidden from a purely structural perspective. Molecular dynamics simulations can uncover these alternative conformations but often at a prohibitively high computational cost. Here we apply our recent statistical mechanics and artificial intelligence-based molecular dynamics framework for enhanced sampling of protein loops. We exemplify the approach through the study of three mutants of the classical test-piece protein T4 lysozyme. We are able to correctly rank these according to the stability of their excited state. By analyzing reaction coordinates, we also obtain crucial insight into why these specific perturbations in sequence space lead to tremendous variations in conformational flexibility. Our framework thus allows an accurate comparison of loop conformation populations with minimal prior human bias and should be directly applicable to a range of macromolecules in biology, chemistry, and beyond.

The effects of mutation and modification on the structure and stability of human lysozyme: A molecular link between carbamylation and atherosclerosis

Amino acid mutations in some proteins such as lysozyme lead to genetically disorder variants and adverse pathogenic consequences. Recently, amino acid modifications were known as a risk factor in many related diseases such as uremia and atherosclerosis, showing the importance of these surface-structure changes. Although the structural consequences of the hereditary proteins have been examined extensively, such effects for the protein modifications are known to a lesser extent. One drawback in the examination of protein modifications is hardness in experimental detection of modifications by techniques such as NMR and crystallography. Molecular modeling and simulation can help to understand such phenomena at the molecular levels. It is more rational that the effects of both mutation and modification can be compared in a single protein model. Here, molecular dynamics simulation is used to compare the effects of a disease-related carbamylation modification and an amyloidogenic mutation (D67H) in human lysozyme as a model protein. The results show that the carbamylation adversely effects on the tertiary structure, leading to the similar unfolding pathway to the hereditary amyloidogenic form. The carbamylation leads to the instability of the overall protein conformation, especially on the β-domain, which is a characteristic of hereditary amyloidosis in human lysozymes. The aggregation behaviors of both modified and mutant lysozyme were examined by molecular docking calculations. The results showed that the partially unfolded lysozyme might form tight protein aggregates upon carbamylation similar to the amyloidogenic variant. Both single and all-residues carbamylations impose serious conformational changes to the tertiary structure of lysozyme. It was obtained that carbamylation of lysozyme strongly effects on the stability of N-terminal β-sheet, which can produce a highly unstable conformation. The results of this study not only show the adverse structural consequences of a disease-associated post-translational modification, but it also may be very helpful to understand the molecular basis for many carbamylation-related diseases such as atherosclerosis in ESRD patients. The results show that non-native post-translational modifications may be as structurally important as hereditary mutations.

Su1904 – Selectin-Inhibition by Doubly Glycosylated Lysozyme Mutants Bearing Long Lactosamine Chains and Terminal Slex Cognated Structures (LYSII/IVT-FUCTVI)

Su1904 – Selectin-Inhibition by Doubly Glycosylated Lysozyme Mutants Bearing Long Lactosamine Chains and Terminal Slex Cognated Structures (LYSII/IVT-FUCTVI)

A protein interaction free energy model based on amino acid residue contributions: Assessment of point mutation stability of T4 lysozyme.

Here we present a model to estimate the interaction free energy contribution of each amino acid residue of a given protein. Protein interaction energy is described in terms of per-residue interaction factors, μ. Multibody interactions are implicitly captured in μ through the combination of amino acid terms (γ) guided by local conformation indices (σ). The model enables construction of an interaction factor heat map for a protein in a given fold, allows prima facie assessment of the degree of residue-residue interaction, and facilitates a qualitative and quantitative evaluation of protein association properties. The model was used to compute thermal stability of T4 bacteriophage lysozyme mutants across seven sites. Qualitative assessment of mutational effects provides a straightforward rationale regarding whether a particular site primarily perturbs native or non-native states, or both. The presented model was found to be in good agreement with experimental mutational data (R 2 = 0.73) and suggests an approach by which to convert structure space into energy space.

A methyl 1H double quantum CPMG experiment to study protein conformational exchange.

Protein conformational changes play crucial roles in enabling function. The Carr-Purcell-Meiboom-Gill (CPMG) experiment forms the basis for studying such dynamics when they involve the interconversion between highly populated and sparsely formed states, the latter having lifetimes ranging from ~ 0.5 to ~ 5ms. Among the suite of experiments that have been developed are those that exploit methyl group probes by recording methyl 1H single quantum (Tugarinov and Kay in J Am Chem Soc 129:9514-9521, 2007) and triple quantum (Yuwen et al. in Angew Chem Int Ed Engl 55:11490-11494, 2016) relaxation dispersion profiles. Here we build upon these by developing a third experiment in which methyl 1H double quantum coherences evolve during a CPMG relaxation element. By fitting single, double, and triple quantum datasets, akin to recording the single quantum dataset at static magnetic fields of Bo, 2Bo and 3Bo, we show that accurate exchange values can be obtained even in cases where exchange rates exceed 10,000s-1. The utility of the double quantum experiment is demonstrated with a pair of cavity mutants of T4 lysozyme (T4L) with ground and excited states interchanged and with exchange rates differing by fourfold (~ 900s-1 and ~ 3600s-1), as well as with a fast-folding domain where the unfolded state lifetime is ~ 80µs.

In silico prediction of novel residues involved in amyloid primary nucleation of human I56T and D67H lysozyme

BackgroundAmyloidogenic proteins are most often associated with neurodegenerative diseases such as Alzheimer’s disease, Parkinson’s disease, and Huntington’s disease, but there are more than two dozen human proteins known to form amyloid fibrils associated with disease. Lysozyme is an antimicrobial protein that is used as a general model to study amyloid fibril formation. Studies aimed at elucidating the process of amyloid formation of lysozyme tend to focus on partial unfolding of the native state due to the relative instability of mutant amyloidogenic variants. While this is well supported, the data presented here suggest the native structure of the variants may also play a role in primary nucleation.ResultsThree-dimensional structural analysis identified lysozyme residues 21, 62, 104, and 122 as displaced in both amyloidogenic variants compared to wild type lysozyme. Residue interaction network (RIN) analysis found greater clustering of residues 112–117 in amyloidogenic variants of lysozyme compared to wild type. An analysis of the most energetically favored predicted dimers and trimers provided further evidence for a role for residues 21, 62, 104, 122, and 112–117 in amyloid formation.ConclusionsThis study used lysozyme as a model to demonstrate the utility of combining 3D structural analysis with RIN analysis for studying the general process of amyloidogenesis. Results indicated that binding of two or more amyloidogenic lysozyme mutants may be involved in amyloid nucleation by placing key residues (21, 62, 104, 122, and 112–117) in proximity before partial unfolding occurs. Identifying residues in the native state that may be involved in amyloid formation could provide novel drug targets to prevent a range of amyloidoses.

Molecular Dynamics Study of Amyloidogenic Mutants of Human Lysozyme

The mutants of human lysozyme are capable of fibril formation implicated in the etiology of familial systemic or renal amyloidosis. A series of 100 ns thermal unfolding molecular dynamics (MD) simulations with WT human lysozyme and its seven amyloidogenic variants (I56T, D67H, F57I, W64R, Y54N, F57I/T70N and T70N/W112R) have been performed at 500 K. The molecular dynamics simulations were performed with GROMACS software (version 5.1) using the CHARMM36m force field. The MD results have been analysed in terms of the parameters characterizing both the global and local protein structure, such as the backbone root mean-square deviation, gyration radius, solvent accessible surface area, the root mean-square fluctuations and the secondary structure content. Depending on the observed effects, the examined variants of human lysozyme have been roughly divided into three groups comprising of mutants with faster (Y54N and F57I/T70N), similar (D67H and I56T) or slower (W64, F57I and T70N/W112R) unfolding rate compared to the wild-type counterpart. The analysis of the protein fluctuational behavior revealed that in most mutants the β-domain displays stronger fluctuations (except the W64R and F57I) and higher flexibility of the C- and D-helices relative to the native lysozyme with the exception of W64R and Y54N which show marked decrease (W64R) or increase (Y54N) in mobility of almost all residues. The analysis of secondary structure evolution provided evidence for higher stability of α-domain compared to β-domain. The results obtained reinforce the idea that mutation-induced global structural destabilization is not the only factor contributing to protein misfolding, the modifications in conformation and dynamics of selected protein regions may also play significant role in amyloid fibril formation.

Analysis of Saturation Recovery Amplitudes to Characterize Conformational Exchange in Spin-Labeled Proteins

Analysis of saturation recovery data from spin-labeled proteins is extended to include the amplitudes in addition to the recovery rates for two-site exchange. It is shown that the recovery amplitudes depend strongly on the exchange rate between states as well as their populations and this dependence provides a simple criterion to identify exchange rates in the 10–1000 kHz range. Analysis of experimental SR relaxation curves via the uniform penalty (UPEN) method allows for reliable identification of single, double, or other multiple-component traces, and global fitting of a set of relaxation curves using both relaxation rates and amplitudes determined from the UPEN fits allows for the estimation of exchange rate in the above domain. The theory is tested on simple model systems, and applied to the determination of conformational exchange rates in spin-labeled mutants of T4 Lysozyme and intestinal fatty acid binding protein. Finally, an example of T 1-weighted spectral editing is provided for systems in the slow exchange limit.

Double-Arm Lanthanide Tags Deliver Narrow Gd3+ -Gd3+ Distance Distributions in Double Electron-Electron Resonance (DEER) Measurements.

Double-arm cyclen-based Gd3+ tags are shown to produce accurate nanometer scale Gd3+ -Gd3+ distance measurements in double electron-electron resonance (DEER) experiments by confining the space accessible to the metal ion. The results show excellent agreement with predictions both for the maximum and width of the measured distance distributions. For distance measurements in proteins, the tags can be attached to two cysteine residues located in positions i and i+4, or i and i+8, of an α-helix. In the latter case, an additional mutation introducing an aspartic acid at position i+4 achieves particularly narrow distribution widths. The concept is demonstrated with cysteine mutants of T4 lysozyme and maltose binding protein. We report the narrowest Gd3+ -Gd3+ distance distributions observed to date for a protein. By limiting the contribution of tag mobility to the distances measured, double-arm Gd3+ tags open new opportunities to study the conformational landscape of proteins in solution with high sensitivity.

Synthesis of 1,2-Azaborines and the Preparation of Their Protein Complexes with T4 Lysozyme Mutants

Synthesis of 1,2-Azaborines and the Preparation of Their Protein Complexes with T4 Lysozyme Mutants

Synthesis of 1,2-Azaborines and the Preparation of Their Protein Complexes with T4 Lysozyme Mutants.

We describe a general synthesis of 1,2-azaborines using standard air-free techniques and protein complex preparation with T4 lysozyme mutants by vapor diffusion. Oxygen- and moisture-sensitive compounds are prepared and isolated under an inert atmosphere (N2) using either a vacuum gas manifold or a glove box. As an example of azaborine synthesis, we demonstrate the synthesis and purification of the volatile N-H-B-ethyl-1,2-azaborine by a five-step sequence involving distillation and column chromatography for the isolation of products. T4 lysozyme mutants L99A and L99A/M102Q are expressed with Escherichia coli RR1 strain. Standard protocols for chemical cell lysis followed by purification using carboxymethyl ion exchange column affords protein of sufficiently high purity for crystallization. Protein crystallization is performed in various concentrations of precipitant at different pH ranges using the hanging drop vapor diffusion method. Complex preparation with the small molecules is carried out by vapor diffusion method under an inert atmosphere. X-ray diffraction analysis of the crystal complex provides unambiguous structural evidence of binding interactions between the protein binding site and 1,2-azaborines.