Comparative analysis of thermal unfolding simulations of RNA recognition motifs (RRMs) of TAR DNA-binding protein 43 (TDP-43)

Computational insights into dynamics and conformational stability of N-acetylmannosamine kinase mutations

The activity of UDP-GlcNAc 2-epimerase/ManNAc kinase (GNE) is essential for the biosynthesis of sialic acid, which is involved in cellular processes in health and diseases. GNE contains an N-terminal epimerase domain and a C-terminal kinase domain (N-acetylmannosamine kinase, MNK). Mutations of the GNE protein led to hypoactivity of the enzyme and cause sialurea or autosomal recessive inclusion body myopathy/Nonaka myopathy. Here, we used all-atom molecular dynamics (MD) simulations to comprehend the folding, dynamics and conformational stability of MNK variants, including the wild type (WT) and three mutants (H677R, V696M and H677R/V696M). The deleterious and destabilizing nature of MNK mutants were predicted using different prediction tools. Results predicted that mutations modulate the stability, flexibility and function of MNK. The effect of mutations on the conformational stability and dynamics of MNK was next studied through the free-energy landscape (FEL), hydrogen-bonds and secondary structure changes. The FEL results show that the mutations interfere with various conformational transitions in both WT and mutants, exposing the structural underpinnings of protein destabilization and unfolding brought on by mutation. We discover that, when compared to the other two mutations, V696M and H677R/V696M, H677R has the most harmful effects. These findings have a strong correlation with published experimental studies that demonstrate how these mutations disrupt MNK activity. Hence, this computational study describes the structural details to unravel the mutant effects at the atomistic resolution and has implications for understanding the GNE's physiological and pathological role.Communicated by Ramaswamy H. Sarma.

Comprehensive mapping of mutations in TDP-43 and α-Synuclein that affect stability and binding

Abnormal aggregation and amyloid inclusions of TAR DNA-binding protein 43 (TDP-43) and α-Synuclein (α-Syn) are frequently co-observed in amyotrophic lateral sclerosis, Parkinson's disease, and Alzheimer's disease. Several reports showed TDP-43 C-terminal domain (CTD) and α-Syn interact with each other and the aggregates of these two proteins colocalized together in different cellular and animal models. Molecular dynamics simulation was conducted to elucidate the stability of the TDP-43 and Syn complex structure. The interfacial mutations in protein complexes changes the stability and binding affinity of the protein that may cause diseases. Here, we have utilized the computational saturation mutagenesis approach including structure-based stability and binding energy calculations to compute the systemic effects of missense mutations of TDP-43 CTD and α-Syn on protein stability and binding affinity. Most of the interfacial mutations of CTD and α-Syn were found to destabilize the protein and reduced the protein binding affinity. The results thus shed light on the functional consequences of missense mutations observed in TDP-43 associated proteinopathies and may provide the mechanisms of co-morbidities involving these two proteins.Communicated by Ramaswamy H. Sarma.

In silico evaluation of the inhibitory potential of nucleocapsid inhibitors of SARS-CoV-2: a binding and energetic perspective

The COVID-19 outbreak brought on by the SARS-CoV-2 virus continued to infect a sizable population worldwide. The SARS-CoV-2 nucleocapsid (N) protein is the most conserved RNA-binding structural protein and is a desirable target because of its involvement in viral transcription and replication. Based on this aspect, this study focused to repurpose antiviral compounds approved or in development for treating COVID-19. The inhibitors chosen are either FDA-approved or are currently being studied in clinical trials against COVID-19. Initially, they were designed to target stress granules and other RNA biology. We have utilized structure-based molecular docking and all-atom molecular dynamics (MD) simulation approach to investigate in detail the binding energy and binding modes of the different anti-N inhibitors to N protein. The result showed that five drugs including Silmitasterib, Ninetanidinb, Ternatin, Luteolin, Fedratinib, PJ34, and Zotatafin were found interacting with RNA binding sites as well as to predicted protein interface with higher binding energy. Overall, drug binding increases the stability of the complex with maximum stability found in the order, Silmitasertib > PJ34 > Zotatatafin. In addition, the frustration changes due to drug binding brings a decrease in local frustration and this decrease is mainly observed in α-helix, β3, β5, and β6 strands and are important for drug binding. Our in-silico data suggest that an effective interaction occurs for some of the tested drugs and prompt their further validation to reduce the rapid outspreading of SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Dissecting how ALS-associated D290V mutation enhances pathogenic aggregation of hnRNPA2286-291 peptides: Dynamics and conformational ensembles

