Validation of the GROMOS force-field parameter set 45A3 against nuclear magnetic resonance data of hen egg lysozyme

The quality of molecular dynamics (MD) simulations of proteins depends critically on the biomolecular force field that is used. Such force fields are defined by force-field parameter sets, which are generally determined and improved through calibration of properties of small molecules against experimental or theoretical data. By application to large molecules such as proteins, a new force-field parameter set can be validated. We report two 3.5 ns molecular dynamics simulations of hen egg white lysozyme in water applying the widely used GROMOS force-field parameter set 43A1 and a new set 45A3. The two MD ensembles are evaluated against NMR spectroscopic data NOE atom–atom distance bounds, ³J_NH and ³J coupling constants, and ¹5N relaxation data. It is shown that the two sets reproduce structural properties about equally well. The 45A3 ensemble fulfills the atom–atom distance bounds derived from NMR spectroscopy slightly less well than the 43A1 ensemble, with most of the NOE distance violations in both ensembles involving residues located in loops or flexible regions of the protein. Convergence patterns are very similar in both simulations atom-positional root-mean-square differences (RMSD) with respect to the X-ray and NMR model structures and NOE inter-proton distances converge within 1.0–1.5 ns while backbone ³J_HN-coupling constants and ¹H– ¹5N order parameters take slightly longer, 1.0–2.0 ns. As expected, side-chain ³J-coupling constants and ¹H– ¹5N order parameters do not reach full convergence for all residues in the time period simulated. This is particularly noticeable for side chains which display rare structural transitions. When comparing each simulation trajectory with an older and a newer set of experimental NOE data on lysozyme, it is found that the newer, larger, set of experimental data agrees as well with each of the simulations. In other words, the experimental data converged towards the theoretical result.     

Molecular dynamics simulation of thionated hen egg white lysozyme

Understanding of the driving forces of protein folding is a complex challenge because different types of interactions play a varying role. To investigate the role of hydrogen bonding involving the backbone, the effect of thio substitutions in a protein, hen egg white lysozyme (HEWL), was investigated through molecular dynamics simulations of native as well as partly (only residues in loops) and fully thionated HEWL using the GROMOS 54A7 force field. The results of the three simulations show that the structural properties of fully thionated HEWL clearly differ from those of the native protein, while for partly thionated HEWL they only changed slightly compared with native HEWL. The analysis of the torsional-angle distributions and hydrogen bonds in the backbone suggests that the α-helical segments of native HEWL tend to show a propensity to convert to 3(10)-helical geometry in fully thionated HEWL. A comparison of the simulated quantities with experimental NMR data such as nuclear overhauser effect (NOE) atom-atom distance bounds and (3)J((H)(N)(H)(α))-couplings measured for native HEWL illustrates that the information content of these quantities with respect to the structural changes induced by thionation of the protein backbone is rather limited.     

An analysis of the Co2+-induced nuclear magnetic resonance perturbations of hen egg white lysozyme

A general methodology is presented for analyzing dipolar shifts induced by paramagnetic ions in the nuclear magnetic resonance (NMR) spectra of ligand molecules. The method is applied to the shift perturbations induced by Co2+ in the spectrum of hen egg white lysozyme. A hypothesis testing scheme is employed to evaluate statistically the relative precision with which the axially symmetric and non-axially symmetric forms of the dipolar shift equation fit the observed data. The assumption of axial symmetry for the magnetic susceptibility tensor of Co2+ is rejected at the confidence level of 99%. Since the results presented here are similar to those reached in our analysis of lanthanide-induced shifts, we suggest that the assumption of axial symmetry may, in general, not hold. Similar conclusions have been reached by other investigators in studies of paramagnetic metal binding to model systems. We have included the three Co2+ coordinates in an eight-parameter fit of the Co2+ shift data. The Co2+ position obtained from this fit is in statistical agreement with the position inferred from x-ray data. Thus, the analysis of shift data may furnish a means for determining the site of metal complexation in macromolecules whose structure has been determined by x-ray crystallography.     

Molecular modelling of membrane sterols with the use of the GROMOS 96 force field

