Thermal Unfolding Pathway of PHD2 Catalytic Domain in Three Different PHD2 Species: Computational Approaches

Prolyl hydroxylase domain 2 containing protein (PHD2) is a key protein in regulation of angiogenesis and metastasis. In normoxic condition, PHD2 triggers the degradation of hypoxia-inducible factor 1 (HIF-1α) that induces the expression of hypoxia response genes. Therefore the correct function of PHD2 would inhibit angiogenesis and consequent metastasis of tumor cells in normoxic condition. PHD2 mutations were reported in some common cancers. However, high levels of HIF-1α protein were observed even in normoxic metastatic tumors with normal expression of wild type PHD2. PHD2 malfunctions due to protein misfolding may be the underlying reason of metastasis and invasion in such cases. In this study, we scrutinize the unfolding pathways of the PHD2 catalytic domain’s possible species and demonstrate the properties of their unfolding states by computational approaches. Our study introduces the possibility of aggregation disaster for the prominent species of PHD2 during its partial unfolding. This may justify PHD2 inability to regulate HIF-1α level in some normoxic tumor types.     

Predicting Atomic Details of the Unfolding Pathway for YibK,a Knotted Protein from the SPOUT Superfamily

Abstract

Several protein structures have been reported to contain intricate knots of the polypeptide backbone but the mechanism of the (un)folding process of knotted proteins remains unknown. The members of the SPOUT superfamily of RNA methyltransferases are some of the most intensely studied systems for investigation of the knot formation and function. YibK (whose biochemical function remains unknown) is the representative protein of the SPOUT superfamily. This protein exhibits a deep trefoil knot at the C-terminus.

We conducted an extensive computational analysis of the unfolding process for the monomeric form of YibK. In order to predict the (un)folding pathway of YibK, we have calculated the order of secondary structure disassembly using UNFOLD, and performed thermal unfolding simulations using classical Molecular Dynamics (MD), as well as simulations employing reduced representation of the peptide chain using either MD with the UNRES method or the Monte Carlo (MC) unfolding with the REFINER method.

Results obtained from all methods used in this work are in qualitative agreement. We found that YibK unfolds through four intermediate states. The trefoil knot in YibK disappears at the end of the unfolding process, long after the protein loses its native topology. We observed that the C-terminus leaves the knotting loop folded into a hairpin-like structure, in agreement with the results of coarse-grained simulation reported earlier. We propose that the folding pathway of YibK corresponds to the reversed sequence of events observed in the unfolding pathway elucidated in this study. Thus, we predict that the knot formation is the slowest part of the YibK folding process.     

All‐atom and coarse‐grained simulations of the forced unfolding pathways of the SNARE complex

The SNARE complex, consisting of three proteins (VAMP2, syntaxin, and SNAP‐25), is thought to drive membrane fusion by assembling into a four‐helix bundle through a zippering process. In support of the above zippering model, a recent single‐molecule optical tweezers experiment by Gao et al. revealed a sequential unzipping of SNARE along VAMP2 in the order of the linker domain → the C‐terminal domain → the N‐terminal domain. To offer detailed structural insights to this unzipping process, we have performed all‐atom and coarse‐grained steered molecular dynamics (sMD) simulations of the forced unfolding pathways of SNARE using different models and force fields. Our findings are summarized as follows: First, the sMD simulations based on either an all‐atom force field (with an implicit solvent model) or a coarse‐grained Go model were unable to capture the forced unfolding pathway of SNARE as observed by Gao et al., which may be attributed to insufficient simulation time and inaccurate force fields. Second, the sMD simulations based on a reparameterized coarse‐grained model (i.e., modified elastic network model) were able to predict a sequential unzipping of SNARE in good agreement with the findings by Gao et al. The key to this success is to reparameterize the intrahelix and interhelix nonbonded force constants against the pair‐wise residue–residue distance fluctuations collected from all‐atom MD simulations of SNARE. Therefore, our finding supports the importance of accurately describing the inherent dynamics/flexibility of SNARE (in the absence of force), in order to correctly simulate its unfolding behaviors under force. This study has established a useful computational framework for future studies of the zippering function of SNARE and its perturbations by point mutations with amino‐acid level of details, and more generally the forced unfolding pathways of other helix bundle proteins. Proteins 2014; 82:1376–1386. © 2014 Wiley Periodicals, Inc.     

