Hydrogen Bond Networks and Hydrophobic Effects in the Amyloid β30–35 Chain in Water: A Molecular Dynamics Study

Hydrodynamic Control of Alzheimer Aβ Fibrillation with Glucosaminic Acid Containing Click-Cyclized β-Bodies

It is well established that the dynamic hydration shell plays a vital role in macromolecular functions such as protein-ligand, protein-protein, protein-DNA, and protein-lipid interactions. Here we investigate how the water modality affects conformational changes, solubility, and motion of fibrillar proteins. The hypothesis is that the introduction of a poly hydroxyl amino acid would increase solvation of the fibril forming peptides, preventing their misfolding and aggregation. For the amyloid β (Aβ) peptide, which is considered to be connected with nervous system diseases, including dementia and cognitive decline in Alzheimer's disease, the formation of β-sheet fibrils always occurs with a conformational change and a decrease in the dynamic hydration shell around Aβ(1-42). We present novel cyclic d-amino acid peptides that effectively inhibit fibrillation through affecting the dynamic hydration shell of Aβ(1-42) in vitro. Using de novo design within the software Molecular Operating Environment (MOE), five different peptides that recognize Alzheimer's fibrils were designed and synthesized. Three of them were cyclic all-d-amino acid peptides incorporating the same polyhydroxy building block derived from d-glucosaminic acid (GA). One peptide was the parent cyclic all d-amino acid inhibitor with no GA incorporated, and another was an all l-amino acid linear fibrillation inhibitor. The GA-containing peptides were found to show significantly improved inhibition of Aβ(1-42) aggregation. The inhibition was dramatically improved by the synergistic application of two GA peptides targeting each end of the growing fibril. The present study may facilitate future developments of intervention strategies for Alzheimer's disease and similar neurodegenerative diseases.

N-glycosylation induced changes in tau protein dynamics reveal its role in tau misfolding and aggregation: A microsecond long molecular dynamics study

Various posttranslational modifications like hyperphosphorylation, O-GlcNAcylation, and acetylation have been attributed to induce the abnormal folding in tau protein. Recent in vitro studies revealed the possible involvement of N-glycosylation of tau protein in the abnormal folding and tau aggregation. Hence, in this study, we performed a microsecond long all atom molecular dynamics simulation to gain insights into the effects of N-glycosylation on Asn-359 residue which forms part of the microtubule binding region. Trajectory analysis of the stimulations coupled with essential dynamics and free energy landscape analysis suggested that tau, in its N-glycosylated form tends to exist in a largely folded conformation having high beta sheet propensity as compared to unmodified tau which exists in a large extended form with very less beta sheet propensity. Residue interaction network analysis of the lowest energy conformations further revealed that Phe378 and Lys353 are the functionally important residues in the peptide which helped in initiating the folding process and Phe378, Lys347, and Lys370 helped to maintain the stability of the protein in the folded state.

Intracellular Aβ42 Aggregation Leads to Cellular Thermogenesis

The aggregation of Aβ42 is a hallmark of Alzheimer’s disease. It is still not known what the biochemical changes are inside a cell which will eventually lead to Aβ42 aggregation. Thermogenesis has been associated with cellular stress, the latter of which may promote aggregation. We perform intracellular thermometry measurements using fluorescent polymeric thermometers to show that Aβ42 aggregation in live cells leads to an increase in cell-averaged temperatures. This rise in temperature is mitigated upon treatment with an aggregation inhibitor of Aβ42 and is independent of mitochondrial damage that can otherwise lead to thermogenesis. With this, we present a diagnostic assay which could be used to screen small-molecule inhibitors to amyloid proteins in physiologically relevant settings. To interpret our experimental observations and motivate the development of future models, we perform classical molecular dynamics of model Aβ peptides to examine the factors that hinder thermal dissipation. We observe that this is controlled by the presence of ions in its surrounding environment, the morphology of the amyloid peptides, and the extent of its hydrogen-bonding interactions with water. We show that aggregation and heat retention by Aβ peptides are favored under intracellular-mimicking ionic conditions, which could potentially promote thermogenesis. The latter will, in turn, trigger further nucleation events that accelerate disease progression.

