Hen Egg-white Lysozyme Protein Research Articles

PH-Dependent HEWL-AuNPs Interactions: Optical Study.

Optical methods (spectroscopy, spectrofluorometry, dynamic light scattering, and refractometry) were used to study the change in the state of hen egg-white lysozyme (HEWL), protein molecules, and gold nanoparticles (AuNPs) in aqueous colloids with changes in pH, and the interaction of protein molecules with nanoparticles was also studied. It was shown that changing pH may be the easiest way to control the protein corona on gold nanoparticles. In a colloid of nanoparticles, both in the presence and absence of protein, aggregation-deaggregation, and in a protein colloid, monomerization-dimerization-aggregation are the main processes when pH is changed. A specific point at pH 7.5, where a transition of the colloidal system from one state to another is observed, has been found using all the optical methods mentioned. It has been shown that gold nanoparticles can stabilize HEWL protein molecules at alkaline pH while maintaining enzymatic activity, which can be used in practice. The data obtained in this manuscript allow for the state of HEWL colloids and gold nanoparticles to be monitored using one or two simple and accessible optical methods.

Hen Egg White Lysozyme (HEWL) Confers Resistance to Verticillium Wilt in Cotton by Inhibiting the Spread of Fungus and Generating ROS Burst.

Verticillium wilt is a soil-borne vascular disease caused by the fungal pathogen Verticillium dahliae. It causes great harm to upland cotton (Gossypium hirsutum) yield and quality. A previous study has shown that Hen egg white lysozyme (HEWL) exerts strong inhibitory activity against V. dahliae in vitro. In the current study, we introduced the HEWL gene into cotton through the Agrobacterium-mediated transformation, and the exogenous HEWL protein was successfully expressed in cotton. Our study revealed that HEWL was able to significantly inhibit the proliferation of V. dahlia in cotton. Consequently, the overexpression of HEWL effectively improved the resistance to Verticillium wilt in transgenic cotton. In addition, ROS accumulation and NO content increased rapidly after the V. dahliae inoculation of plant leaves overexpressing HEWL. In addition, the expression of the PR genes was significantly up-regulated. Taken together, our results suggest that HEWL significantly improves resistance to Verticillium wilt by inhibiting the growth of pathogenic fungus, triggering ROS burst, and activating PR genes expression in cotton.

In Silico Investigation on the Selective Nanotoxicity of Two-Dimensional Materials to Hen Egg White Lysozyme Protein

Inspired by the urge to determine the underlying mechanism of the induction of toxic effects of state-of-the-art two-dimensional (2D) nanomaterials such as graphene and hexagonal boron nitride (h-BN) toward biomolecules, herein we study the interactions between these nanomaterials with a model protein hen egg white lysozyme (HEWL), employing classical molecular dynamics simulations. It is revealed that the protein gets easily adsorbed on both of the 2D materials, with the interaction energy being much higher in the case of h-BN, suggesting a significantly stronger adsorption affinity. The interactions of aromatic amino acid residues such as tyrosine and tryptophan along with few aliphatic residues such as arginine, lysine, and asparagine with the 2D materials are found to be pronounced, and most of the residues taking part in adsorption are nearly the same for both materials. While the secondary structure of HEWL remains nearly unaltered upon adsorption on graphene, h-BN massively perturbs both the α-helix and β-sheet components through the disruption of intraprotein hydrogen bonds, which are in turn sine quo non for the preservation of the structural integrity. It is demonstrated that the disruption of the secondary structure is due to pronounced thermodynamic preference for the adsorption of the constituent amino acid residues on h-BN compared to their spatial disposition within the proteins. The calculated release times from the adsorbed state are found to be orders of magnitude higher in the case of h-BN compared to graphene, and it is unlikely that the protein would get released in accessible time scales unless dislodged by the application of an external force. The present study contributes to the fundamental understanding of the nanotoxicity of emerging 2D materials toward proteins, thereby aiding experimentalists to design biocompatible 2D materials for nano-biomedical usage and device fabrication.

