Energy condensation and dipole alignment in protein dynamics

The possibility that distant biomolecules in a cell interact via electromagnetic (e.m.) radiation was proposed many years ago to explain the high rate of encounters of partners in some enzymatic reactions. The results of two recent experiments designed to test the propensity of protein bovine serum albumin (BSA) to interact via e.m. radiation with other proteins were interpreted in a theoretical framework based on three main assumptions: (i) in order to experience this kind of interaction the protein must be in an out-of-equilibrium state; (ii) in this state there is a condensation of energy in low-frequency vibrational modes; and (iii) the hydration layers of water around the protein sustain the energy condensation. In the present paper we present the results of molecular dynamics simulations of BSA in four states: at equilibrium and out-of-equilibrium in water, and at room and high temperature in vacuum. By comparing physical properties of the system in the four states, our simulations provide a qualitative and quantitative assessment of the three assumptions on which the theoretical framework is based. Our results confirm the assumptions of the theoretical model showing energy condensation at low frequency and electretlike alignment between the protein's and the water's dipoles; they also allow a quantitative estimate of the contribution of the out-of-equilibrium state and of the water to the observed behavior of the protein. In particular, it has been found that in the out-of-equilibrium state the amplitude of the oscillation of the protein's dipole moment greatly increases, thereby enhancing a possible absorption or emission of e.m. radiation. The analysis of BSA's dynamics outlined in the present paper provides a procedure for checking the propensity of a biomolecule to interact via e.m. radiation with its biochemical partners.

Phase-Dependent Adsorption of Myelin Basic Protein to Phosphatidylcholine Lipid Bilayers

The dense packing of opposite cytoplasmic surfaces of the lipid-enriched myelin membrane, responsible for the proper saltatory conduction of nerve impulses through axons, is ensured by the adhesive properties of myelin basic protein (MBP). Although preferentially interacting with negatively charged phosphatidylserine (PS) lipids, as an intrinsically disordered protein, it can easily adapt its shape to its immediate environment and thus adsorb to domains made of zwitterionic phosphatidylcholine (PC) lipids. As the molecular-level interaction pattern between MBP and PC lipid membranes suffers from scarce characterization, an experimental and computational study of multilamellar liposomes (MLVs) composed of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) in the presence of bovine MBP is presented here. Calorimetric and temperature-dependent UV-Vis measurements identified DPPC pretransition temperature (Tp) and calorimetric enthalpy (ΔHcal) as the physicochemical parameters most responsive to the presence of MBP. Besides suggesting an increase in β-sheet fractions of structured MBP segments as DPPC lipids undergo from the gel (20 °C) to the fluid (50 °C) phase, FTIR spectra unraveled the significant contribution of lysine (Lys) residues in the adsorption pattern, especially when DPPC is in the fluid (50 °C) phase. In addition to highlighting the importance of Lys residues in the MBP adsorption on DPPC lipid bilayer, employing salt bridges (SBs) and hydrogen bonds (HBs), MD data suggest the crucial importance of the orientation of MBP with respect to the surface of the DPPC lipid bilayer.

Tau Fibrillation Induced by Heparin or a Lysophospholipid Show Different Initial Oligomer Formation

The protein tau is involved in several neurogenerative diseases such as Alzheimer's Disease, where tau content and fibrillation have been linked to disease progression. Tau colocalizes with phospholipids and glycosaminoglycans in vivo. We investigated if and how tau fibrillation can be induced by two lysophospholipids, namely the zwitterionic 1-myristoyl-2-hydroxy-sn-glycero-3-phosphocholine (LPC) and the anionic 1-myristoyl-2-hydroxy-sn-glycero-3-phospho-(1'-rac-glycerol) (LPG) as well as the glycosaminoglycan heparin. We used a range of biophysical techniques including small-angle X-ray scattering, Thioflavin T fluorescence, and SDS-PAGE, collecting data at various time points to obtain structural information on each phase of the fibrillation. We find that LPC does not induce fibrillation or interact with tau. Low concentrations of LPG induce fibrillation by formation of small hydrophobic clusters with monomeric tau. At higher LPG concentrations, a core-shell complex is formed where tau wraps around LPG micelles with regions extending away from the micelles. For heparin, loosely associated oligomers are formed rapidly with around ten tau molecules. Fibrils formed with either LPG or heparin show similar final cross-section dimensions. Furthermore, SDS-resistant oligomers are observed for both LPG and heparin. Our study demonstrates that tau fibrillation can be induced by two different biologically relevant cofactors leading to structurally different initial states but similar cross-sectional dimensions for the fibrils. Structural information about initial states prior to fibril formation is important both to gain a better understanding of the onset of fibrillation in vivo, and for the development of targeted drugs that can reduce or abolish tau fibrillation.

