Protein surface and core dynamics show concerted hydration-dependent activation

Inelastic neutron scattering and spectroscopy methods to characterize dynamics in colloidal and soft matter systems

Colloidal systems and soft materials are well suited to neutron scattering, and the community has readily adopted elastic scattering techniques to investigate their structure. Due to their unique properties, neutrons may also be used to characterize the dynamics of soft materials over a wide range of length and time scales in situ. Both static structures and an understanding of how molecules move about their equilibrium positions is essential if we are to deliver on the promise of rationally designing soft materials. In this review we introduce the basics of neutron spectroscopy and explore the ways in which inelastic neutron scattering can be used to study colloidal and soft materials. Illustrative examples are chosen that highlight the phenomena suitable for investigation using this suite of techniques.

Structural characterization of an intrinsically disordered protein complex using integrated small-angle neutron scattering and computing

Characterizing structural ensembles of intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) of proteins is essential for studying structure-function relationships. Due to the different neutron scattering lengths of hydrogen and deuterium, selective labeling and contrast matching in small-angle neutron scattering (SANS) becomes an effective tool to study dynamic structures of disordered systems. However, experimental timescales typically capture measurements averaged over multiple conformations, leaving complex SANS data for disentanglement. We hereby demonstrate an integrated method to elucidate the structural ensemble of a complex formed by two IDRs. We use data from both full contrast and contrast matching with residue-specific deuterium labeling SANS experiments, microsecond all-atom molecular dynamics (MD) simulations with four molecular mechanics force fields, and an autoencoder-based deep learning (DL) algorithm. From our combined approach, we show that selective deuteration provides additional information that helps characterize structural ensembles. We find that among the four force fields, a99SB-disp and CHARMM36m show the strongest agreement with SANS and NMR experiments. In addition, our DL algorithm not only complements conventional structural analysis methods but also successfully differentiates NMR and MD structures which are indistinguishable on the free energy surface. Lastly, we present an ensemble that describes experimental SANS and NMR data better than MD ensembles generated by one single force field and reveal three clusters of distinct conformations. Our results demonstrate a new integrated approach for characterizing structural ensembles of IDPs.

Deuteration Aiming for Neutron Scattering

The distinguished feature of neutron as a scattering probe is an isotope effect, especially the large difference in neutron scattering length between hydrogen and deuterium. The difference renders the different visibility between hydrogenated and deuterated proteins. Therefore, the combination of deuterated protein and neutron scattering enables the selective visualization of a target domain in the complex or a target protein in the multi-component system. Despite of this fascinating character, there exist several problems for the general use of this method: difficulty and high cost for protein deuteration, and control and determination of deuteration ratio of the sample. To resolve them, the protocol of protein deuteration techniques is presented in this report. It is strongly expected that this protocol will offer more opportunity for conducting the neutron scattering studies with deuterated proteins.

Structural insights into protein folding, stability and activity using in vivo perdeuteration of hen egg-white lysozyme

This structural and biophysical study exploited a method of perdeuterating hen egg-white lysozyme based on the expression of insoluble protein in Escherichia coli followed by in-column chemical refolding. This allowed detailed comparisons with perdeuterated lysozyme produced in the yeast Pichia pastoris, as well as with unlabelled lysozyme. Both perdeuterated variants exhibit reduced thermal stability and enzymatic activity in comparison with hydrogenated lysozyme. The thermal stability of refolded perdeuterated lysozyme is 4.9°C lower than that of the perdeuterated variant expressed and secreted in yeast and 6.8°C lower than that of the hydrogenated Gallus gallus protein. However, both perdeuterated variants exhibit a comparable activity. Atomic resolution X-ray crystallographic analyses show that the differences in thermal stability and enzymatic function are correlated with refolding and deuteration effects. The hydrogen/deuterium isotope effect causes a decrease in the stability and activity of the perdeuterated analogues; this is believed to occur through a combination of changes to hydrophobicity and protein dynamics. The lower level of thermal stability of the refolded perdeuterated lysozyme is caused by the unrestrained Asn103 peptide-plane flip during the unfolded state, leading to a significant increase in disorder of the Lys97-Gly104 region following subsequent refolding. An ancillary outcome of this study has been the development of an efficient and financially viable protocol that allows stable and active perdeuterated lysozyme to be more easily available for scientific applications.

Dynamics of proteins with different molecular structures under solution condition