DMSO and TMAO-Differences in Interactions in Aqueous Solutions of the K-Peptide.

Interactions between a solvent and their co-solute molecules in solutions of peptides are crucial for their stability and structure. The K-peptide is a synthetic fragment of a larger hen egg white lysozyme protein that is believed to be able to aggregate into amyloid structures. In this study, a complex experimental and theoretical approach is applied to study systems comprising the peptide, water, and two co-solutes: trimethylamide N-oxide (TMAO) or dimethyl sulfoxide (DMSO). Information about their interactions in solutions and on the stability of the K-peptide was obtained by FTIR spectroscopy and differential scanning microcalorimetry. The IR spectra of various osmolyte–water–model-peptide complexes were simulated with the DFT method (B3LYP/6-311++G(d,p)). The FTIR results indicate that both solutes are neutral for the K-peptide in solution. Both co-solutes affect the peptide to different degrees, as seen in the shape of its amide I band, and have different influences on its thermal stability. DFT calculations helped simplify the experimental data for easier interpretation.

Discovery of a tetracyclic indole alkaloid that postpones fibrillation of hen egg white lysozyme protein

Protein aggregation, such as amyloid fibril formation, is molecular hallmark of many neurodegenerative disorders including Alzheimer's, Parkinson's, and Prion disease. Indole alkaloids are well-known as the compounds having the ability to inhibit protein fibrillation. In this study, we experimentally and computationally have investigated the anti-amyloid property of a derivative of a synthesized tetracyclic indole alkaloid (TCIA), possessing capable functional groups. The fibrillation reaction of Hen White Egg Lysozyme (HEWL) was performed in absence and presence of the indole alkaloid. For quantitative analysis, we used Thioflovin T binding assay which showed ~50% reduction in fibril formation in the presence of 20 μM TCIA. Using TEM imaging, we observed a significant morphological change in our model protein in the presence of TCIA. In addition, we exploited FT-IR assay by which Amide I peak's shifting toward lower wavenumber was clearly observed. Using Molecular Docking, the interaction of the inhibitor (TCIA) with the protein's amyloidogenic region was modeled. Also, different biophysical parameters were calculated by Molecular Dynamics (MD) simulation. Various biochemical assays, conformational change, and hydrophobicity exposure of the protein during amyloid formation indicated that the compound assists HEWL to keep its native structure via destabilizing β-sheet structure.

Heterogeneous nanozymatic activity of Hf oxo-clusters embedded in a metal-organic framework towards peptide bond hydrolysis.

Materials with enzyme-like activities and proteolytic potential are emerging as a robust and effective alternative to natural enzymes. Herein, a Hf6O8-based NU-1000 metal organic framework (Hf-MOF) is shown to act as a heterogeneous catalyst for the hydrolysis of peptide bonds under mild conditions. In the presence of Hf-MOF, a glycylglycine model dipeptide was hydrolysed with a rate constant of kobs = 8.33 × 10-7 s-1 (half-life (t1/2) of 231 h) at 60 °C and pD 7.4, which is significantly faster than the uncatalyzed reaction. Other Gly-X peptides (X = Ser, Asp, Ile, Ala, and His) were also smoothly hydrolysed under the same conditions with similar rates, except for the faster reactions observed for Gly-His and Gly-Ser. Moreover, the Hf6O8-based NU-1000 MOF also exhibits a high selectivity in the cleavage of a protein substrate, hen egg white lysozyme (HEWL). Our results suggest that embedding Hf6O8 oxo-clusters is an efficient strategy to conserve the hydrolytic activity while smoothing the strong substrate adsorption previously observed for a discrete Hf oxo-cluster that hindered further development of its proteolytic potential. Furthermore, comparison with isostructural Zr-NU-1000 shows that although the Hf variant afforded the same cleavage pattern towards HEWL but slightly slower reaction rates, it exhibited a larger stability window and a better recyclability profile. The results suggest that these differences originate from the intrinsic differences between HfIV and ZrIV centers, and from the lower surface area of Hf-NU-1000 in comparison to Zr-NU-1000.

