Nonequilibrium all-atom molecular dynamics simulation of the bubble cavitation and application to dissociate amyloid fibrils

Open-boundary molecular dynamics of ultrasound using supramolecular water models

Ultrasound can be used to manipulate protein function and activity, as well as for targeted drug delivery, making it a powerful diagnostic and therapeutic modality with wide applications in sonochemistry, nanotechnology, and engineering. However, a general particle-based approach to ultrasound modeling remains challenging due to the significant disparity between characteristic time scales governing ultrasound propagation. In this study, we use open-boundary molecular dynamics to simulate ultrasound waves in liquid water under ambient conditions by employing supramolecular water models, i.e., the Martini 3, dissipative particle dynamics, and many-body dissipative particle dynamics models. We demonstrate that our approach successfully reproduces the solution of the traveling wave equation and captures the velocity dispersion characteristic of high-frequency ultrasound waves.

Nonequilibrium Molecular Dynamics Method to Generate Poiseuille-Like Flow between Lipid Bilayers

There are various flows inside and outside cells in vivo. Nonequilibrium molecular dynamics (NEMD) simulation is a useful tool for understanding the effects of these flows on the dynamics of biomolecules. We propose an NEMD method to generate a Poiseuille-like flow between lipid bilayers. We extended the conventional equilibrium MD method to produce a flow by adding constant external force terms to the water molecules. Using the Lagrange multiplier method, the center of mass of the lipid bilayer is constrained so that the flow does not sweep away the lipid bilayer, but the individual lipid molecules fluctuate. The temperature of the system is controlled properly in the solution and membrane by using the Nosé-Hoover thermostat. We found that the flow velocity increases linearly as the applied external force term increases. It is possible to estimate the appropriate value of acceleration to generate a flow with an arbitrary velocity using this proportional relation once a single short MD simulation is performed. We also found that the flow between two lipid bilayers is slower than the analytical solution of the Navier-Stokes equations between rigid parallel plates due to the interactions between water molecules and the membrane. This method can be applied not only to a flow on lipid membranes but also to a flow on soft surfaces generally.

Comprehending the Efficacy of Whitlock's Caffeine-Pincered Molecular Tweezer on β-Amyloid Aggregation

Alzheimer's disease (AD) stands as one of the most prevalent neurodegenerative conditions, leading to cognitive impairment, with no cure and preventive measures. Misfolding and aberrant aggregation of amyloid-β (Aβ) peptides are believed to be the underlying cause of AD. These amyloid aggregates culminate in the development of toxic Aβ oligomers and subsequent accumulation of β-amyloid plaques amidst neuronal cells in the brain, marking the hallmarks of AD. Drug development for the potentially curative treatment of Alzheimer's is, therefore, a tremendous challenge for the scientific community. In this study, we investigate the potency of Whitlock's caffeine-armed molecular tweezer in combating the deleterious effects of Aβ aggregation, with special emphasis on the seven residue Aβ16-22 fragment. Extensive all-atom molecular dynamics simulations are conducted to probe the various structural and conformational transitions of the peptides in an aqueous medium in both the presence and absence of tweezers. To explore the specifics of peptide-tweezer interactions, radial distribution functions, contact number calculations, binding free energies, and 2-D kernel density plots depicting the variation of distance-angle between the aromatic planes of the peptide-tweezer pair are computed. The central hydrophobic core, particularly the aromatic Phe residues, is crucial in the development of harmful amyloid oligomers. Notably, all analyses indicate reduced interpeptide interactions in the presence of the tweezer, which is attributed to the tweezer-Phe aromatic interaction. Upon increasing the tweezer concentration, the residues of the peptide are further encased in a hydrophobic environment created by the self-aggregating tweezer cluster, leading to the segregation of the peptide residues. This is further aided by the weakening of interstrand hydrogen bonding between the peptides, thereby impeding their self-aggregation and preventing the formation of neurotoxic β-amyloid. Furthermore, the study also highlights the efficacy of the molecular tweezer in destabilizing preformed amyloid fibrils as well as hindering the aggregation of the full-length Aβ1-42 peptide.

Can local heating and molecular crowders disintegrate amyloid aggregates?