Membrane located sterols determine the structure and function of eucariotic cell membranes. Moreover, they are targets for important antifungal antibiotic amphotericin B. Knowledge about the geometry and dynamics of sterols in the environment of lipidic membranes is necessary to understand their functions. However, due to the dynamic character of the membrane, no experimental data about sterol behaviour on the molecular level is available. Hence molecular modelling simulations could be a source of useful information. The main goal of this paper is to prove the adequacy of the GROMOS 96 force field for molecular simulations of membrane sterols. We focused our attention on the reproduction of characteristic geometrical features observed in the crystal of cholesterol hemiethanolate by molecular dynamics simulations. The results presented clearly indicate that the GROMOS 96 force field can be a useful tool to simulate the highly lipophilic systems. Moreover, interactions responsible for the stability of such systems can also be recognised.     

Recovery of lysozyme from hen egg white by selective magnetic cake filtration

A new concept for the improvement of the downstream processing and purification is the so‐called magnetic separation. By using surface functionalized magnetic substrate particles, which selectively adsorb the target product, it can be directly separated out of a crude bioprocess stream. These methods are already used for analytical purposes, where only small amounts of functionalized particles are necessary. To apply the same concept at a larger scale, effective and economical procedures have to be provided. First, suitable process equipment has to be developed. Second, the magnetic particles have to be manufactured with a stable surface functionalization and long‐term stability for their reuse. Up to now mainly high‐gradient magnetic separation filter devices are applied for selective magnetic separation. They consist of a magnetic matrix in which the magnetic particles are trapped. In this work, a new magnetic filter is introduced that overcomes the capacity limitations of the current high‐gradient magnetic separation technology. The principle is demonstrated by selective recovery of lysozyme from hen egg white. Prior to the separation experiments magnetic beads with a strong acid cation‐exchange surface functionalization are synthesized. The separation procedure is implemented in only one unit operation. With the implementation of the displacement elution sequence lysozyme can be separated out of a hen egg white solution with a purification factor of PF=36 and a purity P=0.83.     

Phosphorus-31 nuclear magnetic resonance studies of the two phosphoserine residues of hen egg white ovalbumin

Ovalbumin contains two phosphoserine residues that give rise to two well-resolved resonances in a 31P NMR spectrum. Ovalbumin samples that have been digested with a variety of phosphatases may give rise to only one phosphoserine resonance, indicating that one of the two phosphorylated sites is relatively inaccessible for phosphatase action. By comparison of the amino acid sequence of the peptide containing the nonsusceptible phosphate to the overall primary structure, we have assigned the resonances observed (pH 8.3) at 5.0 and 4.75 ppm to phosphoserines-68 and -344, respectively. pH titration behavior and susceptibility of the phosphoserine residues to phosphatases indicate that both are located on the surface of the protein. Both residues have a pKa = 6.00-6.04. Analysis of the Hill coefficients measured for the pH titrations and the JPH coupling constants indicate that neither residue interacts with other charged groups on the surface of the protein. Frequency dependence of 31P NMR parameters shows that at higher magnetic field strengths the contribution of chemical shift anisotropy to the line width becomes very significant. We have calculated from the field-dependent terms that phosphoserine-344 is mobile with respect to the protein surface but that phosphoserine-68 is more restricted in its motion. The latter is also involved in a pH-dependent conformational change, since it is shielded from hydrolysis by phosphatases at higher pH. A comparison of the amino acid sequence of the phosphoserine-68 site shows that it has a striking homology to the active-site peptides of a wide variety of hydrolytic enzymes. Moreover, a comparison with the primary sequences of casein suggests that both proteins are phosphorylated by a protein kinase that specifically recognizes a Ser-X-Glu peptide.     

Structure and internal mobility of proteins: a molecular dynamics study of hen egg white lysozyme

