Prion and water: tight and dynamical hydration sites have a key role in structural stability

Experiment and molecular dynamics simulations reveal proanthocyanidin B2 and B3 can inhibit prion aggregation by different mechanisms

Prion diseases are a group of fatal neurodegenerative diseases caused by the misfolding and aggregation of prion protein (PrP), and the inhibition of PrP aggregation is one of the most effective therapeutic strategies. Proanthocyanidin B2 (PB2) and B3 (PB3), the effective natural antioxidants have been evaluated for the inhibition of amyloid-related protein aggregation. Since PrP has similar aggregation mechanism with other amyloid-related proteins, will PB2 and PB3 affect the aggregation of PrP? In this paper, experimental and molecular dynamics (MD) simulation methods were combined to investigate the influence of PB2 and PB3 on PrP aggregation. Thioflavin T assays showed PB2 and PB3 could inhibit PrP aggregation in a concentrate-dependent manner in vitro. To understand the underlying mechanism, we performed 400 ns all-atom MD simulations. The results suggested PB2 could stabilize the α2 C-terminus and the hydrophobic core of protein by stabilizing two important salt bridges R156-E196 and R156-D202, and consequently made global structure of protein more stable. Surprisingly, PB3 could not stabilize PrP, which may inhibit PrP aggregation through a different mechanism. Since dimerization is the first step of aggregation, will PB3 inhibit PrP aggregation by inhibiting the dimerization? To verify our assumption, we then explored the effect of PB3 on protein dimerization by performing 800 ns MD simulations. The results suggested PB3 could reduce the residue contacts and hydrogen bonds between two monomers, preventing dimerization process of PrP. The possible inhibition mechanism of PB2 and PB3 on PrP aggregation could provide useful information for drug development against prion diseases.Communicated by Ramaswamy H. Sarma.

The Advances and Limitations of the Determination and Applications of Water Structure in Molecular Engineering

Water is a key actor of various processes of nature and, therefore, molecular engineering has to take the structural and energetic consequences of hydration into account. While the present review focuses on the target-ligand interactions in drug design, with a focus on biomolecules, these methods and applications can be easily adapted to other fields of the molecular engineering of molecular complexes, including solid hydrates. The review starts with the problems and solutions of the determination of water structures. The experimental approaches and theoretical calculations are summarized, including conceptual classifications. The implementations and applications of water models are featured for the calculation of the binding thermodynamics and computational ligand docking. It is concluded that theoretical approaches not only reproduce or complete experimental water structures, but also provide key information on the contribution of individual water molecules and are indispensable tools in molecular engineering.

Evolutionary Role of Water-Accessible Cavities in Src Homology 2 (SH2) Domains

Protein excited states are fundamental in the understanding of biological function, despite the fact they are hardly observed using traditional biophysical methodologies. Pressure perturbation coupled with nuclear magnetic resonance (NMR) spectroscopy is a powerful physicochemical tool to glance at these low-populated high-energy states on a residue-by-residue basis and underpin mechanistic insights into protein functionalities. Here we performed pressure titrations using NMR spectroscopy and relaxation dispersion experiments to identify the low-lying energetic states of the c-Abl SH2 domain. By showing that the SH2 excited state contains a hydrated hydrophobic cavity, fast-exchange motions, and highly conserved residues facing the water-accessible hole, we discuss the implications of water-protein interactions in SH2 modules achieving high-affinity binding and promiscuous phospho-Tyr peptide recognition.

Advances in the understanding of protein misfolding and aggregation through molecular dynamics simulation

Aberrant protein folding known as protein misfolding is counted as one of the striking factors of neurodegenerative diseases. The extensive range of pathologies caused by protein misfolding, aggregation and subsequent accumulation are mainly classified into either gain of function diseases or loss of function diseases. In order to seek for novel strategies for treatment and diagnosis of neurodegenerative diseases, insights into the mechanism of misfolding and aggregation is essential. A comprehensive knowledge on the factors influencing misfolding and aggregation is required as well. An extensive experimental study on protein aggregation is somewhat challenging due to the insoluble and noncrystalline nature of amyloid fibrils. Thus there has been a growing use of computational approaches including Monte Carlo simulation, docking simulation, molecular dynamics simulation in the study of protein misfolding and aggregation. The review presents a discussion on molecular dynamics simulation alone as to how it has emerged as a promising tool in the understanding of protein misfolding and aggregation in general, detailing upon three different aspects considering four misfold prone proteins in particular. It is noticeable that all four proteins considered in this review i.e prion, superoxide dismutase1, huntingtin and amyloid β are linked to chronic neurodegenerative diseases with debilitating effects. Initially the review elaborates on the factors influencing the misfolding and aggregation. Next, it addresses our current understanding of the amyloid structures and the associated aggregation mechanisms, finally, summarizing the contribution of this computational tool in the search for therapeutic strategies against the respective protein-deposition diseases.

p53 amyloid aggregation in cancer: function, mechanism, and therapy