The present study employs a blend of molecular dynamics simulations and a theoretical model to explore the potential disintegration mechanism of a matured Aβ octamer, aiming to offer a strategy to combat Alzheimer's disease. We investigate local heating and crowding effects on Aβ disintegration by selectively heating key Aβ segments and varying the concentration of sodium dodecyl sulphate (SDS), respectively. Despite initiation of disruption, Aβ aggregates resist complete disintegration during local heating due to rapid thermal energy distribution to the surrounding water. Conversely, although SDS molecules effectively inhibit Aβ aggregation at higher concentration through micelle formation, they fail to completely disintegrate the aggregate due to the exceedingly high energy barrier. To address the sampling challenge posed by the formidable energy barrier, we have performed well-tempered metadynamics simulations. Simulations reveal a multi-step disintegration mechanism for the Aβ octamer, suggesting a probable sequence: octamer → pentamer/hexamer ⇌ tetramer → monomer, with a rate-determining step constituting 45 kJ mol-1 barrier during the octamer to pentamer/hexamer transition. Additionally, we have proposed a novel two-state mean-field model based on Ising spins that offers an insight into the kinetics of the Aβ growth process and external perturbation effects on disintegration. Thus, the current simulation study, coupled with the newly introduced mean-field model, offers an insight into the detailed mechanisms underlying the Aβ aggregation process, guiding potential strategies for effective disintegration of Aβ aggregates.

Disaggregation of amyloid-beta fibrils via natural metabolites using long timescale replica exchange molecular dynamics simulation studies

CONTEXT

Amyloid fibrils are self-assembled fibrous protein aggregates that are associated with several presently incurable diseases such as Alzheimer's. disease that is characterized by the accumulation of amyloid fibrils in the brain, which leads to the formation of plaques and the death of brain cells. Disaggregation of amyloid fibrils is considered a promising approach to cure Alzheimer's disease. The mechanism of amyloid fibril formation is complex and not fully understood, making it difficult to develop drugs that can target the process. Diacetonamine and cystathionine are potential lead compounds to induce disaggregation of amyloid fibrils.

METHODS

In the current research, we have used long timescale molecular simulation studies and replica exchange molecular dynamics (REMD) for 1000 ns (1 μs) to examine the mechanisms by which natural metabolites can disaggregate amyloid-beta fibrils. Molecular docking was carried out using Glide and with prior protein minimization and ligand preparation. We focused on a screening a database of natural metabolites, as potential candidates for disaggregating amyloid fibrils. We used Desmond with OPLS 3e as a force field. MM-GBSA calculations were performed. Blood-brain barrier permeability, SASA, and radius of gyration parameters were calculated.

Protein Corona Formation and Aggregation of Amyloid β 1-40-Coated Gold Nanocolloids

Amyloid fibrillogenesis is a pathogenic protein aggregation process that occurs through a highly ordered process of protein-protein interactions. To better understand the protein-protein interactions involved in amyloid fibril formation, we formed nanogold colloid aggregates by stepwise additions of ∼2 nmol of amyloid β 1-40 peptide (Aβ1-40) at pH ∼3.7 and ∼25 °C. The processes of protein corona formation and building of gold colloid [diameters (d) of 20 and 80 nm] aggregates were confirmed by a red-shift of the surface plasmon resonance (SPR) band, λpeak, as the number of Aβ1-40 peptides [N(Aβ1-40)] increased. The normalized red-shift of λpeak, Δλ, was correlated with the degree of protein aggregation, and this process was approximated as the adsorption isotherm explained by the Langmuir-Freundlich model. As the coverage fraction (θ) was analyzed as a function of ϕ, which is the N(Aβ1-40) per total surface area of nanogold colloids available for adsorption, the parameters for explaining the Langmuir-Freundlich model were in good agreement for both 20 and 80 nm gold, indicating that ϕ could define the stage of the aggregation process. Surface-enhanced Raman scattering (SERS) imaging was conducted at designated values of ϕ and suggested that a protein-gold surface interaction during the initial adsorption stage may be dependent on the nanosize. The 20 nm gold case seems to prefer a relatively smaller contacting section, such as a -C-N or C═C bond, but a plane of the benzene ring may play a significant role for 80 nm gold. Regardless of the size of the particles, the β-sheet and random coil conformations were considered to be used to form gold colloid aggregates. The methodology developed in this study allows for new insights into protein-protein interactions at distinct stages of aggregation.

Perspective for Molecular Dynamics Simulation Studies of Amyloid-β Aggregates

