Recipes for Inducing Cold Denaturation in an Otherwise Stable Protein.

Although cold denaturation is a fundamental phenomenon common to all proteins, it can only be observed in a handful of cases where it occurs at temperatures above the freezing point of water. Understanding the mechanisms that determine cold denaturation and the rules that permit its observation is an important challenge. A way to approach them is to be able to induce cold denaturation in an otherwise stable protein by means of mutations. Here, we studied CyaY, a relatively stable bacterial protein with no detectable cold denaturation and a high melting temperature of 54 °C. We have characterized for years the yeast orthologue of CyaY, Yfh1, a protein that undergoes cold and heat denaturation at 5 and 35 °C, respectively. We demonstrate that, by transferring to CyaY the lessons learnt from Yfh1, we can induce cold denaturation by introducing a restricted number of carefully designed mutations aimed at destabilizing the overall fold and inducing electrostatic frustration. We used molecular dynamics simulations to rationalize our findings and demonstrate the individual effects observed experimentally with the various mutants. Our results constitute the first example of rationally designed cold denaturation and demonstrate the importance of electrostatic frustration on the mechanism of cold denaturation.

Chemical Chaperones Modulate the Formation of Metabolite Assemblies.

The formation of amyloid-like structures by metabolites is associated with several inborn errors of metabolism (IEMs). These structures display most of the biological, chemical and physical properties of protein amyloids. However, the molecular interactions underlying the assembly remain elusive, and so far, no modulating therapeutic agents are available for clinical use. Chemical chaperones are known to inhibit protein and peptide amyloid formation and stabilize misfolded enzymes. Here, we provide an in-depth characterization of the inhibitory effect of osmolytes and hydrophobic chemical chaperones on metabolite assemblies, thus extending their functional repertoire. We applied a combined in vivo-in vitro-in silico approach and show their ability to inhibit metabolite amyloid-induced toxicity and reduce cellular amyloid content in yeast. We further used various biophysical techniques demonstrating direct inhibition of adenine self-assembly and alteration of fibril morphology by chemical chaperones. Using a scaffold-based approach, we analyzed the physiochemical properties of various dimethyl sulfoxide derivatives and their role in inhibiting metabolite self-assembly. Lastly, we employed whole-atom molecular dynamics simulations to elucidate the role of hydrogen bonds in osmolyte inhibition. Our results imply a dual mode of action of chemical chaperones as IEMs therapeutics, that could be implemented in the rational design of novel lead-like molecules.

Crystal Structure of the DFNKF Segment of Human Calcitonin Unveils Aromatic Interactions between Phenylalanines.

Although intensively studied, the high-resolution crystal structure of the peptide DFNKF, the core-segment of human calcitonin, has never been described. Here we report how the use of iodination as a strategy to promote crystallisation and facilitate phase determination, allowed us to solve, for the first time, the single-crystal X-ray structure of a DFNKF derivative. Computational studies suggest that both the iodinated and the wild-type peptides populate very similar conformations. Furthermore, the conformer found in the solid-state structure is one of the most populated in solution, making the crystal structure a reliable model for the peptide in solution. The crystal structure of DFNKF(I) confirms the overall features of the amyloid cross-ß spine and highlights how aromatic-aromatic interactions are important structural factors in the self-assembly of this peptide. A detailed analysis of such interactions is reported.

Rational design of shepherdin, a novel anticancer agent.

Anticancer agents that selectively kill tumor cells and spare normal tissues are urgently needed. Here, we engineered a cell-permeable peptidomimetic, shepherdin, modeled on the binding interface between the molecular chaperone Hsp90 and the antiapoptotic and mitotic regulator, survivin. Shepherdin makes extensive contacts with the ATP pocket of Hsp90, destabilizes its client proteins, and induces massive death of tumor cells by apoptotic and nonapoptotic mechanisms. Conversely, shepherdin does not reduce the viability of normal cells, and does not affect colony formation of purified hematopoietic progenitors. Systemic administration of shepherdin in vivo is well tolerated, and inhibits human tumor growth in mice without toxicity. Shepherdin could provide a potent and selective anticancer agent in humans.

