Mannose 2, 3, 4, 5, 6-O-pentasulfate (MPS): a partial activator of human heparin cofactor II with anticoagulation potential

Thromboembolic diseases are a major cause of mortality in human and the currently available anticoagulants are associated with various drawbacks, therefore the search for anticoagulants that have better safety profile is highly desirable. Compounds that are part of the dietary routine can be modified to possibly increase their anticoagulant potential. We show mannose 2,3,4,5,6-O-pentasulfate (MPS) as a synthetically modified form of mannose that has appreciable anticoagulation properties. An in silico study identified that mannose in sulfated form can bind effectively to the heparin-binding site of antithrombin (ATIII) and heparin cofactor II (HCII). Mannose was sulfated using a simple sulfation strategy-involving triethylamine-sulfur trioxide adduct. HCII and ATIII were purified from human plasma and the binding analysis using fluorometer and isothermal calorimetry showed that MPS binds at a unique site. A thrombin inhibition analysis using the chromogenic substrate showed that MPS partially enhances the activity of HCII. Further an assessment of in vitro blood coagulation assays using human plasma showed that the activated partial thromboplastin time (APTT) and prothrombin time (PT) were prolonged in the presence of MPS. A molecular dynamics simulation analysis of the HCII-MPS complex showed fluctuations in a N-terminal loop and the cofactor binding site of HCII. The results indicate that MPS is a promising lead due to its effect on the in vitro coagulation rate.Communicated by Ramaswamy H. Sarma.

A computational strategy for therapeutic development against superoxide dismutase (SOD1) amyloid formation: effect of polyphenols on the various events in the aggregation pathway

Pathology of superoxide dismutase 1 (SOD1) aggregation is linked to a neurodegenerative disease known as amyotrophic lateral sclerosis (ALS). Without suitable post-translational modifications (PTMs), the protein structure tends to become aggregation-prone. Understanding the role of PTMs and targeting the aggregation-prone SOD1 with small molecules can be used to design a strategy to inhibit its aggregation. Microsecond long molecular dynamics (MD) simulations followed by free energy surface (FES) analyses show that the loss of structure in the apo monomer happens locally and stepwise. Removing the disulfide bond from apoprotein leads to further instability in the zinc-binding loop, giving rise to non-native protein conformations. Further, it was found that these non-native conformations have a higher propensity to form a non-native dimer. We chose three structurally similar polyphenols based on their binding energies and investigated their impact on SOD1 aggregation kinetics. MD simulations of apo-SOD1^SH/corkscrew fibril-polyphenol complexes were also carried out. The effect of polyphenols was seen on fibril elongation as well. Based on the experiments and MD simulation results, it can be inferred that the choice of inhibitors is influenced not only by the binding energy but also by dimer interface stabilization, the proclivity to form non-native dimers, the propensity to break fibrils, and the propensity to decrease the rate of elongation. The polyphenols with 3' and 4' hydroxyl groups are better inhibitors of SOD1 aggregation.

The Role of Superoxide Dismutase 1 in Amyotrophic Lateral Sclerosis: Identification of Signaling Pathways, Regulators, Molecular Interaction Networks, and Biological Functions through Bioinformatics

Mutations in superoxide dismutase 1 (SOD1) result in misfolding and aggregation of the protein, causing neurodegenerative amyotrophic lateral sclerosis (ALS). In recent years, several new SOD1 variants that trigger ALS have been identified, making it increasingly crucial to understand the SOD1 toxicity pathway in ALS. Here we used an integrated bioinformatics approach, including the Ingenuity Pathway Analysis (IPA) tool to analyze signaling pathways, regulators, functions, and network molecules of SOD1 with an emphasis on ALS. IPA toxicity analysis of SOD1 identified superoxide radicals' degradation, apelin adipocyte, ALS, NRF2-mediated oxidative stress response, and sirtuin signaling as the key signaling pathways, while the toxicity of SOD1 is exerted via mitochondrial swelling and oxidative stress. IPA listed CNR1, APLN, BTG2, MAPK, DRAP1, NFE2L2, SNCA, and CG as the upstream regulators of SOD1. IPA further revealed that mutation in SOD1 results in hereditary disorders, including ALS. The exploration of the relationship between SOD1 and ALS using IPA unveiled SOD1-ALS pathway molecules. The gene ontology (GO) analysis of SOD1-ALS pathway molecules with ShinyGO reaffirmed that SOD1 toxicity results in ALS and neurodegeneration. The GO analysis further identified enriched biological processes, molecular functions, and cellular components for SOD1-ALS pathway molecules. The construction of a protein-protein interaction network of SOD1-ALS pathway molecules using STRING and further analysis of that network with Cytoscape identified ACTB followed by TP53, IL6, CASP3, SOD1, IL1B, APP, APOE, and VEGFA as the major network hubs. Taken together, our study provides insight into the molecular underpinning of SOD1's toxicity in ALS.

