The NSP4 T492I mutation increases SARS-CoV-2 infectivity by altering non-structural protein cleavage

The historically dominant SARS-CoV-2 Delta variant and the currently dominant Omicron variants carry a T492I substitution within the non-structural protein 4 (NSP4). Based on in silico analyses, we hypothesized that the T492I mutation increases viral transmissibility and adaptability, which we confirmed with competition experiments in hamster and human airway tissue culture models. Furthermore, we showed that the T492I mutation increases the replication capacity and infectiveness of the virus and improves its ability to evade host immune responses. Mechanistically, the T492I mutation increases the cleavage efficiency of the viral main protease NSP5 by enhancing enzyme-substrate binding, which increases production of nearly all non-structural proteins processed by NSP5. Importantly, the T492I mutation suppresses viral-RNA-associated chemokine production in monocytic macrophages, which may contribute to the attenuated pathogenicity of Omicron variants. Our results highlight the importance of NSP4 adaptation in the evolutionary dynamics of SARS-CoV-2.

Structure-Based Discovery of a Novel Allosteric Inhibitor against Human Dopamine Transporter

Human dopamine transporter (hDAT) regulates the reuptake of extracellular dopamine (DA) and is an essential therapeutic target for central nervous system (CNS) diseases. The allosteric modulation of hDAT has been identified for decades. However, the molecular mechanism underlying the transportation is still elusive, which hinders the rational design of allosteric modulators against hDAT. Here, a systematic structure-based method was performed to explore allosteric sites on hDAT in inward-open (IO) conformation and to screen compounds with allosteric affinity. First, the model of the hDAT structure was constructed based on the recently reported Cryo-EM structure of the human serotonin transporter (hSERT) and Gaussian-accelerated molecular dynamics (GaMD) simulation was further utilized for the identification of intermediate energetic stable states of the transporter. Then, with the potential druggable allosteric site on hDAT in IO conformation, virtual screening of seven enamine chemical libraries (∼440,000 compounds) was processed, resulting in 10 compounds being purchased for in vitro assay and with Z1078601926 discovered to allosterically inhibit hDAT (IC₅₀ = 0.527 [0.284; 0.988] μM) when nomifensine was introduced as an orthosteric ligand. Finally, the synergistic effect underlying the allosteric inhibition of hDAT by Z1078601926 and nomifensine was explored using additional GaMD simulation and postbinding free energy analysis. The hit compound discovered in this work not only provides a good starting point for lead optimization but also demonstrates the usability of the method for the structure-based discovery of novel allosteric modulators of other therapeutic targets.

Proteomics profiles reveal the potential roles of proteins involved in chicken macrophages stimulated by Lipopolysaccharide

Lipopolysaccharide (LPS), a core part of gram-negative bacteria, is crucial for inducing an inflammatory response in living things. In the current study, we used LPS from Salmonella to stimulate chicken macrophages (HD11). Proteomics was used to investigate immune-related proteins and their roles further. Proteomics investigation revealed 31 differential expression proteins (DEPs) after 4 hours of LPS infection. 24 DEPs expressions were up-regulated, while seven were down-regulated. In this investigation, ten DEPs were mainly enriched in S. aureus infection, complement, and coagulation cascades, which were all implicated in the inflammatory response and clearance of foreign pathogens. Notably, complement C3 was shown to be up-regulated in all immune-related pathways, indicating that it is a potential protein in this study. This work contributes to a better understanding and clarification of the processes of Salmonella infection in chickens. It might bring up new possibilities for treating and breeding Salmonella-infected chickens.

Boosting the capacity and stability of a MoO3 cathode via valence regulation and polypyrrole coating for a rechargeable Zn ion battery

