Insight on the interaction between the scorpion toxin blocker Discrepin on potassium voltage-gated channel Kv4.3 by molecular dynamics simulations.

Discrepin is a 38-residue α-toxin extracted from the venom of the Venezuelan scorpion Tityus discrepans, which inhibits ionic transit in the voltage-dependent potassium channels (Kv) of A-type current. The effect of specific residues on the IC50 between Discrepine and Kv4.3, the main component of A-type currents, is known; however, the molecular details of the toxin-channel interaction are not known. In this work, we present interaction models between Discrepin (wt) and two peptide variants (V6K/D20K and K13A) on the pore-forming domain of the Kv4.3 channel obtained from homology, docking, and molecular dynamics modeling techniques. The free energy calculations in these models correspond to the order of the experimentally determined IC50 values. Our studies shed light on the role of the K13 residue as responsible for occluding the Kv4.3 selectivity filter and the importance of the V6K mutation in the approach and stabilization of toxin-channel complex interactions.Communicated by Ramaswamy H. Sarma.

Calcium inhibits penetration of Alzheimer's Aß1 -42 monomers into the membrane.

Calcium ion regulation plays a crucial role in maintaining neuronal functions such as neurotransmitter release and synaptic plasticity. Copper (Cu2+ ) coordination to amyloid-ß (Aß) has accelerated Aß1-42 aggregation that can trigger calcium dysregulation by enhancing the influx of calcium ions by extensive perturbing integrity of the membranes. Aß1-42 aggregation, calcium dysregulation, and membrane damage are Alzheimer disease (AD) implications. To gain a detail of calcium ions' role in the full-length Aß1-42 and Aß1-42 -Cu2+ monomers contact, the cellular membrane before their aggregation to elucidate the neurotoxicity mechanism, we carried out 2.5 µs extensive molecular dynamics simulation (MD) to rigorous explorations of the intriguing feature of the Aß1-42 and Aß1-42 -Cu2+ interaction with the dimyristoylphosphatidylcholine (DMPC) bilayer in the presence of calcium ions. The outcome of the results compared to the same simulations without calcium ions. We surprisingly noted robust binding energies between the Aß1-42 and membrane observed in simulations containing without calcium ions and is two and a half fold lesser in the simulation with calcium ions. Therefore, in the case of the absence of calcium ions, N-terminal residues of Aß1-42 deeply penetrate from the surface to the center of the bilayer; in contrast to calcium ions presence, the N- and C-terminal residues are involved only in surface contacts through binding phosphate moieties. On the other hand, Aß1-42 -Cu2+ actively participated in surface bilayer contacts in the absence of calcium ions. These contacts are prevented by forming a calcium bridge between Aß1-42 -Cu2+ and the DMPC bilayer in the case of calcium ions presence. In a nutshell, Calcium ions do not allow Aß1-42 penetration into the membranes nor contact of Aß1-42 -Cu2+ with the membranes. These pieces of information imply that the calcium ions mediate the membrane perturbation via the monomer interactions but do not damage the membrane; they agree with the western blot experimental results of a higher concentration of calcium ions inhibit the membrane pore formation by Aß peptides.

Imidazole and nitroimidazole derivatives as NADH-fumarate reductase inhibitors: Density functional theory studies, homology modeling, and molecular docking.

Chagas disease is caused by Trypanosoma cruzi. Benznidazole and nifurtimox are drugs used for its therapy; nevertheless, they have collateral effects. NADH-fumarate (FUM) reductase is a potential pharmacological target since it is essential for survival of parasite and is not found in humans. The objectives are to design and characterize the electronic structure of imidazole and nitroimidazole derivatives at DFT-M06-2X level in aqueous solution; also, to model the NADH-FUM reductase and analyze its intermolecular interactions by molecular docking. Quantum-chemical descriptors allowed to select the molecules with the best physicochemical properties and lowest toxicity. A high-quality three-dimensional structure of NADH-FUM reductase was obtained by homology modeling. Water molecules do not have influence in the interaction between FUM and NADH-FUM reductase. The main hydrogen-binding interactions for FUM were identified in NADH, Lys172, and Arg89; while hydrophobic interactions in Phe479, Thr174, Met63. The molecules S3-8, S2-8, and S1-8 could be inhibitors of NADH-FUM reductase.

Membrane fluidity, composition, and charge affect the activity and selectivity of the AMP ascaphin-8.