Similar to neurodegenerative diseases, the concept that tumors are prion like diseases has been proposed in recent years. p53, the most well-known tumor suppressor, has been extensively studied for its expression, mutation, and function in various tumors. Currently, an interesting phenomenon of p53 prion-like aggregation has been found in several tumors, and studies have found that its pathological aggregation may lead to functional alterations and ultimately affect tumor progression. It has been demonstrated that the mechanism of p53 aggregation involves its mutation, domains, isoform, etc. In addition to p53 itself, some other factors, including Zn²⁺ concentration, pH, temperature and chaperone abnormalities, can also contribute to p53 aggregation. Although there are some studies about the mechanism and role of p53 aggregation and amyloidosis in tumors, there still exist some controversies. In this paper, we review the mechanism of p53 amyloid fibril structure and discuss the characteristics and effects of p53 amyloid aggregation, as well as the pathogenic mechanism leading to the occurrence of aggregation in tumors. Finally, we summarize the various inhibitors targeting p53 aggregation and prion-like behavior. In conclusion, a comprehensive understanding of p53 aggregation can expand our understanding of the causes leading its loss of physiological function and that targeting p53 aggregation might be a promising therapeutic strategy for tumor therapy.

Destabilization of polar interactions in the prion protein triggers misfolding and oligomerization

The prion protein (PrP) misfolds and oligomerizes at pH 4 in the presence of physiological salt concentrations. Low pH and salt cause structural perturbations in the monomeric prion protein that lead to misfolding and oligomerization. However, the changes in stability within different regions of the PrP prior to oligomerization are poorly understood. In this study, we have characterized the local stability in PrP at high resolution using amide temperature coefficients (T_C ) measured by nuclear magnetic resonance (NMR) spectroscopy. The local stability of PrP was investigated under native as well as oligomerizing conditions. We have also studied the rapidly oligomerizing PrP variant (Q216R) and the protective PrP variant (A6). We report that at low pH, salt destabilizes PrP at several polar residues, and the hydrogen bonds in helices α2 and α3 are weakened. In addition, salt changes the curvature of the α3 helix, which likely disrupts α2-α3 contacts and leads to oligomerization. These results are corroborated by the T_C values of rapidly oligomerizing Q216R-PrP. The poly-alanine substitution in A6-PrP stabilizes α2, which prevents oligomerization. Altogether, these results highlight the importance of native polar interactions in determining the stability of PrP and reveal the structural disruptions in PrP that lead to misfolding and oligomerization.

Solvent Relaxation NMR: A Tool for Real-Time Monitoring Water Dynamics in Protein Aggregation Landscape

Solvent dynamics strongly induce the fibrillation of an amyloidogenic system. Probing the solvation mechanism is crucial as it enables us to predict different proteins' functionalities, such as the aggregation propensity, structural flexibility, and toxicity. This work shows that a straightforward NMR method in conjunction with phenomenological models gives a global and qualitative picture of water dynamics at different concentrations and temperatures. Here, we study amyloid system Aβ40 and its fragment AV20 (A21-V40) and G37L (mutation at Gly37 → Leu of AV20), having different aggregation and toxic properties. The independent validation of this method is elucidated using all-atom classical MD simulation. These two state-of-the-art techniques are pivotal in linking the effect of solvent environment in the near hydration-shell to their aggregation nature. The time-dependent modulation in solvent dynamics probed with the NMR solvent relaxation method can be further adopted to gain insight into amyloidogenesis and link with their toxicity profiles.

Residue interaction networks of K-Ras protein with water molecules identifies the potential role of switch II and P-loop

The mutant K-Ras with aberrant signaling is the primary cause of several cancers. The proposed study investigated the influence of water molecules in K-Ras crystal structure, where they have a significant function by understanding their residue interaction networks (RINs). We analyzed the RINs of K-Ras with and without water molecules and determined their interaction properties. RINs were developed with the help of StructureViz2 and RINspector; further, the changes in K-Ras backbone flexibility were predicted with the DynaMine. We found that the residues K42, I142, and L159 are the hotspots from water, including the K-Ras-GTP complex with the highest residue centrality analysis (RCA) Z-score. The DynaMine prediction calculated the NMR S² value for the frequently mutated positions G12, G13, and Q61 showing a minor shift in flexibility, which make up the P-Loop and switch II of the K-Ras protein. This flexibility shift can account for changes in conformational activity and the protein's GTPase activity, making it difficult to recognize by the effectors and exchange factors. Taken together, our study helps in understanding the functional importance of the water molecules in K-Ras protein and the impact of mutation that modulate the conformational state of the protein.

Molecular Simulations Guidelines for Biological Nanomaterials: From Peptides to Membranes

In studying biological processes and focusing on the molecular mechanisms at the basis of these, molecular dynamics (MD) simulations have demonstrated to be a very useful tool for the past 50 years. This suite of computational methods calculates the time-dependent evolution of a molecular system using physics-based first principles. In this chapter, we give a brief introduction to the theory and practical use of molecular dynamics simulations, highlighting the different models and algorithms that have been developed to tackle specific problems, with a special focus on classical force fields. Some examples of how simulations have been used in the past will help the reader in discerning their power, limitations, and significance.