Synchrotron X-ray study of intrinsically disordered and polyampholytic Tau 4RS and 4RL under controlled ionic strength

Aggregated and hyperphosphorylated Tau is one of the pathological hallmarks of Alzheimer's disease. Tau is a polyampholytic and intrinsically disordered protein (IDP). In this paper, we present for the first time experimental results on the ionic strength dependence of the radius of gyration (Rg) of human Tau 4RS and 4RL isoforms. Synchrotron X-ray scattering revealed that 4RS Rg is regulated from 65.4 to 58.5 Å and 4RL Rg is regulated from 70.9 to 57.9 Å by varying ionic strength from 0.01 to 0.592 M. The Rg of 4RL Tau is larger than 4RS at lower ionic strength. This result provides an insight into the ion-responsive nature of intrinsically disordered and polyampholytic Tau, and can be implicated to the further study of Tau-Tau and Tau-tubulin intermolecular structure in ionic environments.

Combining Experiments and Simulations to Examine the Temperature-Dependent Behavior of a Disordered Protein

Intrinsically disordered proteins are a class of proteins that lack stable folded conformations and instead adopt a range of conformations that determine their biochemical functions. The temperature-dependent behavior of such disordered proteins is complex and can vary depending on the specific protein and environment. Here, we have used molecular dynamics simulations and previously published experimental data to investigate the temperature-dependent behavior of histatin 5, a 24-residue-long polypeptide. We examined the hypothesis that histatin 5 undergoes a loss of polyproline II (PPII) structure with increasing temperature, leading to more compact conformations. We found that the conformational ensembles generated by the simulations generally agree with small-angle X-ray scattering data for histatin 5, but show some discrepancies with the hydrodynamic radius as probed by pulsed-field gradient NMR spectroscopy, and with the secondary structure information derived from circular dichroism. We attempted to reconcile these differences by reweighting the conformational ensembles against the scattering and NMR data. By doing so, we were in part able to capture the temperature-dependent behavior of histatin 5 and to link the observed decrease in hydrodynamic radius with increasing temperature to a loss of PPII structure. We were, however, unable to achieve agreement with both the scattering and NMR data within experimental errors. We discuss different possible reasons for this including inaccuracies in the force field, differences in conditions of the NMR and scattering experiments, and issues related to the calculation of the hydrodynamic radius from conformational ensembles. Our study highlights the importance of integrating multiple types of experimental data when modeling conformational ensembles of disordered proteins and how environmental factors such as the temperature influence them.

Disordered proteins mitigate the temperature dependence of site-specific binding free energies

Biophysical characterization of protein-protein interactions involving disordered proteins is challenging. A common simplification is to measure the thermodynamics and kinetics of disordered site binding using peptides containing only the minimum residues necessary. We should not assume, however, that these few residues tell the whole story. Son of sevenless, a multidomain signaling protein from Drosophila melanogaster, is critical to the mitogen-activated protein kinase pathway, passing an external signal to Ras, which leads to cellular responses. The disordered 55 kDa C-terminal domain of Son of sevenless is an autoinhibitor that blocks guanidine exchange factor activity. Activation requires another protein, Downstream of receptor kinase (Drk), which contains two Src homology 3 domains. Here, we utilized NMR spectroscopy and isothermal titration calorimetry to quantify the thermodynamics and kinetics of the N-terminal Src homology 3 domain binding to the strongest sites incorporated into the flanking disordered sequences. Comparing these results to those for isolated peptides provides information about how the larger domain affects binding. The affinities of sites on the disordered domain are like those of the peptides at low temperatures but less sensitive to temperature. Our results, combined with observations showing that intrinsically disordered proteins become more compact with increasing temperature, suggest a mechanism for this effect.

