Replica Exchange Molecular Dynamics Simulations Research Articles

Heat-induced structural and chemical changes to a computationally designed miniprotein.

The de novo design of miniprotein inhibitors has recently emerged as a new technology to create proteins that bind with high affinity to specific therapeutic targets. Their size, ease of expression, and apparent high stability makes them excellent candidates for a new class of protein drugs. However, beyond circular dichroism melts and hydrogen/deuterium exchange experiments, little is known about their dynamics, especially at the elevated temperatures they seemingly tolerate quite well. To address that and gain insight for future designs, we have focused on identifying unintended and previously overlooked heat-induced structural and chemical changes in a particularly stable model miniprotein, EHEE_rd2_0005. Nuclear magnetic resonance (NMR) studies suggest the presence of dynamics on multiple time and temperature scales. Transiently elevating the temperature results in spontaneous chemical deamidation visible in the NMR spectra, which we validate using both capillary electrophoresis and mass spectrometry (MS) experiments. High temperatures also result in greatly accelerated intrinsic rates of hydrogen exchange and signal loss in NMR heteronuclear single quantum coherence spectra from local unfolding. These losses are in excellent agreement with both room temperature hydrogen exchange experiments and hydrogen bond disruption in replica exchange molecular dynamics simulations. Our analysis reveals important principles for future miniprotein designs and the potential for high stability to result in long-lived alternate conformational states.

Effect of Pressure on the Conformational Landscape of Human γD-Crystallin from Replica Exchange Molecular Dynamics Simulations.

Human γD-crystallin belongs to a crucial family of proteins known as crystallins located in the fiber cells of the human lens. Since crystallins do not undergo any turnover after birth, they need to possess remarkable thermodynamic stability. However, their sporadic misfolding and aggregation, triggered by environmental perturbations or genetic mutations, constitute the molecular basis of cataracts, which is the primary cause of blindness in the globe according to the World Health Organization. Here, we investigate the impact of high pressure on the conformational landscape of wild-type HγD-crystallin using replica exchange molecular dynamics simulations augmented with principal component analysis. We find pressure to have a modest impact on global measures of protein stability, such as root-mean-square displacement and radius of gyration. Upon projecting our trajectories along the first two principal components from principal component analysis, however, we observe the emergence of distinct free energy basins at high pressures. By screening local order parameters previously shown or hypothesized as markers of HγD-crystallin stability, we establish correlations between a tyrosine-tyrosine aromatic contact within the N-terminal domain and the protein's end-to-end distance with projections along the first and second principal components, respectively. Furthermore, we observe the simultaneous contraction of the hydrophobic core and its intrusion by water molecules. This exploration sheds light on the intricate responses of HγD-crystallin to elevated pressures, offering insights into potential mechanisms underlying its stability and susceptibility to environmental perturbations, crucial for understanding cataract formation.

Computational investigation of the inhibitory interaction of IRF3 and SARS-CoV-2 accessory protein ORF3b

ORF3b is one of the SARS-CoV-2 accessory proteins. Previous experimental study suggested that ORF3b prevents IRF3 translocating to nucleus. However, the biophysical mechanism of ORF3b-IRF3 interaction is elusive. Here, we explored the conformation ensemble of ORF3b using all-atom replica exchange molecular dynamics simulation. Disordered ORF3b has mixed α-helix, β-turn and loop conformers. The potential ORF3b-IRF3 binding modes were searched by docking representative ORF3b conformers with IRF3, and 50 ORF3b-IRF3 complex poses were screened using molecular dynamics simulations ranging from 500 to 1000 ns. We found that ORF3b binds IRF3 predominantly on its CBP binding and phosphorylated pLxIS motifs, with CBP binding site has the highest binding affinity. The ORF3b-IRF3 binding residues are highly conserved in SARS-CoV-2. Our results provided biophysics insights into ORF3b-IRF3 interaction and explained its interferon antagonism mechanism.

In silico study on graphene quantum dots modified with various functional groups inhibiting α‑synuclein dimerization