Insight From Animals Resistant to Prion Diseases: Deciphering the Genotype - Morphotype - Phenotype Code for the Prion Protein

Prion diseases are a group of neurodegenerative diseases endemic in humans and several ruminants caused by the misfolding of native prion protein (PrP) into pathological conformations. Experimental work and the mad-cow epidemic of the 1980s exposed a wide spectrum of animal susceptibility to prion diseases, including a few highly resistant animals: horses, rabbits, pigs, and dogs/canids. The variable susceptibility to disease offers a unique opportunity to uncover the mechanisms governing PrP misfolding, neurotoxicity, and transmission. Previous work indicates that PrP-intrinsic differences (sequence) are the main contributors to disease susceptibility. Several residues have been cited as critical for encoding PrP conformational stability in prion-resistant animals, including D/E159 in dog, S167 in horse, and S174 in rabbit and pig PrP (all according to human numbering). These amino acids alter PrP properties in a variety of assays, but we still do not clearly understand the structural correlates of PrP toxicity. Additional insight can be extracted from comparative structural studies, followed by molecular dynamics simulations of selected mutations, and testing in manipulable animal models. Our working hypothesis is that protective amino acids generate more compact and stable structures in a C-terminal subdomain of the PrP globular domain. We will explore this idea in this review and identify subdomains within the globular domain that may hold the key to unravel how conformational stability and disease susceptibility are encoded in PrP.

Structure and Dynamics of Nanoconfined Water Between Surfactant Monolayers

The properties of nanoconfined water arise in direct response to the properties of the interfaces that confine it. A great deal of research has focused on understanding how and why the physical properties of confined water differ greatly from the bulk. In this work, we have used all-atom molecular dynamics (MD) simulations to provide a detailed description of the structural and dynamical properties of nanoconfined water between two monolayers consisting of an archetypal ionic surfactant, cetrimonium bromide (CTAB, [CH₃(CH₂)₁₅N(CH₃)₃]⁺Br^-). Small differences in the area per surfactant of the monolayers impart a clear effect on the intrinsic density, mobility, and ordering of the interfacial water layer confined by the monolayers. We find that as the area per surfactant within a monolayer decreases, the mobility of the interfacial water molecules decreases in response. As the monolayer packing density decreases, we find that each individual CTAB molecule has a greater effect on the ordering of water molecules in its first hydration shell. In a denser monolayer, we observe that the effect of individual CTAB molecules on the ordering of water molecules is hindered by increased competition between headgroups. Therefore, when two monolayers with different areas per surfactant are used to confine a nanoscale water layer, we observe the emergence of noncentrosymmetry.

Combining molecular dynamics simulations and experimental analyses in protein misfolding

The fold of a protein determines its function and its misfolding can result in loss-of-function defects. In addition, for certain proteins their misfolding can lead to gain-of-function toxicities resulting in protein misfolding diseases such as Alzheimer's, Parkinson's, or the prion diseases. In all of these diseases one or more proteins misfold and aggregate into disease-specific assemblies, often in the form of fibrillar amyloid deposits. Most, if not all, protein misfolding diseases share a fundamental molecular mechanism that governs the misfolding and subsequent aggregation. A wide variety of experimental methods have contributed to our knowledge about misfolded protein aggregates, some of which are briefly described in this review. The misfolding mechanism itself is difficult to investigate, as the necessary timescale and resolution of the misfolding events often lie outside of the observable parameter space. Molecular dynamics simulations fill this gap by virtue of their intrinsic, molecular perspective and the step-by-step iterative process that forms the basis of the simulations. This review focuses on molecular dynamics simulations and how they combine with experimental analyses to provide detailed insights into protein misfolding and the ensuing diseases.

Surface Topography Effects of Globular Biomolecules on Hydration Water

Hydration water serves as a microscopic manifestation of structural stability and functions of biomolecules. To develop bio-nanomaterials in applications, it is important to study how the surface topography and heterogeneity of biomolecules result in their diversity of the hydration dynamics and energetics. We here performed molecular dynamics simulations combined with the steered molecular dynamics and umbrella sampling to investigate the dynamics and escape process associated with the free energy change of water molecules close to a globular biomolecule, i.e., hemoglobin (Hb) and G-quadruplex DNA (GDNA). The residence time, power of long-time tail, and dipole relaxation time were found to display drastic changes within the averaged hydration shell of 3.0-5.0 Å. Compared with bulk water, in the inner hydration shell, the water dipole moment displays a slower relaxation process and is more oriented toward GDNA than toward Hb, forming a hedgehog-like structure when it surrounds GDNA. In particular, a spine water structure is observed in the GDNA narrow groove. The water isotope effect not only prolongs the dynamic time scales of libration motion in the inner hydration shell and the dipole relaxation processes in the bulk but also strengthens the DNA spine water structure. The potential of the mean force profile reflects the integrity of the hydration shell structure and enables us to obtain detailed insights into the structures formed by water, such as the caged H-bond network and the edge bridge structures; it also reveals that local hydration shell free energy (LHSFE) depends on H-bond rupture processes and ranges from 0.2 to 4.2 kcal/mol. Our results demonstrate that the surface topography of a biomolecule influences the integrity of the hydration shell structure and LHSFE. Our studies are able to identify various further applications in the areas of microfluid devices and nano-dewetting on bioinspired surfaces.

