Insight into the early stages of thermal unfolding of peanut agglutinin by molecular dynamics simulations

Peanut agglutinin is a homotetrameric nonglycosylated protein. The protein has a unique open quaternary structure. Molecular dynamics simulations have been employed follow the atomistic details of its unfolding at different temperatures. The early events of the deoligomerization of the protein have been elucidated in the present study. Simulation trajectories of the monomer as well as those of the tetramer have been compared and the tetramer is found to be substantially more stable than its monomeric counterpart. The tetramer shows retention of most of its.. secondary structure but considerable loss of the tertiary structure at high temperature. e generation of a This observation impies the molten globule-like intermediate in the later stages of deoligomerization. The quaternary structure of the protein has weakened to a large extent, but none of the subunits are separated. In addition, the importance of the metal-binding to the stability of the protein structure has also been investigated. Binding of the metal ions not only enhances the local stability of the metal-ion binding loop, but also imparts a global stability to the overall structure. The dynamics of different interfaces vary significantly as probed through interface clusters. The differences are substantially enhanced at higher temperatures. The dynamics and the stability of the interfaces have been captured mainly by cluster analysis, which has provided detailed information on the thermal deoligomerization of the protein.

Unfolding of Plasmodium falciparum Triosephosphate Isomerase in Urea and Guanidinium Chloride: Evidence for a Novel Disulfide Exchange Reaction in a Covalently Cross-Linked Mutant†

The conformational stability of Plasmodium falciparum triosephosphate isomerase (TIMWT) enzyme has been investigated in urea and guanidinium chloride (GdmCl) solutions using circular dichroism, fluorescence, and size-exclusion chromatography. The dimeric enzyme is remarkably stable in urea solutions. It retains considerable secondary, tertiary, and quaternary structure even in 8 M urea. In contrast, the unfolding transition is complete by 2.4 M GdmCl. Although the secondary as well as the tertiary interactions melt before the perturbation of the quaternary structure, these studies imply that the dissociation of the dimer into monomers ultimately leads to the collapse of the structure, suggesting that the interfacial interactions play a major role in determining multimeric protein stability. The C-m(urea)/C-m(GdmCl) ratio (where C-m is the concentration of the denaturant required at the transition midpoint) is unusually high for triosephosphate isomerase as compared to other monomeric and dimeric proteins. A disulfide crosslinked mutant protein (Y74C) engineered to form two disulfide cross-links across the interface (13-74') and (13'-74) is dramatically destablized in urea. The unfolding transition is complete by 6 M urea and involves a novel mechanism of dimer dissociation through intramolecular thiol-disulfide exchange.

Unfolding studies on soybean agglutinin and concanavalin a tetramers: A comparative account

The unfolding pathway of two very similar tetrameric legume lectins soybean agglutinin (SBA) and Concanavalin A ( ConA) were determined using GdnCl-induced denaturation. Both proteins displayed a reversible two-state unfolding mechanism. The analysis of isothermal denaturation data provided values for conformational stability of the two proteins. It was found that the DeltaG of unfolding of SBA was much higher than ConA at all the temperatures at which the experiments were done. ConA had a T-g 18 degreesC less than SBA. The higher conformational stability of SBA in comparison to ConA is largely due to substantial differences in their degrees of subunit interactions. Ionic interactions at the interface of the two proteins especially at the noncanonical interface seem to play a significant role in the observed stability differences between these two proteins. Furthermore, SBA is a glycoprotein with a GlcNac(2)Man(9) chain attached to Asn-75 of each subunit. The sugar chain in SBA lies at the noncanonical interface of the protein, and it is found to interact with the amino acid residues in the adjacent noncanonical interface. These interactions further stabilize SBA with respect to ConA, which is not glycosylated.

Unfolding of Plasmodium falciparum triosephosphate isomerase in urea and guanidinium chloride solutions. Evidence for a novel disulfide exchange reaction in a covalently crosslinked mutant

The conformational stability of Plasmodium falciparum triosephosphate isomerase (TIMWT) enzyme has been investigated in urea and guanidinium chloride (GdmCl) solutions using circular dichroism, fluorescence, and size-exclusion chromatography. The dimeric enzyme is remarkably stable in urea solutions. It retains considerable secondary, tertiary, and quaternary structure even in 8 M urea. In contrast, the unfolding transition is complete by 2.4 M GdmCl. Although the secondary as well as the tertiary interactions melt before the perturbation of the quaternary structure, these studies imply that the dissociation of the dimer into monomers ultimately leads to the collapse of the structure, suggesting that the interfacial interactions play a major role in determining multimeric protein stability. The Cm(urea)/Cm(GdmCl) ratio (where Cm is the concentration of the denaturant required at the transition midpoint) is unusually high for triosephosphate isomerase as compared to other monomeric and dimeric proteins. A disulfide cross-linked mutant protein (Y74C) engineered to form two disulfide cross-links across the interface (13-74‘) and (13‘-74) is dramatically destablized in urea. The unfolding transition is complete by 6 M urea and involves a novel mechanism of dimer dissociation through intramolecular thiol−disulfide exchange.