The cause of Alzheimer's disease is related to aggregates such as oligomers and amyloid fibrils consisting of amyloid-β (Aβ) peptides. Molecular dynamics (MD) simulation studies have been conducted to understand the molecular mechanism of the formation and disruption of Aβ aggregates. In this Perspective, the MD simulation studies are classified into four categories, focusing on the target systems: aggregation of Aβ peptides in bulk solution, Aβ aggregation at the interface, aggregation inhibitor against Aβ peptides, and nonequilibrium MD simulation of Aβ aggregates. MD simulation studies in these categories are first reviewed. Future perspectives in each category are then presented. Finally, the overall perspective is presented on how MD simulations of Aβ aggregates can be utilized for developing Alzheimer's disease treatment.

Dissociation process of polyalanine aggregates by free electron laser irradiation

Polyalanine (polyA) disease-causative proteins with an expansion of alanine repeats can be aggregated. Although curative treatments for polyA diseases have not been explored, the dissociation of polyA aggregates likely reduces the cytotoxicity of polyA. Mid-infrared free electron laser (FEL) successfully dissociated multiple aggregates. However, whether the FEL dissociates polyA aggregates like other aggregates has not been tested. Here, we show that FEL at 6.1 μm experimentally weakened the extent of aggregation of a peptide with 13 alanine repeats (13A), and the irradiated 13A exerted lesser cytotoxicity to neuron-like cells than non-irradiated 13A. Then, we applied molecular dynamics (MD) simulation to follow the dissociation process by FEL. We successfully observed how the intermolecular β-sheet of polyA aggregates was dissociated and separated into monomers with helix structures upon FEL irradiation. After the dissociation by FEL, water molecules inhibited the reformation of polyA aggregates. We recently verified the same dissociation process using FEL-treated amyloid-β aggregates. Thus, a common mechanism underlies the dissociation of different protein aggregates that cause different diseases, polyA disease and Alzheimer's disease. However, MD simulation indicated that polyA aggregates are less easily dissociated than amyloid-β aggregates and require longer laser irradiation due to hydrophobic alanine repeats.

Metal-Organic Framework-Derived Carbon as a Photoacoustic Modulator of Alzheimer's β-Amyloid Aggregate Structure

Photoacoustic materials emit acoustic waves into the surrounding by absorbing photon energy. In an aqueous environment, light-induced acoustic waves form cavitation bubbles by altering the localized pressure to trigger the phase transition of liquid water into vapor. In this study, we report photoacoustic dissociation of beta-amyloid (Aβ) aggregates, a hallmark of Alzheimer's disease, by metal-organic framework-derived carbon (MOFC). MOFC exhibits a near-infrared (NIR) light-responsive photoacoustic characteristic that possesses defect-rich and entangled graphitic layers that generate intense cavitation bubbles by absorbing tissue-penetrable NIR light. According to our video analysis, the photoacoustic cavitation by MOFC occurs within milliseconds in the water, which was controllable by NIR light dose. The photoacoustic cavitation successfully transforms robust, β-sheet-dominant neurotoxic Aβ aggregates into nontoxic debris by changing the asymmetric distribution of water molecules around the Aβ's amino acid residues. This work unveils the therapeutic potential of NIR-triggered photoacoustic cavitation as a modulator of the Aβ aggregate structure.

Key Residue for Aggregation of Amyloid-β Peptides

It is known that oligomers of amyloid-β (Aβ) peptide are associated with Alzheimer's disease. Aβ has two isoforms: Aβ40 and Aβ42. Although the difference between Aβ40 and Aβ42 is only two additional C-terminal residues, Aβ42 aggregates much faster than Aβ40. It is unknown what role the C-terminal two residues play in accelerating aggregation. Since Aβ42 is more toxic than Aβ40, its oligomerization process needs to be clarified. Moreover, clarifying the differences between the oligomerization processes of Aβ40 and Aβ42 is essential to elucidate the key factors of oligomerization. Therefore, to investigate the dimerization process, which is the early oligomerization process, Hamiltonian replica-permutation molecular dynamics simulations were performed for Aβ40 and Aβ42. We identified a key residue, Arg5, for the Aβ42 dimerization. The two additional residues in Aβ42 allow the C-terminus to form contact with Arg5 because of the electrostatic attraction between them, and this contact stabilizes the β-hairpin. This β-hairpin promotes dimer formation through the intermolecular β-bridges. Thus, we examined the effects of amino acid substitutions of Arg5, thereby confirming that the mutations remarkably suppressed the aggregation of Aβ42. Moreover, the mutations of Arg5 suppressed the Aβ40 aggregation. It was found by analyzing the simulations that Arg5 is important for Aβ40 to form intermolecular contacts. Thus, it was clarified that the role of Arg5 in the oligomerization process varies due to the two additional C-terminal residues.

