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.

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.

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.

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.

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.

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.

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.

A protease-resistant Escherichia coli asparaginase with outstanding stability and enhanced anti-leukaemic activity in vitro.

L-Asparaginases (ASNases) have been used as first line drugs for paediatric Acute Lymphoblastic Leukaemia (ALL) treatment for more than 40 years. Both the Escherichia coli (EcAII) and Erwinia chrysanthemi (ErAII) type II ASNases currently used in the clinics are characterized by high in vivo instability, short half-life and the requirement of several administrations to obtain a pharmacologically active concentration. Moreover, they are sensitive to proteases (cathepsin B and asparagine endopeptidase) that are over-expressed by resistant leukaemia lymphoblasts, thereby impairing drug activity and pharmacokinetics. Herein, we present the biochemical, structural and in vitro antiproliferative characterization of a new EcAII variant, N24S. The mutant shows completely preserved asparaginase and glutaminase activities, long-term storage stability, improved thermal parameters, and outstanding resistance to proteases derived from leukaemia cells. Structural analysis demonstrates a modification in the hydrogen bond network related to residue 24, while Normal Mode-based geometric Simulation and Molecular Dynamics predict a general rigidification of the monomer as compared to wild-type. These improved features render N24S a potential alternative treatment to reduce the number of drug administrations in vivo and to successfully address one of the major current challenges of ALL treatment: spontaneous, protease-dependent and immunological inactivation of ASNase.

Novel PARP-1 Inhibitor Scaffolds Disclosed by a Dynamic Structure-Based Pharmacophore Approach.

PARP-1 inhibition has been studied over the last decades for the treatment of various diseases. Despite the fact that several molecules act as PARP-1 inhibitors, a reduced number of compounds are used in clinical practice. To identify new compounds with a discriminatory PARP-1 inhibitory function, explicit-solvent molecular dynamics simulations using different inhibitors bound to the PARP-1 catalytic domain were performed. The representative structures obtained were used to generate structure-based pharmacophores, taking into account the dynamic features of receptor-ligand interactions. Thereafter, a virtual screening of compound databases using the pharmacophore models obtained was performed and the hits retrieved were subjected to molecular docking-based scoring. The drug-like molecules featuring the best ranking were evaluated for their PARP-1 inhibitory activity and IC50 values were calculated for the top scoring docked compounds. Altogether, three new PARP-1 inhibitor chemotypes were identified.

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.

DNA Polymerase Conformational Dynamics and the Role of Fidelity-Conferring Residues: Insights from Computational Simulations.

Herein we investigate the molecular bases of DNA polymerase I conformational dynamics that underlie the replication fidelity of the enzyme. Such fidelity is determined by conformational changes that promote the rejection of incorrect nucleotides before the chemical ligation step. We report a comprehensive atomic resolution study of wild type and mutant enzymes in different bound states and starting from different crystal structures, using extensive molecular dynamics (MD) simulations that cover a total timespan of ~5 ms. The resulting trajectories are examined via a combination of novel methods of internal dynamics and energetics analysis, aimed to reveal the principal molecular determinants for the (de)stabilization of a certain conformational state. Our results show that the presence of fidelity-decreasing mutations or the binding of incorrect nucleotides in ternary complexes tend to favor transitions from closed toward open structures, passing through an ensemble of semi-closed intermediates. The latter ensemble includes the experimentally observed ajar conformation which, consistent with previous experimental observations, emerges as a molecular checkpoint for the selection of the correct nucleotide to incorporate. We discuss the implications of our results for the understanding of the relationships between the structure, dynamics, and function of DNA polymerase I at the atomistic level.

The five-to-six-coordination transition of ferric human serum heme-albumin is allosterically-modulated by ibuprofen and warfarin: a combined XAS and MD study.

Human serum albumin (HSA) is involved physiologically in heme scavenging; in turn, heme-albumin (HSA-heme-Fe) displays globin-like properties. Here, the allosteric effect of ibuprofen and warfarin on the local atomic structure around the ferric heme-Fe (heme-Fe(III)) atom of HSA-heme-Fe (HSA-heme-Fe(III)) has been probed by Fe-K edge X-ray absorption spectroscopy (XAS). The quantitative analysis of the Fe-K edge extended X-ray absorption fine structure (EXAFS) signals and modeling of the near edge (XANES) spectral features demonstrated that warfarin and ibuprofen binding modify the local structure of the heme-Fe(III). Combined XAS data analysis and targeted molecular dynamics (MD) simulations provided atomic resolution insights of protein structural rearrangements required to accommodate the heme-Fe(III) upon ibuprofen and warfarin binding. In the absence of drugs, the heme-Fe(III) atom is penta-coordinated having distorted 4+1 configuration made by the nitrogen atoms of the porphyrin ring and the oxygen phenoxy atom of the Tyr161 residue. MD simulations show that ibuprofen and warfarin association to the secondary fatty acid (FA) binding site 2 (FA2) induces a reorientation of domain I of HSA-heme-Fe(III), this leads to the redirection of the His146 residue providing an additional bond to the heme-Fe(III) atom, providing the 5+1 configuration. The comparison of Fe-K edge XANES spectra calculated using MD structures with those obtained experimentally confirms the reliability of the proposed structural model. As a whole, combining XAS and MD simulations it has been possible to provide a reliable model of the heme-Fe(III) atom coordination state and to understand the complex allosteric transition occurring in HSA-heme-Fe(III) upon ibuprofen and warfarin binding.

Investigating the dynamic aspects of drug-protein recognition through a combination of MD and NMR analyses: implications for the development of protein-protein interaction inhibitors.

In this paper, we investigate the dynamic aspects of the molecular recognition between a small molecule ligand and a flat, exposed protein surface, representing a typical target in the development of protein-protein interaction inhibitors. Specifically, we analyze the complex between the protein Fibroblast Growth Factor 2 (FGF2) and a recently discovered small molecule inhibitor, labeled sm27 for which the binding site and the residues mainly involved in small molecule recognition have been previously characterized. We have approached this problem using microsecond MD simulations and NMR-based characterizations of the dynamics of the apo and holo states of the system. Using direct combination and cross-validation of the results of the two techniques, we select the set of conformational states that best recapitulate the principal dynamic and structural properties of the complex. We then use this information to generate a multi-structure representation of the sm27-FGF2 interaction. We propose this kind of representation and approach as a useful tool in particular for the characterization of systems where the mutual dynamic influence between the interacting partners is expected to play an important role. The results presented can also be used to generate new rules for the rational expansion of the chemical diversity space of FGF2 inhibitors.

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.

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 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.