Aquaphotomics-From Innovative Knowledge to Integrative Platform in Science and Technology

Aquaphotomics is a young scientific discipline based on innovative knowledge of water molecular network, which as an intrinsic part of every aqueous system is being shaped by all of its components and the properties of the environment. With a high capacity for hydrogen bonding, water molecules are extremely sensitive to any changes the system undergoes. In highly aqueous systems-especially biological-water is the most abundant molecule. Minute changes in system elements or surroundings affect multitude of water molecules, causing rearrangements of water molecular network. Using light of various frequencies as a probe, the specifics of water structure can be extracted from the water spectrum, indirectly providing information about all the internal and external elements influencing the system. The water spectral pattern hence becomes an integrative descriptor of the system state. Aquaphotomics and the new knowledge of water originated from the field of near infrared spectroscopy. This technique resulted in significant findings about water structure-function relationships in various systems contributing to a better understanding of basic life phenomena. From this foundation, aquaphotomics started integration with other disciplines into systematized science from which a variety of applications ensued. This review will present the basics of this emerging science and its technological potential.

Applications of in Silico Methods for Design and Development of Drugs Targeting Protein-Protein Interactions

BACKGROUND

Identification of Protein-Protein Interactions (PPIs) is a major challenge in modern molecular biology and biochemistry research, due to the unquestionable role of proteins in cells, biological process and pathological states. Over the past decade, the PPIs have evolved from being considered a highly challenging field of research to being investigated and examined as targets for pharmacological intervention.

OBJECTIVE

Comprehension of protein interactions is crucial to known how proteins come together to build signalling pathways, to carry out their functions, or to cause diseases, when deregulated. Multiplicity and great amount of PPIs structures offer a huge number of new and potential targets for the treatment of different diseases.

METHODS

Computational techniques are becoming predominant in PPIs studies for their effectiveness, flexibility, accuracy and cost. As a matter of fact, there are effective in silico approaches which are able to identify PPIs and PPI site. Such methods for computational target prediction have been developed through molecular descriptors and data-mining procedures.

RESULTS

In this review, we present different types of interactions between protein-protein and the application of in silico methods for design and development of drugs targeting PPIs. We described computational approaches for the identification of possible targets on protein surface and to detect of stimulator/ inhibitor molecules.

CONCLUSION

A deeper study of the most recent bioinformatics methodologies for PPIs studies is vital for a better understanding of protein complexes and for discover new potential PPI modulators in therapeutic intervention.

Heterogeneity and Architecture of Pathological Prion Protein Assemblies: Time to Revisit the Molecular Basis of the Prion Replication Process?

Prions are proteinaceous infectious agents responsible for a range of neurodegenerative diseases in animals and humans. Prion particles are assemblies formed from a misfolded, β-sheet rich, aggregation-prone isoform (PrP^Sc) of the host-encoded cellular prion protein (PrP^C). Prions replicate by recruiting and converting PrP^C into PrP^Sc, by an autocatalytic process. PrP^Sc is a pleiomorphic protein as different conformations can dictate different disease phenotypes in the same host species. This is the basis of the strain phenomenon in prion diseases. Recent experimental evidence suggests further structural heterogeneity in PrP^Sc assemblies within specific prion populations and strains. Still, this diversity is rather seen as a size continuum of assemblies with the same core structure, while analysis of the available experimental data points to the existence of structurally distinct arrangements. The atomic structure of PrP^Sc has not been elucidated so far, making the prion replication process difficult to understand. All currently available models suggest that PrP^Sc assemblies exhibit a PrP^Sc subunit as core constituent, which was recently identified. This review summarizes our current knowledge on prion assembly heterogeneity down to the subunit level and will discuss its importance with regard to the current molecular principles of the prion replication process.

Hyperosmotic stress induces cell-dependent aggregation of α-synuclein

The aggregation of alpha-synuclein (α-syn) is a pathological feature of a number of neurodegenerative conditions, including Parkinson’s disease. Genetic mutations, abnormal protein synthesis, environmental stress, and aging have all been implicated as causative factors in this process. The importance of water in the polymerisation of monomers, however, has largely been overlooked. In the present study, we highlight the role of hyperosmotic stress in inducing human α-syn to aggregate in cells in vitro, through rapid treatment of the cells with three different osmolytes: sugar, salt and alcohol. This effect is cell-dependent and not due to direct protein-osmolyte interaction, and is specific for α-syn when compared to other neurodegeneration-related proteins, such as Tau or Huntingtin. This new property of α-syn not only highlights a unique aspect of its behaviour which may have some relevance for disease states, but may also be useful as a screening test for compounds to inhibit the aggregation of α-syn in vitro.

Probing the early stages of prion protein (PrP) aggregation with atomistic molecular dynamics simulations

