Non-native Hydrogen Bonds Research Articles

Positional dependence of non-native polar mutations on folding of CFTR helical hairpins

Mutations of the cystic fibrosis transmembrane conductance regulator (CFTR) cause CF disease by altering the biosynthesis, maturation, folding and ion conductance of this protein. Our laboratory has focused on expression and structural analysis of the CFTR transmembrane (TM) domains using two-TM segments (i.e., helix–loop–helix constructs) which we term ‘helical hairpins’; these represent the minimal model of tertiary contacts between two helices in a membrane. Previous studies on a library of TM3/4 hairpins of the first CFTR TM domain suggested that introduction of non-native polar residues into TM4 can compromise CFTR function through side chain-side chain H-bonding interactions with native Q207 in TM3 [Choi, M. Y., Cardarelli, L., Therien, A. G., and Deber, C. M. Non-native interhelical hydrogen bonds in the cystic fibrosis transmembrane conductance regulator domain modulated by polar mutations, Biochemistry 43 (2004) 8077–8083]. In the present work, we combine gel shift assays with a series of NMR experiments for comparative structural characterization of the wild type TM3/4 hairpin and its mutants V232D, I231D, Q207N/V232E. Over 95% of the backbone resonances of a 15N, 13C-labelled V232D-TM3/4 construct in the membrane-mimetic environment of perfluorooctanoate (PFO) micelles were successfully assigned, and the presence and boundaries of helical segments within TM3 and TM4 were defined under these conditions. Comparative analysis of 15N and 1H chemical shift variations among HSQC spectra of WT-, V232D-, I231D- and Q207N/V232E-TM3/4 indicated that hairpin conformations vary with the position of a polar mutation (i.e., V232D and I231D vs. WT), but remain similar when hairpins with identically-positioned polar partners are compared (i.e., V232D vs. Q207N-V232E). The overall findings suggest that a polar mutation in a TM helix can potentially distort native interfacial packing determinants in membrane proteins such as CFTR, with consequences that may lead to disease.

Stretching the Immunoglobulin 27 Domain of the Titin Protein: The Dynamic Energy Landscape

Molecular simulations are carried out on the Immunoglobulin 27 domain of the titin protein. The energy landscape is mapped out using an implicit solvent model, and molecular dynamics simulations are run with the solvent explicitly modeled. Stretching a protein is shown to produce a dynamic energy landscape in which the energy minima move in configuration space, change in depth, and are created and destroyed. The connections of these landscape changes to the mechanical unfolding of the Immunoglobulin 27 domain are addressed. Hydrogen bonds break upon stretching by either intrabasin processes associated with the movement of energy minima, or interbasin processes associated with transitions between energy minima. Intrabasin changes are reversible and dominate for flexible interactions, whereas interbasin changes are irreversible and dominate for stiff interactions. The most flexible interactions are Glu-Lys salt bridges, which can act like tethers to bind strands even after all backbone interactions between the strands have been broken. As the protein is stretched, different types of structures become the lowest energy structures, including structures that incorporate nonnative hydrogen bonds. Structures that have flat energy versus elongation profiles become the lowest energy structures at elongations of several Angstroms, and are associated with the unfolding intermediate state observed experimentally.

Complex folding pathways in a simple beta-hairpin.

The determination of the folding mechanisms of proteins is critical to understand the topological change that can propagate Alzheimer and Creutzfeld-Jakobs diseases, among others. The computational community has paid considerable attention to this problem; however, the associated time scale, typically on the order of milliseconds or more, represents a formidable challenge. Ab initio protein folding from long molecular dynamics simulations or ensemble dynamics is not feasible with ordinary computing facilities and new techniques must be introduced. Here we present a detailed study of the folding of a 16-residue beta-hairpin, described by a generic energy model and using the activation-relaxation technique. From a total of 90 trajectories at 300 K, three folding pathways emerge. All involve a simultaneous optimization of the complete hydrophobic and hydrogen bonding interactions. The first two pathways follow closely those observed by previous theoretical studies (folding starting at the turn or by interactions between the termini). The third pathway, never observed by previous all-atom folding, unfolding, and equilibrium simulations, can be described as a reptation move of one strand of the beta-sheet with respect to the other. This reptation move indicates that non-native interactions can play a dominant role in the folding of secondary structures. Furthermore, such a mechanism mediated by non-native hydrogen bonds is not available for study by unfolding and Gō model simulations. The exact folding path followed by a given beta-hairpin is likely to be influenced by its sequence and the solvent conditions. Taken together, these results point to a more complex folding picture than expected for a simple beta-hairpin.

Describing Protein Folding Kinetics by Molecular Dynamics Simulations. 2. Example Applications to Alanine Dipeptide and a β-Hairpin Peptide

In this work we demonstrate the use of a rigorous formalism for the extraction of state-to-state transition functions as a way to study the kinetics of protein folding in the context of a Markov chain. The approach is illustrated by its application to two different systems: a blocked alanine dipeptide in a vacuum and the C-terminal β-hairpin motif from protein G in water. The first system displays some of the desired features of the approach, whereas the second illustrates some of the challenges that must be overcome to apply the method to more complex biomolecular systems. For both example systems, Boltzmann weighted conformations produced by a replica exchange Monte Carlo procedure were used as starting states for kinetic trajectories. The alanine dipeptide displays Markovian behavior in a state space defined with respect to φ−ψ torsion angles. In contrast, Markovian behavior was not observed for the β-hairpin in a state space where all possible native hydrogen bonding patterns were resolved. This may be due to our choice of state definitions or sampling limitations. Furthermore, the use of different criteria for hydrogen bonding results in the apparent observation of different mechanisms from the same underlying data: one set of criteria indicate a zipping type of process, but another indicates more of a collapse followed by almost simultaneous formation of a large number of contacts. Analysis of long-lived states observed during the simulations of the β-hairpin suggests that important aspects of the folding process that are not captured by order parameters in common use include the formation of non-native hydrogen bonds and the degree and nature of salt bridge formation.

High Pressure Refolding of Recombinant Human Growth Hormone from Insoluble Aggregates: STRUCTURAL TRANSFORMATIONS, KINETIC BARRIERS, AND ENERGETICS

Two different types of insoluble, non-native aggregates of recombinant human growth hormone were formed by agitation in buffer or buffer containing 0.75 m guanidine HCl (GdnHCl) and characterized by infrared and second derivative UV spectroscopies. The degree of secondary structural perturbation was greater in the aggregates formed in 0.75 m GdnHCl. Both aggregate types were dissolved and refolded using high hydrostatic pressures in combination with either elevated temperature or non-denaturing levels of guanidine HCl or urea. The effects of a range of temperature, pressure, and chaotrope concentrations were tested and led to optimal conditions that approached 100% yield of native protein. The aggregates formed in 0.75 m GdnHCl required higher concentrations of urea or GdnHCl, or higher temperature or pressure for a yield equivalent to that for aggregates formed in buffer alone. Investigation of the effects of pressure, temperature, and chaotrope on unfolding of rhGH documented that under conditions used for optimal high pressure disaggregation and refolding, the native state is greatly favored thermodynamically (e.g. 25 kJ/mol at 2000 bar and 0.75 m GdnHCl). Dissolution of aggregates under pressure is a kinetically limited process. Comparison of refolding yields in GdnHCl and urea solutions suggest that pressure effects on electrostatic interactions do not dominate pressure effects on disaggregation. We suggest that non-native hydrogen bonds between protein molecules within aggregates of recombinant human growth hormone are responsible for the rate-limiting kinetic barrier in pressure-induced disaggregation.