Non-native Hydrogen Bonds Research Articles

Molecular Details of a Coupled Binding and Folding Reaction between the Amyloid Precursor Protein and a Folded Domain.

Intrinsically disordered regions in proteins often function as binding motifs in protein–protein interactions. The mechanistic aspects and molecular details of such coupled binding and folding reactions, which involve formation of multiple noncovalent bonds, have been broadly studied theoretically, but experimental data are scarce. Here, using a combination of protein semisynthesis to incorporate phosphorylated amino acids, backbone amide-to-ester modifications, side chain substitutions, and binding kinetics, we examined the interaction between the intrinsically disordered motif of amyloid precursor protein (APP) and the phosphotyrosine binding (PTB) domain of Mint2. We show that the interaction is regulated by a self-inhibitory segment of the PTB domain previously termed ARM. The helical ARM linker decreases the association rate constant 30-fold through a fast pre-equilibrium between an open and a closed state. Extensive side chain substitutions combined with kinetic experiments demonstrate that the rate-limiting transition state for the binding reaction is governed by native and non-native hydrophobic interactions and hydrogen bonds. Hydrophobic interactions were found to be particularly important during crossing of the transition state barrier. Furthermore, linear free energy relationships show that the overall coupled binding and folding reaction involves cooperative formation of interactions with roughly 30% native contacts formed at the transition state. Our data support an emerging picture of coupled binding and folding reactions following overall chemical principles similar to those of folding of globular protein domains but with greater malleability of ground and transition states.

Curcumin Inhibits the Primary Nucleation of Amyloid-Beta Peptide: A Molecular Dynamics Study.

The amyloid plaques are a key hallmark of neurodegenerative diseases such as Alzheimer’s disease and Parkinson’s disease. Amyloidogenesis is a complex long-lasting multiphase process starting with the formation of nuclei of amyloid peptides: a process assigned as a primary nucleation. Curcumin (CU) is a well-known inhibitor of the aggregation of amyloid-beta (Aβ) peptides. Even more, CU is able to disintegrate preformed Aβ firbils and amyloid plaques. Here, we simulate by molecular dynamics the primary nucleation process of 12 Aβ peptides and investigate the effects of CU on the process. We found that CU molecules intercalate among the Aβ chains and bind tightly to them by hydrogen bonds, hydrophobic, π–π, and cation–π interactions. In the presence of CU, the Aβ peptides form a primary nucleus of a bigger size. The peptide chains in the nucleus become less flexible and more disordered, and the number of non-native contacts and hydrogen bonds between them decreases. For comparison, the effects of the weaker Aβ inhibitor ferulic acid (FA) on the primary nucleation are also examined. Our study is in good agreement with the observation that taken regularly, CU is able to prevent or at least delay the onset of neurodegenerative disorders.

The molecular mechanism of protein denaturation in supercritical CO2: The role of exposed lysine residues is explored

Although there are many interests to use proteins in supercritical CO2, protein destabilization in this solvent is the major drawback. This is because of the non-natural CO2 environment which imposes essential instability on the protein structure. Due to the limitations of the experimental techniques to obtain in situ molecular levels information about protein microstructure at high pressures of the supercritical state, the molecular level reasons of the protein destabilization remains controversial and relatively unknown in supercritical CO2. While there is a consensus about the probable roles of surface lysine residues in destabilization of proteins in supercritical CO2, no clear mechanism is presented about the way that these amino acids impose instabilities on the protein conformation. Here, using molecular dynamics simulations, we extensively examined the structural behavior of the surface lysine residues of chymotrypsin inhibitor 2 as a model protein in supercritical CO2. In this way, a new molecular mechanism is proposed which states that the inappropriate interactions of polar side chain of lysine residues with the non-polar nature of supercritical CO2 enforce these exposed residues to return from the solvent bulk to the protein exterior and also to form new non-native hydrogen bonds with the acidic residues such as glutamic acid. This induces some surface regions to be very labile and dynamics. These regions are the starting sequence in the process of protein denaturation in supercritical CO2. Denaturation is extended by collapsing the structure due to the formation of multiple non-native interactions between different residues. To evaluate the proposed mechanism, surface exposed lysine residues were replaced to leucine amino acids with the hydrophobic side chains. Interestingly, the protein remains intact by this replacement even at 200 bar in supercritical CO2. However, replacing of lysine residues by negatively charged glutamic acid and also by glycine residues cannot prevent protein denaturation process. These molecular levels explorations not only address some unknown aspects of protein denaturation mechanism in supercritical CO2, but it also can essentially help to design new CO2-tolerant proteins which can be used in popular processes such as biocatalysis in supercritical CO2.

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