Ascaphins are cationic antimicrobial peptides that have been shown to have potential in the treatment of infectious diseases caused by multidrug-resistant pathogens (MDR). However, to date, their principal molecular target and mechanism of action are unknown. Results from peptide prediction software and molecular dynamics simulations confirmed that ascaphin-8 is an alpha-helical peptide. For the first time, the peptide was described as membranotrophic using biophysical approaches including calcein liposome leakage, Laurdan general polarization, and dynamic light scattering. Ascaphin-8's activity and selectivity were modulated by rearranging the spatial distribution of lysine (Var-K5), aspartic acid (Var-D4) residues, or substitution of phenylalanine with tyrosine (Var-Y). The parental peptide and its variants presented high affinity toward the bacterial membrane model (≤2 µM), but lost activity in sterol-enriched membranes (mammal and fungal models, with cholesterol and ergosterol, respectively). The peptide-induced pore size was estimated to be >20 nm in the bacterial model, with no difference among peptides. The same pattern was observed in membrane fluidity (general polarization) assays, where all peptides reduced membrane fluidity of the bacterial model but not in the models containing sterols. The peptides also showed high activity toward MDR bacteria. Moreover, peptide sensitivity of the artificial membrane models compared with pathogenic bacterial isolates were in good agreement.

Computational studies of membrane pore formation induced by Pin2.

Understanding, at the molecular level, the effect of AMPs on biological membranes is of crucial importance given the increasing number of multidrug-resistant bacteria. Being part of an ancient type of innate immunity system, AMPs have emerged as a potential solution for which bacteria have not developed resistance. Traditional antibiotics specifically act on biosynthetic pathways, while AMPs may directly destabilize the lipid membrane, but it is unclear how AMPs affect the membrane's stability. We performed multiscale molecular dynamics simulations to investigate the structural features leading to membrane pores formation on zwitterionic and anionic membranes by the antimicrobial peptide (AMP) Pandinin 2 (Pin2). Some experimental reports propose that Pin2 could form barrel-stave pores, while others suggest that it could form toroidal pores. Since there is no conclusive evidence of which type of pore is formed by Pin2 on bilayers, performing molecular dynamics simulations on these systems could shed some light on whether or not or what type of pore Pin2 forms on model membranes. Our results are focused on a detailed description of the pore formation by Pin2 in POPC and POPE:POPG membranes., which strongly suggest that Pin2 forms a toroidal pore and not a barrel-shaped pore; this type of pore also affects the membrane properties. In the process, a phospholipid remodeling in the POPE:POPG membrane takes place. Moreover, the pores formed by Pin2 indicate that they are selective for the chlorine ion. There are no previous ion selectivity reports for other AMPs with similar physicochemical properties, such as melittin and magainin.Communicated by Ramaswamy H. Sarma.

An Overview of Several Inhibitors for Alzheimer's Disease: Characterization and Failure.

Amyloid beta (Aß) oligomers are the most neurotoxic aggregates causing neuronal death and cognitive damage. A detailed elucidation of the aggregation pathways from oligomers to fibril formation is crucial to develop therapeutic strategies for Alzheimer's disease (AD). Although experimental techniques rely on the measure of time- and space-average properties, they face severe difficulties in the investigation of Aß peptide aggregation due to their intrinsically disorder character. Computer simulation is a tool that allows tracing the molecular motion of molecules; hence it complements Aß experiments, as it allows to explore the binding mechanism between metal ions and Aß oligomers close to the cellular membrane at the atomic resolution. In this context, integrated studies of experiments and computer simulations can assist in mapping the complete pathways of aggregation and toxicity of Aß peptides. Aß oligomers are disordered proteins, and due to a rapid exploration of their intrinsic conformational space in real-time, they are challenging therapeutic targets. Therefore, no good drug candidate could have been identified for clinical use. Our previous investigations identified two small molecules, M30 (2-Octahydroisoquinolin-2(1H)-ylethanamine) and Gabapentin, capable of Aß binding and inhibiting molecular aggregation, synaptotoxicity, intracellular calcium signaling, cellular toxicity and memory losses induced by Aß. Thus, we recommend these molecules as novel candidates to assist anti-AD drug discovery in the near future. This review discusses the most recent research investigations about the Aß dynamics in water, close contact with cell membranes, and several therapeutic strategies to remove plaque formation.