Origin of Selectivity in Protein Hydrolysis by Zr(IV)-Containing Metal Oxides as Artificial Proteases

The origin of selectivity in protein hydrolysis promoted by Zr(IV)-substituted polyoxometalates (POMs) has been investigated through a variety of computational techniques. Initially, we analyzed the reaction mechanism for the observed hydrolysis at the Asn44—Arg45 site in the hen egg-white lysozyme protein (HEWL) by the Zr-substituted Lindqvist anion [W5O18Zr(H2O)(OH)]3– (ZrL) using cluster models obtained from molecular dynamics (MD) simulations, quantum mechanics/molecular mechanics (QM/MM) calculations, and metadynamics simulations. The mechanism characterization shows that the overall activity is governed by the energy cost to reach the transition state for C–N bond cleavage from the reactants, resulting in a calculated overall free-energy barrier of 121 kJ mol–1 that is in excellent agreement with the values derived from the experimental rate constants (113–134 kJ mol–1). In addition, QM/MM metadynamics simulations on the early stages of the mechanism revealed the formation of an exergonic non-covalent POM···protein complex at the protein surface that was stabilized by positively charged amino acids maintained during the Zr coordination to the amide oxygen. For nonreactive related sites containing Arg (Asn113—Arg114, Arg45—Asn46, and Arg21—Gly22,) we found very similar overall barriers within the cluster model approach (124, 124, and 120 kJ mol–1, respectively); however, their nonbonding POM···protein interactions along the simulated coordination of Zr to the amide oxygen were significantly weaker than those for the reactive Asn44—Arg45 site. Thus, for the HEWL protein the selectivity is governed by an enzyme-like recognition of ZrL at the cleavage site that results in an overall acceleration of the reaction rate compared to those at other sites. Conversely for human serum albumin, (HSA) the observed selectivity was not directed by nonbonding POM···protein interactions but instead was controlled by the protein secondary structure. Calculations on several Arg—Leu sites placed in positive patches showed that peptide bonds in an α-helix structure have higher overall free-energy barriers, while for the active Arg114—Leu115 site in a random coil region the C–N cleavage is facilitated by the extended conformation of the protein chain. All in all, this study has identified and evaluated two complementary factors controlling the selectivity in peptide hydrolysis promoted by transition metal-substituted POMs; hydrolysis is disfavored at α-helical regions of the protein, and then specific positively charged patches can trap the POM via electrostatic-type POM···protein interactions and accelerate the reaction.

Microprocessing on Single Protein Crystals Using Femtosecond Pulse Laser.

Proteins with different micropatterns have various applications in biosensing, structural analysis, and other biomedical fields. However, processing of micropatterns on single protein crystals remains a challenge due to the fragility of protein molecules. In this work, we studied femtosecond laser processing on single hen egg white lysozyme protein crystals. Optimized laser parameters were found to achieve micropatterning without cracking of protein crystals. The ablation morphology dependence on the laser fluence and the pulse number was discussed to control the processing results. Under a laser fluence higher than 1 J/cm2, the ablation hole was formed. While multipulses with fluence lower than the ablation threshold were applied, the foaming area was observed due to the denaturation of protein. The numerical simulation shows that the ablation results were influenced by the ionization and energy deposition process. Micropatterns including lines, areas, and microarrays can be processed with a minimum size of 2 μm. Processed patterns on the crystal surface can be used for biosensing microarrays and the enhancement of crystal growth. The microprocessing method proposed in this study has potential applications in different fields including biodevices and biomedicine.

Structure-Activity Relationships for the Affinity of Chaotropic Polyoxometalate Anions towards Proteins.

