Molecular dynamics simulations of native and substrate-bound lysozyme. A study of the average structures and atomic fluctuations

A biophysical toolset to probe the microscopic processes underlying protein aggregation and its inhibition by molecular chaperones

Given the breadth and depth of the scientific contributions of Sir Christopher Dobson, with over 870 publications to date, it is inconceivable to convey in a single review the impact of his work and its legacy. This review therefore primarily focuses on his contributions to the development of strategies for preventing aberrant protein misfolding. The first section of this review highlights his seminal work on the elucidation of the microscopic nucleation processes underlying protein aggregation. Next, we discuss the specific inhibition of these steps by candidate drugs and biologics, with a particular emphasis on the role of molecular chaperones. In the final section, we review how protein aggregation principles can be exploited for the rational design of novel and more potent aggregation inhibitors. These milestones serve as excellent examples of the profound impact of Dobson's seminal work on fundamental science and its translation into drug discovery.

Understanding the Structure-Function Relationship of Lysozyme Resistance in Staphylococcus aureus by Peptidoglycan O-Acetylation Using Molecular Docking, Dynamics, and Lysis Assay

Lysozyme is an important component of the host innate defense system. It cleaves the β-1,4 glycosidic bonds between N-acetylmuramic acid and N-acetylglucosamine of bacterial peptidoglycan and induce bacterial lysis. Staphylococcus aureus (S. aureus), an opportunistic commensal pathogen, is highly resistant to lysozyme, because of the O-acetylation of peptidoglycan by O-acetyl transferase (oatA). To understand the structure-function relationship of lysozyme resistance in S. aureus by peptidoglycan O-acetylation, we adapted an integrated approach to (i) understand the effect of lysozyme on the growth of S. aureus parental and the oatA mutant strain, (ii) study the lysozyme induced lysis of exponentially grown and stationary phase of both the S. aureus parental and oatA mutant strain, (iii) investigate the dynamic interaction mechanism between normal (de-O-acetylated) and O-acetylated peptidoglycan substrate in complex with lysozyme using molecular docking and molecular dynamics simulations, and (iv) quantify lysozyme resistance of S. aureus parental and the oatA mutant in different human biological fluids. The results indicated for the first time that the active site cleft of lysozyme binding with O-acetylated peptidoglycan in S. aureus was sterically hindered and the structural stability was higher for the lysozyme in complex with normal peptidoglycan. This could have conferred reduced survival of the S. aureus oatA mutant in different human biological fluids. Consistent with this computational analysis, the experimental data confirmed decrease in the growth, lysozyme induced lysis, and lysozyme resistance, due to peptidoglycan O-acetylation in S. aureus.

Structure of a low-population intermediate state in the release of an enzyme product

Enzymes can increase the rate of biomolecular reactions by several orders of magnitude. Although the steps of substrate capture and product release are essential in the enzymatic process, complete atomic-level descriptions of these steps are difficult to obtain because of the transient nature of the intermediate conformations, which makes them largely inaccessible to standard structure determination methods. We describe here the determination of the structure of a low-population intermediate in the product release process by human lysozyme through a combination of NMR spectroscopy and molecular dynamics simulations. We validate this structure by rationally designing two mutations, the first engineered to destabilise the intermediate and the second to stabilise it, thus slowing down or speeding up, respectively, product release. These results illustrate how product release by an enzyme can be facilitated by the presence of a metastable intermediate with transient weak interactions between the enzyme and product.

Single molecule recordings of lysozyme activity

Single molecule bioelectronic circuits provide an opportunity to study chemical kinetics and kinetic variability with bond-by-bond resolution. To demonstrate this approach, we examined the catalytic activity of T4 lysozyme processing peptidoglycan substrates. Monitoring a single lysozyme molecule through changes in a circuit's conductance helped elucidate unexplored and previously invisible aspects of lysozyme's catalytic mechanism and demonstrated lysozyme to be a processive enzyme governed by 9 independent time constants. The variation of each time constant with pH or substrate crosslinking provided different insights into catalytic activity and dynamic disorder. Overall, ten lysozyme variants were synthesized and tested in single molecule circuits to dissect the transduction of chemical activity into electronic signals. Measurements show that a single amino acid with the appropriate properties is sufficient for good signal generation, proving that the single molecule circuit technique can be easily extended to other proteins.

Dynamics-stability relationships in apo- and holomyoglobin: a combined neutron scattering and molecular dynamics simulations study

The removal of the heme group from myoglobin (Mb) results in a destabilization of the protein structure. The dynamic basis of the destabilization was followed by comparative measurements on holo- (holo-Mb) and apomyoglobin (apo-Mb). Mean-squared displacements (MSD) and protein resilience on the picosecond-to-nanosecond timescale were measured by elastic incoherent neutron scattering. Differences in thermodynamic parameters, MSD, and resilience were observed for both proteins. The resilience of holo-Mb was significantly lower than that of apo-Mb, indicating entropic stabilization by a higher degree of conformational sampling in the heme-bound folded protein. Molecular dynamics simulations provided site-specific information. Averaged over the whole structure, the molecular dynamics simulations yielded similar MSD and resilience values for the two proteins. The mobility of residues around the heme group in holo-Mb showed a smaller MSD and higher resilience compared to the same residue group in apo-Mb. It is of interest that in holo-Mb, higher MSD values are observed for the residues outside the heme pocket, indicating an entropic contribution to protein stabilization by heme binding, which is in agreement with experimental results.

THz SPECTROSCOPY OF PROTEINS IN WATER: DIRECT ABSORPTION AND CIRCULAR DICHROISM

Many of the functionally relevant collective vibrations of proteins and other biopolymers are expected to occur at terahertz frequencies. Precise absorption measurements combined with careful titration of biopolymers in water have allowed us to directly measure the terahertz absorption spectra associated with these motions, despite the strong background absorption of the solvent. We have also explored the circular dichroism spectroscopy of biomolecules over this same frequency range. Since circular dichroism requires the presence of net chirality in a molecule and chirality is present in nearly all biomaterial, it has the potential to capture the background free spectral features in biopolymers. To undertake these studies we have developed a broad band terahertz spectrometer suitable for both direct absorption and circular dichroism measurements of proteins in water between 0.75 – 3.72 THz. Direct terahertz absorption spectra of prototypical proteins bovine serum albumin (BSA) and hen egg white lysozyme have been documented and are described here. We have also successfully demonstrated the magnetic circular dichroism in semiconductors, and placed an upper bound on the terahertz circular dichroism signature of solvated BSA. In the terahertz frequency range, it appears that circular dichroism signatures are exceedingly small and detection remains a challenge.

