Binding of 12-Crown-4 with Alzheimer’s Aβ40 and Aβ42 Monomers and Its Effect on Their Conformation: Insight from Molecular Dynamics Simulations

Examine stability polyvinyl alcohol-stabilized nanosuspensions to overcome the challenge of poor drug solubility utilizing molecular dynamic simulation

The pharmaceutical industry faces a significant challenge from the low water solubility of nearly 90% of newly developed Active Pharmaceutical Ingredients (APIs). Despite extensive efforts to improve solubility, approximately 40% of these APIs encounter commercialization hurdles, impacting drug efficacy. In this context, a promising strategy will be introduced in which nanosuspensions, particularly polyvinyl alcohol (PVA) as a stabilizer, are applied to increase drug solubility. In this work using molecular dynamics simulations, the nanosuspension of four poorly water-soluble drugs (flurbiprofen, bezafibrate, miconazole, and phenytoin) stabilized with PVA is investigated. The simulation data showed van der Waals energies between polyvinyl alcohol with flurbiprofen and bezafibrate are - 101.12 and - 58.42 kJ/mol, respectively. The results indicate that PVA is an effective stabilizer for these drugs, and superior interactions are obtained with flurbiprofen and bezafibrate. The study also explores the impact of PVA on water molecule diffusion, providing insights into the stability of nanosuspensions. Obtained results also provide valuable insights into hydrogen bond formation, diffusion coefficients, and nanosuspension stability, contributing to the rational design and optimization of pharmaceutical formulations.

Elucidating bis-pyrimidines as new and efficient mushroom tyrosinase inhibitors: synthesis, SAR, kinetics and computational studies

In this study, a series of novel bis-pyrimidine derivatives (1P-8P) were designed, synthesized, characterized, and investigated for their in vitro inhibitory activity against mushroom tyrosinase, an enzyme critical in melanin biosynthesis and implicated in various hyperpigmentation disorders. To the best of our knowledge, the bispyrimidine scaffold has been evaluated for the first time for its tyrosinase inhibitory activity. Their inhibitory activities were assessed, revealing inhibition with IC50 values in the micromolar range. Additionally, this series of compounds were found to inhibit tyrosinase activity in a mixed-type manner, with IC50 values ranging from 12.36 ± 1.24 to 86.67 ± 3.08 μM. To further elucidate the binding interactions, molecular docking simulations were performed, identifying key residues in the active site responsible for binding affinity. Furthermore, molecular dynamics (MD) simulations were conducted to assess the dynamic behavior, stability, and binding affinity of the most potent inhibitor, compound 6P. Quantitative Structure-Activity Relationship (QSAR) models were developed to correlate the structural features of the bis-pyrimidines with their inhibitory activity, providing insights into the structure-activity relationships (SAR) that govern their potency. The experimental and theoretical findings demonstrated excellent agreement. These findings pave the way for the development of novel bis-pyrimidine-based therapeutic agents for treating hyperpigmentation and related conditions.

Investigating the Effects of the POPC-POPG Lipid Bilayer Composition on PAP248-286 Binding Using CG Molecular Dynamics Simulations

PAP248-286 is a fusogenic peptide derived from prostatic acid phosphatase, commonly found in human semen, and is known to mediate HIV fusion with cell membranes. In this study, we performed 120 independent coarse-grained molecular dynamics simulations to investigate the spontaneous binding of PAP248-286 monomers, considering both charged and neutral histidine (His) residues, to membrane bilayers composed of different lipid compositions: 100% POPC, 70% POPC-30% POPG, and 50% POPC-50% POPG. Our simulations revealed that PAP248-286 displayed spontaneous binding to the membrane, with increased binding observed in the presence of anionic lipid POPG. Specifically, in systems containing 30% and 50% POPG lipids, monomer residues, particularly in the systems containing charged histidine (His) residues, exhibited prolonged binding with the membrane. Furthermore, our simulations indicated that PAP248-286 adopted a parallel orientation with the membrane, exposing its positively charged residues to the lipid bilayer. Interestingly, systems containing charged His residues showed a higher lipid occupancy around the peptide. These findings are consistent with previous experimental data, suggesting that PAP248-286 binding is enhanced in membranes with charged His residues, resembling the conditions found in the acidic vaginal pH environment. The results of our study provide further insights into the molecular mechanisms underlying the membrane binding of PAP248-286, contributing to our understanding of its potential role in HIV fusion and infection.