TAR DNA-binding protein 43 (TDP-43) inclusions have been found in Amyotrophic lateral sclerosis (ALS) and several other neurodegenerative diseases. Many studies suggest the involvement of RNA recognition motifs (RRMs) in TDP-43 proteinopathy. To elucidate the structural stability and the unfolding dynamics of RRMs, we have carried out atomistic molecular dynamics simulations at two different temperatures (300 and 500 K). The simulations results indicate that there are distinct structural differences in the unfolding pathway between the two domains and RRM1 unfolds faster than RRM2 in accordance with the lower thermal stability found experimentally. The unfolding behaviors of secondary structures showed that the α-helix was more stable than β-sheet and structural rearrangements of β-sheets results in formation of additional α-helices. At higher temperature, RRM1 exhibit increased overall flexibility and unfolding than RRM2. The temperature-dependent free energy landscapes consist of multiple metastable states stabilized by non-native contacts and hydrogen bonds in RRM2, thus rendering the RRM2 more prone to misfolding. The structural rearrangements of RRM2 could lead to aberrant protein–protein interactions that may account for enhanced aggregation and toxicity of TDP-43. Our analysis, thus identify the structural and thermodynamic characteristics of the RRMs of TDP-43, which will serve to uncover molecular mechanisms and driving forces in TDP-43 misfolding and aggregation.

Temperature-Induced Misfolding in Prion Protein: Evidence of Multiple Partially Disordered States Stabilized by Non-Native Hydrogen Bonds.

The structural basis of pathways of misfolding of a cellular prion (PrPC) into the toxic scrapie form (PrPSC) and identification of possible intermediates (e.g., PrP*) still eludes us. In this work, we have used a cumulative ∼65 μs of replica exchange molecular dynamics simulation data to construct the conformational free energy landscapes and capture the structural and thermodynamic characteristics associated with various stages of the thermal denaturation process in human prion protein. The temperature-dependent free energy surfaces consist of multiple metastable states stabilized by non-native contacts and hydrogen bonds, thus rendering the protein prone to misfolding. We have been able to identify metastable conformational states with high β-content (∼30-40%) and low α-content (∼10-20%) that might be precursors of PrPSC oligomer formation. These conformations also involve participation of the unstructured N-terminal domain, and its role in misfolding has been investigated. All the misfolded or partially unfolded states are quite compact in nature despite having large deviations from the native structure. Although the number of native contacts decreases dramatically at higher temperatures, the radius of gyration and number of intraprotein hydrogen bonds and contacts remain relatively unchanged, leading to stabilization of the misfolded conformations by non-native interactions. Our results are in good agreement with the established view that the C-terminal regions of the second and third helices (H2 and H3, respectively) of mammal prions might be the Achilles heels of their stability, while separation of B1-H1-B2 and H2-H3 domains seems to play a key role, as well.

Elucidating Important Sites and the Mechanism for Amyloid Fibril Formation by Coarse-Grained Molecular Dynamics.

Fibrils formed by the β-amyloid (Aβ) peptide play a central role in the development of Alzheimer's disease. In this study, the principles governing their growth and stability are investigated by analyzing canonical and replica-exchange molecular dynamics trajectories of Aβ(9-40) fibrils. In particular, an unstructured monomer was allowed to interact freely with an Aβ fibril template. Trajectories were generated with the coarse-grained united-residue force field, and one- and two-dimensional free-energy landscapes (FELs) along the backbone virtual-bond angle θ and backbone virtual-bond-dihedral angle γ of each residue and principal components, respectively, were analyzed. Also, thermal unbinding (unfolding) of an Aβ peptide from the fibril template was investigated. These analyses enable us to illustrate the entire process of Aβ fibril elongation and to elucidate the key residues involved in it. Several different pathways were identified during the search for the fibril conformation by the monomer, which finally follows a dock-lock mechanism with two distinct locking stages. However, it was found that the correct binding, with native hydrogen bonds, of the free monomer to the fibril template at both stages is crucial for fibril elongation. In other words, if the monomer is incorrectly bound (with nonnative hydrogen bonds) to the fibril template during the first "docking" stage, it can remain attached to it for a long time before it dissociates and either attempts a different binding or allows another monomer to bind. This finding is consistent with an experimentally observed "stop-and-go" mechanism of fibril growth.

Dynamics of Seeded Aβ40-Fibril Growth from Atomistic Molecular Dynamics Simulations: Kinetic Trapping and Reduced Water Mobility in the Locking Step.