Native Mass Spectrometry Coupled to Spectroscopic Methods to Investigate the Effect of Soybean Isoflavones on Structural Stability and Aggregation of Zinc Deficient and Metal-Free Superoxide Dismutase

The deficiency or wrong combination of metal ions in Cu, Zn-superoxide dismutase (SOD1), is regarded as one of the main factors causing the aggregation of SOD1 and then inducing amyotrophic lateral sclerosis (ALS). A ligands-targets screening process based on native electrospray ionization ion mobility mass spectrometry (ESI-IMS-MS) was established in this study. Four glycosides including daidzin, sophoricoside, glycitin, and genistin were screened out from seven soybean isoflavone compounds and were found to interact with zinc-deficient or metal-free SOD1. The structure and conformation stability of metal-free and zinc-deficient SOD1 and their complexes with the four glycosides was investigated by collision-induced dissociation (CID) and collision-induced unfolding (CIU). The four glycosides could strongly bind to the metal-free and copper recombined SOD1 and enhance the folding stability of these proteins. Additionally, the ThT fluorescence assay showed that these glycosides could inhibit the toxic aggregation of the zinc-deficient or metal-free SOD1. The competitive interaction experiments together with molecular docking indicate that glycitin, which showed the best stabilizing effects, binds with SOD1 between β-sheet 6 and loop IV. In short, this study provides good insight into the relationship between inhibitors and different SOD1s.

Targeted protein degraders march towards the clinic for neurodegenerative diseases

Protein degraders are emerging as potent therapeutic tools to address neurological disorders and many complex diseases. It offered several key advantages, including the doses, drug resistance, and side effects over traditional occupancy-based inhibitors. Translation of chemical degraders into a clinical therapy for neurodegenerative disorders has a well-documented knowledge and resource gap. Researchers strive to develop clinical candidates employing chemical degraders' technologies, including hydrophobic tagging, molecular glues, proteolysis targeting chimeras (PROTACs), specific and nongenetic Inhibitor of Apoptosis Protein (IAP)-dependent protein erasers (SNIPERs), autophagy targeted chimeras, and autophagosome-tethered compounds for targeted degradation of pathological markers in neurodegenerative disease. Herein, we examined the present state of chemical-mediated targeted protein degradation in the quest for medications to treat neurodegenerative diseases. We further identified targeted degraders under clinical development for neurodegenerative diseases summarizing pertinent discoveries guiding the future of degradation therapeutics. We also addressed the necessary pharmacological interventions needed to achieve unprecedented therapeutic efficacy and its associated challenges.

First Principles Calculation of Protein-Protein Dimer Affinities of ALS-Associated SOD1 Mutants

Cu,Zn superoxide dismutase (SOD1) is a 32 kDa homodimer that converts toxic oxygen radicals in neurons to less harmful species. The dimerization of SOD1 is essential to the stability of the protein. Monomerization increases the likelihood of SOD1 misfolding into conformations associated with aggregation, cellular toxicity, and neuronal death in familial amyotrophic lateral sclerosis (fALS). The ubiquity of disease-associated mutations throughout the primary sequence of SOD1 suggests an important role of physicochemical processes, including monomerization of SOD1, in the pathology of the disease. Herein, we use a first-principles statistical mechanics method to systematically calculate the free energy of dimer binding for SOD1 using molecular dynamics, which involves sequentially computing conformational, orientational, and separation distance contributions to the binding free energy. We consider the effects of two ALS-associated mutations in SOD1 protein on dimer stability, A4V and D101N, as well as the role of metal binding and disulfide bond formation. We find that the penalty for dimer formation arising from the conformational entropy of disordered loops in SOD1 is significantly larger than that for other protein-protein interactions previously considered. In the case of the disulfide-reduced protein, this leads to a bound complex whose formation is energetically disfavored. Somewhat surprisingly, the loop free energy penalty upon dimerization is still significant for the holoprotein, despite the increased structural order induced by the bound metal cations. This resulted in a surprisingly modest increase in dimer binding free energy of only about 1.5 kcal/mol upon metalation of the protein, suggesting that the most significant stabilizing effects of metalation are on folding stability rather than dimer binding stability. The mutant A4V has an unstable dimer due to weakened monomer-monomer interactions, which are manifested in the calculation by a separation free energy surface with a lower barrier. The mutant D101N has a stable dimer partially due to an unusually rigid β-barrel in the free monomer. D101N also exhibits anticooperativity in loop folding upon dimerization. These computational calculations are, to our knowledge, the most quantitatively accurate calculations of dimer binding stability in SOD1 to date.