GATA Factors Regulate Inner Ear Development and Midbrain Neurogenesis

The zinc-finger transcription factors GATA2 and GATA3 in vertebrates belong to the six-member family that are essential regulators in the development of various organs. The aim of this study was to gain new information of the roles of GATA2 and GATA3 in inner ear morphogenesis and of the function of GATA2 in neuronal fate specification in the midbrain using genetically modified mouse and chicken embryos as models. A century ago the stepwise process of inner ear epithelial morphogenesis was described, but the molecular players regulating the cellular differentiation of the otic epithelium are still not fully resolved. This study provided novel data on GATA factor roles in several developmental processes during otic development. The expression analysis in chicken suggested that GATA2 and GATA3 possess redundant roles during otic cup and vesicle formation, but complementary cell-type specific functions during vestibular and cochlear morphogenesis. The comparative analysis between mouse and chicken Gata2 and Gata3 expression revealed many conserved aspects, especially during later stages of inner ear development, while the expression was more divergent at early stages. Namely, expression of both Gata genes was initiated earlier in chicken than mouse otic epithelium relative to the morphogenetic stages. Likewise, important differences concerning Gata3 expression in the otic cup epithelium were detected between mouse and chicken, suggesting that distinct molecular mechanisms regulate otic vesicle closure in different vertebrate species. Temporally distinct Gata2 and Gata3 expression was also found during otic ganglion formation in mouse and chicken. Targeted inactivation of Gata3 in mouse embryos caused aberrant morphology of the otic vesicle that in severe cases was disrupted into two parts, a dorsal and a ventral vesicle. Detailed analyses of Gata3 mutant embryos unveiled a crucial role for GATA3 in the initial inner ear morphogenetic event, the invagination of the otic placode. A large-scale comparative expression analysis suggested that GATA3 could control cell adhesion and motility in otic epithelium, which could be important for early morphogenesis. GATA3 was also identified as the first factor to directly regulate Fgf10 expression in the otic epithelium and could thus influence the development of the semicircular ducts. Despite the serious problems in the early inner ear development, the otic sensory fate establishment and some vestibular hair cell differentiation was observable in pharmacologically rescued Gata3-/- embryos. Cochlear sensory differentiation was, however, completely blocked so that no auditory hair cells were detected. In contrast to the early morphogenetic phenotype in Gata3-/- mutants, conditional inactivation of Gata2 in mouse embryos resulted in a relatively late growth defect of the three semicircular ducts. GATA2 was required for the proliferation of the vestibular nonsensory epithelium to support growing of the three ducts. Concurrently, with the role in epithelial semicircular ducts, GATA2 was also required for the mesenchymal cell clearance from the vestibular perilymphatic region between the membranous labyrinth and bony capsule. The gamma-aminobutyric acid-secreting (GABAergic) neurons in the midbrain are clinically relevant since they contribute to fear, anxiety, and addiction regulation. The molecular mechanisms regulating the GABAergic neuronal development, however, are largely unknown. Using tissue-specific mutagenesis in mice, GATA2 was characterized as a critical determinant of the GABAergic neuronal fate in the midbrain. In Gata2-deficient mouse midbrain, GABAergic neurons were not produced, instead the Gata2-mutant cells acquired a glutamatergic neuronal phenotype. Gain-of-function experiments in chicken also revealed that GATA2 was sufficient to induce GABAergic differentiation in the midbrain.

Differential stability of beta-sheets and alpha-helices in beta-lactamase: a high temperature molecular dynamics study of unfolding intermediates

Beta-Lactamase, which catalyzes beta-lactam antibiotics, is prototypical of large alpha/beta proteins with a scaffolding formed by strong noncovalent interactions. Experimentally, the enzyme is well characterized, and intermediates that are slightly less compact and having nearly the same content of secondary structure have been identified in the folding pathway. In the present study, high temperature molecular dynamics simulations have been carried out on the native enzyme in solution. Analysis of these results in terms of root mean square fluctuations in cartesian and [phi, psi] space, backbone dihedral angles and secondary structural hydrogen bonds forms the basis for an investigation of the topology of partially unfolded states of beta-lactamase. A differential stability has been observed for alpha-helices and beta-sheets upon thermal denaturation to putative unfolding intermediates. These observations contribute to an understanding of the folding/unfolding processes of beta-lactamases in particular, and other alpha/beta proteins in general.