The biophysics of disordered proteins from the point of view of single-molecule fluorescence spectroscopy

Intrinsically disordered proteins (IDPs) and regions (IDRs) have emerged as key players across many biological functions and diseases. Differently from structured proteins, disordered proteins lack stable structure and are particularly sensitive to changes in the surrounding environment. Investigation of disordered ensembles requires new approaches and concepts for quantifying conformations, dynamics, and interactions. Here, we provide a short description of the fundamental biophysical properties of disordered proteins as understood through the lens of single-molecule fluorescence observations. Single-molecule Förster resonance energy transfer (FRET) and fluorescence correlation spectroscopy (FCS) provides an extensive and versatile toolbox for quantifying the characteristics of conformational distributions and the dynamics of disordered proteins across many different solution conditions, both in vitro and in living cells.

Negative Thermal Expansion and Disorder-to-Order Collapse of an Intrinsically Disordered Protein under Marginally Denaturing Conditions

Recent work with intrinsically disordered proteins (IDPs) has projected a myriad of their survival instincts based mainly on the total charge content, the abundance of polar residues, and the paucity of hydrophobic amino acids. This work uses a plant IDP AtPP16-1 (Arabidopsis thaliana phloem protein class 16-1), whose solution NMR structure was determined by us recently, to show legitimate negative thermal expansion (NTE) of the native state. The thermal expansion continues to be negative even as the tertiary structure is perturbed by ultralow levels of urea up to 0.4 M. The NTE of these subdenatured states is called apparent NTE because they are prone to undergo conformational changes with temperature. Hydrodynamic shrinkage of the NTE IDP is also observed by dynamic light scattering (DLS) and NMR-measured global rotational correlation time (τ_c). The protein with denatured tertiary structure but otherwise preserved native-state secondary structure collapses to a dynamically rigid state. The data are mainly based on thermal coefficients of chemical shift and nuclear relaxation measured by heteronuclear NMR. The hydrodynamic shrinkage and collapse under marginally varying solvent compositions that may arise from unstable tertiary structure and dynamic disorder of chain segments across the backbone could be a generic property of IDPs.

Enhancing the Biological Activities of Food Protein-Derived Peptides Using Non-Thermal Technologies: A Review

Bioactive peptides (BPs) derived from animal and plant proteins are important food functional ingredients with many promising health-promoting properties. In the food industry, enzymatic hydrolysis is the most common technique employed for the liberation of BPs from proteins in which conventional heat treatment is used as pre-treatment to enhance hydrolytic action. In recent years, application of non-thermal food processing technologies such as ultrasound (US), high-pressure processing (HPP), and pulsed electric field (PEF) as pre-treatment methods has gained considerable research attention owing to the enhancement in yield and bioactivity of resulting peptides. This review provides an overview of bioactivities of peptides obtained from animal and plant proteins and an insight into the impact of US, HPP, and PEF as non-thermal treatment prior to enzymolysis on the generation of food-derived BPs and resulting bioactivities. US, HPP, and PEF were reported to improve antioxidant, angiotensin-converting enzyme (ACE)-inhibitory, antimicrobial, and antidiabetic properties of the food-derived BPs. The primary modes of action are due to conformational changes of food proteins caused by US, HPP, and PEF, improving the susceptibility of proteins to protease cleavage and subsequent proteolysis. However, the use of other non-thermal techniques such as cold plasma, radiofrequency electric field, dense phase carbon dioxide, and oscillating magnetic fields has not been examined in the generation of BPs from food proteins.