Insights into Molecular Mechanisms of EGCG and Apigenin on Disrupting Amyloid-Beta Protofibrils Based on Molecular Dynamics Simulations

The fibrillization and deposition of amyloid-beta (Aβ) protofibrils are one of the important factors leading to Alzheimer's disease. Molecular dynamics simulations can offer information on intermolecular interaction mechanisms between Aβ protofibrils and Aβ fibrillization inhibitors. Here, in this work, we explore the early molecular mechanisms of (-)-epigallocatechin-3-gallate (EGCG) and apigenin on disrupting Aβ42 protofibrils based on molecular simulations. The binding modes of EGCG and apigenin with the Aβ42 protofibril are obtained. Furthermore, we compare the behavioral mechanisms of EGCG and apigenin on disturbing the Aβ42 protofibril. Both EGCG and apigenin are able to decrease the proportion of the β-sheet and bend structures of the Aβ42 protofibril while inducing random coil structures. The results of hydrogen bonds and D23-K28 salt bridges illustrate that EGCG and apigenin have the ability of destabilizing the Aβ42 protofibril. Meanwhile, the van der Waals interactions between the EGCG and Aβ42 protofibril are shown to be larger than that of apigenin with the Aβ42 protofibril, but the electrostatic interactions between apigenin and the Aβ42 protofibril are dominant in the binding affinity. Our findings may help in designing effective drug candidates for disordering the Aβ protofibril and impeding Aβ fibrillization.

Stable Cavitation Interferes with Aβ16-22 Oligomerization

Ultrasound and microbubbles are used for many medical applications nowadays. Scanning ultrasound can remove amyloid-β (Aβ) aggregates in the mouse brain and restores memory in an Alzheimer's disease mouse model. In vitro studies showed that amyloid fibrils are fragmented due to the ultrasound-induced bubble inertial cavitation, and ultrasonic pulses accelerate the depolymerization of Aβ fibrils into monomers at 1 μM of concentration. Under applied ultrasound, microbubbles can be in a stable oscillating state or unstable inertial cavitation state. The latter occurs when ultrasound causes a dramatic change of bubble sizes above a certain acoustic pressure. We have developed and implemented a nonequilibrium molecular dynamics simulation algorithm to the AMBER package, to facilitate the investigation of the molecular mechanism of Aβ oligomerization under stable cavitation. Our results indicated that stable cavitation not only inhibited oligomeric formation, but also prevented the formation of β-rich oligomers. The network analysis of state transitions revealed that stable cavitation altered the oligomerization pathways of Aβ_16-22 peptides. Our simulation tool may be applied to optimize the experimental conditions to achieve the best therapeutical effect.

Molecular Dynamics Simulation Studies on the Aggregation of Amyloid-β Peptides and Their Disaggregation by Ultrasonic Wave and Infrared Laser Irradiation

Alzheimer’s disease is understood to be caused by amyloid fibrils and oligomers formed by aggregated amyloid-β (Aβ) peptides. This review article presents molecular dynamics (MD) simulation studies of Aβ peptides and Aβ fragments on their aggregation, aggregation inhibition, amyloid fibril conformations in equilibrium, and disruption of the amyloid fibril by ultrasonic wave and infrared laser irradiation. In the aggregation of Aβ, a β-hairpin structure promotes the formation of intermolecular β-sheet structures. Aβ peptides tend to exist at hydrophilic/hydrophobic interfaces and form more β-hairpin structures than in bulk water. These facts are the reasons why the aggregation is accelerated at the interface. We also explain how polyphenols, which are attracting attention as aggregation inhibitors of Aβ peptides, interact with Aβ. An MD simulation study of the Aβ amyloid fibrils in equilibrium is also presented: the Aβ amyloid fibril has a different structure at one end from that at the other end. The amyloid fibrils can be destroyed by ultrasonic wave and infrared laser irradiation. The molecular mechanisms of these amyloid fibril disruptions are also explained, particularly focusing on the function of water molecules. Finally, we discuss the prospects for developing treatments for Alzheimer’s disease using MD simulations.