New insights into FAK signaling and localization based on detection of a FAT domain folding intermediate

Mounting evidence suggests that the focal adhesion targeting (FAT) domain, an antiparallel four-helix bundle, exists in alternative conformations that may modulate phosphorylation, ligand binding, and the subcellular localization of focal adhesion kinase (FAK). In order to characterize the conformational dynamics of the FAT domain, we have developed a novel method for reconstructing the folding pathway of the FAT domain by using discrete molecular dynamics (DMD) simulations, with free energy constraints derived from NMR hydrogen exchange data. The DMD simulations detect a folding intermediate, in which a cooperative unfolding event causes helix 1 to lose helical character while separating from the helix bundle. The conformational dynamic features of helix 1 in the intermediate state of the FAT domain are likely to facilitate Y926 phosphorylation, yet interfere with paxillin binding. The presence of this intermediate state in vivo may promote FAK signaling via the ERK/MAPK pathway and by release of FAK from focal adhesions.     

Diffusing and colliding: the atomic level folding/unfolding pathway of a small helical protein

Proteins with ultra-fast folding/unfolding kinetics are excellent candidates for study by molecular dynamics. Here, we describe such simulations of a three helix bundle protein, the engrailed homeodomain (En-HD), which folds via the diffusion-collision model. The unfolding pathway of En-HD was characterized by seven simulations of the protein and 12 simulations of its helical fragments yielding over 1.1 micros of simulation time in water. Various conformational states along the unfolding pathway were identified. There is the compact native-like transition state, a U-shaped helical intermediate and an unfolded state with dynamic helical segments. Each of these states is in good agreement with experimental data. Examining these states as well as the transitions between them, we find the role of long-range tertiary contacts, specifically salt-bridges, important in the folding/unfolding pathway. In the folding direction, charged residues form long-range tertiary contacts before the hydrophobic core is formed. The formation of HII is assisted by a specific salt-bridge and by non-specific (fluctuating) tertiary contacts, which we call contact-assisted helix formation. Salt-bridges persist as the protein approaches the transition state, stabilizing HII until the hydrophobic core is formed. To complement this information, simulations of fragments of En-HD illustrate the helical propensities of the individual segments. By thermal denaturation, HII proved to be the least stable helix, unfolding in less than 450 ps at high temperature. We observed the low helical propensity of C-terminal residues from HIII in fragment simulations which, when compared to En-HD unfolding simulations, link the unraveling of HIII to the initial event that drives the unfolding of En-HD.     

LpxA: A natural nanotube

UDP‐N‐acetylglucosamine 3‐O‐acyltransferase is a protein with a left‐handed parallel β‐helix, which is a natural nanotube. They are associated with unusual high stability. To identify the reason behind the structural stability of β‐helical nanotubular structure, we have performed a total of 4 μs molecular dynamics simulations of the protein in implicit solvent at four different temperatures and monitored the unfolding pathway. The correlation in movement between different regions of the nanotubular structure has been identified from the dynamical cross‐correlation map and contribution of some specific residues towards unfolding transition has been identified by principal component analysis. Difference in stability of the three loop regions has also been characterized. Construction of the unfolding conformational energy landscape identifies the probable intermediates that can appear in the unfolding pathway of the protein. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 845–853, 2010.     

Studying pressure denaturation of a protein by molecular dynamics simulations

