Replica Exchange Molecular Dynamics Simulations Research Articles

Solvent-Dependent Folding Behavior of a Helix-Forming Aryl-Triazole Foldamer.

Aromatic foldamers make up a novel class of bioinspired molecules that display helical conformations and have functions that rely on control over their coil-helix folding preferences. While the folding has been extensively examined by experiment, it has rarely been paired with the types of atomic level insights offered by theory. We present the results of all-atom molecular dynamics (MD) simulations to examine the role of solvent polarity on driving the helical folding behavior of the aryl-triazole foldamer. The temperature-dependent enhanced sampling technique, replica-exchange MD simulations, was employed to understand the folding phenomena. The simulation results show that in a low polarity solvent (dichloromethane), the foldamer prefers to stay in the unfolded state. The unfolded state has four dipolar triazoles (5 D) in their favored all-anti geometries and favoring anti-parallel geometries. However, an increase in solvent polarity using acetonitrile resulted in solvophobic collapse, yielding the helically folded form as the predominant state. The folded helix has an all-syn geometry, with triazoles in parallel arrangements. Intermediate conformations with a mixture of syn and anti arrangements of the triazoles are of lower abundance in both the DCM and MeCN solvents. The chiral handedness of the helix observed experimentally is assigned as left-handed by correlation with computed electronic circular dichroism spectra using time-dependent density functional theory.

Zn2+ Binding Increases Parallel Structure in the Aβ(16-22) Oligomer by Disrupting Salt Bridge in Antiparallel Structure.

The aggregation of monomeric amyloid β protein (Aβ) into oligomers and amyloid plaque in the brain is associated with Alzheimer's disease. The hydrophobic central core Aβ16-22 has been widely studied due to its essential role in the fibrillization of full-length Aβ peptides. Compared to the homogeneous antiparallel structure of Aβ16-22 at the late stage, the early-stage prefibrillar aggregates contain varying proportions of different β structures. In this work, we studied the appearance probabilities of various self-assembly structures of Aβ16-22 and the effects of Zn2+ on these probabilities by replica exchange molecular dynamics simulations. It was found that at room temperature, Aβ16-22 can readily form assembled β-sheet structures in pure water, where a typical antiparallel arrangement dominates (24.8% of all sampled trimer structures). The addition of Zn2+ to the Aβ16-22 solution will dramatically decrease the appearance probability of antiparallel trimer structures to 12.5% by disrupting the formation of the Lys16-Glu22 salt bridge. Meanwhile, the probabilities of hybrid antiparallel/parallel structures increase. Our simulation results not only reveal the competition between antiparallel and parallel structures in the Aβ16-22 oligomers but also show that Zn2+ can affect the oligomer structures. The results also provide insights into the role of metal ions in the self-assembly of short peptides.

The Inhibition Effect of Epigallocatechin-3-Gallate on the Co-Aggregation of Amyloid-β and Human Islet Amyloid Polypeptide Revealed by Replica Exchange Molecular Dynamics Simulations.

Alzheimer's disease and Type 2 diabetes are two epidemiologically linked diseases which are closely associated with the misfolding and aggregation of amyloid proteins amyloid-β (Aβ) and human islet amyloid polypeptide (hIAPP), respectively. The co-aggregation of the two amyloid proteins is regarded as the fundamental molecular mechanism underlying their pathological association. The green tea extract epigallocatechin-3-gallate (EGCG) has been extensively demonstrated to inhibit the amyloid aggregation of Aβ and hIAPP proteins. However, its potential role in amyloid co-aggregation has not been thoroughly investigated. In this study, we employed the enhanced-sampling replica exchange molecular dynamics simulation (REMD) method to investigate the effect of EGCG on the co-aggregation of Aβ and hIAPP. We found that EGCG molecules substantially diminish the β-sheet structures within the amyloid core regions of Aβ and hIAPP in their co-aggregates. Through hydrogen-bond, π-π and cation-π interactions targeting polar and aromatic residues of Aβ and hIAPP, EGCG effectively attenuates both inter-chain and intra-chain interactions within the co-aggregates. All these findings indicated that EGCG can effectively inhibit the co-aggregation of Aβ and hIAPP. Our study expands the potential applications of EGCG as an anti-amyloidosis agent and provides therapeutic options for the pathological association of amyloid misfolding disorders.

Effect of Histidine Behaviors on Structural Properties of Human Islet Amyloid Peptides