Pharmacological enhancement of α-glucosidase by the allosteric chaperone N-acetylcysteine.

Pompe disease (PD) is a metabolic myopathy due to the deficiency of the lysosomal enzyme α-glucosidase (GAA). The only approved treatment for this disorder, enzyme replacement with recombinant human GAA (rhGAA), has shown limited therapeutic efficacy in some PD patients. Pharmacological chaperone therapy (PCT), either alone or in combination with enzyme replacement, has been proposed as an alternative therapeutic strategy. However, the chaperones identified so far also are active site-directed molecules and potential inhibitors of target enzymes. We demonstrated that N-acetylcysteine (NAC) is a novel allosteric chaperone for GAA. NAC improved the stability of rhGAA as a function of pH and temperature without disrupting its catalytic activity. A computational analysis of NAC-GAA interactions confirmed that NAC does not interact with GAA catalytic domain. NAC enhanced the residual activity of mutated GAA in cultured PD fibroblasts and in COS7 cells overexpressing mutated GAA. NAC also enhanced rhGAA efficacy in PD fibroblasts. In cells incubated with NAC and rhGAA, GAA activities were 3.7-8.7-fold higher than those obtained in cells treated with rhGAA alone. In a PD mouse model the combination of NAC and rhGAA resulted in better correction of enzyme activity in liver, heart, diaphragm and gastrocnemia, compared to rhGAA alone.

A Hamiltonian replica exchange molecular dynamics (MD) method for the study of folding, based on the analysis of the stabilization determinants of proteins.

Herein, we present a novel Hamiltonian replica exchange protocol for classical molecular dynamics simulations of protein folding/unfolding. The scheme starts from the analysis of the energy-networks responsible for the stabilization of the folded conformation, by means of the energy-decomposition approach. In this framework, the compact energetic map of the native state is generated by a preliminary short molecular dynamics (MD) simulation of the protein in explicit solvent. This map is simplified by means of an eigenvalue decomposition. The highest components of the eigenvector associated with the lowest eigenvalue indicate which sites, named "hot spots", are likely to be responsible for the stability and correct folding of the protein. In the Hamiltonian replica exchange protocol, we use modified force-field parameters to treat the interparticle non-bonded potentials of the hot spots within the protein and between protein and solvent atoms, leaving unperturbed those relative to all other residues, as well as solvent-solvent interactions. We show that it is possible to reversibly simulate the folding/unfolding behavior of two test proteins, namely Villin HeadPiece HP35 (35 residues) and Protein A (62 residues), using a limited number of replicas. We next discuss possible implications for the study of folding mechanisms via all atom simulations.

Entropically-Driven Co-assembly of l-Histidine and l-Phenylalanine to Form Supramolecular Materials.

Molecular self- and co-assembly allow the formation of diverse and well-defined supramolecular structures with notable physical properties. Among the associating molecules, amino acids are especially attractive due to their inherent biocompatibility and simplicity. The biologically active enantiomer of l-histidine (l-His) plays structural and functional roles in proteins but does not self-assemble to form discrete nanostructures. In order to expand the structural space to include l-His-containing materials, we explored the co-assembly of l-His with all aromatic amino acids, including phenylalanine (Phe), tyrosine (Tyr), and tryptophan (Trp), all in both enantiomeric forms. In contrast to pristine l-His, the combination of this building block with all aromatic amino acids resulted in distinct morphologies including fibers, rods, and flake-like structures. Electrospray ionization mass spectrometry (ESI-MS) indicated the formation of supramolecular co-assemblies in all six combinations, but time-of-flight secondary-ion mass spectrometry (ToF-SIMS) indicated the best seamless co-assembly occurs between l-His and l-Phe while in the other cases, different degrees of phase separation could be observed. Indeed, isothermal titration calorimetry (ITC) suggested the highest affinity between l-His and l-Phe where the formation of co-assembled structures was driven by entropy. In accordance, among all the combinations, the co-assembly of l-His and l-Phe produced single crystals. The structure revealed the formation of a 3D network with nanocavities stabilized by hydrogen bonding between -N (l-His) and -NH (l-Phe). Taken together, using the co-assembly approach we expanded the field of amino acid nanomaterials and showed the ability to obtain discrete supramolecular nanostructures containing l-His based on its specific interactions with l-Phe.