The aggregation of RNA binding proteins, including hnRNPA1/2, TDP-43 and FUS, is heavily implicated in causing or increasing disease risk for a series of neurodegenerative diseases such as amyotrophic lateral sclerosis (ALS). A recent experimental study demonstrated that an ALS-related D290V mutation in the low complexity domain (LCD) of hnRNPA2 can enhance the aggregation propensity of wild type (WT) hnRNPA2_286-291 peptide. However, the underlying molecular mechanisms remain elusive. Herein, we investigated effects of D290V mutation on aggregation dynamics of hnRNPA2_286-291 peptide and the conformational ensemble of hnRNPA2_286-291 oligomers by performing all-atom molecular dynamic and replica-exchange molecular dynamic simulations. Our simulations demonstrate that D290V mutation greatly reduces the dynamics of hnRNPA2_286-291 peptide and that D290V oligomers possess higher compactness and β-sheet content than WT, indicative of mutation-enhanced aggregation capability. Specifically, D290V mutation strengthens inter-peptide hydrophobic, main-chain hydrogen bonding and side-chain aromatic stacking interactions. Those interactions collectively lead to the enhancement of aggregation capability of hnRNPA2_286-291 peptides. Overall, our study provides insights into the dynamics and thermodynamic mechanisms underlying D290V-induced disease-causing aggregation of hnRNPA2_286-291, which could contribute to better understanding of the transitions from reversible condensates to irreversible pathogenic aggregates of hnRNPA2 LCD in ALS-related diseases.

ALS-Linked A315T and A315E Mutations Enhance β-Barrel Formation of the TDP-43307-319 Hexamer: A REST2 Simulation Study

Pathogenic mutations of transactivation response element DNA-binding protein 43 (TDP-43) are closely linked with amyotrophic lateral sclerosis (ALS). It was recently reported that two ALS-linked familial mutants A315T and A315E of TDP-43_307-319 peptides can self-assemble into oligomers including tetramers, hexamers, and octamers, among which hexamers were suggested to form the β-barrel structure. However, due to the transient nature of oligomers, their conformational properties and the atomic mechanisms underlying the β-barrel formation remain largely elusive. Herein, we investigated the hexameric conformational distributions of the wild-type (WT) TDP-43_307-319 fragment and its A315T and A315E mutants by performing all-atom explicit-solvent replica exchange with solute tempering 2 simulations. Our simulations reveal that each peptide can self-assemble into diverse conformations including ordered β-barrels, bilayer β-sheets and/or monolayer β-sheets, and disordered complexes. A315T and A315E mutants display higher propensity to form β-barrel structures than the WT, which provides atomic explanation for their enhanced neurotoxicity reported previously. Detailed interaction analysis shows that A315T and A315E mutations increase inter-molecular interactions. Also, the β-barrel structures formed by the three different peptides are stabilized by distinct inter-peptide side-chain hydrogen bonding, hydrophobic, and aromatic stacking interactions. This study demonstrates the enhanced β-barrel formation of the TDP-43_307-319 hexamer by the pathogenic A315T and A315E mutations and reveals the underlying molecular determinants, which may be helpful for in-depth understanding of the ALS-mutation-induced neurotoxicity of TDP-43 protein.

ALS-associated A315E and A315pT variants exhibit distinct mechanisms in inducing irreversible aggregation of TDP-43312-317 peptides

Amyotrophic lateral sclerosis (ALS) is intensively associated with insoluble aggregates formed by transactivation response element DNA-binding protein 43 (TDP-43) in the cytoplasm of neuron cells. A recent experimental study reported that two ALS-linked familial variants, A315E and A315pT (pT, phosphorylated threonine), can induce irreversible aggregation of the TDP-43 ³¹²NFGAFS³¹⁷ segment (TDP-43_312-317). However, the underlying molecular mechanism remains largely elusive. Here, we investigated the early aggregation process of the wild type (WT) ³¹²NFGAFS³¹⁷ segment and its A315E and A315pT variants by performing multiple microsecond all-atom molecular dynamics simulations. Our simulations show that the two variants display lower fluidity than WT, consistent with their decreased labilities observed in previous denaturation assay experiments. Despite each of the two variants carrying one negative charge, unexpectedly, we find that both A315E mutation and A315pT phosphorylation enhance intermolecular interactions and result in the formation of more compact oligomers. Compared to WT, A315E oligomers possess low β-sheet content but a compact hydrophobic core, while A315pT oligomers have high β-sheet content and large β-sheets. Side chain hydrogen-bonding and hydrophobic interactions as well as N312-E315 salt bridges contribute most to the increased aggregation propensity of the A315E mutant. By contrast, main chain and side chain hydrogen-bonding interactions, side chain hydrophobic and aromatic interactions, are crucial to the enhanced aggregation capability of the A315pT variant. These results indicate that glutamate mutation and phosphorylation at position 315 induce the irreversible aggregation of TDP-43_312-317 peptides through differential mechanisms, which remind us that we should be careful in the investigation of the phosphorylation effect on protein aggregation by using phosphomimetic substitutions. This study provides mechanistic insights into the A315E/A315pT-induced irreversible aggregation of TDP-43_312-317, which may be helpful for the in-depth understanding of ALS-mutation/phosphorylation-associated liquid-to-solid phase transition of TDP-43 protein aggregates.