Replica permutation with solute tempering for molecular dynamics simulation and its application to the dimerization of amyloid-β fragments

We propose the replica permutation with solute tempering (RPST) by combining the replica-permutation method (RPM) and the replica exchange with solute tempering (REST). Temperature permutations are performed among more than two replicas in RPM, whereas temperature exchanges are performed between two replicas in the replica-exchange method (REM). The temperature transition in RPM occurs more efficiently than in REM. In REST, only the temperatures of the solute region, the solute temperatures, are exchanged to reduce the number of replicas compared to REM. Therefore, RPST is expected to be an improved method taking advantage of these methods. For comparison, we applied RPST, REST, RPM, and REM to two amyloid-β(16–22) peptides in explicit water. We calculated the transition ratio and the number of tunneling events in the temperature space and the number of dimerization events of amyloid-β(16–22) peptides. The results indicate that, in RPST, the number of replicas necessary for frequent random walks in the temperature and conformational spaces is reduced compared to the other three methods. In addition, we focused on the dimerization process of amyloid-β(16–22) peptides. The RPST simulation with a relatively small number of replicas shows that the two amyloid-β(16–22) peptides form the intermolecular antiparallel β-bridges due to the hydrophilic side-chain contact between Lys and Glu and hydrophobic side-chain contact between Leu, Val, and Phe, which stabilizes the dimer of the peptides.

All-Atom Molecular Dynamics Simulation Methods for the Aggregation of Protein and Peptides: Replica Exchange/Permutation and Nonequilibrium Simulations

Protein aggregates are associated with more than 40 serious human diseases. To understand the formation mechanism of protein aggregates at atomic level, all-atom molecular dynamics (MD) simulation is a powerful computational tool. In this chapter, we review the all-atom MD simulation methods that are useful for study on the protein aggregation. We first explain conventional MD simulation methods in physical statistical ensembles, such as the canonical and isothermal-isobaric ensembles. We then describe the generalized-ensemble algorithms such as replica-exchange and replica-permutation MD methods. These methods can overcome a difficulty, in which simulations tend to get trapped in local-minimum free-energy states. Finally we explain the nonequilibrium MD method. Some simulation results based on these methods are also presented.

Computer Simulations Aimed at Exploring Protein Aggregation and Dissociation

Protein aggregation can lead to well-defined structures that are functional, but is also the cause of the death of neuron cells in many neurodegenerative diseases. The complexity of the molecular events involved in the aggregation kinetics of amyloid proteins and the transient and heterogeneous characters of all oligomers prevent high-resolution structural experiments. As a result, computer simulations have been used to determine the atomic structures of amyloid proteins at different association stages as well as to understand fibril dissociation. In this chapter, we first review the current computer simulation methods used for aggregation with some atomistic and coarse-grained results aimed at better characterizing the early formed oligomers and amyloid fibril formation. Then we present the applications of non-equilibrium molecular dynamics simulations to comprehend the dissociation of protein assemblies.

Dynamics of Amyloid Formation from Simplified Representation to Atomistic Simulations

Amyloid fibril formation is an intrinsic property of short peptides, non-disease proteins, and proteins associated with neurodegenerative diseases. Aggregates of the Aβ and tau proteins, the α-synuclein protein, and the prion protein are observed in the brain of Alzheimer's, Parkinson's, and prion disease patients, respectively. Due to the transient short-range and long-range interactions of all species and their high aggregation propensities, the conformational ensemble of these devastating proteins, the exception being for the monomeric prion protein, remains elusive by standard structural biology methods in bulk solution and in lipid membranes. To overcome these limitations, an increasing number of simulations using different sampling methods and protein models have been performed. In this chapter, we first review our main contributions to the field of amyloid protein simulations aimed at understanding the early aggregation steps of short linear amyloid peptides, the conformational ensemble of the Aβ40/42 dimers in bulk solution, and the stability of Aβ aggregates in lipid membrane models. Then we focus on our studies on the interactions of amyloid peptides/inhibitors to prevent aggregation, and long amyloid sequences, including new results on a monomeric tau construct.

Role of Water Molecules and Helix Structure Stabilization in the Laser-Induced Disruption of Amyloid Fibrils Observed by Nonequilibrium Molecular Dynamics Simulations