AbstractHistidine behaviors (including tautomeric behaviors and protonation behaviors) have been considered as an origin factor in protein folding and misfolding. However, such effects on human islet amyloid peptide (hIAPP) is still unclear, which are the key points in understanding hIAPP misfolding mechanisms. In current study, to investigate the structural properties of hIAPP peptide with histidine in (δ), (ϵ), and (p) states, three independent replica exchange molecular dynamics simulations were performed. Our results show that β‐sheet structure contents are very low no matter histidine on which state, while high α‐helix contents are obtained in A8‐L16/H18 and G24‐L27/S28 regions. Cluster data suggested that three hIAPP systems are prone to flexible structural features. Further analysis confirmed that the histidine related H‐bonding networks were detected in hIAPP(δ) (N14‐H18) and hIAPP(p) (N14‐H18 and F15‐H18). Moreover, based on overall hIAPP(δ), hIAPP(ϵ), and hIAPP(p) equilibration trajectories, we discovered both sheet and helix structural features in basin of free energy landscape. Our current study contributes to the understanding on conformational changes of hIAPP due to histidine behaviors, which could lead to new insights into histidine behaviors mechanisms.

Dissecting the geometric and hydrophobic constraints of stapled peptides.

Stapled peptides are a promising class of molecules with potential as highly specific probes of protein-protein interactions and as therapeutics. Hydrocarbon stapling affects the peptide properties through the interplay of two factors: enhancing the overall hydrophobicity and constraining the conformational flexibility. By constructing a series of virtual peptides, we study the role of each factor in modulating the structural properties of a hydrocarbon-stapled peptide PM2, which has been shown to enter cells, engage its target Mouse Double Minute 2 (MDM2), and activate p53. Hamiltonian replica exchange molecular dynamics (HREMD) simulations suggest that hydrocarbon stapling favors helical populations of PM2 through a combination of the geometric constraints and the enhanced hydrophobicity of the peptide. To further understand the conformational landscape of the stapled peptides along the binding pathway, we performed HREMD simulations by restraining the peptide at different distances from MDM2. When the peptide approaches MDM2, the binding pocket undergoes dehydration which appears to be greater in the presence of the stapled peptide compared with the linear peptide. In the binding pocket, the helicity of the stapled peptide is increased due to the favorable interactions between the peptide residues as well as the staple and the microenvironment of the binding pocket, contributing to enhanced affinity. The dissection of the multifaceted mechanism of hydrocarbon stapling into individual factors not only deepens fundamental understanding of peptide stapling, but also provides guidelines for the design of new stapled peptides.

Silver nanoparticles alter the dimerization of Aβ42 studied by REMD simulations.

The aggregation of amyloid beta (Aβ) peptides is associated with the development of Alzheimer's disease (AD). However, there has been a growing belief that the oligomerization of Aβ species in different environments has a neurotoxic effect on the patient's brain, causing damage. It is necessary to comprehend the compositions of Aβ oligomers in order to develop medications that may effectively inhibit these neurotoxic forms that affect the nervous system of AD patients. Thus, dissociation or inhibition of Aβ aggregation may be able to prevent AD. To date, the search for traditional agents and biomolecules has largely been unsuccessful. In this context, nanoparticles have emerged as potential candidates to directly inhibit the formation of Aβ oligomers. The oligomerization of the dimeric Aβ peptides with or without the influence of a silver nanoparticle was thus investigated using temperature replica-exchange molecular dynamics (REMD) simulations. The physical insights into the dimeric Aβ oligomerization were clarified by analyzing intermolecular contact maps, the free energy landscape of the dimeric oligomer, secondary structure terms, etc. The difference in obtained metrics between Aβ with or without a silver nanoparticle provides a picture of the influence of silver nanoparticles on the oligomerization process. The underlying mechanisms that are involved in altering Aβ oligomerization will be discussed. The obtained results may play an important role in searching for Aβ inhibitor pathways.

A theoretical study on the activity and selectivity of IDO/TDO inhibitors.