Modeling Solution Behavior of Poly(N-isopropylacrylamide): A Comparison between Water Models

Water is known to play a fundamental role in determining the structure and functionality of macromolecules. The same crucial contribution is also found in the in silico description of polymer aqueous solutions. In this work, we exploit the widely investigated synthetic polymer poly(N-isopropylacrylamide) (PNIPAM) to understand the effect of the adopted water model on its solution behavior and to refine the computational setup. By means of atomistic molecular dynamics simulations, we perform a comparative study of PNIPAM aqueous solution using two advanced water models: TIP4P/2005 and TIP4P/Ice. The conformation and hydration features of an atactic 30-mer at infinite dilution are probed at a range of temperature and pressure suitable to detect the coil-to-globule transition and to map the P–T phase diagram. Although both water models can reproduce the temperature-induced coil-to-globule transition at atmospheric pressure and the polymer hydration enhancement that occurs with increasing pressure, the PNIPAM–TIP4P/Ice solution shows better agreement with experimental findings. This result can be attributed to a stronger interaction of TIP4P/Ice water with both hydrophilic and hydrophobic groups of PNIPAM, as well as to a less favorable contribution of the solvent entropy to the coil-to-globule transition.

Biochemical and biophysical features of disease-associated tau mutants V363A and V363I

The comprehension of pathogenetic mechanisms in tauopathy-associated neurodegenerative diseases can be improved by the knowledge of the biochemical and biophysical features of mutated tau proteins. Here, we used the full-length, wild-type tau, the V363A and V363I mutated species, associated with pathology, and the P301L mutated tau as a benchmark. Using several techniques, including small-angle X-ray scattering, atomic force microscopy, thioflavin T binding, and electrophoretic separation, we compared their course from intrinsically disordered monomers in solution to early-stage recruitment in complexes and then aggregates of increasing size over long periods up to the asymptotic aggregative behavior of full-length tau proteins. We showed that diversity in the kinetics of recruitment and aggregate structure occurs from the beginning and spreads all over their pathway to very large objects. The different extents of conformational changes and types of molecular assemblies among the proteins were also reflected in their in vitro toxicity; this variation could correlate with physiopathology in humans, considering that the P301L mutation is more aggressive than V363A, especially V363I. This study identified the presence of aggregation intermediates and corroborated the oligomeric hypothesis of tauopathies.

Recent Advances in the Label-Free Characterization of Exosomes for Cancer Liquid Biopsy: From Scattering and Spectroscopy to Nanoindentation and Nanodevices

Exosomes (EXOs) are nano-sized vesicles secreted by most cell types. They are abundant in bio-fluids and harbor specific molecular constituents from their parental cells. Due to these characteristics, EXOs have a great potential in cancer diagnostics for liquid biopsy and personalized medicine. Despite this unique potential, EXOs are not yet widely applied in clinical settings, with two main factors hindering their translational process in diagnostics. Firstly, conventional extraction methods are time-consuming, require large sample volumes and expensive equipment, and often do not provide high-purity samples. Secondly, characterization methods have some limitations, because they are often qualitative, need extensive labeling or complex sampling procedures that can induce artifacts. In this context, novel label-free approaches are rapidly emerging, and are holding potential to revolutionize EXO diagnostics. These methods include the use of nanodevices for EXO purification, and vibrational spectroscopies, scattering, and nanoindentation for characterization. In this progress report, we summarize recent key advances in label-free techniques for EXO purification and characterization. We point out that these methods contribute to reducing costs and processing times, provide complementary information compared to the conventional characterization techniques, and enhance flexibility, thus favoring the discovery of novel and unexplored EXO-based biomarkers. In this process, the impact of nanotechnology is systematically highlighted, showing how the effectiveness of these techniques can be enhanced using nanomaterials, such as plasmonic nanoparticles and nanostructured surfaces, which enable the exploitation of advanced physical phenomena occurring at the nanoscale level.

Temperature Dependence of Internal Friction of Peptides