Prions are self-replicating infectious proteinaceous agents whose conformations are capable of forming amyloid-like aggregate fibrils. Here we present molecular dynamics simulations aimed at investigating the aggregation process of the β-rich H2H3 domain of the ovine prion protein (H2H3-OvPrPSc), known to be the portion of prion protein carrying oligomerization activity.

Distribution of tightly and loosely bound water in biological macromolecules and age-related diseases

This mini-review article is focused on publications devoted to the changes in water binding energy and content of bound water in biological tissues during aging processes, when bound water lost from the hydration layer becomes free water. Bound water is released during cataractogenesis. In skin, water bound to proteins and other biomacromolecules becomes more mobile with increasing skin age. Extracellular to intracellular water ratio increases with age and was associated with muscle cell atrophy. Bound water concentration decreases with age in hydrated human bone and can be correlated with the strength and toughness of the bone. Higher fraction of free water in malignant tissues compared to normal tissues was observed. Hydration water mobility is enhanced around tau amyloid fibers. Water plays a decisive role in amyloid formation as entropic driving force. In the natural aging processes dehydration and glycation may be considered as subsequent steps.

Free Energy, Enthalpy and Entropy from Implicit Solvent End-Point Simulations

Free energy is the key quantity to describe the thermodynamics of biological systems. In this perspective we consider the calculation of free energy, enthalpy and entropy from end-point molecular dynamics simulations. Since the enthalpy may be calculated as the ensemble average over equilibrated simulation snapshots the difficulties related to free energy calculation are ultimately related to the calculation of the entropy of the system and in particular of the solvent entropy. In the last two decades implicit solvent models have been used to circumvent the problem and to take into account solvent entropy implicitly in the solvation terms. More recently outstanding advancement in both implicit solvent models and in entropy calculations are making the goal of free energy estimation from end-point simulations more feasible than ever before. We review briefly the basic theory and discuss the advancements in light of practical applications.

P. de Oliveira G. Targeting the Prion-like Aggregation of Mutant p53 to Combat Cancer

Prion-like behavior of several amyloidogenic proteins has been demonstrated in recent years. Despite having functional roles in some cases, irregular aggregation can have devastating consequences. The most commonly known amyloid diseases are Alzheimer's, Parkinson's, and Transmissible Spongiform Encephalopathies (TSEs). The pathophysiology of prion-like diseases involves the structural transformation of wild-type (wt) proteins to transmissible forms that can convert healthy proteins, generating aggregates. The mutant form of tumor suppressor protein, p53, has recently been shown to exhibit prion-like properties. Within the context of p53 aggregation and the search for ways to avert it, this review emphasizes discoveries, approaches, and research from our laboratory and others. Although its standard functions are strongly connected to tumor suppression, p53 mutants and aggregates are involved in cancer progression. p53 aggregates are heterogeneous assemblies composed of amorphous aggregates, oligomers, and amyloid-like fibrils. Evidence of these structures in tumor tissues, the in vitro capability for p53 mutants to coaggregate with wt protein, and the detection of cell-to-cell transmission indicate that cancer has the basic characteristics of prion and prion-like diseases. Various approaches aim to restore p53 functions in cancer. Methods include the use of small-molecule and peptide stabilizers of mutant p53, zinc administration, gene therapy, alkylating and DNA intercalators, and blockage of p53-MDM2 interaction. A primary challenge in developing small-molecule inhibitors of p53 aggregation is the large number of p53 mutations. Another issue is the inability to recover p53 function by dissociating mature fibrils. Consequently, efforts have emerged to target the intermediate species of the aggregation reaction. Φ-value analysis has been used to characterize the kinetics of the early phases of p53 aggregation. Our experiments using high hydrostatic pressure (HHP) and chemical denaturants have helped to clarify excited conformers of p53 that are prone to aggregation. Molecular dynamics (MD) and phasor analysis of single Trp fluorescence signals point toward the presence of preamyloidogenic conformations of p53, which are not observed for p63 or p73. Exploring the features of competent preamyloidogenic states of wt and different p53 mutants may provide a framework for designing personalized drugs for the restoration of p53 function. Protection of backbone hydrogen bonds (BHBs) has been shown to be an important factor for the stability of amyloidogenic proteins and was employed to identify and stabilize the structural defect resulting from the p53 Y220C mutation. Using MD simulations, we compared BHB protection factors between p53 family members to determine the donor-acceptor pairs in p53 that exhibit lower protection. The identification of structurally vulnerable sites in p53 should provide new insights into rational designs that can rapidly be screened using our experimental methodology. Through continued and combined efforts, the outlook is positive for the development of strategies for regulating p53 amyloid transformation.

Molecular Dynamics Simulation Study on the Binding and Stabilization Mechanism of Antiprion Compounds to the "Hot Spot" Region of PrPC