Water plays a crucial role in the formation and destruction of biomolecular structures. The mechanism for destroying biomolecular structures was thought to be an active breaking of hydrogen bonds by water molecules. However, using nonequilibrium molecular dynamics simulations, in which an amyloid-β amyloid fibril was destroyed via infrared free-electron laser (IR-FEL) irradiation, we discovered a new mechanism, in which water molecules disrupt protein aggregates. The intermolecular hydrogen bonds formed by C═O and N-H in the fibril are broken at each pulse of laser irradiation. These bonds spontaneously re-form after the irradiation in many cases. However, when a water molecule happens to enter the gap between C═O and N-H, it inhibits the re-formation of the hydrogen bonds. Such sites become defects in the regularly aligned hydrogen bonds, from which all hydrogen bonds in the intermolecular β-sheet are broken as the fraying spreads. This role of water molecules is entirely different from other known mechanisms. This new mechanism can explain the recent experiments showing that the amyloid fibrils are not destroyed by laser irradiation under dry conditions. Additionally, we found that helix structures form more after the amyloid disruption; this is because the resonance frequency is different in a helix structure. Our findings provide a theoretical basis for the application of IR-FEL to the future treatment of amyloidosis.

Molecular mechanism of ultrasound interaction with a blood brain barrier model

The brain is strictly protected by the blood brain barrier preventing the crossing of therapeutics to treat brain diseases. The high and low intensity focused ultrasound methods have been used to temporarily open the blood brain barrier, facilitating the transport of drugs. The methods are very promising because the opening is transient, localized, and noninvasive. However, the molecular mechanism of the opening is unknown, and this limits the development and application of these methods. With this in mind, we carry out a molecular dynamics simulation study to understand the interaction of ultrasound with the cell membrane and the tight junction. Our minimal blood brain barrier model is composed of two lipid bilayers, mimicking two portions of neighboring cells, connected together by a tight junction formed by a pair of two cis-dimers of the claudin-5 protein. Using an experimental ultrasound frequency of 50 MHz, simulations show that at low intensities, ultrasound does not impact the structure of the cell membranes and tight junction, implying that the direct interaction of ultrasound with the blood brain barrier is not responsible for the experimentally observed opening. At high intensities, the ultrasound pulls the monolayers of individual cell membrane lipid bilayers apart, creating air compartments inside the bilayers. This reduces the free energy barrier for the translocation of drugs across the lipid bilayer and enhances drug permeability. At very high intensities, the two monolayers are largely separated, resulting in cell damage and implying that the blood brain barrier is primarily opened at the experimentally observed damaged areas.

Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis

Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aβ, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.

Promotion and Inhibition of Amyloid-β Peptide Aggregation: Molecular Dynamics Studies

Aggregates of amyloid-β (Aβ) peptides are known to be related to Alzheimer’s disease. Their aggregation is enhanced at hydrophilic–hydrophobic interfaces, such as a cell membrane surface and air-water interface, and is inhibited by polyphenols, such as myricetin and rosmarinic acid. We review molecular dynamics (MD) simulation approaches of a full-length Aβ peptide, Aβ40, and Aβ(16–22) fragments in these environments. Since these peptides have both hydrophilic and hydrophobic amino acid residues, they tend to exist at the interfaces. The high concentration of the peptides accelerates the aggregation there. In addition, Aβ40 forms a β-hairpin structure, and this structure accelerates the aggregation. We also describe the inhibition mechanism of the Aβ(16–22) aggregation by polyphenols. The aggregation of Aβ(16–22) fragments is caused mainly by the electrostatic attraction between charged amino acid residues known as Lys16 and Glu22. Since polyphenols form hydrogen bonds between their hydroxy and carboxyl groups and these charged amino acid residues, they inhibit the aggregation.

Nonequilibrium molecular dynamics simulations of infrared laser-induced dissociation of a tetrameric Aβ42 β-barrel in a neuronal membrane model