Internal friction is a valuable concept to describe the kinetics of proteins. As is well known, internal friction can be modulated by solvent features (such as viscosity). How can internal friction be affected by environmental temperature? The answer to this question is not evident. In the present work, we approach this problem with simulations on two model peptides. The thermodynamics and relaxation kinetics are characterized through long molecular dynamics simulations, with the viscosity modulated by varying the mass of solvent molecules. Based on the extrapolation to zero viscosity together with scaling of the relaxation time scales, we discover that internal friction is almost invariant at various temperatures. Controlled simulations further support the idea that internal friction is independent of environmental temperature. Comparisons between the two model peptides help us to understand the diverse phenomena in experiments.

Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis

Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aβ, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.

Conformation of Myoglobin-Poly(Ethyl Ethylene Phosphate) Conjugates Probed by SANS: Correlation with Polymer Grafting Density and Interaction

One can take advantage of the influence of a polymer conjugated with a protein to control the thermal stability and the deployment of the protein. Here, the structural properties are reported of the protein-polymer conjugate myoglobin (Mb)-poly(ethyl ethylene phosphate) (PEEP) in the native and unfolded conformations, in order to understand the respective roles of the protein and of the polymer size in the stability of the conjugate. The effect is also investigated of the grafting density of the linear biodegradable polyphosphoesters covalently attached to the protein. It is observed that, while the conjugation process at room temperature does not modify the secondary and tertiary structure of the Mb, the unfolding process, as a function of temperature, depends on the grafting density. Small angle neutron scattering reveals that, at room temperature, conjugation does not alter the size of the native protein and that the thickness of the polymer shell around the protein increases as a function of grafting density and of polymer molecular weight. The denatured form of all conjugates is described by an unfolded chain and a correlation length due to the presence of local stiffness.

Kinetic coherence underlies the dynamics of disordered proteins

The dynamics of two proteins of similar size, the globular lysozyme and the intrinsically disordered Huntingtin interacting protein, has been simulated in three states resembling a globule, a pre-molten globule, and a molten globule. A coherence time τ has been defined, measuring the delay in the display of a stochastic behaviour after a perturbation of the system. This time has been computed for two sets of collective variables: the projection of the phase point onto the positions and momenta subspaces (τ_r and τ_p), and the principal components (PCs) of positions q and momenta π produced by a covariance analysis in these subspaces (τ_q and τ_π). In all states τ_p ≈ 3.5τ_r, and τ_π ≈ 3.5τ_q. The coherence times of individual PCs, τ^(l)_q and τ^(l)_π, have also been computed, and τ^(l)_π > τ^(l)_q in all states. The prevalence of τ_p over τ_r, or of τ_π over τ_q, drives the dynamics of the protein over a time range of ≈1–2 ps; moreover, a hidden synchronism appears to raise the momenta subspace's coherence above that of its individual PCs. In the transition of lysozyme to the molten globule the τ^(l)_q decrease but, unexpectedly, the τ^(l)_π increase; after this transition τ_p ≈ 5τ_r and τ_π ≈ 5τ_q. A gain of kinetic coherence accompanies thus the loss of structural coherence caused by the denaturation of the protein in the transition from globule to molten globule. The increase of the τ^(l)_π does not take place in the analogous transition of the Huntingtin protein. These results are compared with those of a similar analysis performed on three pseudo-proteins designed by scrambling the primary sequence of the Huntingtin interacting protein, and on two oligopeptides. The hidden synchronism appears to be a generic property of these polypeptides. The τ^(l)_π spectrum is similar in denaturated and in intrinsically disordered biomolecules; but the gain of kinetic coherence as a result of denaturation seems to be a specific property of the biologically functional lysozyme.

In the phase space of a globular or intrinsically disordered protein, the momenta's dynamics is less chaotic than the coordinates' dynamics. When a protein is denaturated, a gain in kinetic coherence accompanies the loss of structural coherence.

Dynamic structural determinants underlie the neurotoxicity of the N-terminal tau 26-44 peptide in Alzheimer's disease and other human tauopathies