Structural transitions in the prion protein from the cellular form, PrP^C, into the pathological isoform, PrP^Sc, are regarded as the main cause of the transmissible spongiform encephalopathies, also known as prion diseases. Hence, discovering and designing effective antiprion drugs that can inhibit PrP^C to PrP^Sc conversion is regarded as a promising way to cure prion disease. Among several strategies to inhibit PrP^C to PrP^Sc conversion, stabilizing the native PrP^C via specific binding is believed to be one of the valuable approaches and many antiprion compounds have been reported based on this strategy. However, the detailed mechanism to stabilize the native PrP^C is still unknown. As such, to unravel the stabilizing mechanism of these compounds to PrP^C is valuable for the further design and discovery of antiprion compounds. In this study, by molecular dynamics simulation method, we investigated the stabilizing mechanism of several antiprion compounds on PrP^C that were previously reported to have specific binding to the "hot spot" region of PrP^C. Our simulation results reveal that the stabilization mechanism of specific binding compounds can be summarized as (I) to stabilize both the flexible C-terminal of α2 and the hydrophobic core, such as BMD42-29 and GN8; (II) to stabilize the hydrophobic core, such as J1 and GJP49; (III) to stabilize the overall structure of PrP^C by high binding affinity, as NPR-056. In addition, as indicated by the H-bond analysis and decomposition analysis of binding free energy, the residues N159 and Q160 play an important role in the specific binding of the studied compounds and all these compounds interact with PrP^C in a similar way with the key interacting residues L130 in the β1 strand, P158, N159, Q160, etc. in the α1-β2 loop, and H187, T190, T191, etc. in the α2 C-terminus although the compounds have large structural difference. As a whole, our obtained results can provide some insights into the specific binding mechanism of main antiprion compounds to the "hot spot" region of PrP^C at the molecular level and also provide guidance for effective antiprion drug design in the future.

Studying the role of cooperative hydration in stabilizing folded protein states

Understanding and modelling protein folding remains a key scientific and engineering challenge. Two key questions in protein folding are (1) why many proteins adopt a folded state and (2) how these proteins transition from the random coil ensemble to a folded state. In this paper we employ molecular dynamics simulations to address the first of these questions. Computational methods are well-placed to address this issue due to their ability to analyze systems at atomic-level resolution. Traditionally, the stability of folded proteins has been ascribed to the balance of two types of intermolecular interactions: hydrogen-bonding interactions and hydrophobic contacts. In this study, we explore a third type of intermolecular interaction: cooperative hydration of protein surface residues. To achieve this, we consider multiple independent simulations of the villin headpiece domain to quantify the contributions of different interactions to the energy of the native and fully extended states. In addition, we consider whether these findings are robust with respect to the protein forcefield, the water model, and the presence of salt. In all cases, we identify many cooperatively hydrated interactions that are transient but energetically favor the native state. Whilst further work on additional protein structures, forcefields, and water models is necessary, these results suggest a role for cooperative hydration in protein folding that should be explored further. Rational design of cooperative hydration on the protein surface could be a viable strategy for increasing protein stability.

Aggregation tendencies in the p53 family are modulated by backbone hydrogen bonds

The p53 family of proteins is comprised of p53, p63 and p73. Because the p53 DNA binding domain (DBD) is naturally unstable and possesses an amyloidogenic sequence, it is prone to form amyloid fibrils, causing loss of functions. To develop p53 therapies, it is necessary to understand the molecular basis of p53 instability and aggregation. Light scattering, thioflavin T (ThT) and high hydrostatic pressure (HHP) assays showed that p53 DBD aggregates faster and to a greater extent than p63 and p73 DBDs, and was more susceptible to denaturation. The aggregation tendencies of p53, p63, and p73 DBDs were strongly correlated with their thermal stabilities. Molecular Dynamics (MD) simulations indicated specific regions of structural heterogeneity unique to p53, which may be promoted by elevated incidence of exposed backbone hydrogen bonds (BHBs). The results indicate regions of structural vulnerability in the p53 DBD, suggesting new targetable sites for modulating p53 stability and aggregation, a potential approach to cancer therapy.

Pathogenic Mutations within the Disordered Palindromic Region of the Prion Protein Induce Structure Therein and Accelerate the Formation of Misfolded Oligomers

Little is understood about how the intrinsically disordered N-terminal region (NTR) of the prion protein modulates its misfolding and aggregation, which lead to prion disease. In this study, two pathogenic mutations, G113V and A116V, in the palindromic region of the NTR are shown to have no effect on the structure, stability, or dynamics of native mouse prion protein (moPrP) but nevertheless accelerate misfolding and oligomerization. For wild-type moPrP, misfolding and oligomerization appear to occur concurrently, while for both mutant variants, oligomerization is shown to precede misfolding. Kinetic hydrogen-deuterium exchange-mass spectrometry experiments show that sequence segment 89-132 from the NTR becomes structured, albeit weakly, during the oligomerization of both mutant variants. Importantly, this structure formation occurs prior to structural conversion in the C-terminal domain and appears to be the reason that the formation of misfolded oligomers is accelerated by the pathogenic mutations.

Confined Water in Amyloid-Competent Oligomers of the Prion Protein