Crown Ether-Capped Gold Nanoclusters as a Multimodal Platform for Bioimaging

The distinct polarity of biomolecule surfaces plays a pivotal role in their biochemistry and functions as it is involved in numerous processes, such as folding, aggregation, or denaturation. Therefore, there is a need to image both hydrophilic and hydrophobic bio-interfaces with markers of distinct responses to hydrophobic and hydrophilic environments. In this work, we present a synthesis, characterization, and application of ultrasmall gold nanoclusters capped with a 12-crown-4 ligand. The nanoclusters present an amphiphilic character and can be successfully transferred between aqueous and organic solvents and have their physicochemical integrity retained. They can serve as probes for multimodal bioimaging with light (as they emit near-infrared luminescence) and electron microscopy (due to the high electron density of gold). In this work, we used protein superstructures, namely, amyloid spherulites, as a hydrophobic surface model and individual amyloid fibrils with a mixed hydrophobicity profile. Our nanoclusters spontaneously stained densely packed amyloid spherulites as observed under fluorescence microscopy, which is limited for hydrophilic markers. Moreover, our clusters revealed structural features of individual amyloid fibrils at a nanoscale as observed under a transmission electron microscope. We show the potential of crown ether-capped gold nanoclusters in multimodal structural characterization of bio-interfaces where the amphiphilic character of the supramolecular ligand is required.

A coarse-grained molecular dynamics investigation on spontaneous binding of Aβ1-40 fibrils with cholesterol-mixed DPPC bilayers

Alzheimer's disease is the most common form of dementia. Its aetiology is characterized by the misfolding and aggregation of amyloid-β (Aβ) peptides into β-sheet-rich Aβ oligomers/fibrils. Although multiple experimental studies have suggested that Aβ oligomers/fibrils interact with the cell membranes and perturb their structures and dynamics, the molecular mechanism of this interaction is still not fully understood. In the present work, we have performed a total of 120 μs-long simulations to investigate the interaction between trimeric or hexameric Aβ_1-40 fibrils with either a 100% DPPC bilayer, a 70% DPPC-30% cholesterol bilayer or a 50% DPPC-50% cholesterol bilayer. Our simulation data capture the spontaneous binding of the aqueous Aβ_1-40 fibrils with the membranes and show that the central hydrophobic amino acid cluster, the lysine residue adjacent to it and the C-terminal hydrophobic residues are all involved in the process. Moreover, our data show that while the Aβ_1-40 fibril does not bind to the 100% DPPC bilayer, its binding affinity for the membrane increases with the amount of cholesterol. Overall, our data suggest that two clusters of hydrophobic residues and one lysine help Aβ_1-40 fibrils establish stable interactions with a cholesterol-rich DPPC bilayer. These residues are likely to represent potential target regions for the design of inhibitors, thus opening new avenues in structure-based drug design against Aβ oligomer/fibril-membrane interaction.

Early Stages of Misfolding of PAP248-286 at two different pH values: An Insight from Molecular Dynamics Simulations

PAP248-286 peptides, which are highly abundant in human semen, aggregate and form amyloid fibrils that enhance HIV infection. Previous experimental studies have shown that the infection-promoting activity of PAP248-286 begins to increase well before amyloid formation takes place and that pH plays a key role in the enhancement of PAP248-286-related infection. Hence, understanding the early stages of misfolding of the PAP2482-86 peptide is crucial. To this end, we have performed 60 independent MD simulations for a total of 24 µs at two different pH values (4.2 and 7.2). Our data shows that early stages of misfolding of the PAP248-286 peptide is a multistage process and that the first step of the process is a transition from an “I-shaped” structure to a “U-shaped” structure. We further observed that the structure of PAP248-286 at the two different pH values shows significantly different features. At pH 4.2, the peptide has less intra-molecular H-bonds and a reduced α-helical content than at pH 7.2. Moreover, differences in intra-peptide residues contacts are also observed at the two pH values. Finally, free energy landscape analysis shows that there are more local minima in the energy surface of the peptide at pH 7.2 than at pH 4.2. Overall, the present study elucidates the early stages of misfolding of the PAP248-286 peptide at the atomic level, thus possibly opening new avenues in structure-based drug discovery against HIV infection.