The structure and internal motions of the protein hen egg white lysozyme are studied by analysis of simulation and experimental data. A molecular dynamics simulation and an energy minimization of the protein in vacuum have been made and the results compared with high-resolution structures and temperature factors of hen egg white lysozyme in two different crystal forms and of the homologous protein human lysozyme. The structures obtained from molecular dynamics and energy minimization have root-mean-square deviations for backbone atoms of 2.3 Å and 1.1–1.3 Å, respectively, relative to the crystal structures; the different crystal structures have root-mean-square deviations of 0.73–0.81 Å for the backbone atoms. In comparing the backbone dihedral angles, the difference between the dynamics and the crystal structure on which it is based is the same as that between any two crystal structures. The internal fluctuations of atomic positions calculated from the molecular dynamics trajectory agree well with the temperature factors from the three structures. Simulation and crystal results both show that there are large motions for residues involved in exposed turns of the backbone chain, relatively smaller motions for residues involved in the middle of helices or β-sheet structures, and relatively small motions of residues near disulfide bridges. Also, both the simulation and crystal data show that side-chain atoms have larger fluctuations than main-chain atoms. Moreover, the regions that have large deviations among the x-ray crystal structures, which indicates flexibility, are found to have large fluctuations in the simulation.     

Anti-fibrillation propensity of a flavonoid baicalein against the fibrils of hen egg white lysozyme: potential therapeutics for lysozyme amyloidosis

More than 20 human diseases involve the fibrillation of a specific protein/peptide which forms pathological deposits at various sites. Hereditary lysozyme amyloidosis is a systemic disorder which mostly affects liver, spleen and kidney. This conformational disorder is featured by lysozyme fibril formation. In vivo lysozyme fibrillation was simulated under in vitro conditions using a strong denaturant GdHCl at 3 M concentration. Sharp decline in the ANS fluorescence intensity compared to the partially unfolded states, almost 20-fold increase in ThT fluorescence intensity, increase in absorbance at 450 nm suggesting turbidity, negative ellipticity peak in the far-UVCD at 217 nm, red shift of 50 nm compared to the native state in Congo red assay and appearance of a network of long rope-like fibrils in transmission electron microscope (TEM) analysis suggested HEWL fibrillation. Anti-fibrillation potency of baicalein against the preformed fibrils of HEWL was investigated following ThT assay in which there was a dose-dependent decrease in ThT fluorescence intensity compared to the fibrillar state of HEWL with the maximum effect observed at 150-μM baicalein concentration, loss of negative ellipticity peak in the far-UVCD region, dip in the Rayleigh scattering intensity and absorbance at 350 and 450 nm, respectively, together with a reduction in the density of fibrillar structure in TEM imaging. Thus, it could be suggested that baicalein could prove to be a positive therapeutics for hereditary human lysozyme amyloidosis.     

Stabilization of hen egg white lysozyme by a cavity-filling mutation

Stabilization of a protein using cavity-filling strategy has hardly been successful because of unfavorable van der Waals contacts. We succeeded in stabilizing lysozymes by cavity-filling mutations. The mutations were checked by a simple energy minimization in advance. It was shown clearly that the sum of free energy change caused by the hydrophobicity and the cavity size was correlated very well with protein stability. We also considered the aromatic-aromatic interaction. It is reconfirmed that the cavity-filling mutation in a hydrophobic core is a very useful method to stabilize a protein when the mutation candidate is selected carefully.     

Participation of tryptophan 62 in the self-association of hen egg white lysozyme. Application of natural abundance carbon 13 nuclear magnetic resonance spectroscopy.

Self-association of hen egg white lysozyme in solution of 38 degrees) is examined by means of natural abundance 13C nuclear magnetic resonance spectroscopy. The effect of pH on the resonances of the nonprotonated aromatic carbons of 9 mM lysozyme, and the effect of protein concentration (at pH 7) on these resonances, both indicate that self-association significantly affects the chemical shift of Cgamma of Trp-62, but not the chemical shifts of the other nonprotonated aromatic carbons. This result is consistent with the reported participation of Trp-62 in the intermolecular contact (Banerjee, S.K., Pogolotti, A., and Rupley, J.A. (1975) J. Biol. Chem. 250, 8260-8266). Our results indicate that the resonance of Cgamma or Trp-62 is a convenient monitor of lysozyme self-association. The chemical shift of this resonance reflects the extent of aggregation, while the line width yields information about the lifetime of the intermolecular contact. This lifetime is 1 to 2 ms at 38 degrees (9 mM protein, 0.1 M NaCl, pH 7). Our results also indicate that self-association of lysozyme is not accompanied by any general conformational change, and that binding of a lanthanide ion (at the metal ion binding site near the carboxylate groups of ASP-52 AND Glu-35) strongly suppresses self-association.     