Indoleamine 2,3-dioxygenase 1 (IDO) is a tryptophan (Trp) metabolic enzyme along the kynurenine (NFK) pathway. Under pathological conditions, IDO overexpressed by tumor cells causes depletion of tryptophan and the accumulation of metabolic products, which inhibit the local immune response and form immune escape. Therefore, the suppression of IDO activity is one of the strategies for tumor immunotherapy, and drug design for this target has been the focus of research for more than two decades. Apart from IDO, tryptophan dioxygenase (TDO) of the same family can also catalyze the same biochemical reaction in the human body, but it has different tissue distribution and substrate selectivity from IDO. Based on the principle of drug design with high potency and low cross-reactivity to specific targets, in this subject, the activity and selectivity of IDO and TDO toward small molecular inhibitors were studied from the perspective of thermodynamics and kinetics. The aim was to elucidate the structural requirements for achieving favorable biological activity and selectivity of IDO and TDO inhibitors. Specifically, the interactions of inhibitors from eight families with IDO and TDO were initially investigated through molecular docking and molecular dynamics simulations, and the thermodynamic data for binding of inhibitors were predicted by the molecular mechanics/generalized Born surface area (MM/GBSA) method. Secondly, we explored the free energy landscape of JKloops, the kinetic control element of IDO/TDO, using temperature replica exchange molecular dynamics (T-REMD) simulations and elucidated the connection between the rules of IDO/TDO conformational changes and the inhibitor selectivity mechanism. Furthermore, the binding and dissociation processes of the C1 inhibitor (NLG919) were simulated by the adaptive steering molecular dynamics (ASMD) method, which not only addressed the possible stable, metastable, and transition states for C1 inhibitor-IDO/TDO interactions, but also accurately predicted kinetic data for C1 inhibitor binding and dissociation. In conclusion, we have constructed a complete process from enzyme (IDO/TDO) conformational activation to inhibitor binding/dissociation and used the thermodynamic and kinetic data of each link as clues to verify the control mechanism of IDO/TDO on inhibitor selectivity. This is of great significance for us to understand the design principles of tumor immunotherapy drugs and to avoid drug resistance of immunotherapy drugs.

Abnormal Phosphorylation at Ser293 promotes the Oligomerization of Tau Repeat R2: A Computational Study

AbstractBackgroundIn tauopathies such as Alzheimer’s disease (AD), aberrant phosphorylation causes the dissociation of tau protein from microtubules. The dissociated tau then aggregates into sequent forms from soluble oligomers to paired helical filaments (PHF), and insoluble neurofibrillary tangles (NFTs). NFTs is a hallmark of AD, while oligomers are found to be the most toxic form of the tau aggregates. Therefore, understanding tau oligomerization with regard of abnormal phosphorylation is important to therapeutic development of AD. In this study, we investigated the impact of phosphorylated Ser293, one of 40 aberrant phosphorylation sites of full‐length tau protein, on monomeric and dimeric structures of tau repeat R2 peptidesMethodWild‐type tau repeat R2 (wtR2) was taken from an experimental structure of R2‐microtubule complex (PDB code 6CVN). The phosphorylated Ser293 R2 peptide (pS293R2) was obtained by phosphorylating the wtR2 at Ser293 residue. Extensive replica exchange molecular dynamics simulations using all‐atom force field and explicit solvent were carried out for five systems including two monomers (wtR2 and pS293R2) and three dimers (wtR2+wtR2, wtR2+pS293R2, and pS293R2+ pS293R2). The structures of R2 monomers and dimers sampled by the simulations were characterized and compared.ResultFor monomeric conformations, the phosphorylation enhances intramolecular interaction, compactness and β content. For dimeric cases, the phosphorylation increases intermolecular interaction, β‐sheet formation, and turn‐coil structural transition. Interestingly, the phosphorylated residue strongly interacts with cations, resulting in a pS293‐cation‐pS293 bridge in pS293R2+ pS293R2 system.ConclusionOur results suggest that abnormal phosphorylation at Ser293 accelerates the oligomerization of R2 peptides, and consequently it may promote tau aggregation. Furthermore, we also found a pS293‐cation‐pS293 bridge formed, thus uncovered the essential role of cation ions in the oligomerization of the phosphorylated R2 repeat. Our findings suggest that phosphorylation at Ser293 should be taken into consideration when screening the inhibitors of tau oligomerization.

Development of Self-Assembled Biomimetic Nanoscale Collagen-like Peptide-Based Scaffolds for Tissue Engineering: An In Silico and Laboratory Study.