Molecular Insights into the Inhibitory Effect of GV971 Components Derived from Marine Acidic Oligosaccharides against the Conformational Transition of Aβ42 Monomers

GV971 derived from marine acidic oligosaccharides has been used to cure Alzheimer's disease (AD). However, the molecular mechanism of its inhibition of the conformational transition of amyloid β-proteins (Aβ) is still unclear. Herein, molecular dynamics simulations were used to explore the molecular mechanism of the main GV971 components including DiM, TetraM, HexaM, and OctaM to inhibit the conformational conversion of the Aβ42 monomer. It is found that the GV971 components inhibit the conformational transition from α-helix to β-sheet and the hydrophobic collapse of the Aβ42 monomer. In addition, the binding energy analysis implies that both electrostatic and van der Waals interactions are beneficial to the binding of GV971 components to the Aβ42 monomer. Among them, electrostatic interactions occupy the dominant position. Moreover, the GV971 components mainly interact directly with the charged residues D1, R5, K16, and K28 by forming salt bridges and hydrogen bonds, which specifically bind to the N-terminal region of Aβ42.

An α-helix mimetic oligopyridylamide, ADH-31, modulates Aβ42 monomer aggregation and destabilizes protofibril structures: insights from molecular dynamics simulations

Alzheimer's disease (AD), an epidemic growing worldwide due to no effective medical aid available in the market, is a neurological disorder. AD is known to be directly associated with the toxicity of amyloid-β (Aβ) aggregates. In search of potent inhibitors of Aβ aggregation, Hamilton and co-workers reported an α-helix mimetic, ADH-31, which acts as a powerful antagonist of Aβ42 aggregation. To identify the key interactions between protein-ligand complexes and to gain insights into the inhibitory mechanism of ADH-31 against Aβ42 aggregation, molecular dynamics (MD) simulations were performed in the present study. The MD simulations highlighted that ADH-31 showed distinct binding capabilities with residues spanning from the N-terminal to the central hydrophobic core (CHC) region of Aβ42 and restricted the conformational transition of the helix-rich structure of Aβ42 into another form of secondary structures (coil/turn/β-sheet). Hydrophobic contacts, hydrogen bonding and π-π interaction contribute to the strong binding between ADH-31 and Aβ42 monomer. The Dictionary of Secondary Structure of Proteins (DSSP) analysis highlighted that the probability of helical content increases from 38.5% to 50.2% and the turn content reduces from 14.7% to 6.2% with almost complete loss of the β-sheet structure (4.5% to 0%) in the Aβ42 monomer + ADH-31 complex. The per-residue binding free energy analysis demonstrated that Arg5, Tyr10, His14, Gln15, Lys16, Val18, Phe19 and Lys28 residues of Aβ42 are responsible for the favourable binding free energy in Aβ42 monomer + ADH-31 complex, which is consistent with the 2D HSQC NMR of the Aβ42 monomer that depicted a change in the chemical shift of residues spanning from Glu11 to Phe20 in the presence of ADH-31. The MD simulations highlighted the prevention of sampling of amyloidogenic β-strand conformations in Aβ42 trimer in the presence of ADH-31 as well as the ability of ADH-31 to destabilize Aβ42 trimer and protofibril structures. The lower binding affinity between Aβ42 trimer chains in the presence of ADH-31 highlights the destabilization of the Aβ42 trimer structure. Overall, MD results highlighted that ADH-31 inhibited Aβ42 aggregation by constraining Aβ peptides into helical conformation and destabilized Aβ42 trimer as well as protofibril structures. The present study provides a theoretical insight into the atomic level details of the inhibitory mechanism of ADH-31 against Aβ42 aggregation as well as protofibril destabilization and could be implemented in the structure-based drug design of potent therapeutic agents for AD.