The stretching of single poly-ubiquitin molecules: Static versus dynamic disorder in the non-exponential kinetics of chain unfolding

Static disorder has recently been implicated in the non-exponential kinetics of the unfolding of single molecules of poly-ubiquitin under a constant force Kuo, Garcia-Manyes, Li, Barel, Lu, Berne, Urbakh, Klafter, and Fernandez, Proc. Natl. Acad. Sci. U. S. A. 107, 11336 (2010)]. In the present paper, it is suggested that dynamic disorder may provide a plausible, alternative description of the experimental observations. This suggestion is made on the basis of a model in which the barrier to chain unfolding is assumed to be modulated by a control parameter r that evolves in a parabolic potential under the action of fractional Gaussian noise according to a generalized Langevin equation. The treatment of dynamic disorder within this model is pursued using Zwanzig's indirect approach to noise averaging Acc. Chem. Res. 23, 148 (1990)]. In conjunction with a self-consistent closure scheme developed by Wilemski and Fixman J. Chem. Phys. 58, 4009 (1973); ibid. 60, 866 (1974)], this approach eventually leads to an expression for the chain unfolding probability that can be made to fit the corresponding experimental data very closely. (C) 2011 American Institute of Physics.

Structural studies of protein unfolding

Methods for macromolecular structure determination (NMR and crystallography) are now being used to get structural information on partially folded and unfolded states of proteins. These techniques, in combination with proton hydrogen exchange studies are powerful tools to extract information on non-native states of proteins. This review discusses progress In this area of protein folding.

Energy transfer efficiency distributions in polymers in solution during folding and unfolding

Distribution of fluorescence resonance energy transfer (FRET) efficiency between the two ends of a Lennard-Jones polymer chain both at equilibrium and during folding and unfolding has been calculated, for the first time, by Brownian dynamics simulations. The distribution of FRET efficiency becomes bimodal during folding of the extended state subsequent to a temperature quench, with the width of the distribution for the extended state broader than that for the folded state. The reverse process of unfolding subsequent to a upward temperature jump shows different characteristics. The distributions show significant viscosity dependence which can be tested against experiments.

Denaturant mediated unfolding of both native and molten globule states of maltose binding protein are accompanied by large [Delta]Cp's

Maltose binding protein (MBP) is a large, monomeric two domain protein containing 370 amino acids. In the absence of denaturant at neutral pH, the protein is in the native state, while at pH 3.0 it forms a molten globule. The molten globule lacks a tertiary circular dichroism signal but has secondary structure similar to that of the native state. The molten globule binds 8-anilino-1-naphthalene sulfonate (ANS). The unfolding thermodynamics of MBP at both pHs were measured by carrying out a series of isothermal urea melts at temperatures ranging from 274–329 K. At 298 K, values of [Delta]G°, [Delta]Cp, and Cm were 3.1 ± 0.2 kcal mol−1, 5.9 ± 0.8 kcal mol−1 K−1 (15.9 cal (mol-residue)−1 K−1), and 0.8 M, respectively, at pH 3.0 and 14.5 ± 0.4 kcal mol−1, 8.3 ± 0.7 kcal mol−1 K−1 (22.4 kcal (mol-residue)−1 K−1), and 3.3 M, respectively, at pH 7.1. Guanidine hydrochloride denaturation at pH 7.1 gave values of [Delta]G° and [Delta]Cp similar to those obtained with urea. The m values for denaturation are strongly temperature dependent, in contrast to what has been previously observed for small globular proteins. The value of [Delta]Cp per mol-residue for the molten globule is comparable to corresponding values of [Delta]Cp for the unfolding of typical globular proteins and suggests that it is a highly ordered structure, unlike molten globules of many small proteins. The value of [Delta]Cp per mol-residue for the unfolding of the native state is among the highest currently known for any protein.

Ligand-modulated Parallel Mechanical Unfolding Pathways of Maltose-binding Proteins