Many globular proteins unfold when subjected to several kilobars of hydrostatic pressure. This “unfolding‐up‐on‐squeezing” is counter‐intuitive in that one expects mechanical compression of proteins with increasing pressure. Molecular simulations have the potential to provide fundamental understanding of pressure effects on proteins. However, the slow kinetics of unfolding, especially at high pressures, eliminates the possibility of its direct observation by molecular dynamics (MD) simulations. Motivated by experimental results—that pressure denatured states are water‐swollen, and theoretical results—that water transfer into hydrophobic contacts becomes favorable with increasing pressure, we employ a water insertion method to generate unfolded states of the protein Staphylococcal Nuclease (Snase). Structural characteristics of these unfolded states—their water‐swollen nature, retention of secondary structure, and overall compactness—mimic those observed in experiments. Using conformations of folded and unfolded states, we calculate their partial molar volumes in MD simulations and estimate the pressure‐dependent free energy of unfolding. The volume of unfolding of Snase is negative (approximately ?60 mL/mol at 1 bar) and is relatively insensitive to pressure, leading to its unfolding in the pressure range of 1500–2000 bars. Interestingly, once the protein is sufficiently water swollen, the partial molar volume of the protein appears to be insensitive to further conformational expansion or unfolding. Specifically, water‐swollen structures with relatively low radii of gyration have partial molar volume that are similar to that of significantly more unfolded states. We find that the compressibility change on unfolding is negligible, consistent with experiments. We also analyze hydration shell fluctuations to comment on the hydration contributions to protein compressibility. Our study demonstrates the utility of molecular simulations in estimating volumetric properties and pressure stability of proteins, and can be potentially extended for applications to protein complexes and assemblies. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Potentially amyloidogenic conformational intermediates populate the unfolding landscape of transthyretin: Insights from molecular dynamics simulations

Protein aggregation into insoluble fibrillar structures known as amyloid characterizes several neurodegenerative diseases, including Alzheimer's, Huntington's and Creutzfeldt‐Jakob. Transthyretin (TTR), a homotetrameric plasma protein, is known to be the causative agent of amyloid pathologies such as FAP (familial amyloid polyneuropathy), FAC (familial amyloid cardiomiopathy) and SSA (senile systemic amyloidosis). It is generally accepted that TTR tetramer dissociation and monomer partial unfolding precedes amyloid fibril formation. To explore the TTR unfolding landscape and to identify potential intermediate conformations with high tendency for amyloid formation, we have performed molecular dynamics unfolding simulations of WT‐TTR and L55P‐TTR, a highly amyloidogenic TTR variant. Our simulations in explicit water allow the identification of events that clearly discriminate the unfolding behavior of WT and L55P‐TTR. Analysis of the simulation trajectories show that (i) the L55P monomers unfold earlier and to a larger extent than the WT; (ii) the single α‐helix in the TTR monomer completely unfolds in most of the L55P simulations while remain folded in WT simulations; (iii) L55P forms, early in the simulations, aggregation‐prone conformations characterized by full displacement of strands C and D from the main β‐sandwich core of the monomer; (iv) L55P shows, late in the simulations, severe loss of the H‐bond network and consequent destabilization of the CBEF β‐sheet of the β‐sandwich; (v) WT forms aggregation‐compatible conformations only late in the simulations and upon extensive unfolding of the monomer. These results clearly show that, in comparison with WT, L55P‐TTR does present a much higher probability of forming transient conformations compatible with aggregation and amyloid formation.     

FF12MC: A revised AMBER forcefield and new protein simulation protocol