The intrinsically disordered tau protein plays a pivotal role in the pathogenesis of Alzheimer's disease (AD) and other human tauopathies. Abnormal post-translational modifications of tau, such as truncation, are causally involved in the onset/development of these neurodegenerative diseases. In this context, the AD-relevant N-terminal fragment mapping between 26 and 44 amino acids of protein (tau26-44) is interesting, being endowed with potent neurotoxic effects in vitro and in vivo. However, the understanding of the mechanism(s) of tau26-44 toxicity is a challenging task because, similarly to the full-length tau, it does not have a unique 3D structure but exists as dynamic ensemble of conformations. Here we use Atomic Force Spectroscopy, Small Angle X-ray Scattering and Molecular Dynamics simulation to gather structural and functional information on the tau26-44. We highlight the presence, the type and the location of its temporary secondary structures and we unveil the occurrence of relevant transient tertiary conformations that could contribute to tau26-44 toxicity. Data are compared with those obtained on the biologically-inactive, reverse-sequence (tau44-26 peptide) which has the same mass, charge, aminoacidic composition as well as the same overall unfolded character of tau26-44.

Application of native mass spectrometry in studying intrinsically disordered proteins: A special focus on neurodegenerative diseases

Intrinsically disordered proteins (IDPs) are integral part of the proteome, regulating vital biological processes. Such proteins gained further visibility due to their key role in neurodegenerative diseases and cancer. IDPs however, escape structural characterization by traditional biophysical tools owing to their extreme flexibility and heterogeneity. In this review, we discuss the advantages of native mass spectrometry (MS) in analysing the atypical conformational dynamics of IDPs and recent advances made in the field. Especially, MS studies unravelling the conformational facets of IDPs involved in neurodegenerative diseases are highlighted. The limitations and the future promises of native MS while studying IDPs have been discussed.

Human neuroblastoma SH-SY5Y cells treated with okadaic acid express phosphorylated high molecular weight tau-immunoreactive protein species

BACKGROUND

Early stages of Alzheimer's disease (AD) are characterized by high phosphorylation of microtubule-associated protein tau, which may result from the downregulation of protein phosphatases.

NEW METHOD

In order to model phosphatase downregulation and analyze its effect on tau aggregation in vitro, we treated neuroblastoma SH-SY5Y cells with okadaic acid (OA), a protein phosphatase inhibitor, and examined high molecular weight phospho-tau species.

RESULTS AND COMPARISON WITH EXISTING METHODS

OA treatment led to the appearance of heat-stable protein species with apparent molecular weight around 100 kDa, which were immunoreactive to anti-tau antibodies against phosphorylated Ser202 and Ser396. As these high molecular weight tau-immunoreactive proteins (HMW-TIPs) corresponded to the predicted size of two tau monomers, we considered the possibility that they represent phosphorylation-induced tau oligomers. We attempted to dissociate HMW-TIPs by urea and guanidine, as well as by alkaline phosphatase treatment, but HMW-TIPs were stable under all conditions tested. These characteristics resemble properties of certain sodium dodecyl sulfate (SDS)-resistant tau oligomers from AD brains. The absence of HMW-TIPs detection by anti-total tau antibodies Tau46, CP27 and Tau13 may be a consequence of epitope masking and protein truncation. Alternatively, HMW-TIPs may represent previously unreported phosphoproteins cross-reacting with tau.

CONCLUSIONS

Taken together, our data provide a novel characterization of an OA-based cell culture model in which OA induces the appearance of HMW-TIPs. These findings have implications for further studies of tau under the conditions of protein phosphatase downregulation, aiming to explain mechanisms involved in early events leading to AD.

Toward a unified scoring function for native state discrimination and drug-binding pocket recognition

