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.

The Double-Layered Structure of Amyloid-ß Assemblage on GM1-Containing Membranes Catalytically Promotes Fibrillization.

Alzheimer's disease (AD) is associated with progressive accumulation of amyloid-ß (Aß) cross-ß fibrils in the brain. Aß species tightly associated with GM1 ganglioside, a glycosphingolipid abundant in neuronal membranes, promote amyloid fibril formation; therefore, they could be attractive clinical targets. However, the active conformational state of Aß in GM1-containing lipid membranes is still unknown. The present solid-state nuclear magnetic resonance study revealed a nonfibrillar Aß assemblage characterized by a double-layered antiparallel ß-structure specifically formed on GM1 ganglioside clusters. Our data show that this unique assemblage was not transformed into fibrils on GM1-containing membranes but could promote conversion of monomeric Aß into fibrils, suggesting that a solvent-exposed hydrophobic layer provides a catalytic surface evoking Aß fibril formation. Our findings offer structural clues for designing drugs targeting catalytically active Aß conformational species for the development of anti-AD therapeutics.

Ingenuity in performing replica permutation: How to order the state labels for improving sampling efficiency.

In the replica-permutation method, an advanced version of the replica-exchange method, all combinations of replicas and parameters are considered for parameter permutation, and a list of all the combinations is prepared. Here, we report that the temperature transition probability depends on how the list is created, especially in replica permutation with solute tempering (RPST). We found that the transition probabilities decrease at large replica indices when the combinations are sequentially assigned to the state labels as in the originally proposed list. To solve this problem, we propose to modify the list by randomly assigning the combinations to the state labels. We performed molecular dynamics simulations of amyloid-ß(16-22) peptides using RPST with the "randomly assigned" list (RPST-RA) and RPST with the "sequentially assigned" list (RPST-SA). The results show the decreases in the transition probabilities in RPST-SA are eliminated, and the sampling efficiency is improved in RPST-RA.

Inhibition of amyloid-ß(16-22) aggregation by polyphenols using replica permutation with solute tempering molecular dynamics simulation.

Aggregates of amyloid-ß (Aß) peptides are thought to cause Alzheimer's disease. Polyphenolic compounds are known to inhibit Aß aggregation. We applied replica permutation with solute tempering (RPST) to the system of Aß fragments, Aß(16-22), and polyphenols to elucidate the mechanism of inhibition of Aß aggregation. The RPST molecular dynamics simulations were performed for two polyphenols, myricetin (MYC) and rosmarinic acid (ROA). Two Aß fragments were distant, and the number of residues forming the intermolecular ß-sheet was reduced in the presence of MYC and ROA compared with that in the absence of polyphenols. MYC was found to interact with glutamic acid and phenylalanine of Aß fragments. These interactions induce helix structure formation of Aß fragments, making it difficult to form ß-sheet. ROA interacted with glutamic acid and lysine, which reduced the hydrophilic interaction between Aß fragments. These results indicate that these polyphenols inhibit the aggregation of Aß fragments with different mechanisms.

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.

State-of-the-Art Molecular Dynamics Simulation Studies of RNA-Dependent RNA Polymerase of SARS-CoV-2.

Molecular dynamics (MD) simulations are powerful theoretical methods that can reveal biomolecular properties, such as structure, fluctuations, and ligand binding, at the level of atomic detail. In this review article, recent MD simulation studies on these biomolecular properties of the RNA-dependent RNA polymerase (RdRp), which is a multidomain protein, of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are presented. Although the tertiary structures of RdRps in SARS-CoV-2 and SARS-CoV are almost identical, the RNA synthesis activity of RdRp of SARS-CoV is higher than SARS-CoV-2. Recent MD simulations observed a difference in the dynamic properties of the two RdRps, which may cause activity differences. RdRp is also a drug target for Coronavirus disease 2019 (COVID-19). Nucleotide analogs, such as remdesivir and favipiravir, are considered to be taken up by RdRp and inhibit RNA replication. Recent MD simulations revealed the recognition mechanism of RdRp for these drug molecules and adenosine triphosphate (ATP). The ligand-recognition ability of RdRp decreases in the order of remdesivir, favipiravir, and ATP. As a typical recognition process, it was found that several lysine residues of RdRp transfer these ligand molecules to the binding site such as a "bucket brigade." This finding will contribute to understanding the mechanism of the efficient ligand recognition by RdRp. In addition, various simulation studies on the complexes of SARS-CoV-2 RdRp with several nucleotide analogs are reviewed, and the molecular mechanisms by which these compounds inhibit the function of RdRp are discussed. The simulation studies presented in this review will provide useful insights into how nucleotide analogs are recognized by RdRp and inhibit the RNA replication.