The influence of the composition of chaotropic polyoxometalate (POM) anions on their affinity to biological systems was studied by means of atomistic molecular dynamics (MD) simulations. The variations in the affinity to hen egg-white lysozyme (HEWL) were analyzed along two series of POMs whereby the charge or the size and shape of the metal cluster are modified systematically. Our simulations revealed a quadratic relationship between the charge of the POM and its affinity to HEWL as a consequence of the parabolic growth of POM⋅⋅⋅water interaction with the charge. As the charge increases, POMs become less chaotropic (more kosmotropic) increasing the number and the strength of POM-water hydrogen bonds and structuring the solvation shell around the POM. This atomistic description explains the proportionally larger desolvation energies and less protein affinity for highly charged POMs, and consequently, the preference for moderate charge densities (q/M=0.33). Also, our simulations suggest that POM⋅⋅⋅protein interactions are size-specific. The cationic pockets of HEWL protein show a preference for Keggin-like structures, which display the optimal dimensions (≈1 nm). Finally, we developed a quantitative multidimensional model for protein affinity with predictive ability (r2 =0.97; q2 =0.88) using two molecular descriptors that account for the charge density (charge per metal atom ratio; q/M) and the size and shape (shape weighted-volume; VS ).

Interplay between structural parameters and reactivity of Zr6-based MOFs as artificial proteases.

Structural parameters influencing the reactivity of metal–organic frameworks (MOF) are challenging to establish. However, understanding their effect is crucial to further develop their catalytic potential. Here, we uncovered a correlation between reaction kinetics and the morphological structure of MOF-nanozymes using the hydrolysis of a dipeptide under physiological pH as model reaction. Comparison of the activation parameters in the presence of NU-1000 with those observed with MOF-808 revealed the reaction outcome is largely governed by the Zr6 cluster. Additionally, its structural environment completely changes the energy profile of the hydrolysis step, resulting in a higher energy barrier ΔG‡ for NU-1000 due to a much larger ΔS‡ term. The reactivity of NU-1000 towards a hen egg white lysozyme protein under physiological pH was also evaluated, and the results pointed to a selective cleavage at only 3 peptide bonds. This showcases the potential of Zr-MOFs to be developed into heterogeneous catalysts for non-enzymatic but selective transformation of biomolecules, which are crucial for many modern applications in biotechnology and proteomics.

Do carbon nanotubes inhibit or promote amyloid fibrils formation?

Objectives: The purpose of the work was to study the effect of carbon nanotubes on the formation of fibril structures in lysozyme at room temperature under different pH values. Materials and methods: For the preparation of the samples, hen egg-white lysozyme protein (HEWL, Fluka), as well as single-walled (SWCNT, Sigma-Aldrich) and multi-walled (MWCNT, OOO TM “Spetsmash”, Kyiv, Ukraine) carbon nanotubes were used. Used techniques: IR-Fourier Absorption Spectroscopy; confocal microscopy. Results: In this paper, the study of molecular mechanisms of interaction of lysozyme with carbon nanotubes by vibrational spectroscopy was carried out and a conformational analysis of the formed complexes was performed. It is shown that carbon nanotubes can affect the structure of lysozyme even at room temperature and normal pH values, as evidenced by conformational changes in lysozyme due to interaction with carbon nanotubes. Complexes which are formed as a result of such interaction, have characteristic features of amyloid fibrillar structures. It reveals one of possible mechanisms of carbon nanotubes cytotoxicity. On the other hand, such a technique can be introduced to obtain model amyloid fibrils for further study. Conclusion: The method of vibtarional spectroscopy has shown that carbon nanotubes can influence the structure of lysozyme, as it is shown by the conformational analysis of the absorption band Amide I. After the interaction of lysozyme with CNT, an increase in the contribution of antiparallel β-conformation in the structure of lysozyme is observed, and the contribution of the α-helix conformation is reduced, which are characteristic features in the formation of fibrillar structures. The possibility of amyloid fibril formation without the use of high temperatures at different pH values with the interaction of lysozyme and carbon nanotubes, which can be applied as a method for obtaining the model amyloid fibrils, is shown.