Insights into the Molecular Inhibition of the Oncogenic Channel KV10.1 by Globular Toxins.

Inhibition of the expression of the human ether-à-go-go (hEAG1 or hKV10.1) channel is associated with a dramatic reduction in the growth of several cancerous tumors. The modulation of this channel's activity is a promising target for the development of new anticancer drugs. Although some small molecules have shown inhibitory activity against KV10.1, their lack of specificity has prevented their use in humans. In vitro studies have recently identified a limited number of peptide toxins with proven specificity in their hKV10.1 channel inhibitory effect. These peptide toxins have become desirable candidates to use as lead compounds to design more potent and specific hKV10.1 inhibitors. However, the currently available studies lack the atomic resolution needed to characterize the molecular features that favor their binding to hKV10.1. In this work, we present the first attempt to locate the possible hKV10.1 binding sites of the animal peptide toxins APETx4, Aa1a, Ap1a, and k-hefutoxin 1, all of which described as hKV10.1 inhibitors. Our studies incorporated homology modeling to construct a robust three-dimensional (3D) model of hKV10.1, applied protein docking, and multiscale molecular dynamics techniques to reveal in atomic resolution the toxin-channel interactions. Our approach suggests that some peptide toxins bind in the outer vestibule surrounding the pore of hKV10.1; it also identified the channel residues Met397 and Asp398 as possible anchors that stabilize the binding of the evaluated toxins. Finally, a description of the possible mechanism for inhibition and gating is presented.

Estimation of pore dimensions in lipid membranes induced by peptides and other biomolecules: A review.

The cytoplasmic membrane is one of the most frequent cell targets of antimicrobial peptides (AMPs) and other biomolecules. Understanding the mechanism of action of AMPs at the molecular level is of utmost importance for designing of new membrane-specific molecules. In particular, the formation of pores, the structure and size of these pores are of great interest and require nanoscale resolution approaches, therefore, biophysical strategies are essential to achieve an understanding of these processes at this scale. In the case of membrane active peptides, pore formation or general membrane disruption is usually the last step before cell death, and so, pore size is generally directly associated to pore structure and stability and loss of cellular homeostasis, implicated in overall peptide activity. Up to date, there has not been a critical review discussing the methods that can be used specifically for estimating the pore dimensions induced by membrane active peptides. In this review we discuss the scope, relevance and popularity of the different biophysical techniques such as liposome leakage experiments, advanced microscopy, neutron or X-ray scattering, electrophysiological techniques and molecular dynamics studies, all of them useful for determining pore structure and dimension.

Folding profiles of antimicrobial scorpion venom-derived peptides on hydrophobic surfaces: a molecular dynamics study.

Most helical antimicrobial peptides (AMPs) are usually unfolded in aqueous solution; however they acquire their secondary structure in the presence of a hydrophobic environment such as lipid membranes. Being the biological membranes the main target of many AMPs it is necessary to understand their way of action. Pandinin 2 (Pin2) is an alpha-helical AMP isolated from the venom of the African scorpion Pandinus imperator which shows high antimicrobial activity against Gram-positive bacteria and it is less active against Gram-negative bacteria, nevertheless, it has strong hemolytic activity. Its chemically synthesized Pin2GVG analog has low hemolytic activity while keeping its antimicrobial activity. With the aim of exploring the partition and subsequent folding of these peptides, in this work we report the results of extensive molecular dynamics simulations of Pin2 and Pin2GVG peptides in the presence of 2 hydrophobic environments such as dodecyl-phosphocholine (DPC) micelle and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocoline (POPC) membrane. Our results indicate that Pin2 folds in DPC with a 79% of alpha-helical content, which is in agreement with the experimental results, while in POPC it has 62.5% of alpha-helical content. On the other hand, Pin2GVG presents a higher percentage of alpha-helical structure in POPC and a smaller content in DPC when compared with Pin2. These results can help to better choose the starting structures in future molecular dynamics simulations of AMPs, because these peptides can adopt slightly different conformations depending on the hydrophobic environment.Communicated by Ramaswamy H. Sarma.

Biophysical characterization of the insertion of two potent antimicrobial peptides-Pin2 and its variant Pin2[GVG] in biological model membranes.