Endocyclic cleavage in glycosides with 2,3-trans cyclic protecting groups

An endocyclic pathway is proposed as a reaction mechanism for the anomerization from the β (1,2-trans) to the α (1,2-cis) configuration observed in glycosides carrying 2,3-trans cyclic protecting groups. This reaction occurs in the presence of a weak Lewis or Brønsted acid, while endocyclic cleavage (endocleavage) in typical glycosides was observed only when mediated by protic media or strong Lewis acids. To rationalize the behavior of this class of compounds, the reaction mechanism and the promoting factors of the endocleavage are investigated using quantum-mechanical (QM) calculations and experimental studies. We examine anomerization reactions of thioglycosides carrying 2,3-trans cyclic protecting groups, employing boron trifluoride etherate (BF(3)·OEt(2)) as a Lewis acid. The estimated theoretical reactivity, based on a simple model to predict transition state (TS) energies from the strain caused by the fused rings, is very close to the TS energies calculated by the TS search along the C1-C2 bond rotation after the endo C-O bond breaking. Excellent agreement is found between the predicted TS energies and the experimental reactivity ranking. The series of calculations and experiments strongly supports the predominance of the endocyclic rather than the exocyclic mechanism. Furthermore, these investigations suggest that the inner strain is the primary factor enhancing the endocleavage reaction. The effect of the cyclic protecting group in restricting the pyranoside ring to a (4)C(1) conformation, extensively discussed in conjunction with the stereoelectronic effect theory, is shown to be a secondary factor.

Discovering conformational sub-states relevant to protein function

Background

Internal motions enable proteins to explore a range of conformations, even in the vicinity of native state. The role of conformational fluctuations in the designated function of a protein is widely debated. Emerging evidence suggests that sub-groups within the range of conformations (or sub-states) contain properties that may be functionally relevant. However, low populations in these sub-states and the transient nature of conformational transitions between these sub-states present significant challenges for their identification and characterization.

Methods and Findings

To overcome these challenges we have developed a new computational technique, quasi-anharmonic analysis (QAA). QAA utilizes higher-order statistics of protein motions to identify sub-states in the conformational landscape. Further, the focus on anharmonicity allows identification of conformational fluctuations that enable transitions between sub-states. QAA applied to equilibrium simulations of human ubiquitin and T4 lysozyme reveals functionally relevant sub-states and protein motions involved in molecular recognition. In combination with a reaction pathway sampling method, QAA characterizes conformational sub-states associated with cis/trans peptidyl-prolyl isomerization catalyzed by the enzyme cyclophilin A. In these three proteins, QAA allows identification of conformational sub-states, with critical structural and dynamical features relevant to protein function.

Conclusions

Overall, QAA provides a novel framework to intuitively understand the biophysical basis of conformational diversity and its relevance to protein function.

Stochastic molecular dynamics of peanut lectin PNA complex with T-antigen disaccharide

The stochastic boundary molecular dynamics simulation method was applied to investigate the structure of a complex comprised of a tetrameric peanut lectin and the tumour-associated disaccharide, Galbeta1-3GalNAc (T-antigen). Only a small region encompassing the active site was explicitly included in the calculations, but the electrical contribution of most outer atoms was taken into account by adding to the effective potential a term coming from an electrostatic potential grid that was pre-calculated and used to approximate the electrostatic energy and the force at any point in the interacting site. Results of simulating the intermolecular hydrogen bond network agree fairly well with X-ray experiments. An estimation of the direct and water-mediated interaction mean lifetimes and mean water residence times around the T-antigen oxygen atoms was computed over 400 ps. Monitoring the behaviour of water molecules within the active site revealed that there is a constant exchange of water with the bulk, especially in the proximity of ASN41, ASN127 and GLU129. The temporal evolution of the glycosidic linkage was also investigated and compared with simulations of T-antigen in solution. These studies of peanut lectins-sugar complexes clearly emphasize the importance of bound water molecules in generating carbohydrate specificity.

A helix-loop-helix peptide at the upper lip of the active site cleft of lysozyme confers potent antimicrobial activity with membrane permeabilization action

Recently, we have found that partially unfolded lysozyme exerts broad spectrum antimicrobial action in vitro against Gram-negative and Gram-positive bacteria independent of its catalytic activity. In parallel, an internal peptide (residues 98-112) of hen egg white lysozyme, obtained after digestion with clostripain, possessed broad spectrum antimicrobial action in vitro. This internal peptide is part of a helix-loop-helix domain (87-114 sequence of hen lysozyme) located at the upper lip of the active site cleft of lysozyme. The helix-loop-helix (HLH) structures are known motifs commonly found in membrane-active and DNA-binding proteins. To evaluate the contribution of the HLH peptide to the antimicrobial properties of lysozyme, the HLH sequence and its secondary structure derivatives of chicken and human lysozyme were synthesized and tested for antimicrobial activity against several bacterial strains. We found that the full HLH peptide of both chicken and human lysozymes was potently microbicidal against both Gram-positive and Gram-negative bacteria and the fungus Candida albicans. The N-terminal helix of HLH was specifically bactericidal to Gram-positive bacteria, whereas the C-terminal helix was bactericidal to all tested strains. Outer and inner membrane permeabilization studies, as well as measurements of transmembrane electrochemical potentials, provided evidence that HLH peptide and its C-terminal helix domain kill Gram-negative bacteria by crossing the outer membrane via self-promoted uptake and causing damage to the inner membrane through channel formation. The results are discussed in terms of proposed mechanisms for the catalytically independent antimicrobial activity of lysozyme that offer a new strategy for the design of potential antimicrobial drugs in the treatment of infectious diseases.