Molecular dynamics simulations of amyloid-ß peptides in heterogeneous environments.

Alzheimer's disease is thought to be caused by the aggregation of amyloid-ß (Aß) peptides. Their aggregation is accelerated at hydrophilic/hydrophobic interfaces such as the air-water interface and the surface of monosialotetrahexosylganglioside (GM1) clusters on neuronal cell membranes. In this review, we present recent studies of full-length Aß (Aß40) peptides and Aß(16-22) fragments in such heterogeneous environments by molecular dynamics (MD) simulations. These peptides have both hydrophilic and hydrophobic amino-acid residues and tend to exist at the hydrophilic/hydrophobic interface. Therefore, the peptide concentration increases at the interface, which is one of the factors that promote aggregation. Furthermore, it was found that Aß40 forms an α-helix structure and then a ß-hairpin structure at the interface. The ß-hairpin promotes the formation of oligomers with intermolecular ß-sheets. It means that not only the high concentration of Aß40 at the interface but also the structure of Aß40 itself promotes aggregation. In addition, MD simulations of Aß40 on recently-developed GM1-glycan clusters showed that the HHQ (13-15) segment of Aß40 is important for the recognition of GM1-glycan clusters. It was also elucidated that Aß40 forms a helix structure in the C-terminal region on the GM1-glycan cluster. This result suggests that the helix formation, which is the first step in the conformational changes toward pathological aggregation, is initiated at the GM1-glycan moieties rather than at the lipid-ceramide moieties. These studies will enhance the physicochemical understanding of the structural changes of Aß at the heterogeneous interfaces and the mechanism of Alzheimer's disease pathogenesis.

Critical contributions of pre-S1 shoulder and distal TRP box in DAG-activated TRPC6 channel by PIP2 regulation.

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2 or PIP2) regulates the activities of numerous membrane proteins, including diacylglycerol(DAG)-activated TRPC3/6/7 channels. Although PIP2 binding is known to support DAG-activated TRP channel activity, its binding site remains unknown. We screened for PIP2 binding sites within TRPC6 channels through extensive mutagenesis. Using voltage-sensitive phosphatase (DrVSP), we found that Arg437 and Lys442, located in the channel's pre-S1 domain/shoulder, are crucial for interaction with PIP2. To gain structural insights, we conducted computer protein-ligand docking simulations with the pre-S1 domain/shoulder of TRPC6 channels. Further, the functional significance of PIP2 binding to the pre-S1 shoulder was assessed for receptor-operated channel functions, cross-reactivity to DAG activation, and the kinetic model simulation. These results revealed that basic residues in the pre-S1 domain/shoulder play a central role in the regulation of PIP2-dependent gating. In addition, neutralizing mutation of K771 in the distal TRP box reversed the effect of PIP2 depletion from inhibiting to potentiating channel activity. A similar effect was seen in TRPV1 channels, which suggests that TRPC6 possesses a common but robust polarity switch mediating the PIP2-dependent effect. Overall, these mutagenesis studies reveal functional and structural insights for how basic residues and channel segments in TRP channels are controlled through phosphoinositides recognition.

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.

Tardigrade Secretory-Abundant Heat-Soluble Protein Varies Entrance Propensity Depending on the Amino-Acid Sequence.

Secretory-abundant heat-soluble (SAHS) proteins, which constitute a protein family unique to tardigrades, are thought to be essential for anhydrobiosis. Our previous study has revealed that one of the SAHS proteins of Ramazzottius varieornatus (RvSAHS1) has a more flexible entrance than a mammalian fatty-acid-binding protein, which has a crystal structure similar to that of RvSAHS1. Recently, SAHS paralogs that are expressed abundantly and specifically in the early embryos of this tardigrade and Hypsibius exemplaris have been identified. Comparing these amino-acid sequences with that of RvSAHS1, we have found characteristic differences as I113F and D146T. In this study, we investigate I113F and D146T mutants' properties of RvSAHS1 using molecular dynamics simulations and compare the structures and fluctuations of their entrances with those of the wild type. The two mutants exhibit different properties at the entrance of the ß-barrel structure. The I113F mutant tends to close the entrance more than the wild type due to the enhanced hydrophobic network inside the cavity. The D146T mutant, in contrast to the I113F mutant, tends to open the entrance. The mechanism by which this mutation opens the entrance is also discussed. Even though only a single mutation located far from the entrance is added to the wild type, there is a clear difference in the tendency to open and close the ß-barrel entrance. It indicates that the entrance properties of the SAHS protein are sensitive to the amino-acid sequence.

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.

Tardigrade Secretory-Abundant Heat-Soluble Protein Has a Flexible ß-Barrel Structure in Solution and Keeps This Structure in Dehydration.