Filamentous β-amyloid aggregates are crucial for the pathology of Alzheimer's disease. Despite the tremendous biomedical importance, the molecular pathway of growth propagation is not completely understood and remains challenging to investigate by simulations due to the long time scales involved. Here, we apply extensive all-atom molecular dynamics simulations in explicit water to obtain free energy profiles and kinetic information from position-dependent diffusion profiles for three different Aβ9-40-growth processes: fibril elongation by single monomers at the structurally unequal filament tips and association of larger filament fragments. Our approach provides insight into the molecular steps of the kinetic pathway and allows close agreement with experimental binding free energies and macroscopic growth rates. Water plays a decisive role, and solvent entropy is identified as the main driving force for assembly. Fibril growth is disfavored energetically due to cancellation of direct peptide-peptide interactions and solvation effects. The kinetics of growth is consistent with the characteristic dock/lock mechanism, and docking is at least 2 orders of magnitude faster. During initial docking, interactions are mediated by transient non-native hydrogen bonds, which efficiently catch the incoming monomer or fragment already at separations of about 3 nm. In subsequent locking, the dynamics is much slower due to formation of kinetically trapped conformations caused by long-lived non-native hydrogen bonds. Fibril growth additionally requires collective motion of water molecules to create a dry binding interface. Fibril growth is further retarded due to reduced mobility of the involved hydration water, evident from a 2-fold reduction of the diffusion coefficient.

Free energy landscape of activation in a signalling protein at atomic resolution.

The interconversion between inactive and active protein states, traditionally described by two static structures, is at the heart of signaling. However, how folded states interconvert is largely unknown due to the inability to experimentally observe transition pathways. Here we explore the free energy landscape of the bacterial response regulator NtrC by combining computation and NMR, and discover unexpected features underlying efficient signaling. We find that functional states are defined purely in kinetic and not structural terms. The need of a well-defined conformer, crucial to the active state, is absent in the inactive state, which comprises a heterogeneous collection of conformers. The transition between active and inactive states occurs through multiple pathways, facilitated by a number of nonnative transient hydrogen bonds, thus lowering the transition barrier through both entropic and enthalpic contributions. These findings may represent general features for functional conformational transitions within the folded state.

The transition mechanism of DNA overstretching: a microscopic view using molecular dynamics.

The overstretching transition in torsionally unconstrained DNA is studied by means of atomistic molecular dynamics simulations. The free-energy profile as a function of the length of the molecule is determined through the umbrella sampling technique providing both a thermodynamic and a structural characterization of the transition pathway. The zero-force free-energy profile is monotonic but, in accordance with recent experimental evidence, becomes two-state at high forces. A number of experimental results are satisfactorily predicted: (i) the entropic and enthalpic contributions to the free-energy difference between the basic (B) state and the extended (S) state; (ii) the longitudinal extension of the transition state and (iii) the enthalpic contribution to the transition barrier. A structural explanation of the experimental finding that overstretching is a cooperative reaction characterized by elementary units of approximately 22 base pairs is found in the average distance between adenine/thymine-rich regions along the molecule. The overstretched DNA adopts a highly dynamical and structurally disordered double-stranded conformation which is characterized by residual base pairing, formation of non-native intra-strand hydrogen bonds and effective hydrophobic screening of apolar regions.

The trehalose coating effect on the internal protein dynamics

(15)N and (13)C NMR experiments were applied to conduct a comparative study of a cold shock protein (Csp) in two states-lyophilized powder and a protein embedded in a glassy trehalose matrix. Both samples were studied at various levels of rehydration. The experiments used (measuring relaxation rates R(1) and R(1ρ), motionally averaged dipolar couplings and solid state exchange method detecting reorientation of the chemical shift anisotropy tensor) allow obtaining abundant information on the protein structural features and internal motions in a range of correlation times from nanoseconds to seconds. The main results are: (a) the trehalose coating makes the protein structure more native in comparison with the dehydrated lyophilized powder, however, trehalose still cannot remove all non-native hydrogen bonds which are present in a dehydrated protein; (b) trehalose has an appreciable effect on the internal dynamics: the motion of the backbone N-H groups in the nanosecond and microsecond time scales becomes slower while the motional amplitude remains constant; (c) upon adding water to the Csp-trehalose mixture, water molecules accumulate around proteins forming a layer between the protein surface and the trehalose matrix. The protein dynamics become faster, however, not as fast as in the fully hydrated state; (d) the hydration response of dynamics of the NH and CH(CH(2)) groups in a protein is qualitatively different: upon increasing protein hydration, the correlation times of the N-H motions become shorter and the amplitude remains stable, and for CH(CH(2)) groups the motional amplitude increases and the correlation times do not change. This can be explained by a different ability of the NH and CH(CH(2)) groups to form hydrogen bonds.