Specialized to simulate proteins in molecular dynamics (MD) simulations with explicit solvation, FF12MC is a combination of a new protein simulation protocol employing uniformly reduced atomic masses by tenfold and a revised AMBER forcefield FF99 with (i) shortened C? H bonds, (ii) removal of torsions involving a nonperipheral sp³ atom, and (iii) reduced 1–4 interaction scaling factors of torsions ? and ψ. This article reports that in multiple, distinct, independent, unrestricted, unbiased, isobaric–isothermal, and classical MD simulations FF12MC can (i) simulate the experimentally observed flipping between left‐ and right‐handed configurations for C14–C38 of BPTI in solution, (ii) autonomously fold chignolin, CLN025, and Trp‐cage with folding times that agree with the experimental values, (iii) simulate subsequent unfolding and refolding of these miniproteins, and (iv) achieve a robust Z score of 1.33 for refining protein models TMR01, TMR04, and TMR07. By comparison, the latest general‐purpose AMBER forcefield FF14SB locks the C14–C38 bond to the right‐handed configuration in solution under the same protein simulation conditions. Statistical survival analysis shows that FF12MC folds chignolin and CLN025 in isobaric–isothermal MD simulations 2–4 times faster than FF14SB under the same protein simulation conditions. These results suggest that FF12MC may be used for protein simulations to study kinetics and thermodynamics of miniprotein folding as well as protein structure and dynamics. Proteins 2016; 84:1490–1516. © 2016 The Authors Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Early developmental stress negatively affects neuronal recruitment to avian song system nucleus HVC

Adverse environmental conditions can impact the life history trajectory of animals. Adaptive responses enable individuals to cope with unfavorable conditions, but altered metabolism and resource allocation can bear long‐term costs. In songbirds, early developmental stress can cause lifelong changes in learned song, a culturally transmitted trait, and nestlings experiencing developmental stress develop smaller song control nucleus HVCs. We investigated whether nutrition‐related developmental stress impacts neurogenesis in HVC, which may explain how poor nutrition leads to smaller HVC volume. We provided different quality diets (LOW and HIGH) by varying the husks‐to‐seeds ratio to zebra finch families for the first 35 days after the young hatched (PHD). At PHD14–18 and again at nutritional independence (PHD35), juveniles were injected with different cell division markers. To monitor growth, we took body measures at PHD10, 17, and 35. At PHD35 the number of newly recruited neurons in HVC and the rate of proliferation in the adjacent ventricular zone (VZ) were counted. Males raised on the LOW diet for their first weeks of life had significantly fewer new neurons in HVC than males raised on the HIGH diet. At the time when these new HVC neurons were born and labeled in the VZ (PHD17) the birds exposed to the LOW diet had significantly lower body mass. At PHD35 body mass or neuronal proliferation no longer differed. Our study shows that even transitory developmental stress can have negative consequences on the cellular processes underlying the development of neural circuits. © 2015 Wiley Periodicals, Inc. Develop Neurobiol 76: 107–118, 2016     

Folding and unfolding of gammaTIM monomers and dimers

Kinetic simulations of the folding and unfolding of triosephosphate isomerase (TIM) from yeast were conducted using a single monomer gammaTIM polypeptide chain that folds as a monomer and two gammaTIM chains that fold to the native dimer structure. The basic protein model used was a minimalist Gō model using the native structure to determine attractive energies in the protein chain. For each simulation type--monomer unfolding, monomer refolding, dimer unfolding, and dimer refolding--thirty simulations were conducted, successfully capturing each reaction in full. Analysis of the simulations demonstrates four main conclusions. First, all four simulation types have a similar "folding order", i.e., they have similar structures in intermediate stages of folding between the unfolded and folded state. Second, despite this similarity, different intermediate stages are more or less populated in the four different simulations, with 1), no intermediates populated in monomer unfolding; 2), two intermediates populated with beta(2)-beta(4) and beta(1)-beta(5) regions folded in monomer refolding; 3), two intermediates populated with beta(2)-beta(3) and beta(2)-beta(4) regions folded in dimer unfolding; and 4), two intermediates populated with beta(1)-beta(5) and beta(1)-beta(5) + beta(6) + beta(7) + beta(8) regions folded in dimer refolding. Third, simulations demonstrate that dimer binding and unbinding can occur early in the folding process before complete monomer-chain folding. Fourth, excellent agreement is found between the simulations and MPAX (misincorporation proton alkyl exchange) experiments. In total, this agreement demonstrates that the computational Gō model is accurate for gammaTIM and that the energy landscape of gammaTIM appears funneled to the native state.     

