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.     

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.     

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.     

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.     

Highly perturbed pKa values in the unfolded state of hen egg white lysozyme

The majority of pK_a values in protein unfolded states are close to the amino acid model pK_a values, thus reflecting the weak intramolecular interactions present in the unfolded ensemble of most proteins. We have carried out thermal denaturation measurements on the WT and eight mutants of HEWL from pH 1.5 to pH 11.0 to examine the unfolded state pK_a values and the pH dependence of protein stability for this enzyme. The availability of accurate pK_a values for the folded state of HEWL and separate measurements of mutant-induced effects on the folded state pK_a values, allows us to estimate the pK_a values of seven acidic residues in the unfolded state of HEWL. Asp-48 and Asp-66 display pK_a values of 2.9 and 3.1 in our analysis, thus representing the most depressed unfolded state pK_a values observed to date. We observe a strong correlation between the folded state pK_a values and the unfolded state pK_a values of HEWL, thus suggesting that the unfolded state of HEWL possesses a large degree of native state characteristics.     

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

Abstract

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     

Ionic liquids as selectors for controlling the crystallization nucleation of hen egg white lysozyme

Ionic liquids (ILs) have biomaterial applications and are used for protein crystallization. The effect of two imidazolium-based ILs, 1-butyl-3-methylimidazolium chloride ([C₄mim]Cl) and 1,3-dimethylimidazolium iodide ([dmim]I), on the nucleation kinetics of lysozyme was investigated by determining the nucleation induction time, and nucleation parameters were evaluated. The values of interfacial tension calculated for solutions with added 30 g/L ILs [C₄mim]Cl and [dmim]I, and without added ILs were 99.03, 109.7, and 107.3 mJ/m², respectively. Compared with solutions without IL addition, the critical free energy change, size, and molecular number of critical nuclei decreased and the nucleation rate increased after the addition of [C₄mim]Cl. In contrast, the critical free energy change, size, and molecular number of critical nuclei increased and the nucleation rate decreased after the addition of [dmim]I. These new findings provide insights into controlling lysozyme crystallization separation, and present ILs as potentially useful additives for controlling the crystallization of macromolecules.     

1H-NMR study on the structure of lysozyme from guinea hen egg white.

The structure of lysozyme from guinea hen egg white (GEWL), which differs from hen egg white lysozyme (HEWL) by ten amino acid substitutions, was investigated by nuclear magnetic resonance (NMR) spectroscopy. GEWL and HEWL were very similar to each other in their tertiary structure as judged from the profile of 1H-NMR spectra, pH titration, and an N-acetylglucosamine trisaccharide [(GlcNAc)3 binding experiment. However, we have noticed several characteristics which distinguish GEWL from HEWL. The signal of Trp 108 indole N1H of GEWL was shifted upfield by about 0.3 ppm when compared with that of HEWL, and its hydrogen exchange was faster than that of HEWL. The pKa values of Glu 35 estimated from the pH titration curve of Trp 108 indole N1H were different between GEWL and HEWL. From a careful examination of spectral changes caused by (GlcNAc)3 binding, the changes in the chemical shift values of Trp 28 C5H and Asn 59 alpha CH of GEWL were found to be slightly larger than those of HEWL. Ile 55 of HEWL is replaced by valine in GEWL. Such a replacement may affect the neighboring hydrogen bonding between the main chain C = O of Leu 56 and Trp 108 indole N1H, resulting in a change in the microenvironment of the substrate-binding site near Trp 108.