Novel insights into the structural changes induced by disease-associated mutations in TDP-43: a computational approach

Over the last two decades, the pathogenic aggregation of TAR DNA-binding protein 43 (TDP-43) is found to be strongly associated with several fatal neurodegenerative diseases such as amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTD), etc. While the mutations and truncation in TDP-43 protein have been suggested to be responsible for TDP-43 pathogenesis by accelerating the aggregation process, the effects of these mutations on the bio-mechanism of pathological TDP-43 protein remained poorly understood. Investigating this at the molecular level, we formulized an integrated workflow of molecular dynamic simulation and machine learning models (MD-ML). By performing an extensive structural analysis of three disease-related mutations (i.e., I168A, D169G, and I168A-D169G) in the conserved RNA recognition motifs (RRM1) of TDP-43, we observed that the I168A-D169G double mutant delineates the highest packing of the protein inner core as compared to the other mutations, which may indicate more stability and higher chances of pathogenesis. Moreover, through our MD-ML workflow, we identified the biological descriptors of TDP-43 which includes the interacting residue pairs and individual protein residues that influence the stability of the protein and could be experimentally evaluated to develop potential therapeutic strategies.Communicated by Ramaswamy H. Sarma.

Computational Insights of Unfolding of N-Terminal Domain of TDP-43 Reveal the Conformational Heterogeneity in the Unfolding Pathway

TDP-43 proteinopathies is a disease hallmark that characterizes amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD). The N-terminal domain of TDP-43 (NTD) is important to both TDP-43 physiology and TDP-43 proteinopathy. However, its folding and dimerization process is still poorly characterized. In the present study, we have investigated the folding/unfolding of NTD employing all-atom molecular dynamics (MD) simulations in 8 M dimethylsulfoxide (DMSO) at high temperatures. The MD results showed that the unfolding of the NTD at high temperature evolves through the formation of a number of conformational states differing in their stability and free energy. The presence of structurally heterogeneous population of intermediate ensembles was further characterized by the different extents of solvent exposure of Trp80 during unfolding. We suggest that these non-natives unfolded intermediate ensembles may facilitate NTD oligomerization and subsequently TDP-43 oligomerization, which might lead to the formation of irreversible pathological aggregates, characteristics of disease pathogenesis.

Genomics-guided targeting of stress granule proteins G3BP1/2 to inhibit SARS-CoV-2 propagation

SARS-CoV-2 nucleocapsid (N) protein undergoes RNA-induced phase separation (LLPS) and sequesters the host key stress granule (SG) proteins, Ras-GTPase-activating protein SH3-domain-binding protein 1 and 2 (G3BP1 and G3BP2) to inhibit SG formation. This will allow viral packaging and propagation in host cells. Based on a genomic-guided meta-analysis, here we identify upstream regulatory elements modulating the expression of G3BP1 and G3BP2 (collectively called G3BP1/2). Using this strategy, we have identified FOXA1, YY1, SYK, E2F-1, and TGFBR2 as activators and SIN3A, SRF, and AKT-1 as repressors of G3BP1/2 genes. Panels of the activators and repressors were then used to identify drugs that change their gene expression signatures. Two drugs, imatinib, and decitabine have been identified as putative modulators of G3BP1/2 genes and their regulators, suggesting their role as COVID-19 mitigation agents. Molecular docking analysis suggests that both drugs bind to G3BP1/2 with a much higher affinity than the SARS-CoV-2 N protein. This study reports imatinib and decitabine as candidate drugs against N protein and G3BP1/2 protein.

Dynamics of binding interactions of TDP-43 and RNA: An equally weighted multiscale elastic network model study