Short hydrogen bonds enhance nonaromatic protein-related fluorescence

Intrinsic fluorescence of nonaromatic amino acids is a puzzling phenomenon with an enormous potential in biophotonic applications. The physical origins of this effect, however, remain elusive. Herein, we demonstrate how specific hydrogen bond networks can modulate fluorescence. We highlight the key role played by short hydrogen bonds, present in the protein structure, on the ensuing fluorescence. We provide detailed experimental and molecular evidence to explain these unusual nonaromatic optical properties. Our findings should benefit the design of novel optically active biomaterials for applications in biosensing and imaging.

Fluorescence in biological systems is usually associated with the presence of aromatic groups. Here, by employing a combined experimental and computational approach, we show that specific hydrogen bond networks can significantly affect fluorescence. In particular, we reveal that the single amino acid L-glutamine, by undergoing a chemical transformation leading to the formation of a short hydrogen bond, displays optical properties that are significantly enhanced compared with L-glutamine itself. Ab initio molecular dynamics simulations highlight that these short hydrogen bonds prevent the appearance of a conical intersection between the excited and the ground states and thereby significantly decrease nonradiative transition probabilities. Our findings open the door to the design of new photoactive materials with biophotonic applications.

Fibril fragments from the amyloid core of lysozyme: An accelerated molecular dynamics study

Protein aggregation and formation of amyloid fibrils are associated with many diseases and present a ubiquitous problem in protein science. Hen egg white lysozyme (HEWL) can form fibrils both from the full length protein and from its fragments. In the present study, we simulated unfolding of the amyloidogenic fragment of HEWL encompassing residues 49-101 to study the conformational aspects of amyloidogenesis. The accelerated molecular dynamics approach was used to speed up the sampling of the fragment conformers under enhanced temperature. Analysis of conformational transformation and intermediate structures was performed. During the unfolding, the novel short-living and long-living β-structures are formed along with the unstructured random coils. Such β-structure enriched monomers can interact with each other and propagate into fibril-like forms. The stability of oligomers assembled from these monomers was evaluated in the course of MD simulations with explicit water. The residues playing a key role in fibril stabilization were determined. The work provides new insights into the processes occurring at the early stages of amyloid fibril assembly.

Intracellular Thermometry at the Micro-/Nanoscale and its Potential Application to Study Protein Aggregation Related to Neurodegenerative Diseases

Temperature is a fundamental physical parameter that influences biological processes in living cells. Hence, intracellular temperature mapping can be used to derive useful information reflective of thermodynamic properties and cellular behaviour. Herein, existing publications on different thermometry systems, focusing on those that employ fluorescence-based techniques, are reviewed. From developments based on fluorescent proteins and inorganic molecules to metal nanoclusters and fluorescent polymers, the general findings of intracellular measurements from different research groups are discussed. Furthermore, the contradiction of mitochondrial thermogenesis and nuclear-cytoplasmic temperature differences to current thermodynamic understanding are highlighted. Lastly, intracellular thermometry is proposed as a tool to quantify the energy flow and cost associated with amyloid-β42 (Aβ42) aggregation, a hallmark of Alzheimer's disease.

A computational study on the role of water and conformational fluctuations in Hsp90 in response to inhibitors

Molecular chaperone Heat Shock Protein 90 (Hsp90) represents an interesting chemotherapeutic target for cancer treatments as it plays a role in cancer proliferation. Thus, continued effort to identify novel inhibitors of this target is an important task. Drug design using computational approach has gained significant attention in recent years. This work aims to propose docking protocols to re-purpose FDA-approved drugs targeting Hsp90. Sensitivity of results to different docking protocols such apo, holo and receptor ensembles (relaxed complex) structures, the role of water and conformational changes of Hsp90, are described. We show that the protein conformation and water have effects on drug binding. Holo relaxed complex receptors ensembles improves the binding energy of ligands to the protein. We also compare and contrast structural stability of three drugs namely: ezetimibe, pitavastatin and vilazodon in the Hsp90 protein. The results obtained serves as a possible basis towards developing Hsp90 inhibitors.