Conformational switching of the prion protein into the abnormal form involves the formation of (obligatory) molten-oligomers that mature into ordered amyloid fibrils. The role of water in directing the course of amyloid formation remains poorly understood. Here, we show that the mobility of the water molecules within the on-pathway oligomers is highly retarded. The water relaxation time within the oligomers was estimated to be ≈1 ns which is about three orders of magnitude slower than the bulk water and resembles the characteristics of (trapped) nano-confined water. We propose that the coalescence of these obligatory oligomers containing trapped water is entropically favored because of the release of ordered water molecules in the bulk milieu and results in the sequestration of favorable inter-chain amyloid contacts via nucleated conformational conversion. The dynamic role of water in protein aggregation will have much broader implications in a variety of protein misfolding diseases.

Conformational Changes of Trialanine in Water Induced by Vibrational Relaxation of the Amide I Mode

Most of the protein-based diseases are caused by anomalies in the functionality and stability of these molecules. Experimental and theoretical studies of the conformational dynamics of proteins are becoming in this respect essential to understand the origin of these anomalies. However, a description of the conformational dynamics of proteins based on mechano-energetic principles still remains elusive because of the intrinsic high flexibility of the peptide chains, the participation of weak noncovalent interactions, and the role of the ubiquitous water solvent. In this work, the conformational dynamics of trialanine dissolved in water (D2O) is investigated through Molecular Dynamics (MD) simulations combined with instantaneous normal modes (INMs) analysis both at equilibrium and after the vibrational excitation of the C-terminal amide I mode. The conformational equilibrium between α and pPII conformers is found to be altered by the intramolecular relaxation of the amide I mode as a consequence of the different relaxation pathways of each conformer which modify the amount of vibrational energy stored in the torsional motions of the tripeptide, so the α → pPII and pPII → α conversion rates are increased differently. The selectivity of the process comes from the shifts of the vibrational frequencies with the conformational changes that modify the resonance conditions driving the intramolecular energy flows.

Perturbations in inter-domain associations may trigger the onset of pathogenic transformations in PrP(C): insights from atomistic simulations

Conversion of the predominantly α-helical cellular prion protein (PrP(C)) to the misfolded β-sheet enriched Scrapie form (PrP(Sc)) is a critical event in prion pathogenesis. However, the conformational triggers that lead to the isoform conversion (PrP(C) to PrP(Sc)) remain obscure, and conjectures about the role of unusually hydrophilic, short helix H1 of the C-terminal globular domain in the transition are varied. Helix H1 is anchored to helix H3 via a few stabilizing polar interactions. We have employed fully atomistic molecular dynamics simulations to study the effects triggered by a minor perturbation in the network of these non-bonded interactions in PrP(C). The elimination of just one of the key H1-H3 hydrogen bonds led to a cascade of conformational changes that are consistent with those observed in partially unfolded intermediates of PrP(C), with pathogenic mutations and in low pH environments. Our analyses reveal that the perturbation results in the enhanced conformational flexibility of the protein. The resultant enhancement in the dynamics leads to overall increased solvent exposure of the hydrophobic core residues and concomitant disruption of the H1-H3 inter-domain salt bridge network. This study lends credence to the hypothesis that perturbing the cooperativity of the stabilizing interactions in the PrP(C) globular domain can critically affect its dynamics and may lead to structural transitions of pathological relevance.

Stay Wet, Stay Stable? How Internal Water Helps the Stability of Thermophilic Proteins

We present a systematic computational investigation of the internal hydration of a set of homologous proteins of different stability content and molecular complexities. The goal of the study is to verify whether structural water can be part of the molecular mechanisms ensuring enhanced stability in thermophilic enzymes. Our free-energy calculations show that internal hydration in the thermophilic variants is generally more favorable, and that the cumulated effect of wetting multiple sites results in a meaningful contribution to stability. Moreover, thanks to a more effective capability to retain internal water, some thermophilic proteins benefit by a systematic gain from internal wetting up to their optimal working temperature. Our work supports the idea that internal wetting can be viewed as an alternative molecular variable to be tuned for increasing protein stability.

The Volumetric Diversity of Misfolded Prion Protein Oligomers Revealed by Pressure Dissociation

Protein oligomerization has been associated with a wide range of diseases. High pressure approaches offer a powerful tool for deciphering the underlying molecular mechanisms by revealing volume changes associated with the misfolding and assembly reactions. We applied high pressure to induce conformational changes in three distinct β-sheet-rich oligomers of the prion protein PrP, a protein characterized by a variety of infectious quaternary structures that can propagate stably and faithfully and cause diseases with specific phenotypic traits. We show that pressure induces dissociation of the oligomers and leads to a lower volume monomeric PrP state that refolds into the native conformation after pressure release. By measuring the different pressure and temperature sensitivity of the tested PrP oligomers, we demonstrate significantly different void volumes in their quaternary structure. In addition, by focusing on the kinetic and energetic behavior of the pressure-induced dissociation of one specific PrP oligomer, we reveal a large negative activation volume and an increase in both apparent activation enthalpy and entropy. This suggests a transition state ensemble that is less structured and significantly more hydrated than the oligomeric state. Finally, we found that site-specific fluorescent labeling allows monitoring of the transient population of a kinetic intermediate in the dissociation reaction. Our results indicate that defects in atomic packing may deserve consideration as a new factor that influences differences between PrP assemblies and that could be relevant also for explaining the origin of prion strains.