How Does the Mono-Triazole Derivative Modulate Aβ42 Aggregation and Disrupt a Protofibril Structure: Insights from Molecular Dynamics Simulations

Clinical studies have identified that abnormal self-assembly of amyloid-β (Aβ) peptide into toxic fibrillar aggregates is associated with the pathology of Alzheimer's disease (AD). The most acceptable therapeutic approach to stop the progression of AD is to inhibit the formation of β-sheet-rich structures. Recently, we designed and evaluated a series of novel mono-triazole derivatives 4(a-x), where compound 4v was identified as the most potent inhibitor of Aβ₄₂ aggregation and disaggregates preformed Aβ₄₂ fibrils significantly. Moreover, 4v strongly averts the Cu²⁺-induced Aβ₄₂ aggregation and disaggregates the preformed Cu²⁺-induced Aβ₄₂ fibrils, halts the generation of reactive oxygen species, and shows neuroprotective effects in SH-SY5Y cells. However, the underlying molecular mechanism of inhibition of Aβ₄₂ aggregation by 4v and disaggregation of preformed Aβ₄₂ fibrils remains obscure. In this work, molecular dynamics (MD) simulations have been performed to explore the conformational ensemble of the Aβ₄₂ monomer and a pentameric protofibril structure of Aβ₄₂ in the presence of 4v. The MD simulations highlighted that 4v binds preferentially at the central hydrophobic core region of the Aβ₄₂ monomer and chains D and E of the Aβ₄₂ protofibril. The dictionary of secondary structure of proteins analysis indicated that 4v retards the conformational conversion of the helix-rich structure of the Aβ₄₂ monomer into the aggregation-prone β-sheet conformation. The binding free energy calculated by the molecular mechanics Poisson-Boltzmann surface area method revealed an energetically favorable process with ΔG _binding = -44.9 ± 3.3 kcal/mol for the Aβ₄₂ monomer-4v complex. The free energy landscape analysis highlighted that the Aβ₄₂ monomer-4v complex sampled conformations with significantly higher helical contents (35 and 49%) as compared to the Aβ₄₂ monomer alone (17%). Compound 4v displayed hydrogen bonding with Gly37 (chain E) and π-π interactions with Phe19 (chain D) of the Aβ₄₂ protofibril. Further, the per-residue binding free energy analysis also highlighted that Phe19 (chain D) and Gly37 (chain E) of the Aβ₄₂ protofibril showed the maximum contribution in the binding free energy. The decreased binding affinity and residue-residue contacts between chains D and E of the Aβ₄₂ protofibril in the presence of 4v indicate destabilization of the Aβ₄₂ protofibril structure. Overall, the structural information obtained through MD simulations indicated that 4v stabilizes the native helical conformation of the Aβ₄₂ monomer and persuades a destabilization in the protofibril structure of Aβ₄₂. The results of the study will be useful in the rational design of potent inhibitors against amyloid aggregation.

Effect of Piedmont mutation (L34V) on the structure, dynamics, and aggregation of Alzheimer's Aβ40 peptide

The amyloid-β (Aβ) aggregation in the brain has been associated with the development of Alzheimer's disease (AD). The previous studies have reported that Piedmont mutation (L34V) increases the rate of Aβ₄₀ aggregation. However, the underlying molecular mechanism of the effect of L34V mutation on Aβ₄₀ structure, dynamics, and aggregation remains largely unclear. In the present study, molecular dynamics (MD) simulations were performed to elucidate the effect of L34V mutation on the structural changes and conformational dynamics of Aβ₄₀. The secondary structure analysis highlight that L34V mutation enhances Aβ₄₀ self-assembly due to the formation of aggregation-prone β-sheet structure at the C-terminus of Aβ₄₀ monomeric structure. The higher probability of Asp23-Lys28 salt bridge interaction in Aβ₄₀(L34V) leads to aggregation prone β-sheet conformations, which has the potential to increase the fibril formation rate. The free energy landscape (FEL) analysis depict a sampling of coil conformation in the free energy minima of Aβ₄₀, whereas the aggregation-prone β-sheet conformation was observed at the C-terminal region of Aβ₄₀(L34V) in the minimum energy conformations extracted from FEL of Aβ₄₀(L34V). MD simulations, in agreement with experiment, highlight that L34V mutation increases Aβ₄₀ aggregation as the sampling of the aggregation-prone β-sheet conformation substantially increased. Overall, MD simulations provided atomic level details into the increased fibril formation tendency upon L34V mutation and physical insights into the L34V-mediated conformational as well as structural changes in Aβ₄₀.