Low energy optical excitations as an indicator of structural changes initiated at the termini of amyloid proteins

There is a growing body of experimental work showing that protein aggregates associated with amyloid fibrils feature intrinsic fluorescence. In order to understand the microscopic origin of this behavior observed in non-aromatic aggregates of peptides and proteins, we conducted a combined experimental and computational study on the optical properties of amyloid-derived oligopeptides in the near-UV region. We have focused on a few model systems having charged termini (zwitterionic) or acetylated termini. For the zwitterionic system, we were able to simulate the longer tail absorption in the near UV (250-350 nm), supporting the experimental results in terms of excitation spectra. We analyzed the optical excitations responsible for the low-energy absorption and found a large role played by charge-transfer states around the termini. These charge-transfer excitations are very sensitive to the conformation of the peptide and in realistic fibrils may involve inter and intra chain charge reorganization.

Prediction of amyloid aggregation rates by machine learning and feature selection

A novel data-based machine learning algorithm for predicting amyloid aggregation rates is reported in this paper. Based on a highly nonlinear projection from 16 intrinsic features of a protein and 4 extrinsic features of the environment to the protein aggregation rate, a feedforward fully connected neural network (FCN) with one hidden layer is trained on a dataset composed of 21 different kinds of amyloid proteins and tested on 4 rest proteins. FCN shows a much better performance than traditional algorithms, such as multivariable linear regression and support vector regression, with an average accuracy higher than 90%. Furthermore, by the correlation analysis and the principal component analysis, seven key features, folding energy, HP patterns for helix, sheet and helices cross membrane, pH, ionic strength, and protein concentration, are shown to constitute a minimum feature set for characterizing the amyloid aggregation kinetics.

Underestimated Noncovalent Interactions in Protein Data Bank

Noncovalent interactions (NCIs) play essential roles in the structure and function of biomacromolecules. There are various NCIs, e.g., hydrogen bonds (HBs), cation-π and π-π interactions, and ionic bonds, among which HBs are the most widespread and well-studied. By utilizing the ratio of the observed HBs over pseudo HBs (1.0 Å longer than the HB distance criteria without angle constraints), we demonstrated that HBs in both protein-ligand and protein-protein interfaces are overlooked in structures deposited in PDB. After the QM/MM optimization of 12 protein-ligand complexes, we showed that the overlooked HBs could be recovered. With a systematic search in the PDB, we found that the HB number per residue (N_HB/R) in proteins decreases as structural resolution becomes lower, implying that HBs are overlooked even today, regardless of the type of refinement approach used. Similarly, cation-π, π-π, and ionic interactions were found to be significantly lost, manifesting the universal underestimation of various NCIs. Considering the vital role of NCIs, it is important to recover the NCIs to facilitate drug design, to explore protein-protein interaction, and to study protein structure and function.

Computer-based techniques for lead identification and optimization II: Advanced search methods

This paper focuses on advanced computational techniques for identifying and optimizing lead molecules, such as metadynamics and a novel dynamic 3D pharmacophore analysis method called Dynophores. In this paper, the first application of the funnel metadynamics of the Berberine binding to G-quadruplex DNA is depicted, disclosing hints for drug design, in particular clarifying water’s role and suggesting the design of derivatives able to replace the solvent-mediated interactions between ligand and DNA to achieve more potent and selective activity. Secondly, the novel dynamic pharmacophore approach is an extension of the classic 3D pharmacophores, with statistical and sequential information about the conformational flexibility of a molecular system derived from molecular dynamics (MD) simulations.

Screening of Suitable Ionic Liquids as Green Solvents for Extraction of Eicosapentaenoic Acid (EPA) from Microalgae Biomass Using COSMO-RS Model