Transactive response DNA binding protein 43 (TDP-43), an alternative-splicing regulator, can specifically bind long UG-rich RNAs, associated with a range of neurodegenerative diseases. Upon binding RNA, TDP-43 undergoes a large conformational change with two RNA recognition motifs (RRMs) connected by a long linker rearranged, strengthening the binding affinity of TDP-43 with RNA. We extend the equally weighted multiscale elastic network model (ewmENM), including its Gaussian network model (ewmGNM) and Anisotropic network model (ewmANM), with the multiscale effect of interactions considered, to the characterization of the dynamics of binding interactions of TDP-43 and RNA. The results reveal upon RNA binding a loss of flexibility occurs to TDP-43's loop3 segments rich in positively charged residues and C-terminal of high flexibility, suggesting their anchoring RNA, induced fit and conformational adjustment roles in recognizing RNA. Additionally, based on movement coupling analyses, it is found that RNA binding strengthens the interactions among intra-RRM β-sheets and between RRMs partially through the linker's mediating role, which stabilizes RNA binding interface, facilitating RNA binding efficiency. In addition, utilizing our proposed thermodynamic cycle method combined with ewmGNM, we identify the key residues for RNA binding whose perturbations induce a large change in binding free energy. We identify not only the residues important for specific binding, but also the ones critical for the conformational rearrangement between RRMs. Furthermore, molecular dynamics simulations are also performed to validate and further interpret the ENM-based results. The study demonstrates a useful avenue to utilize ewmENM to investigate the protein-RNA interaction dynamics characteristics.

Insights into the aggregation mechanism of RNA recognition motif domains in TDP-43: a theoretical exploration

The transactive response DNA-binding protein 43 (TDP-43) is associated with several diseases such as amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD) due to pathogenic aggregations. In this work, we examined the dimer, tetramer and hexamer models built from the RRM domains of TDP-43 using molecular dynamics simulations in combination with the protein-protein docking. Our results showed that the formations of the dimer models are mainly achieved by the interactions of the RRM1 domains. The parallel β-sheet layers between the RRM1 domains provide most of the binding sites in these oligomer models, and thus play an important role in the aggregation process. The approaching of the parallel β-sheet layers from small oligomer models gradually expand to large ones through the allosteric communication between the α1/α2 helices of the RRM1 domains, which maintains the binding affinities and interactions in the larger oligomer models. Using the repeatable-superimposing method based on the tetramer models, we proposed a new aggregation mechanism of RRM domains in TDP-43, which could well characterize the formation of the large aggregation models with the repeated, helical and rope-like structures. These new insights help to understand the amyloid-like aggregation phenomena of TDP-43 protein in ALS and FTLD diseases.

Identification of natural compounds as potent inhibitors of SARS-CoV-2 main protease using combined docking and molecular dynamics simulations

Coronavirus disease 2019 (COVID-19) has emerged from China and globally affected the entire population through the human-to-human transmission of a newly emerged virus called severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The genome of SARS-CoV-2 encodes several proteins that are essential for multiplication and pathogenesis. The main protease (Mpro or 3CLpro) of SARS-CoV-2 plays a central role in its pathogenesis and thus is considered as an attractive drug target for the drug design and development of small-molecule inhibitors. We have employed an extensive structure-based high-throughput virtual screening to discover potential natural compounds from the ZINC database which could inhibit the Mpro of SARS-CoV-2. Initially, the hits were selected on the basis of their physicochemical and drug-like properties. Subsequently, the PAINS filter, estimation of binding affinities using molecular docking, and interaction analyses were performed to find safe and potential inhibitors of SARS-CoV-2 Mpro. We have identified ZINC02123811 (1-(3-(2,5,9-trimethyl-7-oxo-3-phenyl-7H-furo[3,2-g]chromen-6-yl)propanoyl)piperidine-4-carboxamide), a natural compound bearing appreciable affinity, efficiency, and specificity towards the binding pocket of SARS-CoV-2 Mpro. The identified compound showed a set of drug-like properties and preferentially binds to the active site of SARS-CoV-2 Mpro. All-atom molecular dynamics (MD) simulations were performed to evaluate the conformational dynamics, stability and interaction mechanism of Mpro with ZINC02123811. MD simulation results indicated that Mpro with ZINC02123811 forms a stable complex throughout the trajectory of 100 ns. These findings suggest that ZINC02123811 may be further exploited as a promising scaffold for the development of potential inhibitors of SARS-CoV-2 Mpro to address COVID-19.