Molecular dynamics study of active-site interactions with tetracoordinate transients in acetylcholinesterase and its mutants

The role of active-site residues in the dealkylation reaction in the P(S)C(S) diastereomer of 2-(3,3-dimethylbutyl)methylphosphonofluoridate (soman)-inhibited Torpedo californica acetylcholinesterase (AChE) was investigated by full-scale molecular dynamics simulations using CHARMM: >400 ps equilibration was followed by 150-200 ps production runs with the fully solvated tetracoordinate phosphonate adduct of the wild-type, Trp84Ala and Gly199Gln mutants of AChE. Parallel simulations were carried out with the tetrahedral intermediate formed between serine-200 Ogamma of AChE and acetylcholine. We found that the NepsilonH in histidine H(+)-440 is positioned to protonate the oxygen in choline and thus promote its departure. In contrast, NepsilonH in histidine H(+)-440 is not aligned for a favourable proton transfer to the pinacolyl O to promote dealkylation, but electrostatic stabilization by histidine H(+)-440 of the developing anion on the phosphonate monoester occurs. Destabilizing interactions between residues and the alkyl fragment of the inhibitor enforce methyl migration from Cbeta to Calpha concerted with C-O bond breaking in soman-inhibited AChE. Tryptophan-84, phenyalanine-331 and glutamic acid-199 are within 3.7-3.9 A (1 A=10(-10) m) from a methyl group in Cbeta, 4.5-5.1 A from Cbeta and 4.8-5.8 A from Calpha, and can better stabilize the developing carbenium ion on Cbeta than on Calpha. The Trp84Ala mutation eliminates interactions between the incipient carbenium ion and the indole ring, but also reduces its interactions with phenylalanine-331 and aspartic acid-72. Tyrosine-130 promotes dealkylation by interacting with the indole ring of tryptophan-84. Glutamic acid-443 can influence the orientation of active-site residues through tyrosine-421, tyrosine-442 and histidine-440 in soman-inhibited AChE, and thus facilitate dealkylation.

Molecular dynamics simulations of human rhinovirus and an antiviral compound

The human rhinovirus 14 (HRV14) protomer, with or without the antiviral compound WIN 52084s, was simulated using molecular dynamics and rotational symmetry boundary conditions to model the effect of the entire icosahedral capsid. The protein asymmetrical unit, comprising four capsid proteins (VP1, VP2, VP3, and VP4) and two calcium ions, was solvated both on the exterior and the interior to fill the inside of the capsid. The stability of the simulations of this large system (~800 residues and 6,650 water molecules) is comparable to more conventional globular protein simulations. The influence of the antiviral compound on compressibility and positional fluctuations is reported. The compressibility, estimated from the density fluctuations in the region of the binding pocket, was found to be greater with WIN 52084s bound than without the drug, substantiating previous computations on reduced viral systems. An increase in compressibility correlates with an entropically more favorable system. In contrast to the increase in density fluctuations and compressibility, the positional fluctuations decreased dramatically for the external loops of VP1 and the N-terminus of VP3 when WIN 52084s is bound. Most of these VP1 and VP3 loops are found near the fivefold axis, a region whose mobility was not considered in reduced systems, but can be observed with this simulation of the full viral protomer. Altered loop flexibility is consistent with changes in proteolytic sensitivity observed experimentally. Moreover, decreased flexibility in these intraprotomeric loops is noteworthy since the externalization of VP4, part of VP1, and RNA during the uncoating process is thought to involve areas near the fivefold axis. Both the decrease in positional fluctuations at the fivefold axis and the increase in compressibility near the WIN pocket are discussed in relationship to the antiviral activity of stabilizing the virus against uncoating.

Protein modification during antiviral heat bioprocessing

Heat treatment is routinely used in the preparation of therapeutic protein biopharmaceuticals as a means of viral inactivation. However, in undertaking virucidal heat treatments, a balance must be found between the bioprocessing conditions, virus kill, and the maintenance of protein integrity. In this study, we utilize a simple model protein, hen egg-white lysozyme, to investigate the relationship between antiviral bioprocess conditions (protein formulation and temperature) and the extent and type of protein modification. A variety of industrially relevant wet- and dry-heat treatments were undertaken, using formulations that included sucrose as a thermostabilizing excipient. Although there was no evidence of lysozyme aggregation or crosslinking during any of the heat treatments, using liquid chromatography-electrospray ionization-mass spectroscopy (LC-ESI-MS) and peptide mapping we show that protein modifications do occur with increasingly harsh heat treatment. Modifications were predominantly found after wet-heat treatment, the major covalent modification of lysozyme under these conditions being glycation of Lys(97), by either glucose or fructose derived from hydrolyzed sucrose. The extent of sucrose hydrolysis was itself dependent on both the duration of heat treatment and formulation composition. Advanced glycation end products (AGEs) and additional unidentified products were also present in protein samples subjected to extended heat treatment. AGEs were derived primarily from initial glycation by fructose and not glucose. These findings have implications for the improvement of bioprocesses to ensure protein product quality.

Identification of a glutathione-S-transferase effector domain for inhibition of jun kinase, by molecular dynamics

We have recently found that the glutathione-S-transferase pi-isozyme (GST-pi), a cellular detoxification enzyme, potently and selectively inhibits activation of jun protein by its upstream kinase, jun kinase (JNK). This newly identified regulatory activity of GST-pi is strongly inhibited by a group of agents that inhibit its enzymatic activity. Since loss of enzymatic activity in general does not correlate with loss of regulatory activity, it is likely that inhibitor binding induces changes in the structure of one or more domains of GST that block its interaction with JNK. To identify regions of GST that change conformation on the binding of inhibitors, we have performed molecular dynamics calculations on GST-pi to compute its average structure in the presence and absence of the inhibitor, glutathione sulfonate. Superposition of the two average structures reveals that several regions change local structure depending upon whether the inhibitor is bound or not bound. Two of these regions, residues 36-50 and 194-201, are highly exposed. We have synthesized peptides corresponding to these two segments and find that the 194-201 sequence strongly inhibits the ability of GST-pi to block the in vitro phosphorylation of jun by JNK. These results suggest that this region of GST-pi is critical to its functioning as a newly discovered regulator of signal transduction.

