Structural and molecular insights into tacrine-benzofuran hybrid induced inhibition of amyloid-β peptide aggregation and BACE1 activity

Amyloid-β (Aβ) aggregation and β-amyloid precursor protein cleaving enzyme 1 (BACE1) are the potential therapeutic drug targets for Alzheimer's disease (AD). A recent study highlighted that tacrine-benzofuran hybrid C1 displayed anti-aggregation activity against Aβ42 peptide and inhibit BACE1 activity. However, the inhibition mechanism of C1 against Aβ42 aggregation and BACE1 activity remains unclear. Thus, molecular dynamics (MD) simulations of Aβ42 monomer and BACE1 with and without C1 were performed to inspect the inhibitory mechanism of C1 against Aβ42 aggregation and BACE1 activity. In addition, a ligand-based virtual screening followed by MD simulations was employed to explore potent new small-molecule dual inhibitors of Aβ42 aggregation and BACE1 activity. MD simulations highlighted that C1 promotes the non aggregating helical conformation in Aβ42 and destabilizes D23-K28 salt bridge that plays a vital role in the self-aggregation of Aβ42. C1 displays a favourable binding free energy (-50.7 ± 7.3 kcal/mol) with Aβ42 monomer and preferentially binds to the central hydrophobic core (CHC) residues. MD simulations highlighted that C1 strongly interacted with the BACE1 active site (Asp32 and Asp228) and active pockets. The scrutiny of interatomic distances among key residues of BACE1 highlighted the close flap (non-active) position in BACE1 on the incorporation of C1. The MD simulations explain the observed high inhibitory activity of C1 against Aβ aggregation and BACE1 in the in vitro studies. The ligand-based virtual screening followed by MD simulations identified CHEMBL2019027 (C2) as a promising dual inhibitor of Aβ42 aggregation and BACE1 activity.Communicated by Ramaswamy H. Sarma.

Conformational dynamics of superoxide dismutase (SOD1) in osmolytes: a molecular dynamics simulation study

Amyotrophic lateral sclerosis (ALS) is a progressive neurodegenerative disease caused by the misfolding of Cu, Zn superoxide dismutase (SOD1). Several earlier studies have shown that monomeric apo SOD1 undergoes significant local unfolding dynamics and is the predecessor for aggregation. Here, we have employed atomistic molecular dynamics (MD) simulations to study the structure and dynamics of monomeric apo and holo SOD1 in water, aqueous urea and aqueous urea-TMAO (trimethylamine oxide) solutions. Loop IV (zinc-binding loop) and loop VII (electrostatic loop) of holo SOD1 are considered as functionally important loops as they are responsible for the structural stability of holo SOD1. We found larger local unfolding of loop IV and VII of apo SOD1 as compared to holo SOD1 in water. Urea induced more unfolding in holo SOD1 than apo SOD1, whereas the stabilization of both the form of SOD1 was observed in ternary solution (i.e. water/urea/TMAO solution) but the extent of stabilization was higher in holo SOD1 than apo SOD1. The partially unfolded structures of apo SOD1 in water, urea and holo SOD1 in urea were identified by the exposure of the hydrophobic cores, which are highly dynamic and these may be the initial events of aggregation in SOD1. Our simulation studies support the formation of aggregates by means of the local unfolding of monomeric apo SOD1 as compared to holo SOD1 in water.

Discovery and Exploration of the Efficient Acyclic Dehydration of Hexoses in Dimethyl Sulfoxide/Water

Current gaps in the development of sustainable processes include a lack of strategies to systematically identify and optimize the formation of new products. The dehydration of hexoses to 5-hydroxymethylfurfural (HMF) is a particularly widely studied process. In an attempt to identify a new high-selectivity conversion of glucose, quantitative NMR spectroscopy is used to screen conditions that have been reported to yield high conversions of glucose but low formation of HMF. In this manner, an olefinic six-carbon byproduct is identified. By adding water, selectivity for the compound was nearly tripled relative to previous reports. The detection of high-yielding side reactions in the formation of HMF is remarkable, considering how extensively HMF formation has been studied. High selectivity for the acyclic pathway allows hitherto unobserved intermediates in this pathway to be identified by using in situ NMR spectroscopy. An additional, presumably cyclic, pathway contributes to HMF formation.

Molecular insights into the inhibitory mechanism of bi-functional bis-tryptoline triazole against β-secretase (BACE1) enzyme

The β-site amyloid precursor protein-cleaving enzyme 1 (β-secretase, BACE1) is involved in the formation of amyloid-β (Aβ) peptide that aggregates into soluble oligomers, amyloid fibrils, and plaques responsible for the neurodegeneration in Alzheimer disease (AD). BACE1 is one of the prime therapeutic targets for the design of inhibitors against AD as BACE1 participate in the rate-limiting step in Aβ production. Jiaranaikulwanitch et al. reported bis-tryptoline triazole (BTT) compound as a potent inhibitor against BACE1, Aβ aggregation as well as possessing metal chelation and antioxidant activity. However, the molecular mechanism of BACE1 inhibition by BTT remains unclear. Thus, molecular docking and molecular dynamics (MD) simulations were performed to elucidate the inhibitory mechanism of BTT against BACE1. MD simulations highlight that BTT interact with catalytic aspartic dyad residues (Asp32 and Asp228) and active pocket residues of BACE1. The hydrogen-bond interactions, hydrophobic contacts, and π-π stacking interactions of BTT with flap residues (Val67-Asp77) of BACE1 confine the movement of the flap and help to achieve closed (non-active) conformation. The PCA analysis highlights lower conformational fluctuations for BACE1-BTT complex, which suggests enhanced conformational stability in comparison to apo-BACE1. The results of the present study provide key insights into the underlying inhibitory mechanism of BTT against BACE1 and will be helpful for the rational design of novel inhibitors with enhanced potency against BACE1.

