Backbone Hydrogen Research Articles

Mining protein dynamics from sets of crystal structures using “consensus structures”

In this work, we describe two novel approaches to utilize the dynamic structure information implicitly contained in large crystal structure data sets. The first approach visualizes both consistent as well as variable ligand-induced changes in ligand-bound compared with apo protein crystal structures. For this purpose, information was mined from B-factors and ligand-induced residue displacements in multiple crystal structures, minimizing experimental error and noise. With this approach, the mechanism of action of non-nucleoside reverse transcriptase inhibitors (NNRTIs) as an inseparable combination of distortion of protein dynamics and conformational changes of HIV-1 reverse transcriptase was corroborated (a combination of the previously proposed "molecular arthritis" and "distorted site" mechanisms). The second approach presented here uses "consensus structures" to map common binding features that are present in a set of structures of NNRTI-bound HIV-1 reverse transcriptase. Consensus structures are based on different levels of structural overlap of multiple crystal structures and are used to analyze protein-ligand interactions. The structures are shown to yield information about conserved hydrogen bonding interactions as well as binding-pocket flexibility, shape, and volume. From the consensus structures, a common wild type NNRTI binding pocket emerges. Furthermore, we were able to identify a conserved backbone hydrogen bond acceptor at P236 and a novel hydrophobic subpocket, which are not yet utilized by current drugs. Our methods introduced here reinterpret the atom information and make use of the data variability by using multiple structures, complementing classical 3D structural information of single structures.

Importance of Protein Context on the Unfolding Pathways of β-hairpins

Small β-hairpin peptides as well as three-stranded WW domains have been used as models for the folding of β-hairpins in larger proteins. Previous studies of the FBP28 WW domain proposed that side chain contacts between residues in the strands and not the precise order of backbone hydrogen bond formation guide β-hairpin folding. But how applicable is the folding of model systems, such as FBP28 WW domain, to the folding of β-structure in larger proteins with conventional hydrophobic cores? Here we present multiple unfolding molecular dynamics simulations of three proteins that share a double hairpin motif structurally similar to WW domains: cold shock protein A (CspA), cold shock protein B (CspB) and glucose permease IIA domain.

RNA structural motifs that entail hydrogen bonds involving sugar–phosphate backbone atoms of RNA