Development of biocomposite scaffolds has gained tremendous attention due to their potential for tissue regeneration. However, most scaffolds often contain animal-derived collagen that may elicit an immunological response, necessitating the development of new biomaterials. Herein, we developed a new collagen-like peptide,(Pro-Ala-His)10 (PAH)10, and explored its ability to be utilized as a functional biomaterial by incorporating it with a newly synthesized peptide-based self-assembled gel. The gel was prepared by conjugating a pectin derivative, galataric acid, with a pro-angiogenic peptide (LHYQDLLQLQY) and further functionalized with a cortistatin-derived peptide, (Phe-Trp-Lys-Thr)4 (FWKT)4, and the bio-ionic liquid choline acetate. The self-assembly of (PAH)10 and its interactions with the galactarate-peptide conjugates were examined using replica exchange molecular dynamics (REMD) simulations. Results revealed the formation of a multi-layered scaffold, with enhanced stability at higher temperatures. We then synthesized the scaffold and examined its physicochemical properties and its ability to integrate with aortic smooth muscle cells. The scaffold was further utilized as a bioink for bioprinting to form three-dimensional cell-scaffold matrices. Furthermore, the formation of actin filaments and elongated cell morphology was observed. These results indicate that the (PAH)10 hybrid scaffold provides a suitable environment for cell adhesion, proliferation and growth, making it a potentially valuable biomaterial for tissue engineering.

Structural insights into the co-aggregation of Aβ and tau amyloid core peptides: Revealing potential pathological heterooligomers by simulations

The self-aggregation of amyloid-β (Aβ) and tau proteins are closely implicated in Alzheimer's disease (AD). Recent evidence indicates that Aβ and tau proteins can cross-interact to form co-aggregates, which aggravates the development of AD. However, their transient heterooligomer conformations and co-aggregation molecular mechanisms are largely unknown. Herein, we utilize replica exchange molecular dynamics simulations to investigate the conformational ensembles formed by the central hydrophobic core of Aβ (Aβ16–22) and each of two fibril-nucleating core segments of tau (PHF6* and PHF6). Both PHF6 and PHF6* are found to co-aggregate with Aβ16–22 into β-sheet-rich heterooligomers. Intriguingly, PHF6 and Aβ16–22 peptides formed closed β-barrels, while PHF6* and Aβ16–22 formed open β-barrels, implying their distinct co-aggregation property. Compared to Aβ16–22-PHF6*, Aβ16–22-PHF6 heterooligomers have higher β-sheet content, and contain longer β-strands and larger β-sheets, indicative of stronger co-aggregation ability of PHF6 with Aβ16–22. Further analyses reveal that hydrophobic and π-π stacking interactions between Y310 of PHF6 and Aβ16–22 are crucial for the closed β-barrel/larger β-sheet formation in Aβ16–22-PHF6 heterooligomers. These results highlight the paramount importance of PHF6 fragment, particularly Y310 residue, as a potential target for inhibiting Aβ-tau co-aggregation, which could help for effective therapeutic design in mitigating Aβ-tau co-aggregation related amyloidogenesis.

Structural impairment of p53 C-terminal due to the effect of phosphorylation and acetylation: a study on the interdependence of PTM

The C-terminal of tumor suppressor protein p53 is intrinsically disordered while unbound. This particular segment often shows structural plasticity when bound to other binding partners. The disordered component undergoes a disordered to ordered transition upon recognition. Post-translational modifications (PTMs), namely phosphorylation and acetylation, significantly alter the structural motifs of the segment. Among the various types of PTMs, phosphorylation, and acetylation of p53 at both N- and C- terminals lead to stabilization and activation. It has been noted experimentally that phosphorylation often regulates (enhances or reduces) the acetylation at specific sites. The phosphorylation of Thr377 and Ser378 reduces the acetylation of Lys373 and Lys382. Mutations of Thr377 and Ser378 to neutral Ala enhance and phospho mimic Asp reduce the acetylation of Lys373 and Lys382. Simulations of several single-point and pair-wise mutated systems have been generated to compare how the presence or absence of phosphorylation favors or disfavors the acetylation by thermodynamic and conformational analysis. We are using implicit solvent replica exchange molecular dynamics simulations to get 200 ns well-converged conformational ensembles of each system. Different sets of systems having both single and double PTMs are simulated. The results admit the appreciable change in the secondary structural level upon specific PTM. Also, the residual structure of the unbound p53 with single-point PTM varies significantly with pair-wise modifications. These observations further shed light on the relationship between the interdependencies of the specific PTM sites and the secondary structural levels. Communicated by Ramaswamy H. Sarma