A multiscale approach to characterize the early aggregation steps of the amyloid-forming peptide GNNQQNY from the yeast prion sup-35.

The self-organization of peptides into amyloidogenic oligomers is one of the key events for a wide range of molecular and degenerative diseases. Atomic-resolution characterization of the mechanisms responsible for the aggregation process and the resulting structures is thus a necessary step to improve our understanding of the determinants of these pathologies. To address this issue, we combine the accelerated sampling properties of replica exchange molecular dynamics simulations based on the OPEP coarse-grained potential with the atomic resolution description of interactions provided by all-atom MD simulations, and investigate the oligomerization process of the GNNQQNY for three system sizes: 3-mers, 12-mers and 20-mers. Results for our integrated simulations show a rich variety of structural arrangements for aggregates of all sizes. Elongated fibril-like structures can form transiently in the 20-mer case, but they are not stable and easily interconvert in more globular and disordered forms. Our extensive characterization of the intermediate structures and their physico-chemical determinants points to a high degree of polymorphism for the GNNQQNY sequence that can be reflected at the macroscopic scale. Detailed mechanisms and structures that underlie amyloid aggregation are also provided.

Effects of water miscible organic solvents on the activity and conformation of the Baeyer-Villiger monooxygenases from Thermobifida fusca and Acinetobacter calcoaceticus: a comparative study.

A broader exploitation of enzymes in organic synthesis can be achieved by increasing their tolerance toward organic solvents. In this study, the stability and activity of Baeyer-Villiger monooxygenases from Thermobifida fusca (PAMO) and Acinetobacter sp. (CHMO) in the presence of water miscible organic solvents were compared. PAMO was more stable than CHMO. The concentration of solvent (v/v) at which it halved its activity (C(50) ) was 4- to 16-fold higher than that observed for CHMO. For PAMO, the C(50) varied from 16% to 55% of solvent and followed the destabilizing order methanol < ethanol < 1,4-dioxane < acetonitrile < trifluoroethanol. In the case of CHMO, the maximal C(50) was 7% with methanol and even lower with the other solvents. Therefore, methanol was the most tolerated solvent. In the case of PAMO, methanol induced a significant increase of enzyme activity (up to fivefold), which was optimal at 20% (v/v) solvent. Only minor spectral variations were observed with PAMO in 20% methanol, suggesting that the increase of activity observed in this condition is not due to marked conformational changes. Fluorescence and circular dichroism analyses showed that the lower stability of CHMO toward organic solvent correlates with a more pronounced destructive effect on its secondary and tertiary structure. A possible rationale for the higher stability of PAMO could be inferred from inspection of the PAMO and CHMO (two enzymes of similar size) structure, which revealed a higher (up to twofold) number of ionic bridges in PAMO with respect to CHMO.

The structural intolerance of the PrP alpha-fold for polar substitution of the helix-3 methionines.