Interactions of a multifunctional di-triazole derivative with Alzheimer's Aβ42monomer and Aβ42protofibril: a systematic molecular dynamics study

The molecular dynamics simulations results highlighted that the multi-target-directed ligand6nstabilizes the native α-helix conformation of the Aβ₄₂monomer and induces a sizable destabilization in the Aβ₄₂protofibril structure.

Effects of natural compounds on conformational properties and hairpin formation of amyloid-β42 monomer: docking and molecular dynamics simulation study

The β₄₂ amyloid peptides (Aβ) are identified as a candidate target for Alzheimer's drugs. Phenolic compounds can bind to the Aβ and inhibit amyloid formation. However, the inhibitory mechanism of phenolic compounds remains unclear. In this study, the molecular dynamic simulation and docking program were used to characterize the molecular details of inhibitory mechanism of the phenolic compounds. Our Results show that the phenolic compounds can bind to hydrophobic region in Aβ₄₂ monomer and alter hydrophobic interactions network at Aβ₄₂ which play a key role in β-sheet formation. The cluster analysis and interactions network analysis were used to probe conformational changes in Aβ₄₂. In most populated clusters of Aβ₄₂-phenolic compounds complexes, the sheet structures were not observed or reduced. It seems that the binding of phenolic compounds can induce larger conformational diversity for amyloid peptide and changes conformational properties of amyloid peptide. The phenolic compounds can deform β-Hairpin structure of Aβ by destabilizing salt bridges E22-K28 and D23-K28 which can alter the conformation of Aβ₄₂ in aqueous solution. These findings are in accordance with experimental results, to some extent give a molecular level interpretation for the inhibitory mechanism of phenolic compounds .In addition, this study may add important new details to the inhibitory mechanism of Alzheimer's drug.Communicated by Ramaswamy H. Sarma.

Molecular Simulations Reveal Terminal Group Mediated Stabilization of Helical Conformers in Both Amyloid-β42 and α-Synuclein

The presence of partially structured helices in natively unfolded amyloid-β42 (Aβ42) and α-synuclein (αS) has been shown to accelerate fibrillation in the onset of Alzheimer's and Parkinson's disease, respectively. At the other extreme, folded stable helical conformers have also been reported to resist amyloid formation. Recent studies indicate that amyloidogenic aggregation can be impeded using small molecules that stabilize the α-helical monomers and switch off the neurotoxic pathway. We predict a common intrapeptide route to stabilization based on the plasticity of helical conformations of Aβ42 and αS as assessed through extensive atomistic molecular dynamics (MD) computer simulations (∼36 μs) across ten distinct protein force field and water model combinations. Computed free energies and interaction maps (not obtainable from experiments alone) show that flexible terminal groups (N-terminus of Aβ42 and C-terminus of αS) show a tendency to stabilize folded helical conformations in both peptides via primary hydrophobic interactions with central hydrophobic domains, and secondary salt bridges with other domains. These interactions confer aggregation resistance by decreasing the population of partially structured helices and are absent in control simulations of complete unfolding. Computed helical stability is also significantly reduced in terminal-deleted variants. The models suggest new strategies to tackle neurodegeneration by rationally re-engineering terminal groups to optimize their predicted ability to deactivate helical monomers.

Antimicrobial cell penetrating peptides with bacterial cell specificity: pharmacophore modelling, quantitative structure activity relationship and molecular dynamics simulation