The aim of this study was to investigate the factors that govern the activity and selectivity of two potent antimicrobial peptides (AMPs) using lipid membrane models of bacterial, erythrocyte and fungal cells. These models were used in calcein liposome leakage experiments to explore peptide efficiency. The AMPs (Pin2 and its variant Pin2[GVG]) showed highest affinity towards the bacterial models in the nanomolar range, followed by the erythrocyte and fungal systems. The presence of sterols modulated the variant's selectivity, while the wild type was unaffected. Liposome leakage experiments with Fluorescein Isothiocyanate-dextran (FITC)-dextran conjugates indicated that pore size depended on peptide concentration. Dynamic Light Scattering revealed peptide aggregation in aqueous solution, and that aggregate size was related to activity. The interacting peptides did not alter liposome size, suggesting pore forming activity rather than detergent activity. Atomic Force Microscopy showed differential membrane absorption, being greater in the bacterial model compared to the mammalian model, and pore-like defects were observed. Electrophysiological assays with the Tip-Dip Patch Clamp method provided evidence of changes in the electrical resistance of the membrane. Membrane potential experiments showed that liposomes were also depolarized in the presence of the peptides. Both peptides increased the Laurdan Generalized Polarization of the bacterial model indicating increased viscosity, on the contrary, no effect was observed with the erythrocyte and the fungal models. Peptide membrane insertion and pore formation was corroborated with Langmuir Pressure-Area isotherms and Brewster Angle Microscopy. Finally, molecular dynamics simulations were used to get an insight into the molecular mechanism of action.

Vasoinhibin comprises a three-helix bundle and its antiangiogenic domain is located within the first 79 residues.

Vasoinhibin belongs to a family of angiogenesis inhibitors generated when the fourth α-helix (H4) of the hormone prolactin (PRL) is removed by specific proteolytic cleavage. The antiangiogenic properties are absent in uncleaved PRL, indicating that conformational changes create a new bioactive domain. However, the solution structure of vasoinhibin and the location of its bioactive domain are unknown. Molecular dynamic simulation (MD) showed that the loss of H4 exposes the hydrophobic nucleus of PRL and leads to the compression of the molecule into a three-helix bundle that buries the hydrophobic nucleus again. Compression occurs by the movement of loop 1 (L1) and its interaction with α-helix 1 (H1) generating a new L1 conformation with electrostatic and hydrophobic surfaces distinct from those of PRL, that may correspond to a bioactive domain. Consistent with this model, a recombinant protein containing the first 79 amino acids comprising H1 and L1 of human PRL inhibited the proliferation and migration of endothelial cells and upregulated the vasoinhibin target genes, IL1A and ICAM1. This bioactivity was comparable to that of a conventional vasoinhibin having the 123 residues encompassing H1, L1, Η2, L2, and Η3 of human PRL. These findings extend the vasoinhibin family to smaller proteins and provide important structural information, which will aid in antiangiogenic drug development.

Molecular dynamics simulation of the membrane binding and disruption mechanisms by antimicrobial scorpion venom-derived peptides.

Pandinin 2 (Pin2) is an alpha-helical polycationic peptide, identified and characterized from venom of the African scorpion Pandinus imperator with high antimicrobial activity against Gram-positive bacteria and less active against Gram-negative bacteria, however it has demonstrated strong hemolytic activity against sheep red blood cells. In the chemically synthesized Pin2GVG analog, the GVG motif grants it low hemolytic activity while keeping its antimicrobial activity. In this work, we performed 12 µs all-atom molecular dynamics simulation of the antimicrobial peptides (AMPs) Pin2 and Pin2GVG to explore their adsorption mechanism and the role of their constituent amino acid residues when interacting with pure POPC and pure POPG membrane bilayers. Starting from an α-helical conformation, both AMPs are attracted at different rates to the POPC and POPG bilayer surfaces due to the electrostatic interaction between the positively charged amino acid residues and the charged moieties of the membranes. Since POPG is an anionic membrane, the PAMs adhesion is stronger to the POPG membrane than to the POPC membrane and they are stabilized more rapidly. This study reveals that, before the insertion begins, Pin2 and Pin2GVG remained partially folded in the POPC surface during the first 300 and 600 ns, respectively, while they are mostly unfolded in the POPG surface during most of the simulation time. The unfolded structures provide for a large number of intermolecular hydrogen bonds and stronger electrostatic interactions with the POPG surface. The results show that the aromatic residues at the N-terminus of Pin2 initiate the insertion process in both POPC and POPG bilayers. As for Pin2GVG in POPC the C-terminus residues seem to initiate the insertion process while in POPG this process seems to be slowed down due to a strong electrostatic attraction. The membrane conformational effects upon PAMs binding are measured in terms of the area per lipid and the contact surface area. Several replicas of the systems lead to the same observations.