Effect of Osmolytes on Conformational Behavior of Intrinsically Disordered Protein α-Synuclein

α-Synuclein is an intrinsically disordered protein whose function in a healthy brain is poorly understood. It is genetically and neuropathologically linked to Parkinson's disease (PD). PD is manifested after the accumulation of plaques of α-synuclein aggregates in the brain cells. Aggregates of α-synuclein are very toxic and lead to the disruption of cellular homeostasis and neuronal death. α-Synuclein can also contribute to disease propagation as it may exert noxious effects on neighboring cells. Understanding the mechanism of α-synuclein aggregation will facilitate the problem of dealing with neurodegenerative diseases in general and that of PD in particular. Here, we have used molecular dynamics simulations to investigate the behavior of α-synuclein at various temperatures and in different concentrations of urea and trimethyl amine oxide. The residue region from 61 to 95 of α-synuclein is experimentally known as amyloidogenic. In our study, we have identified some other regions, which also have the propensity to form an aggregate besides this known sequence. Urea being a denaturant interacts more with these regions of α-synuclein through hydrogen bond formation and inhibits the β-sheet formation, whereas trimethyl amine oxide itself does not interact much with the protein and stabilizes the protein by preferentially distributing water molecules on the surface of the protein.

Pyranose ring puckering in aldopentoses, ketohexoses and deoxyaldohexoses. A molecular dynamics study

Conformation of monosaccharides, including the ring shape, has for years been the subject of intensive research. Although d-aldohexopyranoses are the most extensively studied pyranoses, there also exist other groups of saccharides that contain analogous chemical system of the six-membered ring. Here we describe in details the results of the molecular dynamics-based conformational analysis concerning a series of pyranoses, namely: d-aldopentoses, d-ketohexoses as well as deoxy- (d-quinovose, l-fucose, l-rhamnose) and dideoxy- (abequose, paratose, tyvelose, digitoxose) derivatives of aldohexoses. By using the carbohydrate-dedicated GROMOS 56a6_CARBO force field, we determined the conformational properties of both the lactol and hydroxymethyl groups as well as the anomeric populations for all considered compounds. The orientation of the lactol group follows the trend expected on the basis of the exo-anomeric effect for all compounds whereas the conformation of the hydroxymethyl group in d-ketohexoses is represented by the two gauche (with respect to the ring oxygen atom) rotamers. The special emphasis is put on the ring-inversion properties studied in the context of both the full chair-chair inversion and the chair-boat/skew-boat rearrangement. The calculated ring-distortion energies, compared with those obtained for regular d-aldohexopyranoses allowed for estimating the influence of particular substituents on the ring flexibility. Overall, such influence is correlated with the dimension of the substituent and its orientation but is limited to the case of the chair-chair inversion whereas the chair-to-boat/skew-boat rearrangement exhibits roughly the same properties for all pyranoses. For all d-aldopyranoses the α anomers exhibit lower ring-inversion free energies in comparison to the β anomers whereas this trend is inverted in the case of d-ketohexopyranoses.

Scrutiny of the mechanism of small molecule inhibitor preventing conformational transition of amyloid-β42 monomer: insights from molecular dynamics simulations

Alzheimer's disease (AD) is a progressive neurodegenerative disorder that is characterized by loss of intellectual functioning of brain and memory loss. According to amyloid cascade hypothesis, aggregation of amyloid-β₄₂ (Aβ₄₂) peptide can generate toxic oligomers and their accumulation in the brain is responsible for the onset of AD. In spite of carrying out a large number of experimental studies on inhibition of Aβ₄₂ aggregation by small molecules, the detailed inhibitory mechanism remains elusive. In the present study, comparable molecular dynamics (MD) simulations were performed to elucidate the inhibitory mechanism of a sulfonamide inhibitor C1 (2,5-dichloro-N-(4-piperidinophenyl)-3-thiophenesulfonamide), reported for its in vitro and in vivo anti-aggregation activity against Aβ₄₂. MD simulations reveal that C1 stabilizes native α-helix conformation of Aβ₄₂ by interacting with key residues in the central helix region (13-26) with hydrogen bonds and π-π interactions. C1 lowers the solvent-accessible surface area of the central hydrophobic core (CHC), KLVFF (16-20), that confirms burial of hydrophobic residues leading to the dominance of helical conformation in the CHC region. The binding free energy analysis with MM-PBSA demonstrates that Ala2, Phe4, Tyr10, Gln15, Lys16, Leu17, Val18, Phe19, Phe20, Glu22, and Met35 contribute maximum to binding free energy (-43.1 kcal/mol) between C1 and Aβ₄₂ monomer. Overall, MD simulations reveal that C1 inhibits Aβ₄₂ aggregation by stabilizing native helical conformation and inhibiting the formation of aggregation-prone β-sheet conformation. The present results will shed light on the underlying inhibitory mechanism of small molecules that show potential in vitro anti-aggregation activity against Aβ₄₂.