Formation and Effects of Upstream DNA-RNA Base Pairing in Telomerase.

Telomere elongation by telomerase consists of two types of translocation: duplex translocation during each repeat synthesis and template translocation at the end of repeat synthesis. Our replica exchange molecular dynamics simulations show that in addition to the Watson-Crick interactions in the active site, templating RNA can also form base pairs with the upstream regions of DNA, mostly with the second upstream DNA repeat with respect to the 3'-end. At the end of the repeat synthesis, dG10-P442 and dG11-N446 hydrogen bonds form. Then, active-site base pairs dissociate one by one, and the RNA bases reanneal with the complementary base on the upstream DNA repeat. For each dissociating base pair a new one forms, thus conserving the number of base pairs during template translocation.

Minimum Free-Energy Shapes of Ag Nanocrystals: Vacuum vs Solution.

We use two variants of replica-exchange molecular dynamics (MD) simulations, parallel tempering MD and partial replica exchange MD, to probe the minimum free-energy shapes of Ag nanocrystals containing 100-200 atoms in a vacuum, ethylene glycol (EG) solvent, and EG solvent with a PVP polymer containing 100 repeat units. Our simulations reveal a shape intermediate between a Dh and an Ih, a Dh-Ih, that has distinct structural signatures and magic sizes. We find several prominent features associated with entropy: pure FCC nanocrystals are less common than FCC crystals containing stacking faults, and crystals with the minimum potential energy are not always preferred over the range of relevant temperatures. The shapes of the nanocrystals in solution are influenced by the chemical identities of the solution-phase molecules. Comparing Ag nanocrystal shapes in EG to those in an EG+PVP solution, we find more icosahedra in EG and more decahedra in EG+PVP across all of the nanocrystal sizes probed in this study. At certain critical sizes, nanocrystal shapes can change dramatically with the addition and removal of a single atom or with a change in temperature at a fixed size. The information in our study could be useful in efforts to devise processing routes to achieve selective nanocrystal shapes.

Influence of Protonation on the Norepinephrine Inhibiting α-Synuclein 71-82 Oligomerization.

The pathogenesis of Parkinson's disease (PD) is closely linked to the massive presence of Lewy vesicles and Lewy axons in the cytoplasm of neurons, mainly consisting of α-synuclein (αS). Norepinephrine (NE), whose secretion can be increased by exercise, has been demonstrated to prevent the fibrillation of αS and to break down the mature αS fibrils. In this work, we focus on the influence of protonation on the inhibitory ability of NE by using amyloid core fragment αS71-82 as a template. All-atom replica-exchange molecular dynamics simulations (accumulating to 33.6 μs) in explicit water were performed to explore the inhibitory effect of protonated and nonprotonated NE on αS oligomerization. Our results show that NE/NE+ can lead to a significant decrease in β-sheet content with increasing temperature, while isolated αS maintains relatively higher β-sheet conformations until 363 K, implying that both NE and NE+ can lower the critical temperature required for αS fibril decomposition. NE and NE+ also lead to the formation of less compact αS oligomers by preventing the backbone hydrogen bonds and the side-chain packing. The protonation would affect the binding affinity, interaction modes, and binding intensity of NE with αS. Interesting, NE and NE+ have a distinct binding free energy in the electrostatic and solvation terms, which mostly counter each other and produce a weak binding intensity with αS. Our work contributes to a better understanding of the inhibitory mechanism of NE and NE+ on αS oligomerization relevant to PD pathogenesis, which may provide clues for the design of antiamyloid medicine.

Characterizing the transmembrane domains of ADAM10 and BACE1 and the impact of membrane composition

The β-secretase, BACE1, and the α-secretase, ADAM10, are known to competitively cleave amyloid precursor protein (APP) in the amyloid cascades of Alzheimer’s disease. Cleavage of APP by BACE1 produces a 99-residue C-terminal peptide (APP-C99) that is subsequently cleaved by γ-secretase to form amyloid-β (Aβ) protein, whereas cleavage of APP by ADAM10 is nonamyloidogenic. It has been speculated that ADAM10/APP and BACE1/APP interactions are regulated by colocalization within and outside of liquid-ordered membrane domains; however, the mechanism of this regulation and the character of the proteins’ transmembrane domains are not well understood. In this work, we have developed and characterized minimal congener sequences for the transmembrane domains of ADAM10 and BACE1 using a multiscale modeling approach combining both temperature replica exchange and conventional molecular dynamics simulations based on the coarse-grained Martini2.2 and all-atom CHARMM36 force fields. Our results show that membrane composition impacts the character of the transmembrane domains of BACE1 and ADAM10, adding credence to the speculation that membrane domains are involved in the etiology of Alzheimer’s disease.