Ligand recognition properties of the vasopressin V2 receptor studied under QSAR and molecular modeling strategies.

The design of new drugs that target vasopressin 2 receptor (V2R) is of vital importance to develop new therapeutic alternatives to treat diseases such as heart failure, polycystic kidney disease. To get structural insights related to V2R-ligand recognition, we have used a combined approach of docking, molecular dynamics simulations (MD) and quantitative structure-activity relationship (QSAR) to elucidate the detailed interaction of the V2R with 119 of its antagonists. The three-dimensional model of V2R was built by threading methods refining its structure through MD simulations upon which the 119 ligands were subjected to docking studies. The theoretical results show that binding recognition of these ligands on V2R is diverse, but the main pharmacophore (electronic and π-π interactions) is maintained; thus, this information was validated under QSAR results. QSAR studies were performed using MLR analysis followed by ANN analysis to increase the model quality. The final equation was developed by choosing the optimal combination of descriptors after removing the outliers. The applicability domains of the constructed QSAR models were defined using the leverage and standardization approaches. The results suggest that the proposed QSAR models can reliably predict the reproductive toxicity potential of diverse chemicals, and they can be useful tools for screening new chemicals for safety assessment.

Phosphorylation regulates the secondary structure and function of dentin phosphoprotein peptides.

Dentin phosphoprotein (DPP) is the most acidic protein in vertebrates and structurally is classified as an intrinsically disordered protein. Functionally, DPP is related to dentin and bone formation, however the specifics of such association remain unknown. Here, we used atomistic molecular dynamics simulations to screen selected binding domains of DPP onto hydroxyapatite (HA), which is one of its important interacting partners. From these results, we selected a functionally relevant peptide, Ace-SSDSSDSSDSSDSSD-NH2 (named P5) and its phosphorylated form (named P5P), for experimental characterization. SAXS experiments indicated that in solution P5 was disordered, possibly in an extended conformation while P5P displayed more compact globular conformations. Circular dichroism and FTIR confirmed that, either in the presence or absence of Ca2+/HA, P5 adopts a random coil structure, whereas its phosphorylated counterpart, P5P, has a more compact arrangement associated with conformations that display ß-sheet and α-helix motifs when bound to HA. In solution, P5 inhibited HA crystal growth, whereas at similar concentrations, P5P stimulated it. These findings suggest that phosphorylation controls the transient formation of secondary and tertiary structure of DPP peptides, and, most likely of DPP itself, which in turn controls HA growth in solution and possibly HA growth in mineralized tissues.

Simple approach for ranking structure determining residues.

Mutating residues has been a common task in order to study structural properties of the protein of interest. Here, we propose and validate a simple method that allows the identification of structural determinants; i.e., residues essential for preservation of the stability of global structure, regardless of the protein topology. This method evaluates all of the residues in a 3D structure of a given globular protein by ranking them according to their connectivity and movement restrictions without topology constraints. Our results matched up with sequence-based predictors that look up for intrinsically disordered segments, suggesting that protein disorder can also be described with the proposed methodology.

QSAR/QSPR as an application of artificial neural networks.

Quantitative Structure-Activity Relationships (QSARs) and Quantitative Structure-Property Relationships (QSPRs) are mathematical models used to describe and predict a particular activity/property of compounds. On the other hand, the Artificial Neural Network (ANN) is a tool that emulates the human brain to solve very complex problems. The exponential need for new compounds in the drug industry requires alternatives for experimental methods to decrease development time and costs. This is where chemical computational methods have a great relevance, especially QSAR/QSPR-ANN. This chapter shows the importance of QSAR/QSPR-ANN and provides examples of its use.

The role of phosphorylation in dentin phosphoprotein peptide absorption to hydroxyapatite surfaces: a molecular dynamics study.