Molecular dynamics simulations to determine the effect of supercritical carbon dioxide on the structural integrity of hen egg white lysozyme

In this study, various molecular dynamics simulations were conducted to investigate the effect of supercritical carbon dioxide on the structural integrity of hen egg white lysozyme. The analyses of backbone root-mean-square deviation, radius of gyration, and secondary structure stability all show that supercritical CO(2) exhibits the ability to increase the stability of this protein, probably as a result of the solvent with less polarity, where hydrophobic interactions stabilizing the native structure are weakened and simultaneously the local hydrogen bonds are strengthened, resulting in stabilization of the secondary structures. The hydrophobic cores in the alpha- and beta-domains also play an important role in preventing this protein from thermal unfolding. As supercritical CO(2) has been attractive for biomedical applications because of the advantages of mild critical condition, nonflammability, nontoxity, and the purity of the resulting products, the structural stabilizing effect found in this study strongly suggests that it is possible to increase the thermostability of hen egg white lysozyme by pretreatment with supercritical CO(2), leading to better industrial applications of this protein.     

Antibacterial activity of hen egg white lysozyme against Listeria monocytogenes Scott A in foods

Egg white lysozyme killed or prevented growth of Listeria monocytogenes Scott A in several foods. Lysozyme was more active in vegetables than in animal-derived foods that we tested. For maximum activity in certain foods, EDTA was required in addition to lysozyme. Lysozyme with EDTA effectively killed inoculated populations of 10(4) L. monocytogenes per g in fresh corn, fresh green beans, shredded cabbage, shredded lettuce, and carrots during storage at 5 degrees C. Control incubations without lysozyme supported growth of L. monocytogenes to 10(6) to 10(7)/g. Lysozyme had less activity in animal-derived foods, including fresh pork sausage (bratwurst) and Camembert cheese. In bratwurst, lysozyme with EDTA prevented L. monocytogenes from growing for 2 to 3 weeks but did not kill significant numbers of cells and did not prevent eventual growth. The control sausages not containing lysozyme supported rapid and heavy growth, which indicated that lysozyme was bacteriostatic for 2 to 3 weeks in fresh pork sausage. We also prepared Camembert cheese containing 10(4) L. monocytogenes cells per g and investigated the changes during ripening in cheeses supplemented with lysozyme and EDTA. Cheeses with lysozyme by itself or together with EDTA reduced the L. monocytogenes population by approximately 10-fold over the first 3 to 4 weeks of ripening. In the same period, the control cheese wheels without added lysozyme with and without chelator slowly started to grown and eventually reached 10(6) to 10(7) CFU/g after 55 days of ripening.(ABSTRACT TRUNCATED AT 250 WORDS)     

Studies of chemically modified histidine residues of proteins by carbon 13 nuclear magnetic resonance spectroscopy. Reaction of hen egg white lysozyme with iodoacetate.

It is shown that natural abundance 13C NMR spectroscopy can be used to determine the structures and relative amounts of chemically modified forms of a histidine residue of a peptide or protein. The unfractionated product of the reaction of N alpha-acetyl-L-histidine with bromoacetate yields four resonances of nonprotonated aromatic carbons. These resonances are assigned (on a one-to-one basis) to C gamma of the intact amino acid, the two monocarboxymethylated derivatives (at N delta1 and N epsilon2), and the dicarboxymethylated derivative. The effect of pH on the chemical shift of C gamma is characteristic for each of the four species. This property is used to study the carboxymethylation of His-15 of hen egg white lysozyme upon treatment with iodoacetate. With the use of various reaction conditions, His 15 is carboxymethylated in detectable quantities only at N epsilon2. The spectra of the various reaction mixtures indicate which conditions are best for maximizing the yield of this derivative. A comparison of the spectrum of chromatographically pure [N epsilon2-carboxymethylhistidine-15]lysozyme with that of the intact protein indicates that the chemical modification does not significantly affect the conformation of the protein (at least in the regions of all aromatic amino acid residues).     