Phosphorylation Regulation on the Homo-Dimeric Binding of Transactive Response DNA-Binding Protein

The dimerization of transactive response DNA-binding protein of 43 kDa (TDP-43) is crucial for the RNA metabolism, and the higher-order aggregation of TDP-43 would induce several neurodegenerative diseases. The dimerization and aggregation of TDP-43 are regulated by the phosphorylation on its N-terminal domain (NTD). Understanding the regulation mechanism of TDP-43 NTD dimerization is crucial for the preventing of harmful aggregation and the associated diseases. In this study, the dimerization processes of wild-type (WT), phosphorylated S48 (pS48), and phosphomimic S48E mutation (S48E) of TDP-43 NTD are characterized by the enhanced sampling technology. Our results show that the phosphorylation not only shift the conformation population of bound and unbound state of TDP-43 NTD, but also would regulate the dimerization processes, including increase the binding free-energy barrier. The phosphomimic mutation would also shift the conformational space of TDP-43 NTD dimer to the unbound structures; however, the thermodynamic and kinetic properties of the dimerization processes between the phosphorylated and phosphomimic mutant systems are distinct, which reminds us to carefully study the phosphorylation regulation by using the phosphomimic mutations.

Delineating the Aggregation-Prone Hotspot Regions (Peptides) in the Human Cu/Zn Superoxide Dismutase 1

Amyotrophic lateral sclerosis (ALS) is a fatal, incurable neurodegenerative disease described by progressive degeneration of motor neurons. The most common familial form of ALS (fALS) has been associated with mutations in the Cu/Zn superoxide dismutase (SOD1) gene. Mutation-induced misfolding and aggregation of SOD1 is often found in ALS patients. In this work, we probe the aggregation properties of peptides derived from the SOD1. To examine the source of SOD1 aggregation, we have employed a computational algorithm to identify four peptides from the SOD1 protein sequence that aggregates into a fibril. Aided by computational algorithms, we identified four peptides likely involved in SOD1 fibrillization. These four aggregation-prone peptides were 14VQGIINFE21, 30KVWGSIKGL38, 101DSVISLS107, and 147GVIGIAQ153. In addition, the formation of fibril propensities from the identified peptides was investigated through different biophysical techniques. The atomic structures of two fibril-forming peptides from the C-terminal SOD1 showed that the steric zippers formed by 101DSVISLS107 and 147GVIGIAQ153 vary in their arrangement. We also discovered that fALS mutations in the peptide 147GVIGIAQ153 increased the fibril-forming propensity and altered the steric zipper's packing. Thus, our results suggested that the C-terminal peptides of SOD1 have a central role in amyloid formation and might be involved in forming the structural core of SOD1 aggregation observed in vivo.

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.

Structural genomics approach to investigate deleterious impact of nsSNPs in conserved telomere maintenance component 1