Hypothesis: Graphene quantum dots (GQDs) with various functional groups are hypothesized to inhibit the α-synuclein (αS) dimerization, a crucial step in Parkinson’s disease pathogenesis. The potential of differently functionalized GQDs is systematically explored.Experiments: All-atom replica-exchange molecular dynamics simulations (accumulating to 75.6 μs) in explicit water were performed to study the dimerization of the αS non-amyloid component region and the influence of GQDs modified with various functional groups. Conformation ensemble, binding behavior, and free energy analysis were conducted.Findings: All studied GQDs inhibit β-sheet and backbone hydrogen bond formation in αS dimers, leading to looser oligomeric conformations. Charged GQDs severely impede the growth of extended β-sheets by providing extra contact surface. GQD binding primarily disrupts αS inter-peptide interactions through π-π stacking, CH-π interactions, and for charged GQDs, additionally through salt-bridge and hydrogen bonding interactions. GQD-COO– showed the most optimal inhibitory effect, binding mode, and intensity, which holds promise for the development of nanomedicines targeting amyloid aggregation in neurodegenerative diseases.

Molecular dynamics in multidimensional space explains how mutations affect the association path of neomycin to a riboswitch

Riboswitches are messenger RNA (mRNA) fragments binding specific small molecules to regulate gene expression. A synthetic N1 riboswitch, inserted into yeast mRNA controls the translation of a reporter gene in response to neomycin. However, its regulatory activity is sensitive to single-point RNA mutations, even those distant from the neomycin binding site. While the association paths of neomycin to N1 and its variants remain unknown, recent fluorescence kinetic experiments indicate a two-step process driven by conformational selection. This raises the question of which step is affected by mutations. To address this, we performed all-atom two-dimensional replica-exchange molecular dynamics simulations for N1 and U14C, U14C[Formula: see text], U15A, and A17G mutants, ensuring extensive conformational sampling of both RNA and neomycin. The obtained neomycin association and binding paths, along with multidimensional free-energy profiles, revealed a two-step binding mechanism, consisting of conformational selection and induced fit. Neomycin binds to a preformed N1 conformation upon identifying a stable upper stem and U-turn motif in the riboswitch hairpin. However, the positioning of neomycin in the binding site occurs at different RNA-neomycin distances for each mutant, which may explain their different regulatory activities. The subsequent induced fit arises from the interactions of the neomycin's N3 amino group with RNA, causing the G9 backbone to rearrange. In the A17G mutant, the critical C6-A17/G17 stacking forms at a closer RNA-neomycin distance compared to N1. These findings together with estimated binding free energies coincide with experiments and elucidate why the A17G mutation decreases and U15A enhances N1 activity in response to neomycin.

Dissecting the effect of ALS mutation S375G on the conformational properties and aggregation dynamics of TDP-43370-375 fragment

The aggregation of transactive response deoxyribonucleic acid (DNA) binding protein of 43 kDa (TDP-43) into ubiquitin-positive inclusions is closely associated with amyotrophic lateral sclerosis (ALS), frontotemporal lobar degeneration, and chronic traumatic encephalopathy. The 370–375 fragment of TDP-43 (370GNNSYS375, TDP-43370-375), the amyloidogenic hexapeptides, can be prone to forming pathogenic amyloid fibrils with the characteristic of steric zippers. Previous experiments reported the ALS-associated mutation, serine 375 substituted by glycine (S375G) is linked to early onset disease and protein aggregation of TDP-43. Based on this, it is necessary to explore the underlying molecular mechanisms. By utilizing all-atom molecular dynamics (MD) simulations of 102 μs in total, we investigated the impact of S375G mutation on the conformational ensembles and oligomerization dynamics of TDP-43370-375 peptides. Our replica exchange MD simulations show that S375G mutation could promote the unstructured conformation formation and induce peptides to form a loose packed oligomer, thus inhibiting the aggregation of TDP-43370-375. Further analyses suggest that S375G mutation displays a reduction effect on the number of total hydrogen bonds and contacts among TDP-43370-375 peptides. Hydrogen bonding and polar interactions among TDP-43370-375 peptides, as well as Y374-Y374 π-π stacking interaction, are attenuated by S375G mutation. Additional microsecond MD simulations demonstrate that S375G mutation could prohibit the conformational conversion to β-structure-rich aggregates and possess an inhibitory effect on the oligomerization dynamics of TDP-43370-375. This study offers for the first time of molecular insights into the S375G mutation affecting the aggregation of TDP-43370-375 at the atomic level, and may open new avenues in the development of future site-specific mutation therapeutics.

Effect of reaction media on hydrogenolysis of polyethylene plastic waste: Polymer-surface interactions in small alkane/polymer blends