Omega-3 poly unsaturated fatty acids (PUFA) particularly eicosapentaenoic acid (EPA), and docosahexaenoic acid (DHA), have many health benefits including reducing the risk of cancer and cardiovascular disease. Recently, the use of ionic liquids (ILs) in lipid extraction from microalgae provides the potential to overcome common drawbacks, offers several other benefits. To date, very limited researches are available to focus on extracting microalgae lipid and PUFA in particular by using ILs. The objective of current work is to screen the potential ILs that can be applied in EPA extraction. In this study, fast ILs screening was performed with the help of a conductor like screening model for real solvents (COSMO-RS) and the ILs with higher capacity values for use in extraction of EPA were compared. According to the results, the highest capacity for EPA extraction among 352 screened cation/anion combinations belongs to [TMAm][SO₄]. It is expected to achieve a higher yield of EPA once applying this combination as the solvent in the process of extraction. ILs with small anions were observed to have higher capacities, as well possessing higher charge density compared to larger ones, and therefore, they are more preferable for extraction purposes. Moreover, shorter alkyl chain cations are preferred when using imidazolium-based IL, which agrees with experimental data.

Hierarchical lattice models of hydrogen-bond networks in water

We develop a graph-based model of the hydrogen-bond network in water, with a view toward quantitatively modeling the molecular-level correlational structure of the network. The networks formed are studied by the constructing the model on two infinite-dimensional lattices. Our models are built bottom up, based on microscopic information coming from atomistic simulations, and we show that the predictions of the model are consistent with known results from ab initio simulations of liquid water. We show that simple entropic models can predict the correlations and clustering of local-coordination defects around tetrahedral waters observed in the atomistic simulations. We also find that orientational correlations between bonds are longer ranged than density correlations, determine the directional correlations within closed loops, and show that the patterns of water wires within these structures are also consistent with previous atomistic simulations. Our models show the existence of density and compressibility anomalies, as seen in the real liquid, and the phase diagram of these models is consistent with the singularity-free scenario previously proposed by Sastry and coworkers [Phys. Rev. E 53, 6144 (1996)1063-651X10.1103/PhysRevE.53.6144].

The impact of water on the ambivalent behavior and paradoxical phenomenon of the amyloid-β fibril protein

The crucial role of water in amyloid-β(Aβ) fibril proteins is evaluated in several ways including the water's thermodynamic and kinetic solvation effects. As regards the water's character, its hindered-rotation barriers are also considered. The following protein molecules considered here are: the Aβ40 (PDB ID: 2LMN), Aβ42 (PDB ID: 5KK3 and 2NAO) and the double-layered Aβ17-42 fibril. We discuss: (i) extracellular Aβ40 and Aβ42 fibril monomers exhibit an ambivalent propensity to transform into a helical form toward the N-term region and a β-strand-like form near the C-terminal; (ii) interfacial water molecules play a crucial role in protein-protein interactions, as molecular dynamics simulations have shown a significant impact on the protein-protein binding; (iii) it is shown that the spontaneous dimerization process of the Aβ42 fibril protein in water occurs via a two-step nucleation-accommodation mechanism; (iv) MD simulations of the double-layered Aβ17-42 fibril model show that the C↔C interface appears more energetically favorable than the N↔N interface due to large hydrophobic contacts; (v) the water's role in the HET-s prion and in the Aβ fibrillar aggregates; (vi) it was found that the monomer-oligomer equilibrium spontaneously dissociates into stable monomeric species when they are incubated up to 3 μm for a longer time (>1 week) in a physiological buffer.

Direct Evidence of Intrinsic Blue Fluorescence from Oligomeric Interfaces of Human Serum Albumin

The molecular origin behind the concentration-dependent intrinsic blue fluorescence of human serum albumin (HSA) is not known yet. This unusual blue fluorescence is believed to be a characteristic feature of amyloid-like fibrils of protein/peptide and originates due to the delocalization of peptide bond electrons through the extended hydrogen bond networks of cross-β-sheet structure. Herein, by combining the results of spectroscopy, size exclusion chromatography, native gel electrophoresis, and confocal microscopy, we have shown that the intrinsic blue fluorescence of HSA exclusively originates from oligomeric interfaces devoid of any amyloid-like fibrillar structure. Our study suggests that this low energy fluorescence band is not due to any particular residue/sequence, but rather it is a common feature of self-assembled peptide bonds. The present findings of intrinsic blue fluorescence from oligomeric interfaces pave the way for future applications of this unique visual phenomenon for early stage detection of various protein aggregation related human diseases.