Role of Internal Water on Protein Thermal Stability: The Case of Homologous G Domains

In this work, we address the question of whether the enhanced stability of thermophilic proteins has a direct connection with internal hydration. Our model systems are two homologous G domains of different stability: the mesophilic G domain of the elongation factor thermal unstable protein from E. coli and the hyperthermophilic G domain of the EF-1α protein from S. solfataricus. Using molecular dynamics simulation at the microsecond time scale, we show that both proteins host water molecules in internal cavities and that these molecules exchange with the external solution in the nanosecond time scale. The hydration free energy of these sites evaluated via extensive calculations is found to be favorable for both systems, with the hyperthermophilic protein offering a slightly more favorable environment to host water molecules. We estimate that, under ambient conditions, the free energy gain due to internal hydration is about 1.3 kcal/mol in favor of the hyperthermophilic variant. However, we also find that, at the high working temperature of the hyperthermophile, the cavities are rather dehydrated, meaning that under extreme conditions other molecular factors secure the stability of the protein. Interestingly, we detect a clear correlation between the hydration of internal cavities and the protein conformational landscape. The emerging picture is that internal hydration is an effective observable to probe the conformational landscape of proteins. In the specific context of our investigation, the analysis confirms that the hyperthermophilic G domain is characterized by multiple states and it has a more flexible structure than its mesophilic homologue.

Design and application of implicit solvent models in biomolecular simulations

•

Implicit solvent replaces explicit water by a potential of mean force.

•

Popular models are SASA, VOL and Generalized Born.

•

Implicit solvent is used in MD, protein modelling, folding, design, prediction and drug screening.

•

Large-scale simulations allow for parametrisation via force matching.

•

Application to nucleic acids and membranes is challenging.

We review implicit solvent models and their parametrisation by introducing the concepts and recent devlopments of the most popular models with a focus on parametrisation via force matching. An overview of recent applications of the solvation energy term in protein dynamics, modelling, design and prediction is given to illustrate the usability and versatility of implicit solvation in reproducing the physical behaviour of biomolecular systems. Limitations of implicit modes are discussed through the example of more challenging systems like nucleic acids and membranes.

Searching for protein binding sites from Molecular Dynamics simulations and paramagnetic fragment-based NMR studies

Hotspot delineation on protein surfaces represents a fundamental step for targeting protein-protein interfaces. Disruptors of protein-protein interactions can be designed provided that the sterical features of binding pockets, including the transient ones, can be defined. Molecular Dynamics, MD, simulations have been used as a reliable framework for identifying transient pocket openings on the protein surface. Accessible surface area and intramolecular H-bond involvement of protein backbone amides are proposed as descriptors for characterizing binding pocket occurrence and evolution along MD trajectories. TEMPOL induced paramagnetic perturbations on (1)H-(15)N HSQC signals of protein backbone amides have been analyzed as a fragment-based search for surface hotspots, in order to validate MD predicted pockets. This procedure has been applied to CXCL12, a small chemokine responsible for tumor progression and proliferation. From combined analysis of MD data and paramagnetic profiles, two CXCL12 sites suitable for the binding of small molecules were identified. One of these sites is the already well characterized CXCL12 region involved in the binding to CXCR4 receptor. The other one is a transient pocket predicted by Molecular Dynamics simulations, which could not be observed from static analysis of CXCL12 PDB structures. The present results indicate how TEMPOL, instrumental in identifying this transient pocket, can be a powerful tool to delineate minor conformations which can be highly relevant in dynamic discovery of antitumoral drugs.

Structural features of the regulatory ACT domain of phenylalanine hydroxylase

Phenylalanine hydroxylase (PAH) catalyzes the conversion of L-Phe to L-Tyr. Defects in PAH activity, caused by mutations in the human gene, result in the autosomal recessively inherited disease hyperphenylalaninemia. PAH activity is regulated by multiple factors, including phosphorylation and ligand binding. In particular, PAH displays positive cooperativity for L-Phe, which is proposed to bind the enzyme on an allosteric site in the N-terminal regulatory domain (RD), also classified as an ACT domain. This domain is found in several proteins and is able to bind amino acids. We used molecular dynamics simulations to obtain dynamical and structural insights into the isolated RD of PAH. Here we show that the principal motions involve conformational changes leading from an initial open to a final closed domain structure. The global intrinsic motions of the RD are correlated with exposure to solvent of a hydrophobic surface, which corresponds to the ligand binding-site of the ACT domain. Our results strongly suggest a relationship between the Phe-binding function and the overall dynamic behaviour of the enzyme. This relationship may be affected by structure-disturbing mutations. To elucidate the functional implications of the mutations, we investigated the structural effects on the dynamics of the human RD PAH induced by six missense hyperphenylalaninemia-causing mutations, namely p.G46S, p.F39C, p.F39L, p.I65S, p.I65T and p.I65V. These studies showed that the alterations in RD hydrophobic interactions induced by missense mutations could affect the functionality of the whole enzyme.