Superactivity of MOF-808 toward Peptide Bond Hydrolysis.

MOF-808, a Zr(IV)-based metal-organic framework, has been proven to be a very effective heterogeneous catalyst for the hydrolysis of the peptide bond in a wide range of peptides and in hen egg white lysozyme protein. The kinetic experiments with a series of Gly-X dipeptides with varying nature of amino acid side chain have shown that MOF-808 exhibits selectivity depending on the size and chemical nature of the X side chain. Dipeptides with smaller or hydrophilic residues were hydrolyzed faster than those with bulky and hydrophobic residues that lack electron rich functionalities which could engage in favorable intermolecular interactions with the btc linkers. Detailed kinetic studies performed by 1H NMR spectroscopy revealed that the rate of glycylglycine (Gly-Gly) hydrolysis at pD 7.4 and 60 °C was 2.69 × 10-4 s-1 ( t1/2 = 0.72 h), which is more than 4 orders of magnitude faster compared to the uncatalyzed reaction. Importantly, MOF-808 can be recycled several times without significantly compromising the catalytic activity. A detailed quantum-chemical study combined with experimental data allowed to unravel the role of the {Zr6O8} core of MOF-808 in accelerating Gly-Gly hydrolysis. A mechanism for the hydrolysis of Gly-Gly by MOF-808 is proposed in which Gly-Gly binds to two Zr(IV) centers of the {Zr6O8} core via the oxygen atom of the amide group and the N-terminus. The activity of MOF-808 was also demonstrated toward the hydrolysis of hen egg white lysozyme, a protein consisting of 129 amino acids. Selective fragmentation of the protein was observed with 55% yield after 25 h under physiological pH.

Insights into the mechanism of how Morin suppresses amyloid fibrillation of hen egg white lysozyme

This communication describes the inhibitory effect of Morin on the fibrillation of Hen Egg White Lysozyme (HEWL), a generic amyloid-forming model protein. This effect was dose-dependent and stronger than other small molecules we have tested previously. Spectrofluorometric and computational studies support a model suggesting that Morin inhibits amyloid fibril formation of HEWL by binding to the aggregation prone cleft region of the β-domain of HEWL, thereby stabilizing the molecule in its native-like state. Interestingly, transmission electron microscopy observations suggest that, along with increases in Morin concentration, the observed amorphous aggregates became larger and morphologically different. We propose that following occupation of the binding cleft, excess Morin adheres and coats the HEWL protein surface, thereby minimizing the interaction between the protein surface and water molecules.

Improved Relaxation Mode Analysis of a Hen Egg-White Lysozyme Protein

Improved Relaxation Mode Analysis of a Hen Egg-White Lysozyme Protein

Combining weak affinity chromatography, NMR spectroscopy and molecular simulations in carbohydrate–lysozyme interaction studies

By examining the interactions between the protein hen egg-white lysozyme (HEWL) and commercially available and chemically synthesized carbohydrate ligands using a combination of weak affinity chromatography (WAC), NMR spectroscopy and molecular simulations, we report on new affinity data as well as a detailed binding model for the HEWL protein. The equilibrium dissociation constants of the ligands were obtained by WAC but also by NMR spectroscopy, which agreed well. The structures of two HEWL-disaccharide complexes in solution were deduced by NMR spectroscopy using (1)H saturation transfer difference (STD) effects and transferred (1)H,(1)H-NOESY experiments, relaxation-matrix calculations, molecular docking and molecular dynamics simulations. In solution the two disaccharides β-d-Galp-(1→4)-β-D-GlcpNAc-OMe and β-D-GlcpNAc-(1→4)-β-D-GlcpNAc-OMe bind to the B and C sites of HEWL in a syn-conformation at the glycosidic linkage between the two sugar residues. Intermolecular hydrogen bonding and CH/π-interactions form the basis of the protein-ligand complexes in a way characteristic of carbohydrate-protein interactions. Molecular dynamics simulations with explicit water molecules of both the apo-form of the protein and a ligand-protein complex showed structural change compared to a crystal structure of the protein. The flexibility of HEWL as indicated by a residue-based root-mean-square deviation analysis indicated similarities overall, with some residue specific differences, inter alia, for Arg61 that is situated prior to a flexible loop. The Arg61 flexibility was notably larger in the ligand-complexed form of HEWL. N,N'-Diacetylchitobiose has previously been observed to bind to HEWL at the B and C sites in water solution based on (1)H NMR chemical shift changes in the protein whereas the disaccharide binds at either the B and C sites or the C and D sites in different crystal complexes. The present study thus highlights that protein-ligand complexes may vary notably between the solution and solid states, underscoring the importance of targeting the pertinent binding site(s) for inhibition of protein activity and the advantages of combining different techniques in a screening process.