Incoherent quasielastic neutron scattering (iQENS) is a fascinating technique for investigating the internal dynamics of protein. However, low flux of neutron beam, low signal to noise ratio of QENS spectrometers and unavailability of well-established analyzing method have been obstacles for studying internal dynamics under physiological condition (in solution). The recent progress of neutron source and spectrometer provide the fine iQENS profile with high statistics and as well the progress of computational technique enable us to quantitatively reveal the internal dynamic from the obtained iQENS profile. The internal dynamics of two proteins, globular domain protein (GDP) and intrinsically disordered protein (IDP) in solution, were measured with the state-of-the art QENS spectrometer and then revealed with the newly developed analyzing method. It was clarified that the average relaxation rate of IDP was larger than that of GDP and the fraction of mobile H atoms of IDP was also much higher than that of GDP. Combined with the structural analysis and the calculation of solvent accessible surface area of amino acid residue, it was concluded that the internal dynamics were related to the highly solvent exposed amino acid residues depending upon protein's structure.

Role of hydration water in the onset of protein structural dynamics

Proteins are the molecular workhorses in a living organism. Their 3D structures are animated by a multitude of equilibrium fluctuations and specific out-of-equilibrium motions that are required for proteins to be biologically active. When studied as a function of temperature, functionally relevant dynamics are observed at and above the so-called protein dynamical transition (~240 K) in hydrated, but not in dry proteins. In this review we present and discuss the main experimental and computational results that provided evidence for the dynamical transition, with a focus on the role of hydration water dynamics in sustaining functional protein dynamics. The coupling and mutual influence of hydration water dynamics and protein dynamics are discussed and the hypotheses illustrated that have been put forward to explain the physical origin of their onsets.

Dynamics from elastic neutron-scattering via direct measurement of the running time-integral of the van Hove distribution function

We present a new neutron-scattering approach to access the van Hove distribution function directly in the time domain, I(t), which reflects the system dynamics. Currently, I(t) is essentially determined from neutron energy-exchange. Our method consists of the straightforward measurement of the running time-integral of I(t), by computing the portion of scattered neutrons corresponding to species at rest within a time t, (conceptually elastic scattering). Previous attempts failed to recognise this connection. Starting from a theoretical standpoint, a practical realisation is assessed via numerical methods and an instrument simulation.

Aromatic Ring Dynamics, Thermal Activation, and Transient Conformations of a 468 kDa Enzyme by Specific 1H-13C Labeling and Fast Magic-Angle Spinning NMR

Aromatic residues are located at structurally important sites of many proteins. Probing their interactions and dynamics can provide important functional insight but is challenging in large proteins. Here, we introduce approaches to characterize the dynamics of phenylalanine residues using ¹H-detected fast magic-angle spinning (MAS) NMR combined with a tailored isotope-labeling scheme. Our approach yields isolated two-spin systems that are ideally suited for artifact-free dynamics measurements, and allows probing motions effectively without molecular weight limitations. The application to the TET2 enzyme assembly of ∼0.5 MDa size, the currently largest protein assigned by MAS NMR, provides insights into motions occurring on a wide range of time scales (picoseconds to milliseconds). We quantitatively probe ring-flip motions and show the temperature dependence by MAS NMR measurements down to 100 K. Interestingly, favorable line widths are observed down to 100 K, with potential implications for DNP NMR. Furthermore, we report the first ¹³C R_1ρ MAS NMR relaxation-dispersion measurements and detect structural excursions occurring on a microsecond time scale in the entry pore to the catalytic chamber and at a trimer interface that was proposed as the exit pore. We show that the labeling scheme with deuteration at ca. 50 kHz MAS provides superior resolution compared to 100 kHz MAS experiments with protonated, uniformly ¹³C-labeled samples.

Neutron scattering in the biological sciences: progress and prospects

The scattering of neutrons can be used to provide information on the structure and dynamics of biological systems on multiple length and time scales. Pursuant to a National Science Foundation-funded workshop in February 2018, recent developments in this field are reviewed here, as well as future prospects that can be expected given recent advances in sources, instrumentation and computational power and methods. Crystallography, solution scattering, dynamics, membranes, labeling and imaging are examined. For the extraction of maximum information, the incorporation of judicious specific deuterium labeling, the integration of several types of experiment, and interpretation using high-performance computer simulation models are often found to be particularly powerful.

Biological small-angle neutron scattering: recent results and development