Analysis of the relationship between enzyme activity and its internal motion using nuclear magnetic resonance: 15N relaxation studies of wild-type and mutant lysozyme

A mutant lysozyme where R14 and H15 are deleted together has higher activity and a similar binding ability to an inhibitor, trimer of N-acetylglucosamine ((NAG)3), compared with wild-type lysozyme. Since this has been attributed to intrinsic protein dynamic properties, we investigated the relationship between the activity and the internal motions of proteins. Backbone dynamics of the free and the complex forms with the (NAG)3 have been studied by measurement of the 15N T1 and T2 relaxation rates and NOE determinations at 600 MHz. Analysis of the data using the model-free formalism showed that the generalized order parameters (S2) were almost the same in wild-type and mutant lysozyme in unbound state, indicating that the mutation had little effect on the global internal motions. On the other hand, in the presence of (NAG)3, although some signals located around the active site were broadened or decreased in intensity because of strong perturbation by (NAG)3, there were several residues that showed increased or decreased backbone S2 in the complexed lysozymes. A comparison of the internal motions of the wild-type and mutant complexes showed a number of distinct dynamic differences between them. In particular, many residues located at or near active-site regions (turn 1, strand 2, turn 2 and long loop), displayed greater backbone dynamics reflecting the order parameter in mutant complex relative to mutant free. Furthermore, the Rex values at the loop C-D region, which was considered to be important for enzymatic activity, significantly increased. From these results, it was suggested that variations in the dynamics of these regions may play an important role in the enzyme activity.

Molecular dynamics simulation of alpha-lactalbumin and calcium binding c-type lysozyme

Alpha-lactalbumins (LAs) and c-type lysozymes (LYZs) are two classes of proteins which have a 35-40% sequence homology and share a common three dimensional fold but perform different functions. Lysozymes bind and cleave the glycosidic bond linkage in sugars, where as, alpha-lactalbumin does not bind sugar but participates in the synthesis of lactose. Alpha-lactalbumin is a metallo-protein and binds calcium, where as, only a few of the LYZs bind calcium. These proteins consist of two domains, an alpha-helical and a beta-strand domain, separated by a cleft. Calcium is bound at a loop situated at the bottom of the cleft and is important for the structural integrity of the protein. Calcium is an ubiquitous intracellular signal in higher eukaryotes and structural changes induced on calcium binding have been observed in a number of proteins. In the present study, molecular dynamics simulations of equine LYZ and human LA, with and without calcium, were carried out. We detail the differences in the dynamics of equine LYZ and human LA, and discuss it in the light of experimental data already available and relate it to the behavior of the functionally important regions of both the proteins. These simulations bring out the role of calcium in the conformation and dynamics of these metallo-proteins. In the calcium bound LA, the region of the protein around the calcium binding site is not only frozen but the atomic fluctuations are found to increase away from the binding site and peak at the exposed sites of the protein. This channeling of fluctuations away from the metal binding site could serve as a general mechanism by which the effect of metal binding at a site is transduced to other parts of the protein and could play a key role in protein-ligand and/or protein-protein interaction.

Conformational interconversions in peptide beta-turns: analysis of turns in proteins and computational estimates of barriers

The two most important beta-turn features in peptides and proteins are the type I and type II turns, which differ mainly in the orientation of the central peptide unit. Facile conformational interconversion is possible, in principle, by a flip of the central peptide unit. Homologous crystal structures afford an opportunity to structurally characterize both possible conformational states, thus allowing identification of sites that are potentially stereochemically mobile. A representative data set of 250 high-resolution (</=2.0 A), non-homologous protein crystal structures and corresponding variant and homologous entries, obtained from the Brookhaven Protein Data Bank, was examined to identify turns that are assigned different conformational types (type I/type II) in related structures. A total of 55 examples of beta-turns were identified as possible candidates for a stereochemically mobile site. Of the 55 examples, 45 could be classified as a potential site for interconversion between type I and type II beta-turns, while ten correspond to flips from type I' to type II' structures. As a further check, the temperature factors of the central peptide unit carbonyl oxygen atom of the 55 examples were examined. The analysis reveals that the turn assignments are indeed reliable. Examination of the secondary structures at the flanking positions of the flippable beta-turns reveals that seven examples occur in the loop region of beta-hairpins, indicating that the formation of ordered secondary structures on either side of the beta-turn does not preclude local conformational variations. In these beta-turns, Pro (11 examples), Lys (nine examples) and Ser (seven examples) were most often found at the i+1 position. Glycine was found to occur overwhelmingly at position i+2 (28 examples), while Ser (seven examples) and Asn (six examples) were amongst the most frequent residues. Activation energy barriers for the interconversion between type I and type II beta-turns were computed using the peptide models Ac-Pro-Aib-NHMe and Ac-Pro-Gly-NHMe within the framework of the AM1 semi-empirical molecular orbital procedure. In order to have a uniform basis for comparison and to eliminate the distracting influence of the deviation of backbone dihedral angles from that expected for ideal beta-turns, the dihedral angles phii+1 and psii+2 were fixed at the ideal values (phii+1=-60 degrees and psii+2=0 degrees). The other two angles (psii+1 and phii+2) were varied systematically to go from type II to type I beta-turn structures. The computational results suggest that there exists one stereospecific, concerted flip of the central peptide unit involving correlated single bond rotation that can occur with an activation barrier of the order of 3 kcal/mol. The results presented here suggest that conformational variations in beta-turns are observed in protein crystal structures and such changes may be an important dynamic feature in solution.

Energy landscape of a native protein: jumping-among-minima model