Improvement of Structure-Based Potentials for Protein Folding by Native and Nonnative Hydrogen Bonds

Pure Gō models (where every native interaction equally stabilizes the folded state) have widely proved their convenience in the computational investigation of protein folding. However, a chemistry-based description of the real interactions also provides a desirable tune in the analysis of the folding process, and thus some hybrid Gō potentials that combine both aspects have been proposed. Among all the noncovalent interactions that contribute to protein folding, hydrogen bonds are the only ones with a partial covalent character. This feature makes them directional and, thus, more difficult to model as part of the coarse-grained descriptions that are typically employed in Gō models. Thanks to a simplified but rigorous representation of backbone hydrogen bonds that we have recently proposed, we present in this article a combined potential (Gō + backbone hydrogen bond) to study the thermodynamics of protein folding in the frame of very simple simulation models. We show that the explicit inclusion of hydrogen bonds leads to a systematic improvement in the description of protein folding. We discuss a representative set of examples (from two-state folders to downhill proteins, with different types of native structures) that reveal a relevant agreement with experimental data.

Choreographing an enzyme's dance.

While ground state structures combined with chemical tools and enzyme kinetics deliver useful information on possible chemical mechanisms of enzyme catalysis, they do not unravel the finely balanced energy inventory to explain the impressive rate enhancement of enzymes. For this goal, a complete description of enzyme catalysis in the form of an energy landscape is needed. Since the rate of catalysis is determined by the climb over a sequence of energy barriers, we focus here on the critical question of transition pathways. A combination of time-resolved NMR and simulation deliver a glimpse into how proteins can so efficiently move within the ensemble of the native conformations while avoiding unfolding during that journey. The loss of energy due to breakage of native contacts is compensated by non-native transient hydrogen bonds during the transition thereby 'holding on' to the energy until the new native contacts form that define the alternate functional state. The use of kinetic isotope effects (KIE) to study the chemical step show that coordinated atomic fluctuations of the protein component dictate the probability of 'correct' distance and orientation, due to its extreme sensitivity to distance. The examples here stress the point that highly choreographed conformational sampling together with chemical integrity is a prerequisite for efficient enzyme catalysis.

Folding of a Small RNA Hairpin Based on Simulation with Replica Exchange Molecular Dynamics

The folding of a small RNA tetraloop hairpin is studied based on intensive molecular dynamics simulation, aiming to understand the folding mechanism of this small and fast RNA folder. Our results showed that this RNA hairpin has very complicated folding behavior in spite of its small size. It is found that the folding transition has low cooperativity. Instead of a two-state folding, four major states are observed, including the native state, the intermediate, the unfolded state, and the misfolded state. The misfolded state is mainly stabilized by the non-native hydrogen bonds, and is more compact. Two potential folding pathways, in which two basepairs formed with different order, are observed, and the pathway with the inboard basepair formed before the terminal one is much more favorable, and dominates the folding of the RNA hairpin.

Roles of non-native hydrogen-bonding interaction in helix-coil transition of a single polypeptide as revealed by comparison between Gō-like and non-Gō models

To explore the role of non-native interactions in the helix-coil transition, a detailed comparison between a Gō-like model and a non-Gō model has been performed via lattice Monte Carlo simulations. Only native hydrogen bonding interactions occur in the Gō-like model, and the non-native ones with sequence interval more than 4 is also included into the non-Gō model. Some significant differences between the results from those two models have been found. The non-native hydrogen bonds were found most populated at temperature around the helix-coil transition. The rearrangement of non-native hydrogen bonds into native ones in the formation of alpha-helix leads to the increase of susceptibility of chain conformation, and even two peaks of susceptibility of radius of gyration versus temperature exist in the case of non-Gō model for a non-short peptide, while just a single peak exists in the case of Gō model for a single polypeptide chain with various chain lengths. The non-native hydrogen bonds have complicated the temperature-dependence of Zimm-Bragg nucleation constant. The increase of relative probability of non-native hydrogen bonding for long polypeptide chains leads to non-monotonous chain length effect on the transition temperature.

Protein mechanical unfolding: Importance of non-native interactions