Purification and characterization of angiotensin I-converting enzyme inhibitory peptides from enzymatic hydrolysate of hen egg white lysozyme

The aims of this study were to explore new functional bioactivity of hen egg white lysozyme (HEWL) protein hydrolysate and to identify potential angiotensin I-converting enzyme (ACE) inhibitory peptides from HEWL. The protein was hydrolyzed by gastrointestinal enzymes, including pepsin, α-chymotrypsin and trypsin, under simulated physiological conditions, and the hydrolysate exhibited potent ACE inhibitory activity with an IC50 value of 15.6±1.4μg/mL. Three novel ACE inhibitory peptides, specifically, Met-Lys-Arg, Arg-Gly-Tyr and Val-Ala-Trp with IC50 values of 25.7±0.2, 61.9±0.1 and 2.86±0.08μM, respectively, were isolated from the hydrolysate of HEWL in two stages of reverse-phase high-performance liquid chromatography (RP-HPLC), and their sequences were analyzed by UPLC–MS/MS. Kinetics studies suggested that the purified peptides were all competitive inhibitors; preincubation experiments implied that the peptide Arg-Gly-Tyr was a true substrate and that the other two peptides were true ACE inhibitors. These results indicate that HEWL may become an effective source of ACE inhibitory peptides. This study reveals the antihypertensive bioactivity of HEWL protein hydrolysate and provides further evidence that eggs are an excellent source of health-promoting food.

Presentation of Type B Peptide–MHC Complexes from Hen Egg White Lysozyme by TLR Ligands and Type I IFNs Independent of H2-DM Regulation

In APCs, presentation by MHC II molecules of the chemically dominant peptide from the protein hen egg white lysozyme (HEL) generates different conformational isomers of the peptide-MHC II complexes (pMHC). Type B pMHCs are formed in early endosomes from exogenous peptides in the absence of H2-DM, whereas in contrast, type A pMHC complexes are formed from HEL protein in late vesicles after editing by H2-DM. Thus, H2-DM edits off the more unstable pMHC complexes, which are not presented from HEL. In this study, we show that type B pMHC complexes were presented from HEL protein only after stimulation of dendritic cells (DC) with TLR ligands or type I IFN. Type I IFN contributed to most TLR ligand-induced type B pMHC generation, as presentation decreased in DC lacking the receptor for type I IFNs (IFNAR1(-/-)). In contrast, presentation of type A pMHC from HEL and from peptide was minimally affected by TLR ligands. The relative effectiveness of CD8α(+) DC or CD8α(-) DC in presenting type B pMHC complexes varied depending on the TLR ligand used. The mechanisms of generation of type B pMHC from HEL protein with TLR stimulation did not involve H2-DM or release of peptides. DC from H2-DM-deficient mice in the presence of TLR ligands presented type B pMHC. Such DC showed a slight enhancement of HEL catabolism, but peptide release was not evident. Thus, TLR ligands and type I IFN alter the pathways of presentation by MHC II molecules of DC such that type B pMHCs are generated from protein Ag.

Key Residues that Play a Critical Role in Urea-Induced Lysozyme Unfolding