We have investigated energy landscape of human lysozyme in its native state by using principal component analysis and a model, jumping-among-minima (JAM) model. These analyses are applied to 1 nsec molecular dynamics trajectory of the protein in water. An assumption embodied in the JAM model allows us to divide protein motions into intra-substate and inter-substate motions. By examining intra-substate motions, it is shown that energy surfaces of individual conformational substates are nearly harmonic and mutually similar. As a result of principal component analysis and JAM model analysis, protein motions are shown to consist of three types of collective modes, multiply hierarchical modes, singly hierarchical modes, and harmonic modes. Multiply hierarchical modes, the number of which accounts only for 0.5% of all modes, dominate contributions to total mean-square atomic fluctuation. Inter-substate motions are observed only in a small-dimensional subspace spanned by the axes of multiplyhierarchical and singly hierarchical modes. Inter-substate motions have two notable time components: faster component seen within 200 psec and slower component. The former involves transitions among the conformational substates of the low-level hierarchy, whereas the latter involves transitions of the higher level substates observed along the first four multiply hierarchical modes. We also discuss dependence of the subspace, which contains conformational substates, on time duration of simulation.

Protein conformational changes determined by matrix-assisted laser desorption mass spectrometry

It is shown that simultaneously to the unfolding of hen egg white lysozyme and horse heart cytochrome c the sequential conformational changes and molten globule states can be detected by the combination of proteolysis and matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS). This is demonstrated by the differences among the products and the time courses of native lysozyme as well as those unfolded in 1 and 3 M guanidine hydrochloride (GuHCl) when they were proteolyzed by proteinase K and analyzed by MALDI-MS. Due to the absence of disulfide bonds in the cytochrome c molecule, it is more sensitive to the disturbance of the denaturant. The partially unfolded state as detected at low concentrations of guanidine hydrochloride in our experiment resemble the molten globule state. One of the unique properties of the method described herein is to measure directly the peptide fragment liberated from proteolysis of the protein. It allows the identification of the sensitive sites susceptible to denaturation, which are subsequently cleaved by proteinase K proteolysis.

Locally accessible conformations of proteins: multiple molecular dynamics simulations of crambin

Multiple molecular dynamics (MD) simulations of crambin with different initial atomic velocities are used to sample conformations in the vicinity of the native structure. Individual trajectories of length up to 5 ns sample only a fraction of the conformational distribution generated by ten independent 120 ps trajectories at 300 K. The backbone atom conformational space distribution is analyzed using principal components analysis (PCA). Four different major conformational regions are found. In general, a trajectory samples only one region and few transitions between the regions are observed. Consequently, the averages of structural and dynamic properties over the ten trajectories differ significantly from those obtained from individual trajectories. The nature of the conformational sampling has important consequences for the utilization of MD simulations for a wide range of problems, such as comparisons with X-ray or NMR data. The overall average structure is significantly closer to the X-ray structure than any of the individual trajectory average structures. The high frequency (less than 10 ps) atomic fluctuations from the ten trajectories tend to be similar, but the lower frequency (100 ps) motions are different. To improve conformational sampling in molecular dynamics simulations of proteins, as in nucleic acids, multiple trajectories with different initial conditions should be used rather than a single long trajectory.

Domain motions in dihydrofolate reductase: a molecular dynamics study

Molecular dynamics simulations have been carried out on the enzyme dihydrofolate reductase from Lactobacillus casei complexed with methotrexate, NADPH and 264 crystallographic water molecules. Analysis of correlations in atomic fluctuations reveal the presence of highly correlated motion (correlation coefficient > 0.6) in the region between residues 30 to 35 and 85 to 90 leading to the identification of two domains, an "adenosine-binding domain" and a "large domain", which rotate by 3 to 4 degrees with respect to each other. The strongest correlation (> 0.6) within the large domain involves a coupling between the motions of the "teen-loop", and the spatially contiguous loops linking beta 6-beta 7 and beta 7-beta 8. Moreover, there is a significant correlation (approximately 0.5) between the adenosine fragment of NADPH and the pteridine and p-aminobenzoyl fragments of methotrexate, which are separated by approximately 17 A, and is lost on removal of "rigid-body" motion from the original trajectory. This provides support for the idea that the relative motion of the two domains is a means by which the occupation of the binding site for the adenosine end of the coenzyme can affect methotrexate binding and vice versa. Quasiharmonic vibrational analysis of the trajectory reveals that the overall dynamics of the system are governed by domain motions whose contributions are dominant at low frequencies. In addition, different low-frequency modes are responsible for separately coupling the adenosine-binding site and parts of methotrexate.

Lysozyme: a model enzyme in protein crystallography

The review concentrates on the crystal structure results from several protein crystallography laboratories on three different lysozymes, the type-c lysozymes such as hen egg-white lysozyme (HEWL), the type-g lysozyme, such as goose egg-white lysozyme (GEWL), and the lysozyme from T4 bacteriophage (T4L). The crystallographic studies on HEWL in several different crystal forms have shown that the lysozyme molecule is relatively rigid with the residues of the active site Glu35 and Asp52 adopting almost identical conformations in all structures and species variants. The NMR results also confirm the presence of a similar conformation of HEWL in solution. All three enzymes, HEWL, GEWL and T4L are composed of two domains, one that is predominantly alpha-helical and a smaller domain that is mainly beta-sheet in nature. The general acid/general base residue in each lysozyme (Glu35 in HEWL, Glu73 in GEWL and Glu11 in T4L) is contributed by the larger alpha-helical domain. The beta-sheet domains of HEWL and T4L contribute an aspartate to their respective active sites, which is likely involved in electrostatic stabilization of the oxycarbonium ion intermediate of the site D sugar on the hydrolytic pathway of oligosaccharides. There is no analogous aspartate carboxylate group in GEWL although minor conformational changes could position one or other of Asp86 or Asp97 for such a stabilization role. The binding of substrate analogues, transition state mimics and oligosaccharide products of hydrolysis to HEWL, GEWL and T4L have contributed greatly to our understanding of sugar binding to proteins. The observed subtle conformational differences of the free vs bound forms of these enzymes are best described by a narrowing of the active site clefts in the presence of the inhibitors. Details of the binding interactions of those residues lining the oligosaccharide binding clefts of the three-enzymes HEWL, GEWL and T4L with the sugar residues in sites A, B, C and D are presented and discussed. Oligosaccharides of (GlcNAc)n and alternating MurNAc-GlcNAc-MurNAc have been bound to these three enzymes and the structures determined at high resolution. These binding studies have contributed greatly to our understanding of the catalytic mechanism of the lysozyme glycosidase activity. The currently accepted view of this mechanism is presented and discussed in this review.