Secretory-abundant heat-soluble (SAHS) proteins are unique heat-soluble proteins of Tardigrada and are believed to play an essential role in anhydrobiosis, a latent state of life induced by desiccation. To investigate the dynamic properties, molecular dynamics (MD) simulations of a SAHS protein, RvSAHS1, were performed in solution and under dehydrating conditions. For comparison purposes, MD simulations of a human liver-type fatty-acid binding protein (LFABP) were performed in solution. Furthermore, high-speed atomic force microscopy observations were conducted to ascertain the results of the MD simulations. Three properties of RvSAHS1 were found as follows. (1) The entrance region of RvSAHS1 is more flexible and can be more extensive in solutions compared with that of a human LFABP because there is no salt bridge between the ßD and ßE strands. (2) The intrinsically disordered domain in the N-terminal region significantly fluctuates and can form an amphiphilic α-helix. (3) The size of the entrance region gets smaller along with dehydration, keeping the ß-barrel structure. Overall, the obtained results provide atomic-level dynamics of SAHS proteins.

"Bucket brigade" using lysine residues in RNA-dependent RNA polymerase of SARS-CoV-2.

The RNA-dependent RNA polymerase (RdRp) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a promising drug target for coronavirus disease 2019 (COVID-19) because it plays the most important role in the replication of the RNA genome. Nucleotide analogs such as remdesivir and favipiravir are thought to interfere with the RNA replication by RdRp. More specifically, they are expected to compete with nucleoside triphosphates, such as ATP. However, the process in which these drug molecules and nucleoside triphosphates are taken up by RdRp remains unknown. In this study, we performed all-atom molecular dynamics simulations to clarify the recognition mechanism of RdRp for these drug molecules and ATP that were at a distance. The ligand recognition ability of RdRp decreased in the order of remdesivir, favipiravir, and ATP. We also identified six recognition paths. Three of them were commonly found in all ligands, and the remaining three paths were ligand-dependent ones. In the common two paths, it was observed that the multiple lysine residues of RdRp carried the ligands to the binding site like a "bucket brigade." In the remaining common path, the ligands directly reached the binding site. Our findings contribute to the understanding of the efficient ligand recognition by RdRp at the atomic level.

Dynamic properties of SARS-CoV and SARS-CoV-2 RNA-dependent RNA polymerases studied by molecular dynamics simulations.

One of the promising drug targets against COVID-19 is an RNA-dependent RNA polymerase (RdRp) of SARS-CoV-2. The tertiary structures of the SARS-CoV-2 and SARS-CoV RdRps are almost the same. However, the RNA-synthesizing activity of the SARS-CoV RdRp is higher than that of the SARS-CoV-2 RdRp. We performed molecular dynamics simulations and found differences in their dynamic properties. In the SARS-CoV RdRp, motifs A-G, which form the active site, are up to 63% closer to each other. We also observed cooperative domain motion in the SARS-CoV RdRp. Such dynamic differences may cause the activity differences between the two RdRps.

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.

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.

Structural basis for promiscuous action of monoterpenes on TRP channels.

Monoterpenes are major constituents of plant-derived essential oils and have long been widely used for therapeutic and cosmetic applications. The monoterpenes menthol and camphor are agonists or antagonists for several TRP channels such as TRPM8, TRPV1, TRPV3 and TRPA1. However, which regions within TRPV1 and TRPV3 confer sensitivity to monoterpenes or other synthesized chemicals such as 2-APB are unclear. In this study we identified conserved arginine and glycine residues in the linker between S4 and S5 that are related to the action of these chemicals and validated these findings in molecular dynamics simulations. The involvement of these amino acids differed between TRPV3 and TRPV1 for chemical-induced and heat-evoked activation. These findings provide the basis for characterization of physiological function and biophysical properties of ion channels.

Involvement of pore helix in voltage-dependent inactivation of TRPM5 channel.

The transient receptor potential melastatin 5 (TRPM5) channel is a monovalent-permeable cation channel that is activated by intracellular Ca2+. Expression of TRPM5 has been shown in taste cells, pancreas, brainstem and olfactory epithelium, and this channel is thought to be involved in controlling membrane potentials. In whole-cell patch-clamp recordings, TRPM5 exhibited voltage-dependent inactivation at negative membrane potentials and time constant of voltage-dependent inactivation of TRPM5 did not depend on the intracellular Ca2+ concentrations between 100 and 500 nM. Alanine substitution at Y913 and I916 in the pore helix of TRPM5 increased time constant of voltage-dependent inactivation. Meanwhile, voltage-dependent inactivation was reduced in TRPM5 mutants having glycine substitution at L901, Y913, Q915 and I916 in the pore helix. From these results, we conclude that the pore helix in the outer pore loop might play a role in voltage-dependent inactivation of TRPM5.