The polymer reaction media and its properties can be altered by recycling a fraction of liquid products or adding alkane solvents. Less clear is whether this strategy affects hydrogenolysis. Herein, we investigated the effect of short-chain alkanes Cn consisting of n carbons (n=8, 16, and 32) on the upcycling of high-density polyethylene (HDPE) plastic waste to lubricant-range products over Ru/TiO2 catalysts by multiscale simulations and experiments. First, we trained a force field for polymer/surface interactions on a Ru22 nanoparticle (NP) supported on TiO2. Using replica exchange molecular dynamics simulations, we studied the effect of small hydrocarbons on the adsorption of a surrogate polymer, C142, on the catalyst. We found segregation of long chains (C142) at the catalyst surface due to the enthalpy gained by adsorbing more C-C bonds of the long chains, compensating for entropic losses upon adsorption. Short-chain molecules decrease the adsorbed carbons of long chains on the Ru NP due to blocking Ru active sites. Compared to the bulk chains, competitive adsorption results in a broader, heavy-tailed distribution of end-to-end distance of adsorbed chains. Our experiments demonstrated that catalyst activity declines significantly beyond simple dilution due to changes in polymer adsorption, and tuning the reaction media by creating suitable blends impacts hydrogenolysis. Density distributions for a 50:50%wt mixture of PP and PE show that PE chains are segregated at the surface, so they are prone to C-C bond breaking much faster than PP chains. H/D exchange experiments show preferential deuteration of PE, while CH3 groups of PP remain undeuterated. This may be explained by the preferential sorption of PE over PP, leading to specific distribution in the polymer blend.

Molecular Insights into the Differential Effects of Acetylation on the Aggregation of Tau Microtubule-Binding Repeats.

Aggregation of tau protein into intracellular fibrillary inclusions is characterized as the hallmark of tauopathies, including Alzheimer's disease and chronic traumatic encephalopathy. The microtubule-binding (MTB) domain of tau, containing either three or four repeats with sequence similarities, plays an important role in determining tau's aggregation. Previous studies have reported that abnormal acetylation of lysine residues displays a distinct effect on the formation of pathological tau aggregates. However, the underlying molecular mechanism remains mostly elusive. In this study, we performed extensive replica exchange molecular dynamics (REMD) simulations of 144 μs in total to systematically investigate the dimerization of four tau MTB repeats and explore the impacts of Lys280 (K280) or Lys321 (K321) acetylation on the conformational ensembles of the R2 or R3 dimer. Our results show that R3 is the most prone to aggregation among the four repeats, followed by R2 and R4, while R1 displays the weakest aggregation propensity with a disordered structure. Acetylation of K280 could promote the aggregation of R2 peptides by increasing the formation of β-sheet structures and strengthening the interchain interaction. However, K321 acetylation decreases the β-sheet content of the R3 dimer, reduces the ability of R3 peptides to form long β-strands, and promotes the stable helix structure formation. The salt bridge and Y310-Y310 π-π stacking interactions of the R3 dimer are greatly weakened by K321 acetylation, resulting in the inhibition of dimerization. This study uncovers the structural ensembles of tau MTB repeats and provides mechanistic insights into the influences of acetylation on tau aggregation, which may deepen the understanding of the pathogenesis of tauopathies.

Ab initio electron-lattice downfolding: Potential energy landscapes, anharmonicity, and molecular dynamics in charge density wave materials

The interplay of electronic and nuclear degrees of freedom presents an outstanding problem in condensed matter physics and chemistry. Computational challenges arise especially for large systems, long time scales, in nonequilibrium, or in systems with strong correlations. In this work, we show how downfolding approaches facilitate complexity reduction on the electronic side and thereby boost the simulation of electronic properties and nuclear motion-in particular molecular dynamics (MD) simulations. Three different downfolding strategies based on constraining, unscreening, and combinations thereof are benchmarked against full density functional calculations for selected charge density wave (CDW) systems, namely 1H-TaS_22, 1T-TiSe_22, 1H-NbS_22, and a one-dimensional carbon chain. We find that the downfolded models can reproduce potential energy surfaces on supercells accurately and facilitate computational speedup in MD simulations by about five orders of magnitude in comparison to purely ab initio calculations. For monolayer 1H-TaS_22 we report classical and path integral replica exchange MD simulations, revealing the impact of thermal and quantum fluctuations on the CDW transition.

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

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. 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.

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.

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