Simulations of lysozyme: internal motions and the reaction mechanism

Simulations of lysozyme by molecular dynamics have greatly increased our understanding of this enzyme. It has been shown how the internal motions are related to the structural elements (helices, sheets and loops) of the molecule. Comparisons of the motions in the free and substrate-bound form reveal that most are similar but that there are significant differences. Comparisons of the theoretical results with X-ray and nuclear magnetic resource data show good agreement. The hinge-bending motion, which opens and closes the binding site cleft between the two domains, is correlated with substrate binding. From analysis of simulations with bound substrate, an alternative mechanism for oligoglycoside hydrolysis was proposed. It involves cleavage of an endocyclic C-O bond, instead of the exocyclic cleavage proposed in the standard mechanism. Both mechanisms have been demonstrated in solution, but it is still unclear which is prevalent in lysozyme.

Zinc binding in proteins and solution: a simple but accurate nonbonded representation

Force field parameters that use a combination of Lennard-Jones and electrostatic interactions are developed for divalent zinc and tested in solution and protein simulations. It is shown that the parameter set gives free energies of solution in good agreement with experiment. Molecular dynamics simulations of carboxypeptidase A and carbonic anhydrase are performed with these zinc parameters and the CHARMM 22 beta all-atom parameter set. The structural results are as accurate as those obtained in published simulations that use specifically bonded models for the zinc ion and the AMBER force field. The inclusion of longer-range electrostatic interactions by use of the Extended Electrostatics model is found to improve the equilibrium conformation of the active site It is concluded that the present parameter set, which permits different coordination geometries and ligand exchange for the zinc ion, can be employed effectively for both solution and protein simulations of zinc-containing systems.

Structures of partridge egg-white lysozyme with and without tri-N-acetylchitotriose inhibitor at 1.9 A resolution

The three-dimensional structures of native partridge egg-white lysozyme (PEWL) and PEWL complexed with tri-N-acetylchitotriose inhibitor have been determined crystallographically and refined at 1.9 A resolution. Crystals of native and complexed protein are isomorphous and have space group and cell dimensions that are identical to those of hen egg-white lysozyme (HEWL) under similar crystallization conditions. Full occupancy of the trisaccharide in the inhibitor complex has allowed definitive modeling and refinement of all three sugar residues, located at subsites A, B, and C in the PEWL active site. A comparison has been made with HEWL/inhibitor complexes in which coordinates were either not refined (Blake CCF, et al., 1967, Proc R Soc B 167:378-388) or were refined at partial occupancy (Cheetham JC, Artymiuk PJ, Phillips DC, 1992, J Mol Biol 224:613-628). Although the loop comprising residues 70-75 is located on the surface of the protein and not near the active site, it appears to be affected indirectly by trisaccharide binding such that the loop shifts toward the active site and becomes relatively immobilized. The source of this loop movement appears to be the anchoring of Trp62, located in the active site cleft, as it forms a hydrogen bond with O6 of the N-acetylglucosamine at site C. Good electron density for the trisaccharide in the PEWL complex structure shows that Asp 101 is involved in hydrogen bonding interactions with the terminal sugar residue.

Protein folding dynamics: the diffusion-collision model and experimental data

The diffusion-collision model of protein folding is assessed. A description is given of the qualitative aspects and quantitative results of the diffusion-collision model and their relation to available experimental data. We consider alternative mechanisms for folding and point out their relationship to the diffusion-collision model. We show that the diffusion-collision model is supported by a growing body of experimental and theoretical evidence, and we outline future directions for developing the model and its applications.

A structural basis for the interaction of urea with lysozyme

The effect of urea on the crystal structure of hen egg-white lysozyme has been investigated using X-ray crystallography. High resolution structures have been determined from crystals grown in the presence of 0, 0.7, 2, 3, 4, and 5 M urea and from crystals soaked in 9 M urea. All the forms are essentially isomorphous with the native type II crystals, and the derived structures exhibit excellent geometry and RMS differences from ideality in bond distances and angles. Comparison of the urea complex structures with the native enzyme (type II form, at 1.5 A resolution) indicates that the effect of urea is minimal over the concentration range studied. The mean difference in backbone conformation between the native enzyme and its urea complexes varies from 0.18 to 0.49 A. Conformational changes are limited to flexible surface loops (Thr 69-Asn 74, Ser 100-Asn 103), the active site loop (Asn 59-Cys 80), and the C-terminus (Cys 127-Leu 129). Urea molecules are bound to distinct sites on the surface of the protein. One molecule is bound to the active site cleft's C subsite, at all concentrations, in a fashion analogous to that of the N-acetyl substituent of substrate and inhibitor sugars normally bound to this site. Occupation of this subsite by urea alone does not appear to induce the conformational changes associated with inhibitor binding.

Thermal expansion of hen egg-white lysozyme. Comparison of the 1.9 A resolution structures of the tetragonal form of the enzyme at 100 K and 298 K

The average structural and dynamic properties of tetragonal hen egg-white lysozyme have been compared, in structures refined at 1.9 A resolution, using data collected at 100 K and 298 K. The molecule expands by 1.8% over this temperature range with the expansion occurring primarily in its small sub-atomic-sized spaces in an anisotropic manner. Hen egg-white lysozyme consists of two domains: domain 1 (residues 40 to 88) is composed primarily of beta-sheet and is observed to expand by only 0.3%; domain 2 (residues 1 to 39 and 89 to 129) is chiefly alpha-helix and is observed to expand by 2.2%. This is consistent with previous observations that proteins composed primarily of alpha-helix expand more with temperature than do those composed primarily of beta-sheet. The largest movement in the molecule is undergone by the two domains of the structure that move further apart as the temperature is raised. This motion is not a cleft opening but rather consists of a tilt by 2.3 degrees of domain 1 away from domain 2. Within the individual domains the largest movement is undergone by loop T1 of domain 2, consisting of residues 17 to 23. This loop moves in the opposite direction to the rest of the molecule as the temperature is raised. Average temperature factors for the room-temperature and low-temperature structures are 15.2 A2 and 8.1 A2, respectively, when all protein atoms are considered, while these values are 14.0 A2 and 7.8 A2, when only main-chain atoms (N, C alpha, C) are taken into account. An examination of the main-chain averaged B-factor per residue shows that residues involved in intermolecular protein-protein contacts, with symmetry-related molecules, have somewhat lower B-factors than the average and undergo smaller than average changes in B-factor as the temperature is lowered.