Current research has shown cell-penetrating peptides and antimicrobial peptides (AMPs) as probable vectors for use in drug delivery and as novel antibiotics. It has been reported that the higher the therapeutic index (TI) the higher would be the bacterial cell penetrating ability. To the best of our knowledge, no in-silico study has been performed to determine bacterial cell specificity of the antimicrobial cell penetrating peptides (aCPP's) based on their TI. The aim of this study was to develop a quantitative structure activity relationship (QSAR) model, which can estimate antimicrobial potential and cell-penetrating ability of aCPPs against S. aureus, to confirm the relationship between the TI and aCPPs and to identify specific descriptors responsible for aCPPs penetrating ability. Molecular dynamics (MD) simulation was also performed to confirm the membrane insertion of the most active aCPPs obtained from the QSAR study. The most appropriate pharmacophore was identified to predict the aCPP's activity. The statistical results confirmed the validity of the model. The QSAR model was successful in identifying the optimal aCPP with high activity prediction and provided insights into the structural requirements to correlate their TI to cell penetrating ability. MD simulation of the best aCPP with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer confirmed its interaction with the membrane and the C-terminal residues of the aCPP played a key role in membrane penetration. The strategy of combining QSAR and molecular dynamics, allowed for optimal estimation of ligand-target interaction and confirmed the importance of Trp and Lys in interacting with the POPC bilayer. Communicated by Ramaswamy H. Sarma.

Unraveling the mechanism of l-gulonate-3-dehydrogenase inhibition by ascorbic acid: Insights from molecular modeling

l-Gulonate dehydrogenase (GuDH) is a crucial enzyme in the non-phosphorylated sugar metabolism or glucuronate-xylulose (GX) pathway. Some naturally occurring compounds inhibit GuDH. Ascorbic acid is one of such inhibitors for GuDH. However, the exact mechanism by which ascorbic acid inhibits GuDH is still unknown. In this study, we try to investigate GuDH inhibition using computational approaches by generating a model for buffalo GuDH. We used this model to perform blind dockings of ascorbic acid to GuDH. Some docked conformations of ascorbic acid bind near Asp39 and have steric clashes with crystal structure conformation of NADH. To assess the dynamic stability of the GuDH-ascorbic acid complex, we performed six molecular dynamics simulations for GuDH, three each in its free form and in complex with ascorbic acid for 50 ns, to obtain 300 ns of trajectories in total. During the simulations, ascorbic acid interacted with several residues nearby Asp39. As Asp39 is an important residue for NADH binding and specificity, the interaction of ascorbic acid near Asp39 hinders further NADH binding and ultimately affects the enzymatic functioning of GuDH. In this study, we analyze these interactions between ascorbic acid and GuDH. Our analysis reveals novel details on the mechanism of GuDH inhibition by ascorbic acid.

Formulation and Molecular Dynamics Simulations of a Fusidic Acid Nanosuspension for Simultaneously Enhancing Solubility and Antibacterial Activity

The aim of the present study was to formulate a nanosuspension (FA-NS) of fusidic acid (FA) to enhance its aqueous solubility and antibacterial activity. The nanosuspension was characterized using various in vitro, in silico, and in vivo techniques. The size, polydispersity index, and zeta potential of the optimized FA-NS were 265 ± 2.25 nm, 0.158 ± 0.026, and -16.9 ± 0.794 mV, respectively. The molecular dynamics simulation of FA and Poloxamer-188 showed an interaction and binding energy of -74.42 kJ/mol and -49.764 ± 1.298 kJ/mol, respectively, with van der Waals interactions playing a major role in the spontaneous binding. There was an 8-fold increase in the solubility of FA in a nanosuspension compared to the bare drug. The MTT assays showed a cell viability of 75-100% confirming the nontoxic nature of FA-NS. In vitro antibacterial activity revealed a 16- and 18-fold enhanced activity against Staphylococcus aureus (SA) and methicillin-resistant SA (MRSA), respectively, when compared to bare FA. Flowcytometry showed that MRSA cells treated with FA-NS had almost twice the percentage of dead bacteria in the population, despite having an 8-fold lower MIC in comparison to the bare drug. The in vivo skin-infected mice showed a 76-fold reduction in the MRSA load for the FA-NS treated group compared to that of the bare FA. These results show that the nanosuspension of antibiotics can enhance their solubility and antibacterial activity simultaneously.