Protein folding and unfolding are complex phenomena, and it is accepted that multidomain proteins generally follow multiple pathways. Maltose-binding protein (MBP) is a large (a two-domain, 370-amino acid residue) bacterial periplasmic protein involved in maltose uptake. Despite the large size, it has been shown to exhibit an apparent two-state equilibrium unfolding in bulk experiments. Single-molecule studies can uncover rare events that are masked by averaging in bulk studies. Here, we use single-molecule force spectroscopy to study the mechanical unfolding pathways of MBP and its precursor protein (preMBP) in the presence and absence of ligands. Our results show that MBP exhibits kinetic partitioning on mechanical stretching and unfolds via two parallel pathways: one of them involves a mechanically stable intermediate (path I) whereas the other is devoid of it (path II). The apoMBP unfolds via path I in 62% of the mechanical unfolding events, and the remaining 38% follow path II. In the case of maltose-bound MBP, the protein unfolds via the intermediate in 79% of the cases, the remaining 21% via path II. Similarly, on binding to maltotriose, a ligand whose binding strength with the polyprotein is similar to that of maltose, the occurrence of the intermediate is comparable (82% via path I) with that of maltose. The precursor protein preMBP also shows a similar behavior upon mechanical unfolding. The percentages of molecules unfolding via path I are 53% in the apo form and 68% and 72% upon binding to maltose and maltotriose, respectively, for preMBP. These observations demonstrate that ligand binding can modulate the mechanical unfolding pathways of proteins by a kinetic partitioning mechanism. This could be a general mechanism in the unfolding of other large two-domain ligand-binding proteins of the bacterial periplasmic space.

Chemical unfolding of chicken villin headpiece in aqueous dimethyl sulfoxide solution: cosolvent concentration dependence, pathway, and microscopic mechanism

Unfolding of a protein often proceeds through partial unfolded intermediate states (PUIS). PUIS have been detected in several experimental and simulation studies. However, complete analyses of transitions between different PUIS and the unfolding trajectory are sparse. To understand such dynamical processes, we study chemical unfolding of a small protein, chicken villin head piece (HP-36), in aqueous dimethyl sulfoxide (DMSO) solution. We carry out molecular dynamics simulations at various solution compositions under ambient conditions. In each concentration, the initial step of unfolding involves separation of two adjacent native contacts, between phenyl alanine residues (11-18 and 7-18). This first step induces, under appropriate conditions, subsequent separation among other hydrophobic contacts, signifying a high degree of cooperativity in the unfolding process. The observed sequence of structural changes in HP-36 on increasing DMSO concentration and the observed sequence of PUIS, are in approximate agreement with earlier simulation results (in pure water) and experimental observations on unfolding of HP-36. Peculiar to water-DMSO mixture, an intervening structural transformation (around 15% of DMSO) in the binary mixture solvent retards the progression of unfolding as composition is increased. This is reflected in a remarkable nonmonotonic composition dependence of RMSD, radius of gyration and the fraction of native contacts. At 30% mole fraction of DMSO, we find the extended randomly coiled structure of the unfolded protein. The molecular mechanism of DMSO induced unfolding process is attributed to the initial preferential solvation of the hydrophobic side chain atoms through the methyl groups of DMSO, followed by the hydrogen bonding of the oxygen atom of DMSO to the exposed backbone NH groups of HP-36.

Solvent Sensitivity of Protein Unfolding: Dynamical Study of Chicken Villin Headpiece Subdomain in Water Ethanol Binary Mixture

We carry out a series of long atomistic molecular dynamics simulations to study the unfolding of a small protein, chicken villin headpiece (HP-36), in water-ethanol (EtOH) binary mixture. The prime objective of this work is to explore the sensitivity of protein unfolding dynamics toward increasing concentration of the cosolvent and unravel essential features of intermediates formed in search of a dynamical pathway toward unfolding. In water ethanol binary mixtures, HP-36 is found to unfold partially, under ambient conditions, that otherwise requires temperature as high as similar to 600 K to denature in pure aqueous solvent. However, an interesting course of pathway is observed to be followed in the process, guided by the formation of unique intermediates. The first step of unfolding is essentially the separation of the cluster formed by three hydrophobic (phenylalanine) residues, namely, Phe-7, Phe-11, and Phe-18, which constitute the hydrophobic core, thereby initiating melting of helix-2 of the protein. The initial steps are similar to temperature-induced unfolding as well as chemical unfolding using DMSO as cosolvent. Subsequent unfolding steps follow a unique path. As water-ethanol shows composition-dependent anomalies, so do the details of unfolding dynamics. With an increase in cosolvent concentration, different partially unfolded intermediates are found to be formed. This is reflected in a remarkable nonmonotonic composition dependence of several order parameters, including fraction of native contacts and protein-solvent interaction energy. The emergence of such partially unfolded states can be attributed to the preferential solvation of the hydrophobic residues by the ethyl groups of ethanol. We further quantify the local dynamics of unfolding by using a Marcus-type theory.