Small-angle neutron scattering (SANS) has increasingly been used by the structural biology community in recent years to obtain low-resolution information on solubilized biomacromolecular complexes in solution. In combination with deuterium labelling and solvent-contrast variation (H2O/D2O exchange), SANS provides unique information on individual components in large heterogeneous complexes that is perfectly complementary to the structural restraints provided by crystallography, nuclear magnetic resonance and electron microscopy. Typical systems studied include multi-protein or protein–DNA/RNA complexes and solubilized membrane proteins. The internal features of these systems are less accessible to the more broadly used small-angle X-ray scattering (SAXS) technique owing to a limited range of intra-complex and solvent electron-density variation. Here, the progress and developments of biological applications of SANS in the past decade are reviewed. The review covers scientific results from selected biological systems, including protein–protein complexes, protein–RNA/DNA complexes and membrane proteins. Moreover, an overview of recent developments in instruments, sample environment, deuterium labelling and software is presented. Finally, the perspectives for biological SANS in the context of integrated structural biology approaches are discussed.

Dynamic Neutron Scattering by Biological Systems

Dynamic neutron scattering directly probes motions in biological systems on femtosecond to microsecond timescales. When combined with molecular dynamics simulation and normal mode analysis, detailed descriptions of the forms and frequencies of motions can be derived. We examine vibrations in proteins, the temperature dependence of protein motions, and concepts describing the rich variety of motions detectable using neutrons in biological systems at physiological temperatures. New techniques for deriving information on collective motions using coherent scattering are also reviewed.

Influence of cosolvents, self-crowding, temperature and pressure on the sub-nanosecond dynamics and folding stability of lysozyme

We studied the effects of temperature and hydrostatic pressure on the dynamical properties and folding stability of highly concentrated lysozyme solutions in the absence and presence of the osmolytes trimethylamine-N-oxide (TMAO) and urea. Elastic incoherent neutron scattering (EINS) was applied to determine the mean-squared displacement (MSD) of the protein's hydrogen atoms to yield insights into the effects of these cosolvents on the averaged sub-nanosecond dynamics in the pressure range from ambient up to 4000 bar. To evaluate the additional effect of self-crowding, two protein concentrations (80 and 160 mg mL^-1) were used. We observed a distinct effect of TMAO on the internal hydrogen dynamics, namely a reduced mobility. Urea, on the other hand, revealed no marked effect and consequently, no counteracting effect in an urea-TMAO mixture was observed. Different from the less concentrated protein solution, no significant effect of pressure on the MSD was observed for 160 mg mL^-1 lysozyme. The EINS experiments were complemented by Fourier-transform infrared (FTIR) spectroscopy measurements, which led to additional insights into the folding stability of lysozyme under the various environmental conditions. We observed a stabilization of the protein in the presence of the compatible osmolyte TMAO and a destabilization in the presence of urea against temperature and pressure for both protein concentrations. Additionally, we noticed a slight destabilizing effect upon self-crowding at very high protein concentration (160 mg mL^-1), which is attributable to transient destabilizing intermolecular interactions. Furthermore, a pressure-temperature diagram could be obtained for lysozyme at these high protein concentrations that mimics densely packed intracellular conditions.

Ergodicity breaking of iron displacement in heme proteins

We present a model of the dynamical transition of atomic displacements in proteins. Increased mean-square displacement at higher temperatures is caused by the softening of the force constant for atomic/molecular displacements by electrostatic and van der Waals forces from the protein-water thermal bath. Displacement softening passes through a nonergodic dynamical transition when the relaxation time of the force-force correlation function enters, with increasing temperature, the instrumental observation window. Two crossover temperatures are identified. The lower crossover, presently connected to the glass transition, is related to the dynamical unfreezing of rotations of water molecules within nanodomains polarized by charged surface residues of the protein. The higher crossover temperature, usually assigned to the dynamical transition, marks the onset of water translations. All crossovers are ergodicity breaking transitions depending on the corresponding observation windows. Allowing stretched exponential relaxation of the protein-water thermal bath significantly improves the theory-experiment agreement when applied to solid protein samples studied by Mössbauer spectroscopy.

Solvent-Driven Dynamical Crossover in the Phenylalanine Side-Chain from the Hydrophobic Core of Amyloid Fibrils Detected by 2H NMR Relaxation