Efficient and accurate calculation of proline cis/trans isomerization free energies from Hamiltonian replica exchange molecular dynamics simulations

Proline cis/trans isomerization plays an important role in many biological processes but occurs on time scales not accessible to brute-force molecular dynamics (MD) simulations. We have designed a new Hamiltonian replica exchange scheme, ω-bias potential replica exchange molecular dynamics (ωBP-REMD), to efficiently and accurately calculate proline cis/trans isomerization free energies. ωBP-REMD is applied to various proline-containing tripeptides and a biologically important proline residue in the N2-domain of the gene-3-protein of phage fd in the wildtype and mutant variants of the protein. Excellent cis/trans transition rates are obtained. Reweighting of the sampled probability distribution along the peptide bond dihedral angle allows construction of the corresponding free-energy profile and calculation of the cis/trans isomerization free energy with high statistical precision. Very good agreement with experimental data is obtained. ωBP-REMD outperforms standard umbrella sampling in terms of convergence and agreement with experiment and strongly reduces perturbation of the local structure near the proline residue.

Structural Space of the Duffy Antigen/Receptor for Chemokines' Intrinsically Disordered Ectodomain 1 Explored by Temperature Replica-Exchange Molecular Dynamics Simulations.

Plasmodium vivax malaria affects 14 million people each year. Its invasion requires interactions between the parasitic Duffy-binding protein (PvDBP) and the N-terminal extracellular domain (ECD1) of the host's Duffy antigen/receptor for chemokines (DARC). ECD1 is highly flexible and intrinsically disordered, therefore it can adopt different conformations. We computationally modeled the challenging ECD1 local structure. With T-REMD simulations, we sampled its dynamic behavior and collected its most representative conformations. Our results suggest that most of the DARC ECD1 domain remains in a disordered state during the simulated time. Globular local conformations are found in the analyzed local free-energy minima. These globular conformations share an α-helix spanning residues Ser18 to Ser29 and in many cases they comprise an antiparallel β-sheet, whose β-strands are formed around residues Leu10 and Ala49. The formation of a parallel β-sheet is almost negligible. So far, progress in understanding the mechanisms forming the basis of the P. vivax malaria infection of reticulocytes has been hampered by experimental difficulties, along with a lack of DARC structural information. Our collection of the most probable ECD1 structural conformations will help to advance modeling of the DARC structure and to explore DARC-ECD1 interactions with a range of physiological and pathological ligands.

Deciphering the Inhibitory Mechanism of Naphthoquinone-Dopamine on the Aggregation of Tau Core Fragments PHF6* and PHF6.

The formation of neurofibrillary tangles by abnormal aggregation of tau protein is considered to be an important pathological characteristic of tauopathies, including Alzheimer's disease and chronic traumatic encephalopathy. Two hexapeptides 275VQIINK280 and 306VQIVYK311 in the microtubule binding region, named PHF6* and PHF6, are known to be aggregation-prone and responsible for tau fibrillization. Previous experiments reported that naphthoquinone-dopamine (NQDA) could effectively inhibit the aggregation of PHF6* and PHF6 and disrupt the fibrillar aggregates into nontoxic species, displaying a dual effect on the amyloid aggregation. However, the underlying molecular mechanism remains mostly elusive. Herein, we performed all-atom molecular dynamics (MD) simulations for 114 μs in total to systematically investigate the impacts of NQDA on the oligomerization of PHF6* and PHF6. The conformational ensembles of PHF6* and PHF6 peptides generated by replica exchange MD simulations show that NQDA could effectively prevent the hydrogen bond formation, reduce the ability of peptides to self-assemble into long β-strand and large β-sheets, and induce peptides to form a loosely packed and coil-rich oligomer. The interaction analysis shows that the binding of NQDA to PHF6* is mainly through hydrophobic interactions with residue I277 and hydrogen bonding interactions with Q276; for the PHF6 peptides, NQDA displays a strong π-π stacking interaction with residue Y310, thus impeding the Y310-Y310 π-π stacking and I308-Y310 CH-π interactions. The DA group of NQDA displays a stronger cation-π interaction than the NQ group, while the NQ group exhibits a stronger π-π stacking interaction. MD simulations demonstrate that NQDA prevents the conformational conversion to β-sheet-rich aggregates and displays an inhibitory effect on the oligomerization dynamics of PHF6* and PHF6. Our results provide a complete picture of inhibitory mechanisms of NQDA on PHF6* and PHF6 oligomerization, which may pave the way for designing drug candidates for the treatment of tauopathies.