In this paper, we have developed a simple sensitivity score, based on the relative population of solvent molecules near each residue, to analyze the detailed motions of both urea and water around the hen egg-white lysozyme protein (W62G mutant) during its early stage of urea-induced unfolding for a better understanding of the atomic picture of the chemical denaturation process. Our simulation and analysis show that some hydrophobic core residues can keep dry from water for tens of nanoseconds in 8 M urea, while their contacts with urea increase significantly at the same time, forming a molten dry-globule-like state. Also, different from previously proposed actions that urea molecules preferentially absorb onto charged residues, our analysis shows that the noncharged residues, rather than the charged ones, attract more urea molecules in their surroundings (acting as attractants for urea), which is consistent with our earlier findings that urea molecules preferentially bind to protein through their stronger dispersion interactions than water. Once the initial adsorption surrounding the protein surface is accomplished, the further intrusion is found to be facilitated by a group of key residues, including Leu8, Met12, Val29, and Ala95, which play a critical role in the formation of the dry-globule structure. These hydrophobic dry residues form a local contact map which excludes the intrusion of water but accommodates the presence of urea due to their stronger binding to protein during this swelling process, thus maintaining an interesting transient dry-globule state.

Nonnative Electrostatic Interactions Can Modulate Protein Folding: Molecular Dynamics with a Grain of Salt

In recent years, a growing number of protein folding studies have focused on the unfolded state, which is now recognized as playing a major role in the folding process. Some of these studies show that interactions occurring in the unfolded state can significantly affect the stability and kinetics of the protein folding reaction. In this study, we modeled the effect of electrostatic interactions, both native and nonnative, on the folding of three protein systems that underwent selective charge neutralization or reversal or complete charge suppression. In the case of the N-terminal L9 protein domain, our results directly attribute the increase in thermodynamic stability to destabilization of the unfolded ensemble, reaffirming the experimental observations. These results provide a deeper structural insight into the ensemble of the unfolded state and predict a new mutation site for increased protein stability. In the second case, charge reversal mutations of RNase Sa affected protein stability, with the destabilizing mutations being less destabilizing at higher salt concentrations, indicating the formation of charge–charge interactions in the unfolded state. In the N-terminal L9 and RNase Sa systems, changes in electrostatic interactions in the unfolded state that cause an increase in free energy had an overall compaction effect that suggests a decrease in entropy. In the third case, in which we compared the β-lactalbumin and hen egg-white lysozyme protein homologues, we successfully eliminated differences between the folding kinetics of the two systems by suppressing electrostatic interactions, supporting previously reported findings. Our coarse-grained molecular dynamics study not only reproduces experimentally reported findings but also provides a detailed molecular understanding of the elusive unfolded-state ensemble and how charge–charge interactions can modulate the biophysical characteristics of folding.

Molecular Dynamics Simulations of Hen Egg White Lysozyme Adsorption at a Charged Solid Surface

Hen egg white lysozyme (HEWL) adsorption on negatively charged, hydrophilic surfaces has been investigated using atomistic molecular dynamics. Analysis of six 20 ns trajectories performed at pH 7 and ionic strength 0.02 M (NaCl) reveals that conformational alterations are required for HEWL adsorption, and that upon adsorption the protein loses some alpha-helical content. Simulations for a few different initial orientations show that the HEWL protein adsorbs on a flat surface with an angle between the protein long axis and the surface of about 45 degrees . The main adsorption site is located on the N,C-terminal part of the HEWL surface; the major role is played by Lys1, Arg5, Arg14, and Arg128. Adsorption is not found with contrary orientations. Two additional 20 ns trajectories calculated with 0.5 M ionic strength suggest that the main force governing adsorption is electrostatic attraction between parts of the protein and the surface. A trajectory obtained for the protein situated inside a cubic box built from the charged surfaces shows that the adsorption pattern is different for flat and nonflat surfaces, and in particular, adsorption on the nonflat surface requires tertiary structure alterations and partial unfolding. The observed trends are consistent with both experimental and previous computational studies.