Simultaneous Inhibition of SARS-CoV-2 Entry Pathways by Cyclosporine

The COVID-19 pandemic caused by SARS-CoV-2 represents a global public health emergency. The entry of SARS-CoV-2 into host cells requires the activation of its spike protein by host cell proteases. The serine protease, TMPRSS2, and cysteine proteases, Cathepsins B/L, activate spike protein and enable SARS-CoV-2 entry to the host cell through two completely different and independent pathways. Therefore, inhibiting either TMPRSS2 or cathepsin B/L may not sufficiently block the virus entry. We here hypothesized that simultaneous targeting of both the entry pathways would be more efficient to block the virus entry rather than targeting the entry pathways individually. To this end, we utilized the network-based drug repurposing analyses to identify the possible common drugs that can target both the entry pathways. This study, for the first time, reports the molecules like cyclosporine, calcitriol, and estradiol as candidate drugs with the binding ability to the host proteases, TMPRSS2, and cathepsin B/L. Next, we analyzed drug-gene and gene-gene interaction networks using 332 human targets of SARS-CoV-2 proteins. The network results indicate that, out of 332 human proteins, cyclosporine interacts with 216 (65%) proteins. Furthermore, we performed molecular docking and all-atom molecular dynamics (MD) simulations to explore the binding of drug with TMPRSS2 and cathepsin L. The molecular docking and MD simulation results showed strong and stable binding of cyclosporine A (CsA) with TMPRSS2 and CTSL genes. The above results indicate cyclosporine as a potential drug molecule, as apart from interacting with SARS-CoV-2 entry receptors, it also interacts with most of SARS-CoV-2 target host genes; thus it could potentially interfere with functions of SARS-CoV-2 proteins in human cells. We here also suggest that these antiviral drugs alone or in combination can simultaneously target both the entry pathways and thus can be considered as a potential treatment option for COVID-19.

Tethering-induced destabilization and ATP-binding for tandem RRM domains of ALS-causing TDP-43 and hnRNPA1

TDP-43 and hnRNPA1 contain tandemly-tethered RNA-recognition-motif (RRM) domains, which not only functionally bind an array of nucleic acids, but also participate in aggregation/fibrillation, a pathological hallmark of various human diseases including amyotrophic lateral sclerosis (ALS), frontotemporal dementia (FTD), alzheimer's disease (AD) and Multisystem proteinopathy (MSP). Here, by DSF, NMR and MD simulations we systematically characterized stability, ATP-binding and conformational dynamics of TDP-43 and hnRNPA1 RRM domains in both tethered and isolated forms. The results reveal three key findings: (1) upon tethering TDP-43 RRM domains become dramatically coupled and destabilized with Tm reduced to only 49 °C. (2) ATP specifically binds TDP-43 and hnRNPA1 RRM domains, in which ATP occupies the similar pockets within the conserved nucleic-acid-binding surfaces, with the affinity slightly higher to the tethered than isolated forms. (3) MD simulations indicate that the tethered RRM domains of TDP-43 and hnRNPA1 have higher conformational dynamics than the isolated forms. Two RRM domains become coupled as shown by NMR characterization and analysis of inter-domain correlation motions. The study explains the long-standing puzzle that the tethered TDP-43 RRM1–RRM2 is particularly prone to aggregation/fibrillation, and underscores the general role of ATP in inhibiting aggregation/fibrillation of RRM-containing proteins. The results also rationalize the observation that the risk of aggregation-causing diseases increases with aging.

Insights into the Conserved Regulatory Mechanisms of Human and Yeast Aging

Aging represents a significant biological process having strong associations with cancer, diabetes, and neurodegenerative and cardiovascular disorders, which leads to progressive loss of cellular functions and viability. Astonishingly, age-related disorders share several genetic and molecular mechanisms with the normal aging process. Over the last three decades, budding yeast Saccharomyces cerevisiae has emerged as a powerful yet simple model organism for aging research. Genetic approaches using yeast RLS have led to the identification of hundreds of genes impacting lifespan in higher eukaryotes. Numerous interventions to extend yeast lifespan showed an analogous outcome in multi-cellular eukaryotes like fruit flies, nematodes, rodents, and humans. We collected and analyzed a multitude of observations from published literature and provide the contribution of yeast in the understanding of aging hallmarks most applicable to humans. Here, we discuss key pathways and molecular mechanisms that underpin the evolutionarily conserved aging process and summarize the current understanding and clinical applicability of its trajectories. Gathering critical information on aging biology would pave the way for future investigation targeted at the discovery of aging interventions.