Aromatic residues are important markers of dynamical changes in proteins' hydrophobic cores. In this work we investigated the dynamics of the F19 side-chain in the core of amyloid fibrils across a wide temperature range of 300 to 140 K. We utilized solid-state 2H NMR relaxation to demonstrate the presence of a solvent-driven dynamical crossover between different motional regimes, often also referred to as the dynamical transition. In particular, the dynamics are dominated by small-angle fluctuations at low temperatures and by π-flips of the aromatic ring at high temperatures. The crossover temperature is more than 43 degrees lower for the hydrated state of the fibrils compared to the dry state, indicating that interactions with water facilitate π-flips. Further, crossover temperatures are shown to be very sensitive to polymorphic states of the fibrils, such as the 2-fold and 3-fold symmetric morphologies of the wild-type protein as well as D23N mutant protofibrils. We speculate that these differences can be attributed, at least partially, to enhanced interactions with water in the 3-fold polymorph, which has been shown to have a water-accessible cavity. Combined with previous studies of methyl group dynamics, the results highlight the presence of multiple dynamics modes in the core of the fibrils, which was originally believed to be quite rigid.

Static solid-state 2H NMR methods in studies of protein side-chain dynamics

In this review, we discuss the experimental static deuteron NMR techniques and computational approaches most useful for the investigation of side-chain dynamics in protein systems. Focus is placed on the interpretation of line shape and relaxation data within the framework of motional modeling. We consider both jump and diffusion models and apply them to uncover glassy behaviors, conformational exchange and dynamical transitions in proteins. Applications are chosen from globular and membrane proteins, amyloid fibrils, peptide adsorbed on surfaces and proteins specific to connective tissues.

Solid state deuterium NMR study of LKα14 peptide aggregation in biosilica

In nature, organisms including diatoms, radiolaria, and marine sponges use proteins, long chain polyamines, and other organic molecules to regulate the assembly of complex silica-based structures. Here, the authors investigate structural features of small peptides, designed to mimic the silicifying activities of larger proteins found in natural systems. LKα14 (Ac-LKKLLKLLKKLLKL-C), an amphiphilic lysine/leucine repeat peptide with an α-helical secondary structure at polar/apolar interfaces, coprecipitates with silica to form nanospheres. Previous 13C magic angle spinning studies suggest that the tetrameric peptide bundles that LKα14 is known to form in solution may persist in the silica-complexed form, and may also function as catalysts and templates for silica formation. To further investigate LKα14 aggregation in silica, deuterium solid-state nuclear magnetic resonance (2H ssNMR) was used to establish how leucine side-chain dynamics differ in solid LKα14 peptides isolated from aqueous solution, from phosphate-buffered solution, and in the silica-precipitated states. Modeling the 2H ssNMR line shapes probed the mechanisms of peptide preaggregation and silica coprecipitation. The resulting NMR data indicates that the peptide bundles in silica preserve the hydrophobic interior that they display in the hydrated solid state. However, NMR data also indicate free motion of the leucine residues in silica, a condition that may result from structural deformation of the aggregates arising from interactions between the surface lysine side chains and the surrounding silica matrix.

Fast Motions of Key Methyl Groups in Amyloid-β Fibrils

Amyloid-β (Aβ) peptide is the major component of plaques found in Alzheimer's disease patients. Using solid-state 2H NMR relaxation performed on selectively deuterated methyl groups, we probed the dynamics in the threefold symmetric and twofold symmetric polymorphs of native Aβ as well as the protofibrils of the D23N mutant. Specifically, we investigated the methyl groups of two leucine residues that belong to the hydrophobic core (L17 and L34) as well as M35 residues belonging to the hydrophobic interface between the cross-β subunits, which has been previously found to be water-accessible. Relaxation measurements performed over 310-140 K and two magnetic field strengths provide insights into conformational variability within and between polymorphs. Core packing variations within a single polymorph are similar to what is observed for globular proteins for the core residues, whereas M35 exhibits a larger degree of variability. M35 site is also shown to undergo a solvent-dependent dynamical transition in which slower amplitude motions of methyl axes are activated at high temperature. The motions, modeled as a diffusion of methyl axis, have activation energy by a factor of 2.7 larger in the twofold compared with the threefold polymorph, whereas D23N protofibrils display a value similar to the threefold polymorph. This suggests enhanced flexibility of the hydrophobic interface in the threefold polymorph. This difference is only observed in the hydrated state and is absent in the dry fibrils, highlighting the role of solvent at the cavity. In contrast, the dynamic behavior of the core is hydration-independent.

Atomistic details of protein dynamics and the role of hydration water

BACKGROUND

The importance of protein dynamics for their biological activity is now well recognized. Different experimental and computational techniques have been employed to study protein dynamics, hierarchy of different processes and the coupling between protein and hydration water dynamics. Yet, understanding the atomistic details of protein dynamics and the role of hydration water remains rather limited.