Experimental studies have reported that the amyloid-β proteins can form pores in cell membranes, and this could be one possible source of toxicity in Alzheimer's disease. Dissociation of these pores could therefore be a potential therapeutic approach. It is known that high photon density free-electron laser experiments and laser-induced nonequilibrium molecular dynamics simulations (NEMD) can dissociate amyloid fibrils at specific frequencies in vitro. Our question is whether NEMD simulations can dissociate amyloid pores in a bilayer mimicking a neuronal membrane, and as an example, we select a tetrameric Aβ42 β-barrel. Our simulations shows that the resonance between the laser field and the amide I vibrational mode of the barrel destabilises all intramolecular and intermolecular hydrogen bonds of Aβ42 and converts the β-barrel to a random/coil disordered oligomer. Starting from this disordered oligomer, extensive standard MD simulations shows sampling of disordered Aβ42 states without any increase of β-sheet and reports that the orientational order of lipids is minimally disturbed. Interestingly, the frequency to be employed to dissociate this beta-barrel is specific to the amino acid sequence. Taken together with our previous simulation results, this study indicates that infrared laser irradiation can dissociate amyloid fibrils and oligomers in bulk solution and in a membrane environment without affecting the surrounding molecules, offering therefore a promising way to retard the progression of AD.

Disaggregation of Amyloid-β Plaques by a Local Electric Field Generated by a Vertical Nanowire Electrode Array

The aggregation and accumulation of amyloid-β (Aβ) peptides is a characteristic pathology for Alzheimer's disease (AD). Although noninvasive therapies involving stimulation by electric field (EF) have been reported, the efficiency of Aβ disaggregation needs to be further improved for this strategy to be used in clinical settings. In this study, we show that an electrode based on a vertical nanowire electrode array (VNEA) is far more superior to a typical flat-type electrode in disaggregating Aβ plaques. The enhanced disaggregation efficiency of VNEA is due to the formation of high-strength local EF between the nanowires, as verified by in silico and empirical evidence. Compared with those of the flat electrode, the simulation data revealed that 19.8-fold and 8.8-fold higher EFs are generated above and between the nanowires, respectively. Moreover, empirical cyclic voltammetry data demonstrated that VNEA had a 2.7-fold higher charge capacity than the flat electrode; this is associated with the higher surface area of VNEA. The conformational transition of Aβ peptides between the β-sheet and α-helix could be sensitively monitored in real time by the newly designed in situ circular dichroism instrument. This highly efficient EF-configuration of VNEA will lower the stimulation power for disaggregating the Aβ plaques, compared to that of other existing field-mediated modulation systems. Considering the complementary metal-oxide-semiconductor-compatibility and biocompatible strength of the EF for perturbing the Aβ aggregation, our study could pave the way for the potential use of electric stimulation devices for in vivo therapeutic application as well as scientific studies for AD.

Molecular dynamics simulation of soundwave propagation in a simple fluid

A molecular dynamics (MD) simulation was performed to study the propagation of soundwaves in a fluid. Soundwaves are generated by a sinusoidally oscillating wall and annihilated by a locally applied Langevin thermostat near the opposite wall. The waveform changes from sinusoidal to sawtooth with increasing wave amplitude. For low-frequency sounds, the simulation results show a very good agreement with Burgers's equation without any fitting parameters. In contrast, for high-frequency sounds, significant deviations are obtained because of acoustic streaming. The speed of sound can be directly determined from the Fourier transform of a waveform with high accuracy. Although obtaining the attenuation rate directly from the simulation results is difficult because of the nonlinear effects of the wave amplitude, it can be estimated via Burgers's equation. The results demonstrate that MD simulations are a useful tool for the quantitative analysis of soundwaves.

Infrared Laser-Induced Amyloid Fibril Dissociation: A Joint Experimental/Theoretical Study on the GNNQQNY Peptide

Neurodegenerative diseases are usually characterized by plaques made of well-ordered aggregates of distinct amyloid proteins. Dissociating these very stable amyloid plaques is a critical clinical issue. In this study, we present a joint mid-infrared free electron laser experiment/nonequilibrium molecular dynamics simulation to understand the dissociation process of a representative example GNNQQNY fibril. By tuning the laser frequency to the amide I band of the fibril, the resonance takes place and dissociation is occurred. With the calculated and observed wide-angle X-ray scattering profiles and secondary structures before and after laser irradiation being identical, we can propose a dissociation mechanism with high confidence from our simulations. We find that dissociation starts in the core of the fibrils by fragmenting the intermolecular hydrogen bonds and separating the peptides and then propagates to the fibril extremities leading to the formation of unstructured expanded oligomers. We suggest that this should be a generic mechanism of the laser-induced dissociation of amyloid fibrils.

Structures of the intrinsically disordered Aβ, tau and α-synuclein proteins in aqueous solution from computer simulations