Molybdenum trioxide (MoO₃) is emerging as a hugely competitive cathode material for aqueous zinc ion batteries (ZIBs) for its high theoretical capacity and electrochemical activity. Nevertheless, owing to its undesirable electronic transport capability and poor structural stability, the practical capacity and cycling performance of MoO₃ are yet unsatisfactory, which greatly blocks its commercial use. In this work, we report an effective approach to first synthesise nanosized MoO_3-x materials to provide more active specific surface areas, while improving the capacity and cycle life of MoO₃ by introducing low valence Mo and coated polypyrrole (PPy). MoO₃ nanoparticles with low-valence-state Mo and PPy coating (denoted as MoO_3-x@PPy) are synthesized via a solvothermal method and subsequent electrodeposition process. The as-prepared MoO_3-x@PPy cathode delivers a high reversible capacity of 212.4 mA h g^-1 at 1 A g^-1 with good cycling life (more than 75% capacity retention after 500 cycles). In contrast, the original commercial MoO₃ sample only obtains a capacity of 99.3 mA h g^-1 at 1 A g^-1, and a cycling stability of 10% capacity retention over 500 cycles. Additionally, the fabricated Zn//MoO_3-x@PPy battery obtains a maximum energy density of 233.6 W h kg^-1 and a power density of 11.2 kW kg^-1. Our results provide an efficient and practical approach to enhance commercial MoO₃ materials as high-performance cathodes for AZIBs.

Multi-state Model-Based Identification of Cryptic Allosteric Sites on Human Serotonin Transporter