Delineating the conformational dynamics of intermediate structures on the unfolding pathway of β-lactoglobulin in aqueous urea and dimethyl sulfoxide

The funnel shaped energy landscape model of the protein folding suggests that progression of folding proceeds through multiple pathways, having the multiple intermediates which leads to multidimensional free-energy surface. Herein, we applied all-atom MD simulation to conduct a comparative study on the structure of β-lactoglobulin (β-LgA) in aqueous mixture of 8 M urea and 8 M dimethyl sulfoxide (DMSO), at different temperatures. The cumulative results of multiple simulations suggest a common unfolding pathway of β-LgA, occurred through the stable and meta-stable intermediates (I), in both urea and DMSO. However, the free-energy landscape (FEL) analyses show that the structural transitions of I-states are energetically different. In urea, FEL shows distinct ensemble of intermediates, I1 and I2, separated by the energy barrier of ∼3.0 kcal mol^-1. Similarly, we find the population of two distinct I1 and I2 states in DMSO, however, the I1 appeared transiently around ∼30-35 ns and is short-lived. But, the I2 ensemble is observed structurally compact and long-lived (∼50-150 ns) as compared to unfolding in urea. Furthermore, the I1 and I2 are separated through a high energy barrier of ∼6.0 kcal mol^-1. Thus, our results provide the structural insights of intermediates which essentially bear the signature of a different unfolding pathway of β-LgA in urea and DMSO.Abbreviationsβ-LgAβ-lactoglobulinDMSOdimethyl sulfoxideFELfree-energy landscapeGdmClguanidinium chlorideIintermediate stateMGmolten globule statePMEparticle mesh EwaldQfraction of native contactsRMSDroot mean square deviationRMSFroot mean square fluctuationR_gradius of gyrationSASAsolvent Accessible Surface AreascSASAthe side chain SASATrptryptophanCommunicated by Ramaswamy H. Sarma.

Unfolding Transitions of Peripheral Subunit Binding Domains Show Cooperative Behavior

Characterization of native, intermediate, and denatured states is crucial for understanding the factors influencing the stability of proteins. We have carried out molecular dynamics simulations to study the unfolding of three peripheral subunit binding domains (PSBDs): E. coli BBL, Bacillus stearothermophilus E3BD, and human hbSBD, at three different temperatures: 300, 330, and 400 K, and in the presence of two solvents: water and 5 M guanidinium hydrochloride (GndCl) solution. These proteins share similar folds, with two parallel helices, maintained via a hydrophobic core comprising residues from their interconnecting loop. BBL is more sensitive to thermal and chemical denaturation in comparison to hbSBD, and E3BD is the most stable of all of the three proteins. The effect of temperature on the stability of these proteins is more pronounced in "water-only" simulations compared to that in the presence of guanidium hydrochloride in high concentrations. Our results show cooperative unfolding transitions of these proteins, which are triggered by an initial melting of the C-terminal helix H2. The consequent loss of interhelical interactions or native contacts, as observed, leads to the subsequent melting of the N-terminal helix H1.

Elucidation of the structural stability and dynamics of heterogeneous intermediate ensembles in unfolding pathway of the N-terminal domain of TDP-43

The N-terminal domain of the RNA binding protein TDP-43 (NTD) is essential to both physiology and proteinopathy; however, elucidation of its folding/unfolding still remains a major quest. In this study, we have investigated the biophysical behavior of intermediate ensembles employing all-atom molecular dynamics simulations in 8 M urea accelerated with high temperatures to achieve unfolded states in a confined computation time. The cumulative results of the 2.75 μs simulations show that unfolding of the NTD at 350 K evolves through different stable and meta-stable intermediate states. The free-energy landscape reveals two meta-stable intermediates (IN and IU) stabilized by non-native interactions, which are largely hydrophilic and highly energetically frustrated. A single buried tryptophan residue, W80, undergoes solvent exposure to different extents during unfolding; this suggests a structurally heterogeneous population of intermediate ensembles. Furthermore, the structure properties of the IN state show a resemblance to the molten globule (MG) state with most of the secondary structures intact. The unfolding of the NTD is initiated by the loss of β-strands, and the unfolded (U) states exhibit a population of non-native α-helices. These non-native unfolded intermediate ensembles may mediate protein oligomerization, leading to the formation of pathological, irreversible aggregates, characteristics of disease pathogenesis.