Binding of viral nuclear localization signal peptides to importin-α nuclear transport protein

Using all-atom replica-exchange molecular dynamics simulations, we mapped the mechanisms of binding of the nuclear localization signal (NLS) sequence from Venezuelan equine encephalitis virus (VEEV) capsid protein to importin-α (impα) transport protein. Our objective was to identify the VEEV NLS sequence fragment that confers native, experimentally resolved binding to impα as well as to study associated binding energetics and conformational ensembles. The two selected VEEV NLS peptide fragments, KKPK and KKPKKE, show strikingly different binding mechanisms. The minNLS peptide KKPK binds non-natively and nonspecifically by adopting five diverse conformational clusters with low similarity to the x-ray structure 3VE6 of NLS-impα complex. Despite the prevalence of non-native interactions, the minNLS peptide still largely binds to the impα major NLS binding site. In contrast, the coreNLS peptide KKPKKE binds specifically and natively, adopting a largely homogeneous binding ensemble with a dominant, highly native-like conformational cluster. The coreNLS peptide retains most of native binding interactions, including π-cation contacts and a tryptophan cage. While KKPK binding is governed by a complex multistate free energy landscape featuring transitions between multiple binding poses, the coreNLS peptide free energy map is simple, exhibiting a single dominant native-like bound basin. We argue that the origin of the coreNLS peptide binding specificity is several electrostatic interactions formed by the two C-terminal amino acids, Lys10 and Glu11, with impα. The coreNLS sequence is then sufficient for native binding, but none of the amino acids flanking minNLS, including Lys10 and Glu11, are strictly necessary for the native pose. Our analyses indicate that the VEEV coreNLS sequence is virtually unique among human and viral proteins interacting with impα making it a potential target for VEEV-specific inhibitors.

De Novo Design of a β-Helix Tau Protein Scaffold: An Oligomer-Selective Vaccine Immunogen Candidate for Alzheimer's Disease.

Tau pathology is associated with many neurodegenerative disorders, including Alzheimer's disease (AD), where the spatio-temporal pattern of tau neurofibrillary tangles strongly correlates with disease progression, which motivates therapeutics selective for misfolded tau. Here, we introduce a new avidity-enhanced, multi-epitope approach for protein-misfolding immunogen design, which is predicted to mimic the conformational state of an exposed epitope in toxic tau oligomers. A predicted oligomer-selective tau epitope 343KLDFK347 was scaffolded by designing a β-helix structure that incorporated multiple instances of the 16-residue tau fragment 339VKSEKLDFKDRVQSKI354. Large-scale conformational ensemble analyses involving Jensen-Shannon Divergence and the embedding depth showed that the multi-epitope scaffolding approach, employed in designing the β-helix scaffold, was predicted to better discriminate toxic tau oligomers than other "monovalent" strategies utilizing a single instance of an epitope for vaccine immunogen design. Using Rosetta, 10,000 sequences were designed and screened for the linker portions of the β-helix scaffold, along with a C-terminal stabilizing α-helix that interacts with the linkers, to optimize the folded structure and stability of the scaffold. Structures were ranked by energy, and the lowest 1% (82 unique sequences) were verified using AlphaFold. Several selection criteria involving AlphaFold are implemented to obtain a lead-designed sequence. The structure was further predicted to have free energetic stability by using Hamiltonian replica exchange molecular dynamics (MD) simulations. The synthesized β-helix scaffold showed direct binding in surface plasmon resonance (SPR) experiments to several antibodies that were raised to the structured epitope using a designed cyclic peptide. Moreover, the strength of binding of these antibodies to in vitro tau oligomers correlated with the strength of binding to the β-helix construct, suggesting that the construct presents an oligomer-like conformation and may thus constitute an effective oligomer-selective immunogen.