SCOOP OF REVIEW

Based on overview of neutron scattering, molecular dynamic simulations, NMR and dielectric spectroscopy results we present a general picture of protein dynamics covering time scales from faster than ps to microseconds and the influence of hydration water on different relaxation processes.

MAJOR CONCLUSIONS

Internal protein dynamics spread over a wide time range from faster than picosecond to longer than microseconds. We suggest that the structural relaxation in hydrated proteins appears on the microsecond time scale, while faster processes present mostly motion of side groups and some domains. Hydration water plays a crucial role in protein dynamics on all time scales. It controls the coupled protein-hydration water relaxation on 10-100ps time scale. This process defines the friction for slower protein dynamics. Analysis suggests that changes in amount of hydration water affect not only general friction, but also influence significantly the protein's energy landscape.

GENERAL SIGNIFICANCE

The proposed atomistic picture of protein dynamics provides deeper understanding of various relaxation processes and their hierarchy, similarity and differences between various biological macromolecules, including proteins, DNA and RNA. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo".

Dynamics of Hydrophobic Core Phenylalanine Residues Probed by Solid-State Deuteron NMR

We conducted a detailed investigation of the dynamics of two phenylalanine side chains in the hydrophobic core of the villin headpiece subdomain protein (HP36) in the hydrated powder state over the 298-80 K temperature range. Our main tools were static deuteron NMR measurements of longitudinal relaxation and line shapes supplemented with computational modeling. The temperature dependence of the relaxation times reveals the presence of two main mechanisms that can be attributed to the ring-flips, dominating at high temperatures, and small-angle fluctuations, dominating at low temperatures. The relaxation is nonexponential at all temperatures with the extent of nonexponentiality increasing from higher to lower temperatures. This behavior suggests a distribution of conformers with unique values of activation energies. The central values of the activation energies for the ring-flipping motions are among the smallest reported for aromatic residues in peptides and proteins and point to a very mobile hydrophobic core. The analysis of the widths of the distributions, in combination with the earlier results on the dynamics of flanking methyl groups (Vugmeyster et al. J. Phys. Chem. B 2013, 117, 6129-6137), suggests that the hydrophobic core undergoes slow concerted fluctuations. There is a pronounced effect of dehydration on the ring-flipping motions, which shifts the distribution toward more rigid conformers. The crossover temperature between the regions of dominance of the small-angle fluctuations and ring-flips shifts from 195 K in the hydrated protein to 278 K in the dry one. This result points to the role of solvent in softening the core and highlights aromatic residues as markers of the protein dynamical transitions.

The unusual internal motion of the villin headpiece subdomain

The thermostable 36-residue subdomain of the villin headpiece (HP36) is the smallest known cooperatively folding protein. Although the folding and internal dynamics of HP36 and close variants have been extensively studied, there has not been a comprehensive investigation of side-chain motion in this protein. Here, the fast motion of methyl-bearing amino acid side chains is explored over a range of temperatures using site-resolved solution nuclear magnetic resonance deuterium relaxation. The squared generalized order parameters of methyl groups extensively spatially segregate according to motional classes. This has not been observed before in any protein studied using this methodology. The class segregation is preserved from 275 to 305 K. Motions detected in Helix 3 suggest a fast timescale of conformational heterogeneity that has not been previously observed but is consistent with a range of folding and dynamics studies. Finally, a comparison between the order parameters in solution with previous results based on solid-state nuclear magnetic resonance deuterium line shape analysis of HP36 in partially hydrated powders shows a clear disagreement for half of the sites. This result has significant implications for the interpretation of data derived from a variety of approaches that rely on partially hydrated protein samples.

Molecular Dynamics Simulations of a Powder Model of the Intrinsically Disordered Protein Tau

The tau protein, whose aggregates are involved in Alzheimer's disease, is an intrinsically disordered protein (IDP) that regulates microtubule activity in neurons. An IDP lacks a single, well-defined structure and, rather, constantly exchanges among multiple conformations. In order to study IDP dynamics, the combination of experimental techniques, such as neutron scattering, and computational techniques, such as molecular dynamics (MD) simulations, is a powerful approach. Amorphous hydrated powder samples have been very useful for studying protein internal dynamics experimentally, e.g., using neutron scattering. Thus, there is demand for realistic in silico models of hydrated protein powders. Here we present an MD simulation analysis of a powder hydrated at 0.4 g water/g protein of the IDP tau in the temperature range 20-300 K. By comparing with neutron scattering data, we identify the protein-water interface as the predominant feature determining IDP dynamics. The so-called protein dynamical transition is shown to be attenuated, but not suppressed, in the parts of the protein that are not exposed to the solvent. In addition, we find similarities in the mean-squared displacements of the core of a globular protein and "dry" clusters formed by the IDP in hydrated powders. Thus, the ps to ns dynamics of proteins in hydrated powders originate mainly from those residues in contact with solvent. We propose that by measuring the dynamics of protein assemblies, such as aggregates, one might assess qualitatively their state of hydration.