Serotonin transporter (SERT) plays a fundamental role in taking the synaptic cleft serotonin back to the presynaptic neuron. The discovery of allosteric SERT modulators represents the next-generation medication for psychiatric disorders such as depression. Here, based on the cryo-EM structures of ibogaine in complex with SERT in distinct conformations, the multiple functional structures of the transporter bound to serotonin, including outward-open (OO_holo), outward-occluded (OC_holo), and inward-open (IO_holo and IO_holo'), were carefully characterized by induced-fit docking Gaussian-accelerated molecular dynamics (IFD-GaMD) simulation and the free-energy landscape analysis. Further MM/GBSA binding free energy, per-residue contribution, and molecular interaction fingerprint calculations revealed the interaction variations of serotonin with SERT in functional structures, which confirmed the allostery of SERT during serotonin reuptake. Moreover, five unique cryptic allosteric sites, which are druggable and capable of targeting by small molecules, were identified on the characterized multistate structures. These results provide structural and energetic information for the molecular mechanism of serotonin reuptake and will provide opportunities for the development of novel therapeutics based on the identified new allosteric sites on SERT.

Sodium Ion Pre-Intercalation of δ-MnO2 Nanosheets for High Energy Density Aqueous Zinc-Ion Batteries

With the merits of low cost, environmental friendliness and rich resources, manganese dioxide is considered to be a promising cathode material for aqueous zinc-ion batteries (AZIBs). However, its low ion diffusion and structural instability greatly limit its practical application. Hence, we developed an ion pre-intercalation strategy based on a simple water bath method to grow in situ δ-MnO₂ nanosheets on flexible carbon cloth substrate (MnO₂), while pre-intercalated Na⁺ in the interlayer of δ-MnO₂ nanosheets (Na-MnO₂), which effectively enlarges the layer spacing and enhances the conductivity of Na-MnO₂. The prepared Na-MnO₂//Zn battery obtained a fairly high capacity of 251 mAh g^-1 at a current density of 2 A g^-1, a satisfactory cycle life (62.5% of its initial capacity after 500 cycles) and favorable rate capability (96 mAh g^-1 at 8 A g^-1). Furthermore, this study revealed that the pre-intercalation engineering of alkaline cations is an effective method to boost the properties of δ-MnO₂ zinc storage and provides new insights into the construction of high energy density flexible electrodes.

Unbiased molecular dynamics simulation of a first-in-class small molecule inhibitor binds to oncostatin M

Small molecule inhibitors (SMIs) targeting oncostatin M (OSM) signaling pathway represent new therapeutics to combat cancer, inflammatory bowel disease (IBD) and CNS disease. Recently, the first-in-class SMI named SMI-10B that target OSM and block its interaction with receptor (OSMR) were reported. However, the binding pocket and interaction mode of the compound on OSM remain poorly understood, which hampering the rational design of SMIs that target OSM. Here, using SMI-10B as a probe, the multiple pockets on OSM for small molecules binding were extensively explored by unbiased molecular dynamics (MD) simulations. Then, the near-native structure of the complex was identified by molecular mechanics generalized Born surface area (MM/GBSA) binding energy funnel. Moreover, the binding stabilities of the protein-ligand complexes in near- and non-native conformations were verified by additional independent MD runs and absolute free energy perturbation (FEP) calculation. In summary, the unique feature of SMI-10B spontaneously binds to OSM characterized here not only provide detailed information for understanding the molecular mechanism of SMI-10B binding to OSM, but also will facilitate the rational design of novel and more potent SMIs to block OSM signaling.

Unconventional secretion of unglycosylated ORF8 is critical for the cytokine storm during SARS-CoV-2 infection

Coronavirus disease 2019 is a respiratory infectious disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Evidence on the pathogenesis of SARS-CoV-2 is accumulating rapidly. In addition to structural proteins such as Spike and Envelope, the functional roles of non-structural and accessory proteins in regulating viral life cycle and host immune responses remain to be understood. Here, we show that open reading frame 8 (ORF8) acts as messenger for inter-cellular communication between alveolar epithelial cells and macrophages during SARS-CoV-2 infection. Mechanistically, ORF8 is a secretory protein that can be secreted by infected epithelial cells via both conventional and unconventional secretory pathways. Conventionally secreted ORF8 is glycosylated and loses the ability to recognize interleukin 17 receptor A of macrophages, possibly due to the steric hindrance imposed by N-glycosylation at Asn78. However, unconventionally secreted ORF8 does not undergo glycosylation without experiencing the ER-Golgi trafficking, thereby activating the downstream NF-κB signaling pathway and facilitating a burst of cytokine release. Furthermore, we show that ORF8 deletion in SARS-CoV-2 attenuates inflammation and yields less lung lesions in hamsters. Our data collectively highlights a role of ORF8 protein in the development of cytokine storms during SARS-CoV-2 infection.

The Novel lncRNA RP9P Promotes Colorectal Cancer Progression by Modulating miR-133a-3p/FOXQ1 Axis

Background

The long non-coding RNA (lncRNA) RP9 pseudogene (RP9P) is a pseudogene-derived lncRNA that has never been reported in cancer, and its function underlying tumorigenesis in colorectal cancer (CRC) remains unknown.

Methods

RP9P and miR-133a-3p were filtered through bioinformatics analysis. The level of RP9P, miR-133a-3p, and FOXQ1 in CRC cell lines was detected by real-time PCR. Cell Counting Kit-8 and flow cytometric analyses were used to detect cell proliferation and apoptosis, respectively. Interactions between RP9P, miR-133a-3p, and FOXQ1 were confirmed by a dual-luciferase reporter assay.

Results

RP9P was overexpressed in CRC compared to normal control tissues and cells. Knockdown of RP9P inhibited CRC cell viability. RP9P directly interacted with miR-133a-3p, and miR-133a-3p downregulation abrogated the tumor-suppressing effect of RP9P knockdown. miR-133a-3p directly targeted FOXQ, which was positively regulated by RP9P. RP9P knockdown decreased FOXQ1 expression levels in CRC cells by directly targeting miR-133a-3p via a sponge mechanism. In addition, in vivo experiments in a xenograft model revealed that downregulated RP9P expression inhibited CRC cell tumorigenesis.

Conclusion

RP9P promotes colorectal cancer progression by regulating the miR-133a-3p/FOXQ1 axis.

Physicochemical and Biological Insights Into the Molecular Interactions Between Extracellular DNA and Exopolysaccharides in Myxococcus xanthus Biofilms

Extracellular DNA (eDNA) is a critical component in the extracellular matrix (ECM) of bacterial biofilms, while little is known about the mechanisms underlying how eDNA integrates into the ECM through potential macromolecular interactions. Myxococcus xanthus biofilm was employed as a suitable model for the investigation due to the co-distribution of eDNA and exopolysaccharides (EPS) owing to their direct interactions in the ECM. DNA is able to combine with M. xanthus EPS to form a macromolecular conjugate, which is dominated by the electrostatic forces participating in the polymer-polymer interactions. Without intercalation binding, DNA-EPS interactions exhibit a certain degree of reversibility. Acting as a strong extracellular framework during biofilm formation process, the eDNA-EPS complex not only facilitates the initial cell adhesion and subsequent establishment of ECM architecture, but also renders cells within biofilms stress resistances that are relevant to the survival of M. xanthus in some hostile environments. Furthermore, the EPS protects the conjugated DNA from the degradation by nucleic acid hydrolases, which leads to the continuous and stable existence of eDNA in the native ECM of M. xanthus biofilms. These results will shed light on developing prevention and treatment strategies against biofilm-related risks.

Increased odds ratio for erectile dysfunction in COVID-19 patients

PURPOSE

Erectile dysfunction and COVID-19 share similar risk factors, including vascular disruption of integrity, cytokine release, cardiovascular disease, diabetes and obesity. The aim of this study was to investigate the association between erectile dysfunction and COVID-19 patients.

METHODS

Odds ratio for erectile dysfunction in patients with a history of COVID-19 with and without comorbidities were calculated using a patients' registry platform i2b2. ICD-10 diagnoses codes were accessed for queries and data were analyzed using logistic regression.

RESULTS

Patients with COVID-19 were 3.3 times more likely to have erectile dysfunction with 95% CI (2.8, 3.8). The association became stronger with odds ratio 4.8 (95% CI (4.1, 5.7)) after adjusting for age groups. The odds ratio remained the same after adjusting for smoking status with 3.5 (95% CI (3.0, 4.1)). After adjusting for race, COVID-19 patients were 2.6 (95% CI (2.2, 3.1)) times more likely to have erectile dysfunction. The odds ratio were 1.6, 1.8, 1.9 and 2.3 after adjusting for respiratory disease, obesity, circulatory disease and diabetes, respectively.

CONCLUSION

COVID-19 and erectile dysfunction are strongly associated even after adjustment for known risk factors and demographics.

Molecular Mechanism for the Allosteric Inhibition of the Human Serotonin Transporter by Antidepressant Escitalopram

Human serotine transporter (hSERT) is one of the most influential drug targets, and its allosteric modulators (e.g., escitalopram) have emerged to be the next-generation medication for psychiatric disorders. However, the molecular mechanism underlying the allosteric modulation of hSERT is still elusive. Here, the simulation strategies of conventional (cMD) and steered (SMD) molecular dynamics were applied to investigate this molecular mechanism from distinct perspectives. First, cMD simulations revealed that escitalopram's binding to hSERT's allosteric site simultaneously enhanced its binding to the orthosteric site. Then, SMD simulation identified that the occupation of hSERT's allosteric site by escitalopram could also block its dissociation from the orthosteric site. Finally, by comparing the simulated structures of two hSERT-escitalopram complexes with and without allosteric modulation, a new conformational coupling between an extracellular (Arg104-Glu494) and an intracellular (Lys490-Glu494) salt bridge was identified. In summary, this study explored the mechanism underlying the allosteric modulation of hSERT by collectively applying two MD simulation strategies, which could facilitate our understanding of the allosteric modulations of not only hSERT but also other clinically important therapeutic targets.

Optimized, visible light-induced crosslinkable hybrid gelatin/hyaluronic acid scaffold promotes complete spinal cord injury repair

Severe microenvironmental changes after spinal cord injury (SCI) present serious challenges in neural regeneration and tissue repair. Gelatin (GL)- and hyaluronic acid (HA)-based hydrogels are attractive scaffolds because they are major components of the extracellular matrix and can provide a favorable adjustable microenvironment for neurogenesis and motor function recovery. In this study, three-dimensional hybrid GL/HA hydrogel scaffolds were prepared and optimized. The hybrid hydrogels could undergo in-situ gelation and fit the defects perfectly via visible light- induced crosslinking in the complete SCI rats. We found that the transplantation of the hybrid hydrogel scaffold significantly reduced the inflammatory responses and suppressed glial scar formation in an HA concentration-dependent manner. Moreover, the hybrid hydrogel with GL/HA ratios less than 8/2 effectively promoted endogenous neural stem cell migration and neurogenesis, as well as improved neuron maturation and axonal regeneration. The results showed locomotor function improved 60 days after transplantation, thus suggesting that GL/HA hydrogels can be considered as a promising scaffold for complete SCI repair.

Nucleocapsid mutations R203K/G204R increase the infectivity, fitness, and virulence of SARS-CoV-2

Previous work found that the co-occurring mutations R203K/G204R on the SARS-CoV-2 nucleocapsid (N) protein are increasing in frequency among emerging variants of concern or interest. Through a combination of in silico analyses, this study demonstrates that R203K/G204R are adaptive, while large-scale phylogenetic analyses indicate that R203K/G204R associate with the emergence of the high-transmissibility SARS-CoV-2 lineage B.1.1.7. Competition experiments suggest that the 203K/204R variants possess a replication advantage over the preceding R203/G204 variants, possibly related to ribonucleocapsid (RNP) assembly. Moreover, the 203K/204R virus shows increased infectivity in human lung cells and hamsters. Accordingly, we observe a positive association between increased COVID-19 severity and sample frequency of 203K/204R. Our work suggests that the 203K/204R mutations contribute to the increased transmission and virulence of select SARS-CoV-2 variants. In addition to mutations in the spike protein, mutations in the nucleocapsid protein are important for viral spreading during the pandemic.

SYNBIP: synthetic binding proteins for research, diagnosis and therapy

The success of protein engineering and design has extensively expanded the protein space, which presents a promising strategy for creating next-generation proteins of diverse functions. Among these proteins, the synthetic binding proteins (SBPs) are smaller, more stable, less immunogenic, and better of tissue penetration than others, which make the SBP-related data attracting extensive interest from worldwide scientists. However, no database has been developed to systematically provide the valuable information of SBPs yet. In this study, a database named ‘Synthetic Binding Proteins for Research, Diagnosis, and Therapy (SYNBIP)’ was thus introduced. This database is unique in (a) comprehensively describing thousands of SBPs from the perspectives of scaffolds, biophysical & functional properties, etc.; (b) panoramically illustrating the binding targets & the broad application of each SBP and (c) enabling a similarity search against the sequences of all SBPs and their binding targets. Since SBP is a human-made protein that has not been found in nature, the discovery of novel SBPs relied heavily on experimental protein engineering and could be greatly facilitated by in-silico studies (such as AI and computational modeling). Thus, the data provided in SYNBIP could lay a solid foundation for the future development of novel SBPs. The SYNBIP is accessible without login requirement at both official (https://idrblab.org/synbip/) and mirror (http://synbip.idrblab.net/) sites.

[Quantitative analysis of differential proteins in liver tissues of patients with non-alcoholic steatohepatitis using iTRAQ technology]

OBJECTIVE

To screen differentially expressed proteins (DSPs) in the liver tissues of patients with nonalcoholic steatohepatitis (NASH) using proteomic technologies to identify potential therapeutic targets of NASH.

METHODS

Liver tissue specimens were obtained from 3 patients with pathologically confirmed NASH and 3 normal control subjects. The total proteins were extracted from the specimens, and iTRAQ reagent was used to label the peptides for liquid chromatography tandem mass spectrometry (LC-MS/MS) detection. The DSPs were identified by comparing the data against UniProt protein database using Mascot2.3.02 software and were annotated and enriched using GO database; KEGG database was used for enrichment of the pathways involving these proteins. Real-time fluorescent quantitative PCR (qPCR) was performed to detect the mRNA expressions of the significant DSPs in NASH.

RESULTS

By the criteria that a DSP has >1.2 or < 0.8 fold difference between NASH group and the control group and with P < 0.05 as the threshold, a total of 648 significant DSPs in NASH were identified, including 246 up-regulated and 402 down-regulated proteins. GO functional enrichment analysis showed that the DSPs were involved mainly in small molecule metabolism, organic acid metabolism, oxygen acid metabolism and other biological processes, and were enriched in KEGG pathways including the metabolic pathways, complement coagulation cascades, and ribosomes. Among the 25 DEPs with a fold difference >2.0 or < 0.5 (P < 0.05), 6 proteins showed consistent results between qPCR verification and proteomic analysis, including 5 down-regulated proteins: Jumonji protein (JARID2), Lebasillinlike protein (LCA5L), synaptophysin 1 (SYN1) and collagen α-1 (XIII) chain (COL13A1), FYVE, RhoGEF and PH domain protein 5 (FGD5), and 1 upregulated protein glutathione S-transferase Mu 4 (GSTM4).

CONCLUSION

We identified 648 DEPs inthe liver tissue of patients NASH using iTRAQ technology and bioinformatics methods, and among them JARID2, SYN1, COL13A1, FGD5, and GSTM4 may serve as the key target proteins of NASH.

Structure-Based Discovery of Novel Nonpeptide Inhibitors Targeting SARS-CoV-2 Mpro

The continual spread of novel coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), posing a severe threat to the health worldwide. The main protease (Mpro, alias 3CLpro) of SARS-CoV-2 is a crucial enzyme for the maturation of viral particles and is a very attractive target for designing drugs to treat COVID-19. Here, we propose a multiple conformation-based virtual screening strategy to discover inhibitors that can target SARS-CoV-2 Mpro. Based on this strategy, nine Mpro structures and a protein mimetics library with 8960 commercially available compounds were prepared to carry out ensemble docking for the first time. Five of the nine structures are apo forms presented in different conformations, whereas the other four structures are holo forms complexed with different ligands. The surface plasmon resonance assay revealed that 6 out of 49 compounds had the ability to bind to SARS-CoV-2 Mpro. The fluorescence resonance energy transfer experiment showed that the biochemical half-maximal inhibitory concentration (IC50) values of the six compounds could hamper Mpro activities ranged from 0.69 ± 0.05 to 2.05 ± 0.92 μM. Evaluation of antiviral activity using the cell-based assay indicated that two compounds (Z1244904919 and Z1759961356) could strongly inhibit the cytopathic effect and reduce replication of the living virus in Vero E6 cells with the half-maximal effective concentrations (EC50) of 4.98 ± 1.83 and 8.52 ± 0.92 μM, respectively. The mechanism of the action for the two inhibitors were further elucidated at the molecular level by molecular dynamics simulation and subsequent binding free energy analysis. As a result, the discovered noncovalent reversible inhibitors with novel scaffolds are promising antiviral drug candidates, which may be used to develop the treatment of COVID-19.

ANPELA: analysis and performance assessment of the label-free quantification workflow for metaproteomic studies

Label-free quantification (LFQ) with a specific and sequentially integrated workflow of acquisition technique, quantification tool and processing method has emerged as the popular technique employed in metaproteomic research to provide a comprehensive landscape of the adaptive response of microbes to external stimuli and their interactions with other organisms or host cells. The performance of a specific LFQ workflow is highly dependent on the studied data. Hence, it is essential to discover the most appropriate one for a specific data set. However, it is challenging to perform such discovery due to the large number of possible workflows and the multifaceted nature of the evaluation criteria. Herein, a web server ANPELA (https://idrblab.org/anpela/) was developed and validated as the first tool enabling performance assessment of whole LFQ workflow (collective assessment by five well-established criteria with distinct underlying theories), and it enabled the identification of the optimal LFQ workflow(s) by a comprehensive performance ranking. ANPELA not only automatically detects the diverse formats of data generated by all quantification tools but also provides the most complete set of processing methods among the available web servers and stand-alone tools. Systematic validation using metaproteomic benchmarks revealed ANPELA's capabilities in 1 discovering well-performing workflow(s), (2) enabling assessment from multiple perspectives and (3) validating LFQ accuracy using spiked proteins. ANPELA has a unique ability to evaluate the performance of whole LFQ workflow and enables the discovery of the optimal LFQs by the comprehensive performance ranking of all 560 workflows. Therefore, it has great potential for applications in metaproteomic and other studies requiring LFQ techniques, as many features are shared among proteomic studies.

Copper (Cu2+) ion-induced misfolding of tau protein R3 peptide revealed by enhanced molecular dynamics simulation

Tau misfolding plays a significant role in some neurodegenerative diseases such as Alzheimer's disease (AD). It is intrinsically disordered and highly soluble under normal physiological conditions. While the protein will aggregate to form paired helical filaments (PHFs) under copper homeostasis at pathological conditions, which is the main substance of neurofibrillary tangles (NFTs) in the brain of AD patients. However, the molecular mechanism underlying the copper (Cu2+) ion-induced tau misfolding is not fully understood. In this study, using the 1/2 third repeat fragment (R3 peptide) of tau protein (residues 318-335: VTSKCGSLGNIHHKPGGG) as a model, a Gaussian accelerated molecular dynamics (GaMD) method followed by efficient trajectory analysis was carried out to investigate the influences of Cu2+ on the tau about the protein fold and the free energy landscape along the simulation. The two-dimensional potential of mean force (PMF) profiles obtained from reweighting of the GaMD simulations as well as clustering analysis revealed the Cu2+ ion induced α-helix fold R3 peptide located at the low-energy wells of free energy map, which is in agreement with the reported experimental result. In contrast, there is no α-helix fold of R3 peptide that appeared during the GaMD simulation without Cu2+ ion existing. Furthermore, the definition of secondary structure of protein (DSSP) analysis indicated that the R3 peptide with Cu2+ ion forms a stable structure of the helix (Lys321-His330 interval of the peptide) at between 400 and 500 ns. Therefore, the structures and free energy profiles from GaMD simulations proposed that Cu2+ triggers the aggregation of R3 peptide into toxic PHFs through a stable α-helix fold form.