Conserved telomere maintenance component 1 (CTC1) is an important component of the CST (CTC1-STN1-TEN1) complex, involved in maintaining the stability of telomeric DNA. Several non-synonymous single-nucleotide polymorphisms (nsSNPs) in CTC1 have been reported to cause Coats plus syndrome and Dyskeratosis congenital diseases. Here, we have performed sequence and structure analyses of nsSNPs of CTC1 using state-of-the-art computational methods. The structure-based study focuses on the C-terminal OB-fold region of CTC1. There are 11 pathogenic mutations identified, and detailed structural analyses were performed. These mutations cause a significant disruption of noncovalent interactions, which may be a possible reason for CTC1 instability and consequent diseases. To see the impact of such mutations on the protein conformation, all-atom molecular dynamics (MD) simulations of CTC1-wild-type (WT) and two of the selected mutations, R806C and R806L for 200 ns, were carried out. A significant conformational change in the structure of the R806C mutant was observed. This study provides a valuable direction to understand the molecular basis of CTC1 dysfunction in disease progression, including Coats plus syndrome.

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.

Targeting SARS-CoV-2 nucleocapsid oligomerization: Insights from molecular docking and molecular dynamics simulations

The outbreak of COVID-19 caused by SARS-CoV-2 virus continually led to infect a large population worldwide. Currently, there is no specific viral protein-targeted therapeutics. The Nucleocapsid (N) protein of the SARS-CoV-2 virus is necessary for viral RNA replication and transcription. The C-terminal domain of N protein (CTD) involves in the self-assembly of N protein into a filament that is packaged into new virions. In this study, the CTD (PDB ID: 6WJI) was targeted for the identification of possible inhibitors of oligomerization of N protein. Herein, multiple computational approaches were employed to explore the potential mechanisms of binding and inhibitor activity of five antiviral drugs toward CTD. The five anti-N drugs studied in this work are 4E1RCat, Silmitasertib, TMCB, Sapanisertib, and Rapamycin. Among the five drugs, 4E1RCat displayed highest binding affinity (-10.95 kcal/mol), followed by rapamycin (-8.91 kcal/mol), silmitasertib (-7.89 kcal/mol), TMCB (-7.05 kcal/mol), and sapanisertib (-6.14 kcal/mol). Subsequently, stability and dynamics of the protein-drug complex were examined with molecular dynamics (MD) simulations. Overall, drug binding increases the stability of the complex with maximum stability observed in the case of 4E1RCat. The CTD-drug complex systems behave differently in terms of the free energy landscape and showed differences in population distribution. Overall, the MD simulation parameters like RMSD, RMSF, Rg, hydrogen bonds analysis, PCA, FEL, and DCCM analysis indicated that 4E1RCat and TMCB complexes were more stable as compared to silmitasertib and sapanisertib and thus could act as effective drug compounds against CTD.Communicated by Ramaswamy H. Sarma.

Trifluoroethanol Partially Unfolds G93A SOD1 Leading to Protein Aggregation: A Study by Native Mass Spectrometry and FPOP Protein Footprinting

Misfolding of Cu, Zn superoxide dismutase (SOD1) variants may lead to protein aggregation and ultimately amyotrophic lateral sclerosis (ALS). The mechanism and protein conformational changes during this process are complex and remain unclear. To study SOD1 variant aggregation at the molecular level and in solution, we chemically induced aggregation of a mutant variant (G93A SOD1) with trifluoroethanol (TFE) and used both native mass spectrometry (MS) to analyze the intact protein and fast photochemical oxidation of proteins (FPOP) to characterize the structural changes induced by TFE. We found partially unfolded G93A SOD1 monomers prior to oligomerization and identified regions of the N-terminus, C-terminus, and strands β5, β6 accountable for the partial unfolding. We propose that exposure of hydrophobic interfaces of these unstructured regions serves as a precursor to aggregation. Our results provide a possible mechanism and molecular basis for ALS-linked SOD1 misfolding and aggregation.

Targeting virus-host interaction by novel pyrimidine derivative: an in silico approach towards discovery of potential drug against COVID-19