The conversion of the cellular prion protein (PrP(C)) into its disease-associated form (PrP(Sc)) involves a major conformational change and the accumulation of sulfoxidized methionines. Computational and synthetic approaches have shown that this change in the polarity of M206 and M213 impacts the C-terminal domain native alpha-fold allowing the flexibility required for the structural conversion. To test the effect in the full-length molecule with site-specificity, we have generated M-to-S mutations. Molecular dynamics simulations show that the replacement indeed perturbs the native state. When this mutation is placed at the conserved methionines of HaPrP(23-231), only substitutions at the Helix-3 impair the alpha-fold, stabilizing a non-native state with perturbed secondary structure, loss of native tertiary contacts, increased surface hydrophobicity, reduced thermal stability and an enhanced tendency to aggregate into protofibrillar polymers. Our work supports that M206 and M213 function as alpha-fold gatekeepers and suggests that their redox state regulate misfolding routes.

Understanding ligand-based modulation of the Hsp90 molecular chaperone dynamics at atomic resolution.

Molecular switching and ligand-based modulation of the 90-kDa heat-shock protein (Hsp90) chaperone activity may ultimately facilitate conformational coupling to the ATPase cycle along with activation and recruitment of the broad range of client proteins. We present an atomic resolution analysis of the Hsp90 N-terminal domain (NTD) binding energy landscape by simulating protein dynamics with a range of binding partners. We show that the activity of the molecular chaperone may be linked to (i) local folding-unfolding transitions and conformational switching of the "active site lid" upon binding and (ii) differences in the underlying protein dynamics as a function of the binding partner. This study suggests that structural plasticity of the Hsp90 NTD can be exploited by the molecular chaperone machinery to modulate enhanced structural rigidity during ATP binding and increased protein flexibility as a consequence of the inhibitor binding. The present study agrees with the experimental structural data and provides a plausible molecular model for understanding mechanisms of modulation of molecular chaperone activities by binding partners.

Revealing Escherichia coli type II L-asparaginase active site flexible loop in its open, ligand-free conformation.

Since 1993, when the structure of Escherichia coli type II L-asparaginase (EcAII) in complex with L-aspartate was firstly reported, many structures of the wild type and mutated enzyme have been deposited in the Protein Data Bank. None of them report the full structure of the monomer in its ligand-free, open conformation, mainly because of the high dynamic and flexibility of the active site flexible loop. Here we report for the first time the structure of EcAII wild type in its open conformation comprising, for at least one protomer, clear electron density for the active site flexible loop (PDB ID: 6YZI). The structural element is highly mobile and it is transposed onto the rigid part of the active site upon substrate binding to allow completion of the enzyme catalytic center, thanks to key residues that serve as hinges and anchoring points. In the substrate binding pocket, several highly conserved water molecules are coordinated by residues involved in substrate binding, comprising two water molecules very likely involved in the enzyme catalytic process. We also describe, by molecular dynamics simulations, how the transposition of the loop, besides providing the proximity of residues needed for catalysis, causes a general stabilization of the protein.

Simple Model of Protein Energetics To Identify Ab Initio Folding Transitions from All-Atom MD Simulations of Proteins.

A fundamental requirement to predict the native conformation, address questions of sequence design and optimization, and gain insights into the folding mechanisms of proteins lies in the definition of an unbiased reaction coordinate that reports on the folding state without the need to compare it to reference values, which might be unavailable for new (designed) sequences. Here, we introduce such a reaction coordinate, which does not depend on previous structural knowledge of the native state but relies solely on the energy partition within the protein: the spectral gap of the pair nonbonded energy matrix (ENergy Gap, ENG). This quantity can be simply calculated along unbiased MD trajectories. We show that upon folding the gap increases significantly, while its fluctuations are reduced to a minimum. This is consistently observed for a diverse set of systems and trajectories. Our approach allows one to promptly identify residues that belong to the folding core as well as residues involved in non-native contacts that need to be disrupted to guide polypeptides to the folded state. The energy gap and fluctuations criteria are then used to develop an automatic detection system which allows us to extract and analyze folding transitions from a generic MD trajectory. We speculate that our method can be used to detect conformational ensembles in dynamic and intrinsically disordered proteins, revealing potential preorganization for binding.