Molecular dynamics of the H-ras gene-encoded p21 protein; identification of flexible regions and possible effector domains

We previously reported a complete computer-based three-dimensional structure for residues 1-171 of the Gly 12-containing ras-gene-encoded p21 protein complexed with GDP. This structure was subsequently shown to closely agree with a high-resolution x-ray crystallographic structure of p21. In this communication, we report a molecular dynamics stimulation of the modelled structure in an explicit shell of water molecules to identify domains within the protein that are unusually flexible. These domains represent regions which are most likely to undergo important conformational changes when the protein is activated by binding to GTP or by oncogenic amino acid substitutions such as Val for Gly 12. The starting structure was surrounded with water molecules, temperature-equilibrated and then followed over a 100 ps trajectory during which time the energy converged after about 50 ps. Regions of the protein that were found to have the largest coordinate fluctuations involved residues 12-16, 30-35, 40-52, 60-73, 85-89, 101-109, 119-123, and 127-131. Many of these sequences with high flexibility have been implicated in the functioning of this protein. Since the overall largest fluctuations were observed for residues 101-106 and 119-123, p21 peptides containing these residues (96-110 and 115-126) were synthesized and were found to inhibit strongly the effects of oncogenic p21 protein in an oocyte maturation assay. These results indicate that the flexible p21 sequences may constitute critical functional domains of the activated protein and that this general approach may be useful for identification of important functional domains in proteins.

Influence of protein flexibility on the redox potential of rubredoxin: energy minimization studies

A theoretical investigation of the protein contribution to the redox potential of the iron-sulfur protein rubredoxin is presented. Structures of the oxidized and reduced forms of the protein were obtained by energy minimizing the oxidized crystal structure of Clostridium pasteurianum rubredoxin with appropriate charges and parameters. By including 102 crystal waters, structures close to the original crystal structure were obtained (rms difference of 1.16 A), even with extensive minimization, thus allowing accurate calculations of comparative energies. Our calculations indicate an energy change of about -60 kcal/mol (2.58 eV) in the protein alone upon reduction. This energy change was due to both the change in charge of the redox site and the subsequent relaxation of the protein. An energy minimization procedure for the relaxation gives rms differences between the oxidized and reduced states of about 0.2 A. The changes were small and occurred in both the backbone and sidechain mainly near the Fe-S center but contributed about -16 kcal/mol (0.69 eV) to the total protein contribution. Although the neglect of certain effects such as electronic polarization may make the relaxation energies calculated an upper limit, the results indicate that protein relaxation contributes substantially to the redox potential.

1H-NMR analysis of turkey egg-white lysozyme and comparison with hen egg-white lysozyme

The complete main chain and approximately 75% of the side chain 1H-NMR assignments of the 129-residue protein, turkey egg-white lysozyme, are presented. NOE data, hydrogen-exchange rates, chemical shifts and coupling constants are reported and are indicative of a structure in solution that is essentially identical to that of the homologous hen egg-white lysozyme. The NH-alpha CH coupling constants of turkey lysozyme are compared to torsion-angle data from three crystal structures of the protein and the results are interpreted in terms of crystal-structure resolution and refinement.

Molecular dynamics of sickle and normal hemoglobins

Molecular dynamics (MD) simulations have been carried out for 62.5 ps on crystal structures of deoxy sickle cell hemoglobin (HbS) and normal deoxy hemoglobin (HbA) using the CHARMM MD algorithm, with a time step of 0.001 ps. In the trajectory analysis of the 12.5-62.5 (50 ps) simulation, oscillations of the radius of gyration and solvent-accessible surface area were calculated. HbS exhibited a general contraction during the simulation, while HbA exhibited a nearly constant size. The average deviations of simulated structures from the starting structures were found to be 1.8 A for HbA and 2.3 A for HbS. The average rms amplitudes of atomic motions (atomic flexibility) were about 0.7 A HbA and about 1.0 A for HbS. The amplitudes of backbone motion correlate well with temperature factors derived from x-ray crystallography. A comparison of flexibility between the alpha- and beta-chains in both HbA and HbS indicates that the beta-chains generally exhibited greater flexibility than the alpha-chains, and that the HbS beta-chains exhibit greater flexibility in the N-terminal and D- and F-helix regions than do those of HbA. The average amplitude of backbone torsional oscillations was about 9 degrees, a value comparable with that of other simulations, with enhanced torsional oscillation occurring primarily at the ends of helices or in loop regions between helices. Comparison of atomic flexibility and torsional oscillation results suggests that the increased beta-chain flexibility results from relatively concerted motions of secondary structure elements. The increased flexibility may play an important role in HbS polymerization.(ABSTRACT TRUNCATED AT 250 WORDS)

A 500 ps molecular dynamics simulation study of interleukin-1 beta in water. Correlation with nuclear magnetic resonance spectroscopy and crystallography

We report the results of a 500 ps molecular dynamics simulation of the cytokine interleukin-1 beta, a protein of 153 amino acids, immersed in a sphere of 3783 bulk water molecules with a radius of 33 A. The simulation reproduces the amplitudes of the fast librational motions of the backbone N-H bonds determined from 15N nuclear magnetic relaxation data, as well as the crystallographic B-factors. Moreover, this study suggests a molecular picture of the nature of the slow internal motions that have been inferred from nuclear magnetic resonance relaxation experiments. These experiments indicated that, in addition to fast motions common to all residues, 32 surface residues exhibit slow motions on the 400 ps to 5 ns time-scale. While the present simulation is not sufficiently long to provide a quantitative description of events on this time-scale, it is long enough to observe several large amplitude transitions that are likely candidates for these slow motions. Specifically, in many of these 32 residues, the N-H groups are hydrogen bonded and infrequent dihedral transitions cause the N-H vectors to jump between states with well-defined orientations. It is shown that the time course of the angular reorientational correlation functions of these residues calculated from the trajectory is a reflection of the random times at which these infrequent jumps happen to have occurred. Thus, while the rate of these transitions cannot be quantified, the simulated decay of these correlation functions is completely consistent with the physical picture in which the N-H vectors, in addition to fast librational motion, undergo large amplitude jumps between conformations stabilized by hydrogen bonds.