Self‐guided Langevin dynamics study of regulatory interactions in NtrC

Multiple self‐guided Langevin dynamics (SGLD) simulations were performed to examine structural and dynamical properties of the receiver domain of nitrogen regulatory protein C (NtrC^r). SGLD and MD simulations of the phosphorylated active form structure suggest a mostly stable but broad structural ensemble of this protein. The finite difference Poisson–Boltzmann calculations of the pK_a values of the active site residues suggest an increase in the pK_a of His‐84 on phosphorylation of Asp‐54. In SGLD simulations of the phosphorylated active form with charged His‐84, the average position of the regulatory helix α4 is found closer to the starting structure than in simulations with the neutral His‐84. To model the transition pathway, the phosphate group was removed from the simulations. After 7 ns of simulations, the regulatory helix α4 was found approximately halfway between positions in the NMR structures of the active and inactive forms. Removal of the phosphate group stimulated loss of helix α4, suggesting that the pathway of conformational transition may involve partial unfolding mechanism. The study illustrates the potential utility of the SGLD method in studies of the coupling between ligand binding and conformational transitions. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Structural and Biophysical Characterization of the Interactions between Calmodulin and the Pleckstrin Homology Domain of Akt

The translocation of Akt, a serine/threonine kinase, to the plasma membrane is a critical step in the Akt activation pathway. It is established that membrane binding of Akt is mediated by direct interactions between its pleckstrin homology domain (PHD) and phosphatidylinositol 3,4,5-trisphosphate (PI(3,4,5)P₃). There is now evidence that Akt activation in many breast cancer cells is also modulated by the calcium-binding protein, calmodulin (CaM). Upon EGF stimulation of breast cancer cells, CaM co-localizes with Akt at the plasma membrane to enhance activation. However, the molecular details of Akt(PHD) interaction with CaM are not known. In this study, we employed NMR, biochemical, and biophysical techniques to characterize CaM binding to Akt(PHD). Our data show that CaM forms a tight complex with the PHD of Akt (dissociation constant = 100 nm). The interaction between CaM and Akt(PHD) is enthalpically driven, and the affinity is greatly dependent on salt concentration, indicating that electrostatic interactions are important for binding. The CaM-binding interface in Akt(PHD) was mapped to two loops adjacent to the PI(3,4,5)P₃ binding site, which represents a rare CaM-binding motif and suggests a synergistic relationship between CaM and PI(3,4,5)P₃ upon Akt activation. Elucidation of the mechanism by which Akt interacts with CaM will help in understanding the activation mechanism, which may provide insights for new potential targets to control the pathophysiological processes of cell survival.     

Polar or Apolar—The Role of Polarity for Urea-Induced Protein Denaturation

Urea-induced protein denaturation is widely used to study protein folding and stability; however, the molecular mechanism and driving forces of this process are not yet fully understood. In particular, it is unclear whether either hydrophobic or polar interactions between urea molecules and residues at the protein surface drive denaturation. To address this question, here, many molecular dynamics simulations totalling ca. 7 µs of the CI2 protein in aqueous solution served to perform a computational thought experiment, in which we varied the polarity of urea. For apolar driving forces, hypopolar urea should show increased denaturation power; for polar driving forces, hyperpolar urea should be the stronger denaturant. Indeed, protein unfolding was observed in all simulations with decreased urea polarity. Hyperpolar urea, in contrast, turned out to stabilize the native state. Moreover, the differential interaction preferences between urea and the 20 amino acids turned out to be enhanced for hypopolar urea and suppressed (or even inverted) for hyperpolar urea. These results strongly suggest that apolar urea–protein interactions, and not polar interactions, are the dominant driving force for denaturation. Further, the observed interactions provide a detailed picture of the underlying molecular driving forces. Our simulations finally allowed characterization of CI2 unfolding pathways. Unfolding proceeds sequentially with alternating loss of secondary or tertiary structure. After the transition state, unfolding pathways show large structural heterogeneity.     