The Answer Lies in the Energy: How Simple Atomistic Molecular Dynamics Simulations May Hold the Key to Epitope Prediction on the Fully Glycosylated SARS-CoV-2 Spike Protein.

SARS-CoV-2 is a health threat with dire socioeconomical consequences. As the crucial mediator of infection, the viral glycosylated spike protein (S) has attracted the most attention and is at the center of efforts to develop therapeutics and diagnostics. Herein, we use an original decomposition approach to identify energetically uncoupled substructures as antibody binding sites on the fully glycosylated S. Crucially, all that is required are unbiased MD simulations; no prior knowledge of binding properties or ad hoc parameter combinations is needed. Our results are validated by experimentally confirmed structures of S in complex with anti- or nanobodies. We identify poorly coupled subdomains that are poised to host (several) epitopes and potentially involved in large functional conformational transitions. Moreover, we detect two distinct behaviors for glycans: those with stronger energetic coupling are structurally relevant and protect underlying peptidic epitopes, and those with weaker coupling could themselves be prone to antibody recognition.

Mechanisms of Metabolite Amyloid Formation: Computational Studies for Drug Design against Metabolic Disorders.

Ordered self-organization of polypeptides into fibrillar assemblies has been associated with a number of pathological conditions linked to degenerative diseases. Recent experimental observations have demonstrated that even small-molecule metabolites can aggregate into supramolecular arrangements with structural and functional properties reminiscent of peptide-based amyloids. The molecular determinants of such mechanisms, however, are not clear yet. Herein, we examine the process of formation of ordered aggregates by adenine in aqueous solution by molecular dynamics simulations. We also investigate the effects of an inhibiting polyphenol, namely, epigallocatechin gallate (EGCG), on this mechanism. We show that, while adenine alone is able to form extended amyloid-like oligomers, EGCG interferes with the supramolecular organization process. Interestingly, acetylsalicylic acid is shown not to interfere with ordered aggregation, consistent with experiments. The results of these mechanistic studies indicate the main pharmacophoric determinants that a drug-like inhibitor should possess to effectively interfere with metabolite amyloid formation.

Investigating the mechanism of peptide aggregation: insights from mixed monte carlo-molecular dynamics simulations.

The early stages of peptide aggregation are currently not accessible by experimental techniques at atomic resolution. In this article, we address this problem through the application of a mixed simulation scheme in which a preliminary coarse-grained Monte Carlo analysis of the free-energy landscape is used to identify representative conformations of the aggregates and subsequent all-atom molecular dynamics simulations are used to analyze in detail possible pathways for the stabilization of oligomers. This protocol was applied to systems consisting of multiple copies of the model peptide GNNQQNY, whose detailed structures in the aggregated state have been recently solved in another study. The analysis of the various trajectories provides dynamical and structural insight into the details of aggregation. In particular, the simulations suggest a hierarchical mechanism characterized by the initial formation of stable parallel beta-sheet dimers and identify the formation of the polar zipper motif as a fundamental feature for the stabilization of initial oligomers. Simulation results are consistent with experimentally derived observations and provide an atomically detailed view of the putative initial stages of fibril formation.

Computational studies of the structure, dynamics and native content of amyloid-like fibrils of ribonuclease A.