Protein surface topology-probing by selective chemical modification and mass spectrometric peptide mapping

Aminoacetylation of lysine residues and the modification of arginine by 1,2-cyclohexanedione to N7,N8-(dihydroxy-1,2-cyclohexylidene)arginine were used for probing the surface topology of hen-eggwhite lysozyme as a model protein. The molecular identification of lysine and arginine modification sites was provided by molecular weight determinations of modified and unmodified tryptic peptide mixtures (peptide mapping) using 252Cf plasma desorption mass spectrometry. At conditions of limited chemical modification, mass-spectrometric peptide-mapping analyses of lysozyme derivatives enabled the direct assignment of relative reactivities of lysine and arginine residues at different reaction times and reagent concentrations. The relative reactivities of lysine residues showed a direct correlation with their surface accessibilities from x-ray structure data. For the reaction with 1,2-cyclohexanedione, a selective modification at Arg-5, -125, -112, and -73 was identified, and an inverse correlation of relative reactivities with the surface accessibility ratios of the N7- and the N8-guanidino functions was obtained. By examination of the x-ray structural data of lysozyme, this selective modification was attributed to intramolecular catalysis because of the presence of neighboring proton acceptor groups, such as the Asp-119 carboxylate group for Arg-125 and the Trp-123 and Arg-125 carbonyl groups for Arg-5.

Internal motional averaging and three-dimensional structure determination by nuclear magnetic resonance

Dynamic averaging effects from internal motions on interproton distances estimated from nuclear Overhauser effects (NOE) are determined by using a molecular dynamics simulation of lysozyme. Generalized order parameters measuring angular averaging and radial averaging parameters are calculated. The product of these two parameters describes the full averaging effects on cross-relaxation. Analysis of 2778 non-methyl NOE interactions from the protein interior and surface indicates that distances estimated by assuming a rigid molecule have less than 10% error for 89% of the NOE interactions. However, analysis of 1854 methyl interactions found that only 68% of the distances estimated from cross-relaxation rates would have less than 10% error. Qualitative evaluation of distances according to strong, medium and weak NOE intensities, when used to define only the upper bound for interproton separation, would misassign less than 1% of the distance constraints because of motional averaging. Internal motions do not obscure the identification of secondary structure, although some instances of significant averaging effects were found for interactions in alpha-helical regions. Interresidue NOEs for amino acids more than three residues apart in the primary sequence are more extensively averaged than intraresidue or short-range interresidue NOEs. Intraresidue interactions exhibit a greater degree of angular averaging than those involving interresidue proton pairs. An internal motion does not equally affect all NOE interactions for a particular proton. Thus, incorporation of averaging parameters in nuclear magnetic resonance structure determination procedures must be made on a proton-pair-wise basis. On the basis of the motional averaging results, particular fixed-distance proton pairs in proteins are suggested for use as distance references. A small percentage of NOE pairs localized to three regions of the protein exhibit extreme averaging effects from internal motions. The regions and types of motions involved are described.

Refinement of an enzyme complex with inhibitor bound at partial occupancy. Hen egg-white lysozyme and tri-N-acetylchitotriose at 1.75 A resolution

The structure of the tri-N-acetylchitotriose inhibitor complex of hen egg-white lysozyme has been refined at 1.75 A resolution, using data collected from a complex crystal with ligand bound at less than full occupancy. To determine the exact value of the inhibitor occupancy, a model comprising unliganded and sugar-bound protein molecules was generated and refined against the 1.75 A data, using a modified version of the Hendrickson & Konnert least-squares procedure. The crystallographic R-factor for the model was found to fall to a minimum at 55% bound sugar. Conventional refinement assuming unit occupancy was found to yield incorrect thermal and positional parameters. Application of the same refinement procedures to an earlier 2.0 A data set, collected independently on different complex crystals by Blake et al. gave less consistent results than the 1.75 A refinement. From an analysis of the high resolution structure a detailed picture of the protein-carbohydrate interactions in the non-productive complex has emerged, together with the conformation and mobility changes that accompany ligand binding. The specificity of interaction between the protein and inhibitor, bound in subsites A to C of the active site, is seen to be generated primarily by an extensive network of hydrogen bonds, both to the protein itself and to bound solvent molecules. The latter also play an important role in maintaining the structural integrity of the active site cleft in the apo-protein.

A method for determining the positions of polar hydrogens added to a protein structure that maximizes protein hydrogen bonding

An automated method for the optimal placement of polar hydrogens in a protein structure is described. This method treats the polar, side chain hydrogens of lysine, serine, threonine, and tyrosine and the amino terminus of a protein. The program, called NETWORK, divides the potential hydrogen-bonding pairs of a protein into groups of interacting donors and acceptors. A search is conducted on each of the local groups to find an arrangement which forms the most hydrogen bonds. If two or more arrangements have the same number of hydrogen bonds, the arrangement with the shortest set of hydrogen bonds is selected. The polar hydrogens of the histidyl side chain are specifically treated, and the ionization state of this residue is allowed to change, if this change results in additional hydrogen bonds for the local group. The program will accept Protein Data Bank as well as Biosym-format coordinate files. Input and output routines can be easily modified to accept other coordinate file formats. The predictions from this method are compared to known hydrogen positions for bovine pancreatic trypsin inhibitor, insulin, RNase-A, and trypsin for which the neutron diffraction structures have been determined. The usefulness of this program is further demonstrated by a comparison of molecular dynamics simulations for the enzyme cytochrome P-450cam with and without using NETWORK.