The growing number of high-resolution crystal structures of large RNA molecules provides much information for understanding the principles of structural organization of these complex molecules. Several in-depth analyses of nucleobase-centered RNA structural motifs and backbone conformations have been published based on this information, including a systematic classification of base pairs by Leontis and Westhof. However, hydrogen bonds involving sugar-phosphate backbone atoms of RNA have not been analyzed systematically until recently, although such hydrogen bonds appear to be common both in local and tertiary interactions. Here we review some backbone structural motifs discussed in the literature and analyze a set of eight high-resolution multi-domain RNA structures. The analyzed RNAs are highly structured: among 5372 nucleotides in this set, 89% are involved in at least one "long-range" RNA-RNA hydrogen bond, i.e., hydrogen bonds between atoms in the same residue or sequential residues are ignored. These long-range hydrogen bonds frequently use backbone atoms as hydrogen bond acceptors, i.e., OP1, OP2, O2', O3', O4', or O5', or as a donor (2'OH). A surprisingly large number of such hydrogen bonds are found, considering that neither single-stranded nor double-stranded regions will contain such hydrogen bonds unless additional interactions with other residues exist. Among 8327 long-range hydrogen bonds found in this set of structures, 2811, or about one-third, are hydrogen bonds entailing RNA backbone atoms; they involve 39% of all nucleotides in the structures. The majority of them (2111) are hydrogen bonds entailing ribose hydroxyl groups, which can be used either as a donor or an acceptor; they constitute 25% of all hydrogen bonds and involve 31% of all nucleotides. The phosphate oxygens OP1 or OP2 are used as hydrogen bond acceptors in 12% of all nucleotides, and the ribose ring oxygen O4' and phosphodiester oxygens O3' and O5' are used in 4%, 4%, and 1% of all nucleotides, respectively. Distributions of geometric parameters and some examples of such hydrogen bonds are presented in this report. A novel motif involving backbone hydrogen bonds, the ribose-phosphate zipper, is also identified.

Structure of the N-terminal Regulatory Domain of a Plant NADPH Oxidase and Its Functional Implications

Plant NADPH oxidases (Rboh, for respiratory burst oxidase homolog) produce reactive oxygen species that are key regulators of various cellular events including plant innate immunity. Rbohs possess a highly conserved cytoplasmic N-terminal region containing two EF-hand motifs that regulate Rboh activity. Rice (Oryza sativa) RbohB (OsRbohB) is regulated by the direct binding of a small GTPase (Rac1) to this regulatory region as well as by Ca(2+) binding to the EF-hands. Here, we present the atomic structure of the N-terminal region of OsRbohB. The structure reveals that OsRbohB forms a unique dimer stabilized by swapping the EF-hand motifs. We identified two additional EF-hand-like motifs that were not predicted from sequence data so far. These EF-hand-like motifs together with the swapped EF-hands form a structure similar to that found in calcineurin B. We observed conformational changes mediated by Ca(2+) binding to only one EF-hand. Structure-based in vitro pulldown assays and NMR titration experiments defined the OsRac1 binding interface within the coiled-coil region created by swapping the EF-hands. In addition, we demonstrate a direct intramolecular interaction between the N and C terminus, and that the complete N-terminal cytoplasmic region is required for this interaction. The structural features and intramolecular interactions characterized here might be common elements shared by Rbohs that contribute to the regulation of reactive oxygen species production.

Allosteric Mechanism in Neurophysin

Neurophysins are small two- domain proteins that exhibit ligand-facilitated dimerization and its reciprocal-stronger binding by dimer than by monomer. We have previously reported evidence of conformational differences between unliganded monomer and dimer, but had not established their identity, their origin, or how they affect binding affinity.In the present study we have determined the relative structures of a neurophysin monomer and its dimer by NMR spectroscopy. The structures show that a number of neurophysin residues that interact with the essential tyrosine of the ligand peptide move from locations in the unliganded monomer that are unfavorable for binding to positions in the unliganded dimer that more closely resemble those of the liganded state, thereby contributing to the stronger binding by dimer. We have also identified the origin of the dimerization-induced conformational change; i.e., the same residues that form the principal inter-domain bonds in the monomer via side-chain interactions form the principal inter-subunit bonds in the dimer via backbone hydrogen bonding. Steric crowding upon dimerization leads to a modification of inter-domain interactions, a resultant change in the relative orientation of the two domains and associated changes in intra-domain interactions. The principal residues responsible for these interactions are the Cys residues of the 28-34 disulfide bridge (amino domain) and 74-79 disulfide bridge (carboxyl domain). Strong inter-domain van der Waals interactions between sulfur atoms appear to play important roles in determining the relative domain orientation.(Supported by grant GM-17528 from NIH).

FoldAmyloid: a method of prediction of amyloidogenic regions from protein sequence

Amyloidogenic regions in polypeptide chains are very important because such regions are responsible for amyloid formation and aggregation. It is useful to be able to predict positions of amyloidogenic regions in protein chains. Two characteristics (expected probability of hydrogen bonds formation and expected packing density of residues) have been introduced by us to detect amyloidogenic regions in a protein sequence. We demonstrate that regions with high expected probability of the formation of backbone-backbone hydrogen bonds as well as regions with high expected packing density are mostly responsible for the formation of amyloid fibrils. Our method (FoldAmyloid) has been tested on a dataset of 407 peptides (144 amyloidogenic and 263 non-amyloidogenic peptides) and has shown good performance in predicting a peptide status: amyloidogenic or non-amyloidogenic. The prediction based on the expected packing density classified correctly 75% of amyloidogenic peptides and 74% of non-amyloidogenic ones. Two variants (averaging by donors and by acceptors) of prediction based on the probability of formation of backbone-backbone hydrogen bonds gave a comparable efficiency. With a hybrid-scale constructed by merging the above three scales, our method is correct for 80% of amyloidogenic peptides and for 72% of non-amyloidogenic ones. Prediction of amyloidogenic regions in proteins where positions of amyloidogenic regions are known from experimental data has also been done. In the proteins, our method correctly finds 10 out of 11 amyloidogenic regions. The FoldAmyloid server is available at http://antares.protres.ru/fold-amyloid/.

Simulation Studies on the Stabilities of Aggregates Formed by Fibril-Forming Segments of α-Synuclein

We performed molecular dynamics simulations for various oligomers with different β-sheet conformations consisting of α-Synuclein 71–82 residues using an all atom force field and explicit water model. Tetramers of antiparallel β-sheet are shown to be stable, whereas parallel sheets are highly unstable due to the repulsive interactions between bulky and polar side chains as well as the weaker backbone hydrogen bonds. We also investigated the stabilities of double antiparallel β-sheets stacked with asymmetric and symmetric geometries. Our results show that this 12 amino acid residue peptide can form stable β-sheet conformers at 320K and higher temperatures. The backbone hydrogen bonds in β-sheet and the steric packing between hydrophobic side chains between β-sheets are shown to give conformational stabilities.

Structure of Escherichia coli Hfq bound to polyriboadenylate RNA

Hfq is a small, highly abundant hexameric protein that is found in many bacteria and plays a critical role in mRNA expression and RNA stability. As an "RNA chaperone," Hfq binds AU-rich sequences and facilitates the trans annealing of small RNAs (sRNAs) to their target mRNAs, typically resulting in the down-regulation of gene expression. Hfq also plays a key role in bacterial RNA decay by binding tightly to polyadenylate [poly(A)] tracts. The structural mechanism by which Hfq recognizes and binds poly(A) is unknown. Here, we report the crystal structure of Escherichia coli Hfq bound to the poly(A) RNA, A(15). The structure reveals a unique RNA binding mechanism. Unlike uridine-containing sequences, which bind to the "proximal" face, the poly(A) tract binds to the "distal" face of Hfq using 6 tripartite binding motifs. Each motif consists of an adenosine specificity site (A site), which is effected by peptide backbone hydrogen bonds, a purine nucleotide selectivity site (R site), and a sequence-nondiscriminating RNA entrance/exit site (E site). The resulting implication that Hfq can bind poly(A-R-N) triplets, where R is a purine nucleotide and N is any nucleotide, was confirmed by binding studies. Indeed, Hfq bound to the oligoribonucleotides (AGG)(8), (AGC)(8), and the shorter (A-R-N)(4) sequence, AACAACAAGAAG, with nanomolar affinities. The abundance of (A-R-N)(4) and (A-R-N)(5) triplet repeats in the E. coli genome suggests additional RNA targets for Hfq. Further, the structure provides insight into Hfq-mediated sRNA-mRNA annealing and the role of Hfq in RNA decay.

Peptide T exhibits a well‐defined structure in fluorinated solvent at low temperature

The structure of Peptide T was determined by solution NMR spectroscopy, under strong structure-inducing conditions: 40% hexafluoro-2-propanol aqueous solution at 5 degrees C. Under these conditions it was possible to detect medium-range NOEs for the first time for this peptide. This allowed a much better-defined structure to be determined for Peptide T in comparison with earlier NMR and computational studies. Peptide structures consistent with the experimental restraints were generated using a restrained MD simulation with a full empirical force field. Residues 4-8 of Peptide T take on a well-defined structure with a heavy atom RMSD of 0.78 A. The structure is stabilized by hydrogen bonding to side-chain oxygen atoms of Thr 4 and Thr 8, as well as backbone hydrogen bonding between residues 5 and 7 that forms this region into a classic gamma-turn.

Fibronectin forms the most extensible biological fibers displaying switchable force-exposed cryptic binding sites

Rather than maximizing toughness, as needed for silk and muscle titin fibers to withstand external impact, the much softer extracellular matrix fibers made from fibronectin (Fn) can be stretched by cell generated forces and display extraordinary extensibility. We show that Fn fibers can be extended more than 8-fold (>700% strain) before 50% of the fibers break. The Young's modulus of single fibers, given by the highly nonlinear slope of the stress-strain curve, changes orders of magnitude, up to MPa. Although many other materials plastically deform before they rupture, evidence is provided that the reversible breakage of force-bearing backbone hydrogen bonds enables the large strain. When tension is released, the nano-sized Fn domains first contract in the crowded environment of fibers within seconds into random coil conformations (molten globule states), before the force-bearing hydrogen bond networks that stabilize the domain's secondary structures are reestablished within minutes (double exponential). The exposure of cryptic binding sites on Fn type III modules increases steeply upon stretching. Thus fiber extension steadily up-regulates fiber rigidity and cryptic epitope exposure, both of which are known to differentially alter cell behavior. Finally, since stress-strain relationships cannot directly be measured in native extracellular matrix (ECM), the stress-strain curves were correlated with stretch-induced alterations of intramolecular fluorescence resonance energy transfer (FRET) obtained from trace amounts of Fn probes (mechanical strain sensors) that can be incorporated into native ECM. Physiological implications of the extraordinary extensibility of Fn fibers and contraction kinetics are discussed.

Understanding hydrogen-bond patterns in proteins using network motifs

Protein structures can be viewed as networks of contacts (edges) between amino-acid residues (nodes). Here we dissect proteins into sub-graphs consisting of six nodes and their corresponding edges, with an edge being either a backbone hydrogen bond (H-bond) or a covalent interaction. Six thousand three hundred and twenty-two such sub-graphs were found in a large non-redundant dataset of high-resolution structures, from which 35 occur much more frequently than in a random model. Many of these significant sub-graphs (also called network motifs) correspond to sub-structures of alpha helices and beta-sheets, as expected. However, others correspond to more exotic sub-structures such as 3(10) helix, Schellman motif and motifs that were not defined previously. This topological characterization of patterns is very useful for producing a detailed differences map to compare protein structures. Here we analyzed in details the differences between NMR, molecular dynamics (MD) simulations and X-ray structures for Lysozyme, SH3 and the lambda repressor. In these cases, the same structures solved by NMR and simulated by MD showed small but consistent differences in their motif composition from the crystal structures, despite a very small root mean square deviation (RMSD) between them. This may be due to differences in the pair-wise energy functions used and the dynamic nature of these proteins. A web-based tool to calculate network motifs is available at http://bioinfo.weizmann.ac.il/protmot/.

Cis−Trans Proline Isomerization Effects on Collagen Triple-Helix Stability Are Limited

We investigated the effect of restricting cis-trans proline isomerization on collagen triple-helix stability. The Pro residues at the Xaa and Yaa positions of an (Xaa-Yaa-Gly) triplet were replaced by a Pro-trans-Pro alkene isostere in the host-guest peptide, H-(Pro-Pro-Gly)(10)-OH. The resulting alkene isostere peptide had a T(m) value 53.6 degrees C lower than that of the control peptide. The Pro-trans-Pro alkene isostere peptide had a T(m) value 3.9 degrees C higher than that of the previously reported Pro-trans-Gly alkene isostere peptide that did not involve cis-trans Pro isomerization (Jenkins, C. L.; Vasbinder, M. M.; Miller, S. J.; Raines, R. T. Org. Lett. 2005, 7, 2619-22). Thus, single cis-trans proline amide isomerization alone has limited contribution to the overall stability of the collagen triple helix. Since collagen has a high content of imino acid residues, the cumulative effects of cis-trans isomerization may be quite significant. The peptide containing the Pro-trans-Pro isostere was significantly less stable than the previously reported Gly-trans-Pro alkene isostere peptide that retained the backbone interchain hydrogen bond (Dai, N.; Wang, X. J.; Etzkorn, F. A. J. Am. Chem. Soc. 2008, 130, 5396-5397), which confirms that direct interchain backbone hydrogen bonding is a major force for stabilizing the collagen triple helix.

Induced β-Barrel Formation of the Alzheimer's A β25–35 Oligomers on Carbon Nanotube Surfaces: Implication for Amyloid Fibril Inhibition

Recent experimental studies show that carbon nanotubes impact the aggregation process of proteins associated with neurodegenerative diseases. However, the details of molecular interactions between proteins and carbon nanotubes are still not well understood. In this study, we investigate the initial adsorption features and dynamics of the Alzheimer's amyloid- β peptide spanning residues 25–35 (A β25–35) on a single-walled carbon nanotube (SWNT) surface using fully atomic molecular dynamics simulations (MD) in explicit solvent. The initial configurations of the A β25–35 peptides consist of two preformed bilayer β-sheets, each with four or five β-strands in parallel or mixed antiparallel-parallel orientations. Our simulations show, for what we believe is the first time, that two disjointed A β25–35 β-sheets with mixed antiparallel-parallel strands can assemble into β-barrels wrapping the SWNT. In contrast, both simulations of A β25–35 without SWNT, and simulations of SWNT−A β25–35 with purely parallel β-strands, lead to disordered aggregates. We find that A β25–35 β-barrel formation involves at least two steps: i), curving of the A β25–35 β-sheets as a result of strong hydrophobic interactions with carbon nanotube concomitantly with dehydration of the SWNT-peptide interface; and ii), intersheet backbone hydrogen bond formation with fluctuating intrasheet hydrogen bonds. Detailed analysis of the conversion shows that β-barrel formation on SWNT surface results from the interplay of dehydration and peptide-SWNT/peptide-peptide interactions. Implications of our results on amyloid fibril inhibition are discussed.

Dimerization of Fibril-forming Segments of α-Synuclein

We have performed replica-exchange molecular dynamics (REMD) simulations on the dimer formation of fibril- forming segments of α-Synuclein (residues 71 - 82) using implicit solvation models with two kinds of force fields- AMBER parm99SB and parm96. We observed spontaneous formation of dimers from the extensive simulations, demonstrating the self-aggregating and fibril forming properties of the peptides. Secondary structure profile and clustering analysis showed that dimers with antiparallel β-sheet conformations, stabilized by well-defined hydrogen boding, are major species corresponding to global free energy minimum. Parallel dimers with partial β-sheets are found to be off-pathway intermediates. The relative instability of the parallel arrangements is due to the repulsive interactions between bulky and polar side chains as well as weaker backbone hydrogen bonds.

On the remarkable mechanostability of scaffoldins and the mechanical clamp motif

Protein mechanostability is a fundamental biological property that can only be measured by single-molecule manipulation techniques. Such studies have unveiled a variety of highly mechanostable modules (mainly of the Ig-like, beta-sandwich type) in modular proteins subjected to mechanical stress from the cytoskeleton and the metazoan cell-cell interface. Their mechanostability is often attributed to a "mechanical clamp" of secondary structure (a patch of backbone hydrogen bonds) fastening their ends. Here we investigate the nanomechanics of scaffoldins, an important family of scaffolding proteins that assembles a variety of cellulases into the so-called cellulosome, a microbial extracellular nanomachine for cellulose adhesion and degradation. These proteins anchor the microbial cell to cellulose substrates, which makes their connecting region likely to be subjected to mechanical stress. By using single-molecule force spectroscopy based on atomic force microscopy, polyprotein engineering, and computer simulations, here we show that the cohesin I modules from the connecting region of cellulosome scaffoldins are the most robust mechanical proteins studied experimentally or predicted from the entire Protein Data Bank. The mechanostability of the cohesin modules studied correlates well with their mechanical kinetic stability but not with their thermal stability, and it is well predicted by computer simulations, even coarse-grained. This extraordinary mechanical stability is attributed to 2 mechanical clamps in tandem. Our findings provide the current upper limit of protein mechanostability and establish shear mechanical clamps as a general structural/functional motif widespread in proteins putatively subjected to mechanical stress. These data have important implications for the scaffoldin physiology and for protein design in biotechnology and nanotechnology.

Tyrosine kinase inhibition: Ligand binding and conformational change in c-Kit and c-Abl

The conformational flexibility exhibited by protein kinases poses an enormous challenge to the design of cancer therapeutics. Additionally the high degree of structural conservation within the kinase superfamily often leads to inhibitors that exhibit little selectivity and substantial cross reactivity. This work investigates the conformational changes that accompany the binding of Gleevec, or imatinib mesylate, to the tyrosine kinases c-Kit and c-Abl. Our analysis is that this fit is driven, at least in part, by the need to exclude water from solvent-exposed backbone hydrogen bonds. Both experimental and molecular modeling studies of the active state inhibitor of the tyrosine kinase c-Abl indicate that solvent exclusion also plays a role in this system.

Structural Basis of the Catalytic Mechanism Operating in Open-Closed Conformers of Lipocalin Type Prostaglandin D Synthase

Lipocalin type prostaglandin D synthase (L-PGDS) is a multifunctional protein acting as a somnogen (PGD2)-producing enzyme, an extracellular transporter of various lipophilic ligands, and an amyloid-β chaperone in human cerebrospinal fluid. In this study, we determined the crystal structures of two different conformers of mouse L-PGDS, one with an open cavity of the β-barrel and the other with a closed cavity due to the movement of the flexible E-F loop. The upper compartment of the central large cavity contains the catalytically essential Cys65 residue and its network of hydrogen bonds with the polar residues Ser45, Thr67, and Ser81, whereas the lower compartment is composed of hydrophobic amino acid residues that are highly conserved among other lipocalins. SH titration analysis combined with site-directed mutagenesis revealed that the Cys65 residue is activated by its interaction with Ser45 and Thr67 and that the S45A/T67A/S81A mutant showed less than 10% of the L-PGDS activity. The conformational change between the open and closed states of the cavity indicates that the mobile calyx contributes to the multiligand binding ability of L-PGDS.

Enhancement of beta-sheet assembly by cooperative hydrogen bonds potential

The roughness of energy landscapes is a major obstacle to protein structure prediction, since it forces conformational searches to spend much time struggling to escape numerous traps. Specifically, beta-sheet formation is prone to stray, since many possible combinations of hydrogen bonds are dead ends in terms of beta-sheet assembly. It has been shown that cooperative terms for backbone hydrogen bonds ease this problem by augmenting hydrogen bond patterns that are consistent with beta sheets. Here, we present a novel cooperative hydrogen-bond term that is both effective in promoting beta sheets and computationally efficient. In addition, the new term is differentiable and operates on all-atom protein models. Energy optimization of poly-alanine chains under the new term led to significantly more beta-sheet content than optimization under a non-cooperative term. Furthermore, the optimized structure included very few non-native patterns. The new term is implemented within the MESHI package and is freely available at http://cs.bgu.ac.il/ approximately meshi.

New angle‐dependent potential energy function for backbone–backbone hydrogen bond in protein–protein interactions

Backbone-backbone hydrogen bonds (BBHBs) are one of the most abundant interactions at the interface of protein-protein complex. Here, we propose an angle-dependent potential energy function for BBHB based on density functional theory (DFT) calculations and the operation of a genetic algorithm to find the optimal parameters in the potential energy function. The angular part of the energy function is assumed to be the product of the power series of sine and cosine functions with respect to the two angles associated with BBHB. Two radial functions are taken into account in this study: Morse and Leonard-Jones 12-10 potential functions. Of these two functions under consideration, the former is found to be more accurate than the latter in terms of predicting the binding energies obtained from DFT calculations. The new HB potential function also compares well with the knowledge-based potential derived by applying Boltzmann statistics for a variety of protein-protein complexes in protein data bank.

A Single α‐Helical Turn Stabilized by Replacement of an Internal Hydrogen Bond with a Covalent Ethylene Bridge

Walking the carbon bridge to helix stabilization: α-Helix mimetics have been synthesized that contain an ethylene linkage in place of an (i,i+4) hydrogen bond (see structure). The smallest unit possible, a modified pentapeptide, shows the characteristic features of a helix (see CD spectrum). Incorporation of this carbon bridge into a peptide represents the first example of replacement of an internal backbone hydrogen bond.