Antibacterial activity of hen egg white lysozyme against Listeria monocytogenes Scott A in foods.

Egg white lysozyme killed or prevented growth of Listeria monocytogenes Scott A in several foods. Lysozyme was more active in vegetables than in animal-derived foods that we tested. For maximum activity in certain foods, EDTA was required in addition to lysozyme. Lysozyme with EDTA effectively killed inoculated populations of 10(4) L. monocytogenes per g in fresh corn, fresh green beans, shredded cabbage, shredded lettuce, and carrots during storage at 5 degrees C. Control incubations without lysozyme supported growth of L. monocytogenes to 10(6) to 10(7)/g. Lysozyme had less activity in animal-derived foods, including fresh pork sausage (bratwurst) and Camembert cheese. In bratwurst, lysozyme with EDTA prevented L. monocytogenes from growing for 2 to 3 weeks but did not kill significant numbers of cells and did not prevent eventual growth. The control sausages not containing lysozyme supported rapid and heavy growth, which indicated that lysozyme was bacteriostatic for 2 to 3 weeks in fresh pork sausage. We also prepared Camembert cheese containing 10(4) L. monocytogenes cells per g and investigated the changes during ripening in cheeses supplemented with lysozyme and EDTA. Cheeses with lysozyme by itself or together with EDTA reduced the L. monocytogenes population by approximately 10-fold over the first 3 to 4 weeks of ripening. In the same period, the control cheese wheels without added lysozyme with and without chelator slowly started to grown and eventually reached 10(6) to 10(7) CFU/g after 55 days of ripening.(ABSTRACT TRUNCATED AT 250 WORDS)     

Importance of the CMAP correction to the CHARMM22 protein force field: dynamics of hen lysozyme

The recently developed CMAP correction to the CHARMM22 force field (C22) is evaluated from 25 ns molecular dynamics simulations on hen lysozyme. Substantial deviations from experimental backbone root mean-square fluctuations and N-H NMR order parameters obtained in the C22 trajectories (especially in the loops) are eliminated by the CMAP correction. Thus, the C22/CMAP force field yields improved dynamical and structural properties of proteins in molecular dynamics simulations.     

L-Arginine increases the solubility of unfolded species of hen egg white lysozyme

L-Arginine (L-Arg) has been widely used as an enhancer of protein renaturation. The mechanism behind its action is still not fully understood. Using hen egg white lysozyme as a model protein, we present data that clearly demonstrate the suppression of the aggregation of denatured protein by L-Arg. By chemical modification of free cysteines, a series of unfolded lysozyme species were obtained that served as models for unfolded and intermediate states during the process of oxidative refolding. An increased equilibrium solubility of unfolded species and intermediates in the presence of L-Arg seems to be its major mechanism of action.     

The dynamical response of hen egg white lysozyme to the binding of a carbohydrate ligand

It has become clear that the binding of small and large ligands to proteins can invoke significant changes in side chain and main chain motion in the fast picosecond to nanosecond timescale. Recently, the use of a "dynamical proxy" has indicated that changes in these motions often reflect significant changes in conformational entropy. These entropic contributions are sometimes of the same order as the total entropy of binding. Thus, it is important to understand the connections amongst motion between the manifold of states accessible to the native state of proteins, the corresponding entropy, and how this impacts the overall energetics of protein function. The interaction of proteins with carbohydrate ligands is central to a range of biological functions. Here, we examine a classic carbohydrate interaction with an enzyme: the binding of wild-type hen egg white lysozyme (HEWL) to the natural, competitive inhibitor chitotriose. Using NMR relaxation experiments, backbone amide and side chain methyl axial order parameters were obtained across apo and chitotriose-bound HEWL. Upon binding, changes in the apparent amplitude of picosecond to nanosecond main chain and side chain motions are seen across the protein. Indeed, binding of chitotriose renders a large contiguous fraction of HEWL effectively completely rigid. Changes in methyl flexibility are most pronounced closest to the binding site, but average to only a small overall change in the dynamics across the protein. The corresponding change in conformational entropy is unfavorable and estimated to be a significant fraction of the total binding entropy.