Activation of the prolyl hydroxylase oxygen-sensor results in induction of GLUT1, heme oxygenase-1, and nitric-oxide synthase proteins and confers protection from metabolic inhibition to cardiomyocytes

Recently an oxygen-sensing/transducing mechanism has been identified as a family of O2-dependent prolyl hydroxylase domain-containing enzymes (PHD). In normoxia, PHD hydroxylates a specific proline residue that directs the degradation of constitutively synthesized hypoxia-inducible factor-1alpha. During hypoxia, the cessation of hydroxylation of this proline results in less degradation and thus increases hypoxia-inducible factor-1alpha protein levels. In this study we have examined the consequences of activating the PHD oxygen-sensing pathway in cultured neonatal myocytes using ethyl-3,4 dihydroxybenzoate and dimethyloxalylglycine, inhibitors that, similar to hypoxia, inhibit this family of O2-dependent PHD enzymes. Increased glucose uptake and enhanced glycolytic metabolism are classical cellular responses to hypoxia. Ethyl-3,4 dihydroxybenzoate treatment of cardiomyocyte cultures for 24 h increased [3H]deoxy-4-glucose uptake concurrent with an induction of GLUT1 protein. In addition, ethyl-3,4 dihydroxybenzoate, dimethyloxalylglycine, and hypoxia treatments were found to induce protein levels of nitricoxide synthase-2 and heme oxygenase-1, two important cardioregulatory proteins whose expression in response to hypoxic conditions is poorly understood. In conjunction with these changes in gene expression, activation of the PHD oxygen-sensing mechanism was found to preserve myocyte viability in the face of metabolic inhibition with cyanide and 2-deoxyglucose. These results point to a key role for the PHD pathway in the phenotypic changes that are observed in a hypoxic myocyte and may suggest a strategy to pharmacologically induce protection in heart.     

Thermal Stabilization of Dihydrofolate Reductase Using Monte Carlo Unfolding Simulations and Its Functional Consequences

Design of proteins with desired thermal properties is important for scientific and biotechnological applications. Here we developed a theoretical approach to predict the effect of mutations on protein stability from non-equilibrium unfolding simulations. We establish a relative measure based on apparent simulated melting temperatures that is independent of simulation length and, under certain assumptions, proportional to equilibrium stability, and we justify this theoretical development with extensive simulations and experimental data. Using our new method based on all-atom Monte-Carlo unfolding simulations, we carried out a saturating mutagenesis of Dihydrofolate Reductase (DHFR), a key target of antibiotics and chemotherapeutic drugs. The method predicted more than 500 stabilizing mutations, several of which were selected for detailed computational and experimental analysis. We find a highly significant correlation of r = 0.65–0.68 between predicted and experimentally determined melting temperatures and unfolding denaturant concentrations for WT DHFR and 42 mutants. The correlation between energy of the native state and experimental denaturation temperature was much weaker, indicating the important role of entropy in protein stability. The most stabilizing point mutation was D27F, which is located in the active site of the protein, rendering it inactive. However for the rest of mutations outside of the active site we observed a weak yet statistically significant positive correlation between thermal stability and catalytic activity indicating the lack of a stability-activity tradeoff for DHFR. By combining stabilizing mutations predicted by our method, we created a highly stable catalytically active E. coli DHFR mutant with measured denaturation temperature 7.2°C higher than WT. Prediction results for DHFR and several other proteins indicate that computational approaches based on unfolding simulations are useful as a general technique to discover stabilizing mutations.