Intrinsically disordered proteins (IDPs) play many biological roles in the human proteome ranging from vesicular transport, signal transduction to neurodegenerative diseases. The Aβ and tau proteins, and the α-synuclein protein, key players in Alzheimer's and Parkinson's diseases, respectively are fully disordered at the monomer level. The structural heterogeneity of the monomeric and oligomeric states and the high self-assembly propensity of these three IDPs have precluded experimental structural determination. Simulations have been used to determine the atomic structures of these IDPs. In this article, we review recent computer models to capture the equilibrium ensemble of Aβ, tau and α-synuclein proteins at different association steps in aqueous solution and present new results of the PEP-FOLD framework on α-synuclein monomer.

Nonequilibrium atomistic molecular dynamics simulation of tubular nanomotor propelled by bubble propulsion

We develop a molecular nanoscaled model for tubular motors propelled by bubble propulsion. The motor is modeled by a carbon nanotube, and the bubble is represented by a particle interacting with water by a time-dependent potential. Effects of liquid viscosity, fuel concentration, geometry, and size of the tube on the performance of the motor are effectively encoded into two parameters: time scales of the bubble expansion and bubble formation. Our results are qualitatively consistent with experimental data of much larger motors. Simulations suggest that (i) the displacement of the tube is optimized if two time scales are as short as possible, (ii) the compromise between the performance and fuel consumption is achieved if the bubble formation time is shorter than the velocity correlation time of the tube, (iii) the motor efficiency is higher with slow expansion, short formation of the bubble than fast growth but long formation time, and (iv) the tube is propelled by strong forces on the order of mN, reaching high speeds up to ∼60 m/s. Our simulation may be useful for refining and encouraging future experimental work on nanomotors having the size of a few nanometers. The tiny size and high speed motors could have great potential applications in real life.

Interaction mechanism between the focused ultrasound and lipid membrane at the molecular level

Focused ultrasound (FUS) has a wide range of medical applications. Nowadays, the diagnostic and therapeutic ultrasound procedures are routinely used; effects of ultrasound on biological systems at the molecular level are, however, not fully understood. Experimental results on the interaction of the cell membrane, a simplest but important system component, with ultrasound are controversial. Molecular dynamics (MD) simulations could provide valuable insights, but there is no single study on the mechanism of the FUS induced structural changes in cell membranes. With this in mind, we develop a simple method to include FUS into a standard MD simulation. Adopting the 1,2-dioleoyl-sn-glycero-3-phosphocholine lipid membrane as a representative model described by the MARTINI coarse-grained force field, and using experimental values of the ultrasound frequency and intensity, we show that the heat and bubble cavitation are not the primary direct mechanisms that cause structural changes in the membrane. The spatial pressure gradients between the focused and free regions and between the parallel and perpendicular directions to the membrane are the origin of the mechanism. These gradients force lipids to move out of the focused region, forming a lipid flow along the membrane diagonal. Lipids in the free region move in the opposite direction due to the conservation of the total momentum. These opposite motions create wrinkles along the membrane diagonal at low FUS intensities and tear up the membrane at high FUS intensities. Once the membrane is torn up, it is not easy to reform. The implication of our findings in the FUS-induced drug delivery is discussed in some detail.

Effect of Shock-Induced Cavitation Bubble Collapse on the damage in the Simulated Perineuronal Net of the Brain

The purpose of this study is to conduct modeling and simulation to understand the effect of shock-induced mechanical loading, in the form of cavitation bubble collapse, on damage to the brain's perineuronal nets (PNNs). It is known that high-energy implosion due to cavitation collapse is responsible for corrosion or surface damage in many mechanical devices. In this case, cavitation refers to the bubble created by pressure drop. The presence of a similar damage mechanism in biophysical systems has long being suspected but not well-explored. In this paper, we use reactive molecular dynamics (MD) to simulate the scenario of a shock wave induced cavitation collapse within the perineuronal net (PNN), which is the near-neuron domain of a brain's extracellular matrix (ECM). Our model is focused on the damage in hyaluronan (HA), which is the main structural component of PNN. We have investigated the roles of cavitation bubble location, shockwave intensity and the size of a cavitation bubble on the structural evolution of PNN. Simulation results show that the localized supersonic water hammer created by an asymmetrical bubble collapse may break the hyaluronan. As such, the current study advances current knowledge and understanding of the connection between PNN damage and neurodegenerative disorders.