Dentin phosphoprotein (DPP) is a protein expressed mainly in dentin and to a lesser extent in bone. DPP has a disordered structure, rich in glutamic acid, aspartic acid and phosphorylated serine/threonine residues. It has a high capacity for binding to calcium ions and to hydroxyapatite (HA) crystal surfaces. We used molecular dynamics (MD) simulations as a method for virtually screening interactions between DPP motifs and HA. The goal was to determine which motifs are absorbed to HA surfaces. For these simulations, we considered five peptides from the human DPP sequence. All-atom MD simulations were performed using GROMACS, the peptides were oriented parallel to the {100} HA crystal surface, the distance between the HA and the peptide was 3 nm. The system was simulated for 20 ns. Preliminary results show that for the unphosphorylated peptides, the acidic amino acids present an electrostatic attraction where their side chains are oriented towards HA. This attraction, however, is slow to facilitate bulk transport to the crystal surface. On the other hand, the phosphorylated (PP) peptides are rapidly absorbed on the surface of the HA with their centers of mass closer to the HA surface. More importantly, the root mean square fluctuation (RMSF) indicates that the average structures of the phosphorylated peptides are very inflexible and elongate, while that of the unphosphorylated peptides are flexible. Radius of gyration (Rg) analysis showed the compactness of un-phosphorylated peptides is lower than phosphorylated peptides. Phosphorylation of the DPP peptides is necessary for binding to HA surfaces.

Exploring the ligand recognition properties of the human vasopressin V1a receptor using QSAR and molecular modeling studies.

Vaptans are compounds that act as non-peptide vasopressin receptor antagonists. These compounds have diverse chemical structures. In this study, we used a combined approach of protein folding, molecular dynamics simulations, docking, and quantitative structure-activity relationship (QSAR) to elucidate the detailed interaction of the vasopressin receptor V1a (V1aR) with some of its blockers (134). QSAR studies were performed using MLR analysis and were gathered into one group to perform an artificial neural network (ANN) analysis. For each molecule, 1481 molecular descriptors were calculated. Additionally, 15 quantum chemical descriptors were calculated. The final equation was developed by choosing the optimal combination of descriptors after removing the outliers. Molecular modeling enabled us to obtain a reliable tridimensional model of V1aR. The docking results indicated that the great majority of ligands reach the binding site under π-π, π-cation, and hydrophobic interactions. The QSAR studies demonstrated that the heteroatoms N and O are important for ligand recognition, which could explain the structural diversity of ligands that reach V1aR.

A beta-complex statistical four body contact potential combined with a hydrogen bond statistical potential recognizes the correct native structure from protein decoy sets.

We present a new four-body knowledge-based potential for recognizing the native state of proteins from their misfolded states. This potential was extracted from a large set of protein structures determined by X-ray crystallography using BetaMol, a software based on the recent theory of the beta-complex (ß-complex) and quasi-triangulation of the Voronoi diagram of spheres. This geometric construct reflects the size difference among atoms in their full Euclidean metric; property not accounted for in a typical 3D Delaunay triangulation. The ability of this potential to identify the native conformation over a large set of decoys was evaluated. Experiments show that this potential outperforms a potential constructed with a classical Delaunay triangulation in decoy discrimination tests. The addition of a statistical hydrogen bond potential to our four-body potential allows a significant improvement in the decoy discrimination, in such a way that we are able to predict successfully the native structure in 90% of cases.

Low-temperature molecular dynamics simulations of horse heart cytochrome c and comparison with inelastic neutron scattering data.

Molecular dynamics (MD) simulation combined with inelastic neutron scattering can provide information about the thermal dynamics of proteins, especially the low-frequency vibrational modes responsible for large movement of some parts of protein molecules. We performed several 30-ns MD simulations of cytochrome c (Cyt c) in a water box for temperatures ranging from 110 to 300 K and compared the results with those from experimental inelastic neutron scattering. The low-frequency vibrational modes were obtained via dynamic structure factors, S(Q, ω), obtained both from inelastic neutron scattering experiments and calculated from MD simulations for Cyt c in the same range of temperatures. The well known thermal transition in structural movements of Cyt c is clearly seen in MD simulations; it is, however, confined to unstructured fragments of loops Ω1 and Ω2; movement of structured loop Ω3 and both helical ends of the protein is resistant to thermal disturbance. Calculated and experimental S(Q, ω) plots are in qualitative agreement for low temperatures whereas above 200 K a boson peak vanishes from the calculated plots. This may be a result of loss of crystal structure by the protein-water system compared with the protein crystal.