The entire human population over the globe is currently facing appalling conditions due to the spread of infection from coronavirus disease-2019 (COVID-19). The spike glycoprotein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) present on the surface of the virion mediates the virus entry into the host cells and therefore is targeted by several scientific groups as a novel drug target site. The spike glycoprotein binds to the human angiotensin-converting enzyme-2 (hACE2) cell surface receptor abundantly expressed in lung tissues, and this binding phenomenon is a primary determinant of cell tropism and pathogenesis. The binding and internalization of the virus is the primary and most crucial step in the process of infection, and therefore the molecules targeting the inhibition of this process certainly hold a significant therapeutic value. Thus, we systematically applied the computational techniques to identify the plausible inhibitor from a chosen set of well characterized diaryl pyrimidine analogues which may disrupt interfacial interaction of spike glycoprotein (S) at the surface of hACE2. Using molecular docking, molecular dynamics (MD) simulation and binding free energy calculation, we have identified AP-NP (2-(2-amino-5-(naphthalen-2-yl)pyrimidin-4-yl)phenol), AP-3-OMe-Ph (2-(2-amino-5-(3-methoxyphenyl)pyrimidin-4-yl)phenol) and AP-4-Me-Ph (2-(2-amino-5-(p-tolyl) pyrimidin-4-yl)phenol) from a group of diaryl pyrimidine derivatives which appears to bind at the interface of the hACE2-S complex with low binding free energy. Thus, pyrimidine derivative AP-NP may be explored as an effective inhibitor for hACE2-S complex. Furthermore, in vitro and in vivo studies will strengthen the use of these inhibitors as suitable drug candidates against SARS-COV-2.

Communicated by Ramaswamy H. Sarma

Identification of potential inhibitors against SARS-CoV-2 by targeting proteins responsible for envelope formation and virion assembly using docking based virtual screening, and pharmacokinetics approaches

WHO has declared the outbreak of COVID-19 as a public health emergency of international concern. The ever-growing new cases have called for an urgent emergency for specific anti-COVID-19 drugs. Three structural proteins (Membrane, Envelope and Nucleocapsid protein) play an essential role in the assembly and formation of the infectious virion particles. Thus, the present study was designed to identify potential drug candidates from the unique collection of 548 anti-viral compounds (natural and synthetic anti-viral), which target SARS-CoV-2 structural proteins. High-end molecular docking analysis was performed to characterize the binding affinity of the selected drugs-the ligand, with the SARS-CoV-2 structural proteins, while high-level Simulation studies analyzed the stability of drug-protein interactions. The present study identified rutin, a bioflavonoid and the antibiotic, doxycycline, as the most potent inhibitor of SARS-CoV-2 envelope protein. Caffeic acid and ferulic acid were found to inhibit SARS-CoV-2 membrane protein while the anti-viral agent's simeprevir and grazoprevir showed a high binding affinity for nucleocapsid protein. All these compounds not only showed excellent pharmacokinetic properties, absorption, metabolism, minimal toxicity and bioavailability but were also remain stabilized at the active site of proteins during the MD simulation. Thus, the identified lead compounds may act as potential molecules for the development of effective drugs against SARS-CoV-2 by inhibiting the envelope formation, virion assembly and viral pathogenesis.

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.

T54R mutation destabilizes the dimer of superoxide dismutase 1T54R by inducing steric clashes at the dimer interface

Mutations cause abnormalities in protein structure, function and oligomerization. Different mutations in the superoxide dismutase 1 (SOD1) protein cause its misfolding, loss of dimerization and aggravate its aggregation in the amyotrophic lateral sclerosis disease. In this study, we report the mechanistic details of how a threonine-to-arginine mutation at the 54^th position (T54R) of SOD1 results in destabilization of the dimer interface of SOD1^T54R. Using computational and experimental methods, we show that the T54R mutation increases fluctuation of the mutation-harboring loop (R54-loop) of SOD1^T54R. Fluctuation of this loop causes steric clashes that involve arginine-54 (R54) and other residues of SOD1^T54R, resulting in loss of inter-subunit contacts at the dimer interface. Since the T54 residue-containing loop is necessary for the dimerization of wild-type SOD1, fluctuation of the R54-loop, steric clashes involving R54 and loss of inter-subunit contacts give rise to the loss of SOD1^T54R dimer stability. This correlates to energetically unfavorable tethering of the monomers of SOD1^T54R. The outcome is gradual splitting of SOD1^T54R dimers into monomers, thereby exposing the previously buried hydrophobic interface residues to the aqueous environment. This event finally leads to aggregation of SOD1^T54R. T54R mutation has no effect in altering the relative positions of copper and zinc ion binding residues of SOD1^T54R. The native SOD1 structure is stable, and there is no destabilizing effect at its dimer interface. Overall, our study reveals the intricate mechanism of T54R mutation-associated destabilization of the dimer of the SOD1^T54R protein.

T54R mutation destabilizes the dimer of SOD1^T54R.

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.