Protein dynamics: from rattling in a cage to structural relaxation

We present an overview of protein dynamics based mostly on results of neutron scattering, dielectric relaxation spectroscopy and molecular dynamics simulations. We identify several major classes of protein motions on the time scale from faster than picoseconds to several microseconds, and discuss the coupling of these processes to solvent dynamics. Our analysis suggests that the microsecond backbone relaxation process might be the main structural relaxation of the protein that defines its glass transition temperature, while faster processes present some localized secondary relaxations. Based on the overview, we formulate a general picture of protein dynamics and discuss the challenges in this field.

Intrinsically disordered proteins--relation to general model expressing the active role of the water environment

This work discusses the role of unstructured polypeptide chain fragments in shaping the protein's hydrophobic core. Based on the "fuzzy oil drop" model, which assumes an idealized distribution of hydrophobicity density described by the 3D Gaussian, we can determine which fragments make up the core and pinpoint residues whose location conflicts with theoretical predictions. We show that the structural influence of the water environment determines the positions of disordered fragments, leading to the formation of a hydrophobic core overlaid by a hydrophilic mantle. This phenomenon is further described by studying selected proteins which are known to be unstable and contain intrinsically disordered fragments. Their properties are established quantitatively, explaining the causative relation between the protein's structure and function and facilitating further comparative analyses of various structural models.

Effect of subdomain interactions on methyl group dynamics in the hydrophobic core of villin headpiece protein

Thermostable villin headpiece protein (HP67) consists of the N-terminal subdomain (residues 10-41) and the autonomously folding C-terminal subdomain (residues 42-76) which pack against each other to form a structure with a unified hydrophobic core. The X-ray structures of the isolated C-terminal subdomain (HP36) and its counterpart in HP67 are very similar for the hydrophobic core residues. However, fine rearrangements of the free energy landscape are expected to occur because of the interactions between the two subdomains. We detect and characterize these changes by comparing the µs-ms time scale dynamics of the methyl-bearing side chains in isolated HP36 and in HP67. Specifically, we probe three hydrophobic side chains at the interface of the two subdomains (L42, V50, and L75) as well as at two residues far from the interface (L61 and L69). Solid-state deuteron NMR techniques are combined with computational modeling for the detailed characterization of motional modes in terms of their kinetic and thermodynamic parameters. The effect of interdomain interactions on side chain dynamics is seen for all residues but L75. Thus, changes in dynamics because of subdomain interactions are not confined to the site of perturbation. One of the main results is a two- to threefold increase in the value of the activation energies for the rotameric mode of motions in HP67 compared with HP36. Detailed analysis of configurational entropies and heat capacities complement the kinetic view of the degree of the disorder in the folded state.

Hydrophobic core in domains of immunoglobulin-like fold

This work analyzes proteins which contain an immunoglobulin fold, focusing on their hydrophobic core structure. The "fuzzy oil drop" model was used to measure the regularity of hydrophobicity distribution in globular domains belonging to proteins which exhibit the above-mentioned fold. Light-chain IgG domains are found to frequently contain regular hydrophobic cores, unlike the corresponding heavy-chain domains. Enzymes and DNA binding proteins present in the data-set are found to exhibit poor accordance with the hydrophobic core model.

Zaccai neutron resilience and site-specific hydration dynamics in a globular protein

A discussion is presented of contributions of the Zaccai group to the understanding of flexibility in biological macromolecules using dynamic neutron scattering. The concept of resilience as introduced by Zaccai is discussed and investigated using molecular dynamics simulation on camphor-bound cytochrome P450. The resilience of hydrophilic residues is found to be more strongly affected by hydration than that of hydrophobic counterparts. The hydration-induced softening of protein propagates from the surface into the dry core. Moreover, buried hydrophilic residues behave more like those exposed on the protein surface, and are different from their hydrophobic counterparts.