The characterization at atomic resolution of amyloid-like protein aggregates is one of the fundamental problems of modern biology. In particular, the question whether native-like domains are retained or completely refolded in the amyloid state and the identification of possible mechanisms for macromolecular ordered aggregation represent major unresolved puzzles. To address these issues, in this article we examine the stability, dynamics, and conservation of native-like properties of several models of a previously designed amyloid-like fibril of RNase A (Sambashivan et al., Nature 2005; 437:266-269). Through the use of molecular dynamics (MD) simulations, we have provided molecular-level insights into the role of different parts of the sequence on the stability of fibrils, the collective properties of supramolecular complexes, and the presence of native-like conformations and dynamics in supramolecular aggregates. We have been able to show that within the fibrils the three-dimensional globular domain-swapped units preserve the conformational, dynamical, and hydration properties typical of the monomeric state, providing a rationalization for the experimentally observed catalytic activity of fibrils. The nativeness of the globular domains is not affected by the amyloidogenic stretches, which determine the molecular recognition process underlying aggregation through the formation of a stable steric zipper motif. Moreover, through the study of the hydration features of a single sheet model, we have been able to show that polyglutamine stretches of the domain-swapped ribonuclease tend to minimize the interaction with water in favor of sidechain-sidechain interactions, shedding light on the factors leading to the supramolecular assembly of beta-sheet layers into dry steric zippers.

Molecular dynamics simulations of proteins and peptides: from folding to drug design.

Computer simulations of proteins, lipids and nucleic acids at equilibrium have become essentially routine. However, the fact remains that complete sampling of conformational space continues to be a bottle-neck in the field. The challenge for the future is to overcome such problems and use computational approaches to understand recognition and spontaneous self-organization in biomolecular systems (folding, aggregation and assembly of complexes), processes that cannot be directly observed experimentally. In this review, examples illustrating the extent to which simulations can be used to understand these phenomena in biomolecular systems will be presented along with examples of methodological developments to increase our physical understanding of the processes. The study cases will cover the problems of peptide-receptor recognition and the use of the information obtained for the design of new non-peptidic ligands; the study of the folding mechanism of small proteins and finally the study of the initial stages of peptide self-aggregation.

How the Ligand-Induced Reorganization of Protein Internal Energies Is Coupled to Conformational Events.

Here, we introduce a novel computational method to identify the protein substructures most likely to support the functionally oriented structural deformations that occur upon ligand-binding. To this aim, we study the modulation of protein energetics along the trajectory of a molecular dynamics simulation of different proteins in the presence and in the absence of their respective ligands, namely, human FGF, human second PDZ from human PTP1E/PTPL1, and the N terminal domain of human Hsp90. The method is based on the idea that a subset of protein residues (hotspots) may initiate the global response via the disassembly and reassembly of interactions, which is reflected in the modulation of the overall protein energetics. To identify structural hotspots and dynamic states linked to the onset of functionally relevant conformational transitions, we define an energy profile to monitor the protein energetics, based on a previously introduced approach that highlights the essential nonbonded couplings among all residues. The energy profiles are calculated along the trajectory to yield a time-dependent evolution, and their relative population in the presence and absence of the ligand is evaluated by means of a clustering procedure. It is found that interconversion between clusters, as well as their population and the density of specific energy profiles in the vicinity of structural transitions, provides specific information on the impact of the ligand in driving the protein conformational response. This analysis also highlights the hotspot residues that are most responsive to the presence of the ligand. Importantly, identified hotspots are in agreement with experimental evidence in the three considered systems. We propose that this approach can be generally used in the prediction of "allosteric hotspots" and ligand-induced conformational responses, as well as to select conformations more likely to support functional transitions (e.g., in the framework of adaptive sampling approaches).

Are Amyloid Fibrils RNA-Traps? A Molecular Dynamics Perspective.

The self-assembly of proteins and peptides into amyloids is a key feature of an increasing number of diseases. Amyloid fibrils display a unique surface reactivity endowing the sequestration of molecules such as MicroRNAs, which can be the active moiety of the toxic action. To test this hypothesis we studied the recognition between a model RNA and two different steric zipper spines using molecular dynamics simulations. We found that the interaction occurs and displays peptide-sequence dependence. Interestingly, interactions with polar zipper surfaces such as the formed by SNQNNF are more stable and favor the formation of ß-barrel like complexes resembling the structures of toxic oligomers. These sequence-structure-recognition relationships of the two different assemblies may be exploited for the design of compounds targeting the fibers or competing with RNA-amyloid attachment.