Protein folding and receptor-ligand recognition are fundamental processes for any living organism. Although folding and ligand recognition are based on the same chemistry, the existing empirical scoring functions target just one problem: predicting the correct fold or the correct binding pose. We here introduce a statistical potential which considers moieties as fundamental units. The scoring function is able to deal with both folding and ligand pocket recognition problems with a performance comparable to the scoring functions specifically tailored for one of the two tasks. We foresee that the capability of the new scoring function to tackle both problems in a unified framework will be a key to deal with the induced fit phenomena, in which a target protein changes significantly its conformation upon binding. Moreover, the new scoring function might be useful in docking protocols towards intrinsically disordered proteins, whose flexibility cannot be handled with the available docking software.

Molecular Dynamics Simulations Combined with Nuclear Magnetic Resonance and/or Small-Angle X-ray Scattering Data for Characterizing Intrinsically Disordered Protein Conformational Ensembles

The concept of intrinsically disordered proteins (IDPs) has emerged relatively slowly, but over the past 20 years, it has become an intense research area in structural biology. Indeed, because of their considerable flexibility and structural heterogeneity, the determination of IDP conformational ensemble is particularly challenging and often requires a combination of experimental measurements and computational approaches. With the improved accuracy of all-atom force fields and the increasing computing performances, molecular dynamics (MD) simulations have become more and more reliable to generate realistic conformational ensembles. And the combination of MD simulations with experimental approaches, such as nuclear magnetic resonance (NMR) and/or small-angle X-ray scattering (SAXS) allows one to converge toward a more accurate and exhaustive description of IDP structures. In this Review, we discuss the state of the art of MD simulations of IDP conformational ensembles, with a special focus on studies that back-calculated and directly compared theoretical and experimental NMR or SAXS observables, such as chemical shifts (CS), ³J-couplings (³Jc), residual dipolar couplings (RDC), or SAXS intensities. We organize the review in three parts. In the first section, we discuss the studies which used NMR and/or SAXS data to test and validate the development of force fields or enhanced sampling techniques. In the second part, we explore different methods for the refinement of MD-derived structural ensembles, such as NMR or SAXS data-restrained MD simulations or ensemble reweighting to better fit experiments. Finally, we survey some recent studies combining MD simulations with NMR and/or SAXS measurements to investigate the relationship between IDP conformational ensemble and biological activity, as well as their implication in human diseases. From this review, we noticed that quite a few studies compared MD-generated conformational ensembles with both NMR and SAXS measurements to validate IDP structures at both local and global levels. Yet, beside the IDP propensity to form local secondary structures, their dynamic extension or compactness also appears important for their activity. Thus, we believe that a close synergy between MD simulations, NMR, and SAXS experiments would be greatly appropriate to address the challenges of characterizing the disordered structures of proteins and their complexes, relative to their biological functions.

Temperature-induced collapse of a disordered peptide observed by three sampling methods in molecular dynamics simulations

The conformational ensembles of a disordered peptide, polyglutamine Q15, over a wide temperature range were sampled using multiple replicates of conventional molecular dynamics (cMD) simulations as well as two enhanced sampling methods, temperature replica exchange (TREMD) and replica exchange with solute tempering (REST). The radius of gyration, asphericity, secondary structure, and hydrogen bonding patterns were used for the comparison of the sampling methods. Overall, the three sampling methods generated similar conformational ensembles, with progressive collapse at higher temperatures. Although accumulating the longest simulation time (90 μs), cMD at room temperature missed a small subspace that was sampled by both TREMD and REST. This subspace was high in α-helical content and separated from the main conformational space by an energy barrier. REST used less simulation time than TREMD (36 μs versus 42 μs), and this gap is expected to widen significantly for larger disordered proteins. We conclude that REST is the method of choice for conformational sampling of intrinsically disordered proteins.

Comment on “Thermal compaction of the intrinsically disordered protein tau: entropic, structural, and hydrophobic factors” by A. Battisti, G. Ciasca, A. Grottesi and A. Tenenbaum, Phys. Chem. Chem. Phys., 2017, 19, 8435

Chain compaction is favoured on raising the temperature because the entropy gain of water molecules due to the decrease in solvent-excluded volume increases with temperature.