Mechanical unfolding of the fourth domain of Distyostelium discoideum filamin (DDFLN4) was studied by all-atom molecular dynamics simulations, using the GROMOS96 force field 43a1 and the simple point charge explicit water solvent. Our study reveals an important role of non-native interactions in the unfolding process. Namely, the existence of a peak centered at the end-to-end extension DeltaR approximately 22 nm in the force-extension curve is associated with breaking of non-native hydrogen bonds. Such a peak has been observed in experiments but not in Go models, where non-native interactions are neglected. We predict that an additional peak occurs at DeltaR approximately 2 nm using not only GROMOS96 force field 43a1 but also Amber 94 and OPLS force fields. This result would stimulate further experimental studies on elastic properties of DDFLN4.

Non-Native Interactions Are Critical for Mechanical Strength in PKD Domains

SummaryExperimental observation has led to the commonly held view that native state protein topology is the principle determinant of mechanical strength. However, the PKD domains of polycystin-1 challenge this assumption: they are stronger than predicted from their native structure. Molecular dynamics simulations suggest that force induces rearrangement to an intermediate structure, with nonnative hydrogen bonds, that resists unfolding. Here we test this hypothesis directly by introducing mutations designed to prevent formation of these nonnative interactions. We find that these mutations, which only moderately destabilize the native state, reduce the mechanical stability dramatically. The results demonstrate that nonnative interactions impart significant mechanical stability, necessary for the mechanosensor function of polycystin-1. Remarkably, such nonnative interactions result from force-induced conformational change: the PKD domain is strengthened by the application of force.

During Transitions Proteins Make Fleeting Bonds

How do proteins efficiently and precisely shift from one conformation to another? Gardino et al. (2009) show that transient hydrogen bonds are critical to the conformational transition of the nitrogen regulatory protein NtrC between its native state and its active state.

Transient Non-native Hydrogen Bonds Promote Activation of a Signaling Protein

Phosphorylation is a common mechanism for activating proteins within signaling pathways. Yet, the molecular transitions between the inactive and active conformational states are poorly understood. Here we quantitatively characterize the free-energy landscape of activation of a signaling protein, nitrogen regulatory protein C (NtrC), by connecting functional protein dynamics of phosphorylation-dependent activation to protein folding and show that only a rarely populated, pre-existing active conformation is energetically stabilized by phosphorylation. Using nuclear magnetic resonance (NMR) dynamics, we test an atomic scale pathway for the complex conformational transition, inferred from molecular dynamics simulations (Lei et al., 2009). The data show that the loss of native stabilizing contacts during activation is compensated by non-native transient atomic interactions during the transition. The results unravel atomistic details of native-state protein energy landscapes by expanding the knowledge about ground states to transition landscapes.

Segmented Transition Pathway of the Signaling Protein Nitrogen Regulatory Protein C

Recent advances in experimental methods provide increasing evidence that proteins sample the conformational substates that are important for function in the absence of their ligands. An example is the receiver domain of nitrogen regulatory protein C, a member of the phosphorylation-mediated signaling family of “two-component systems.” The receiver domain of nitrogen regulatory protein C samples both inactive conformation and the active conformation before phosphorylation. Here we determine a possible pathway of interconversion between the active state and the inactive state by targeted molecular dynamics simulations and quasi-harmonic analysis; these methods are used because the experimental conversion rate is in the high microsecond range, longer than those that are easily accessible to atomistic molecular dynamics simulations. The calculated pathway is found to be composed of four consecutive stages described by different progress variables. The lowest quasi-harmonic principal components from unbiased molecular dynamics simulations on the active state correspond to the first stage, but not to the subsequent stages of the transition. The targeted molecular dynamics pathway suggests that several transient nonnative hydrogen bonds may facilitate the transition.

Understanding β-Hairpin Formation: Computational Studies for Three Different Hairpins

We have studied the folding mechanism of <TEX>$\beta$</TEX> -hairpins in the proteins 1GB1, 3AIT and 1A2P by conducting unfolding simulations at moderately high temperatures. The analysis of trajectories obtained from molecular dynamics simulations in explicit aqueous solution suggests that the positions of the hydrophobic core residues lead to subtle differences in the details of folding dynamics. However, the folding of three different hairpins can be explained by a unified mechanism that is a blend of the hydrogen-bond-centric and the hydrophobiccentric models. The initial stage of <TEX>$\beta$</TEX>-hairpin folding involves various partially folded intermediate structures which are stabilized by both the van der Waals interactions of hydrophobic core residues and the electrostatic interactions of non-native hydrogen bonds. The native structure is obtained by forming native contacts in the final tune-up process. Depending on the relative positions of the hydrophobic residues, the actual mechanism of hairpi n folding may or may not exhibit well-defined intermediates.