The cis-Pro touch-turn: a rare motif preferred at functional sites

A Joint Experimental and Theoretical Study on the Structural and Spectroscopic Properties of the Piv-Pro-d-Ser-NHMe Peptide

In this paper, we investigate the secondary structure of the Piv-Pro-d-Ser-NHMe peptide by means of nuclear magnetic resonance (NMR) and electronic circular dichroism (ECD) experiments, in conjunction with theoretical simulations based on molecular dynamics and time-dependent density functional theory calculations including polarizable embedding to account for solvent effects. The various experimental and theoretical protocols are assessed and validated, and are shown to provide a consistent description of the turn structure adopted by this peptide in solution. In addition, a simple fitting procedure is proposed to make the simulated and experimental ECD almost perfectly match. This full methodology is finally tested on another small peptide, enlightening its efficiency and robustness.

Characterization of β-turns by electronic circular dichroism spectroscopy: a coupled molecular dynamics and time-dependent density functional theory computational study

Electronic circular dichroism is one of the most used spectroscopic techniques for peptide and protein structural characterization. However, while valuable experimental spectra exist for α-helix, β-sheet and random coil secondary structures, previous studies showed important discrepancies for β-turns, limiting their use as a reference for structural studies. In this paper, we simulated circular dichroism spectra for the best-characterized β-turns in peptides, namely types I, II, I' and II'. In particular, by combining classical molecular dynamics simulations and state-of-the-art quantum time-dependent density functional theory (with the polarizable embedding multiscale model) computations, two common electronic circular dichroism patterns were found for couples of β-turn types (namely, type I/type II' and type II/type I'), at first for a minimal di-peptide model (Ace-Ala-Ala-NHMe), but also for all sequences tested with non-aromatic residues in the central positions. On the other hand, as expected, aromatic substitution causes important perturbations to the previously found ECD patterns. Finally, by applying suitable approximations, these patterns were subsequently rationalized based on the exciton chirality rule. All these results provide useful predictions and pave the way for a possible experimental characterization of β-turns based on circular dichroism spectroscopy.

Common Mode of Remodeling AAA ATPases p97/CDC48 by Their Disassembling Cofactors ASPL/PUX1

The hexameric ring structure of the type II AAA+ ATPases is considered as stable and permanent. Recently, the UBX domain-containing cofactors Arabidopsis thaliana PUX1 and human alveolar soft part sarcoma locus (ASPL) were reported to bind and disassemble the cognate AAA+ ATPases AtCDC48 and human p97. Here, we present two crystal structures related to these complexes: a truncated AtCDC48 (AtCDC48-ND1) and a hybrid complex containing human p97-ND1 and the UBX domain of plant PUX1 (p97-ND1:PUX1-UBX). These structures reveal close similarity between the human and plant AAA+ ATPases, but also highlight differences between disassembling and non-disassembling AAA+ ATPase cofactors. Based on an AtCDC48 disassembly assay with PUX1 and known crystal structures of the p97-bound human cofactor ASPL, we propose a general ATPase disassembly model. Thus, our structural and biophysical investigations provide detailed insight into the mechanism of AAA+ ATPase disassembly by UBX domain cofactors and suggest a general mode of regulating the cellular activity of these molecular machines.

Unconstrained peptoid tetramer exhibits a predominant conformation in aqueous solution

Conformational control in peptoids, N-substituted glycines, is crucial for the design and synthesis of biologically-active compounds and atomically-defined nanomaterials. While there are a growing number of structural studies in solution, most have been performed with conformationally-constrained short sequences (e.g., sterically-hindered sidechains or macrocyclization). Thus, the inherent degree of heterogeneity of unconstrained peptoids in solution remains largely unstudied. Here, we explored the folding landscape of a series of simple peptoid tetramers in aqueous solution by NMR spectroscopy. By incorporating specific ¹³ C-probes into the backbone using bromoacetic acid-2-¹³ C as a submonomer, we developed a new technique for sequential backbone assignment of peptoids based on the 1,n-Adequate pulse sequence. Unexpectedly, two of the tetramers, containing an N-(2-aminoethyl)glycine residue (Nae), had preferred conformations. NMR and molecular dynamics studies on one of the tetramers showed that the preferred conformer (52%) had a trans-cis-trans configuration about the three amide bonds. Moreover, >80% of the ensemble contained a cis amide bond at the central amide. The backbone dihedral angles observed fall directly within the expected minima in the peptoid Ramachandran plot. Analysis of this compound against similar peptoid analogs suggests that the commonly used Nae monomer plays a key role in the stabilization of peptoid structure via a side-chain-to-main-chain interaction. This discovery may offer a simple, synthetically high-yielding approach to control peptoid structure, and suggests that peptoids have strong intrinsic conformational preferences in solution. These findings should facilitate the predictive design of folded peptoid structures, and accelerate application in areas ranging from drug discovery to biomimetic nanoscience.

Broad Analysis of Vicinal Disulfides: Occurrences, Conformations with Cis or with Trans Peptides, and Functional Roles Including Sugar Binding

Vicinal disulfides between sequence-adjacent cysteine residues are very rare and rather startling structural features which play a variety of functional roles. Typically discussed as an isolated curiosity, they have never received a general treatment covering both cis and trans forms. Enabled by the growing database of high-resolution structures, required deposition of diffraction data, and improved methods for discriminating reliable from dubious cases, we identify and describe distinct protein families with reliably genuine examples of cis or trans vicinal disulfides and discuss their conformations, conservation, and functions. No cis-trans interconversions and only one case of catalytic redox function are seen. Some vicinal disulfides are essential to large, functionally coupled motions, whereas most form the centers of tightly packed internal regions. Their most widespread biological role is providing a rigid hydrophobic contact surface under the undecorated side of a sugar or multiring ligand, contributing an important aspect of binding specificity.

Conformational preferences of peptide-peptoid hybrid oligomers

Peptomers are oligomeric molecules composed of both α-amino acids and N-substituted glycine monomers, thus creating a hybrid of peptide and peptoid units. Peptomers have been used in several applications such as antimicrobials, protease inhibitors, and antibody mimics. Despite the considerable promise of peptomers as chemically diverse molecular scaffolds, we know little about their conformational tendencies. This lack of knowledge limits the ability to implement computational approaches for peptomer design. Here we computationally evaluate the local structural propensities of the peptide-peptoid linkage. We find some general similarities between the peptide residue conformational preferences and the Ramachandran distribution of residues that precede proline in folded protein structures. However, there are notable differences. For example, several β-turn motifs are disallowed when the i+2 residue is also a peptoid monomer. Significantly, the lowest energy geometry, when dispersion forces are accounted for, corresponds to a "cis-Pro touch-turn" conformation, an unusual turn motif that has been observed at protein catalytic centers and binding sites. The peptomer touch-turn thus represents a useful design element for the construction of folded oligomers capable of molecular recognition and as modules in the assembly of structurally complex peptoid-protein hybrid macromolecules.

Identifying functionally important cis-peptide containing segments in proteins and their utility in molecular function annotation

Cis-peptide embedded segments are rare in proteins but often highlight their important role in molecular function when they do occur. The high evolutionary conservation of these segments illustrates this observation almost universally, although no attempt has been made to systematically use this information for the purpose of function annotation. In the present study, we demonstrate how geometric clustering and level-specific Gene Ontology molecular-function terms (also known as annotations) can be used in a statistically significant manner to identify cis-embedded segments in a protein linked to its molecular function. The present study identifies novel cis-peptide fragments, which are subsequently used for fragment-based function annotation. Annotation recall benchmarks interpreted using the receiver-operator characteristic plot returned an area-under-curve > 0.9, corroborating the utility of the annotation method. In addition, we identified cis-peptide fragments occurring in conjunction with functionally important trans-peptide fragments, providing additional insights into molecular function. We further illustrate the applicability of our method in function annotation where homology-based annotation transfer is not possible. The findings of the present study add to the repertoire of function annotation approaches and also facilitate engineering, design and allied studies around the cis-peptide neighborhood of proteins.

Atomic-accuracy prediction of protein loop structures through an RNA-inspired Ansatz

Consistently predicting biopolymer structure at atomic resolution from sequence alone remains a difficult problem, even for small sub-segments of large proteins. Such loop prediction challenges, which arise frequently in comparative modeling and protein design, can become intractable as loop lengths exceed 10 residues and if surrounding side-chain conformations are erased. Current approaches, such as the protein local optimization protocol or kinematic inversion closure (KIC) Monte Carlo, involve stages that coarse-grain proteins, simplifying modeling but precluding a systematic search of all-atom configurations. This article introduces an alternative modeling strategy based on a ‘stepwise ansatz’, recently developed for RNA modeling, which posits that any realistic all-atom molecular conformation can be built up by residue-by-residue stepwise enumeration. When harnessed to a dynamic-programming-like recursion in the Rosetta framework, the resulting stepwise assembly (SWA) protocol enables enumerative sampling of a 12 residue loop at a significant but achievable cost of thousands of CPU-hours. In a previously established benchmark, SWA recovers crystallographic conformations with sub-Angstrom accuracy for 19 of 20 loops, compared to 14 of 20 by KIC modeling with a comparable expenditure of computational power. Furthermore, SWA gives high accuracy results on an additional set of 15 loops highlighted in the biological literature for their irregularity or unusual length. Successes include cis-Pro touch turns, loops that pass through tunnels of other side-chains, and loops of lengths up to 24 residues. Remaining problem cases are traced to inaccuracies in the Rosetta all-atom energy function. In five additional blind tests, SWA achieves sub-Angstrom accuracy models, including the first such success in a protein/RNA binding interface, the YbxF/kink-turn interaction in the fourth ‘RNA-puzzle’ competition. These results establish all-atom enumeration as an unusually systematic approach to ab initio protein structure modeling that can leverage high performance computing and physically realistic energy functions to more consistently achieve atomic accuracy.

Doing molecular biophysics: finding, naming, and picturing signal within complexity

A macromolecular structure, as measured data or as a list of coordinates or even on-screen as a full atomic model, is an extremely complex and confusing object. The underlying rules of how it folds, moves, and interacts as a biological entity are even less evident or intuitive to the human mind. To do science on such molecules, or to relate them usefully to higher levels of biology, we need to start with a natural history that names their features in meaningful ways and with multiple representations (visual or algebraic) that show some aspect of their organizing principles. The two of us have jointly enjoyed a highly varied and engrossing career in biophysical research over nearly 50 years. Our frequent changes of emphasis are tied together by two threads: first, by finding the right names, visualizations, and methods to help both ourselves and others to better understand the 3D structures of protein and RNA molecules, and second, by redefining the boundary between signal and noise for complex data, in both directions-sometimes identifying and promoting real signal up out of what seemed just noise, and sometimes demoting apparent signal into noise or systematic error. Here we relate parts of our scientific and personal lives, including ups and downs, influences, anecdotes, and guiding principles such as the title theme.

Investigating diproline segments in proteins: occurrences, conformation and classification

The covalent linkage between the side-chain and the backbone nitrogen atom of proline leads to the formation of the five-membered pyrrolidine ring and hence restriction of the backbone torsional angle ϕ to values of -60 °± 30° for the L-proline. Diproline segments constitute a chain fragment with considerably reduced conformational choices. In the current study, the conformational states for the diproline segment (( L) Pro-( L) Pro) found in proteins has been investigated with an emphasis on the cis and trans states for the Pro-Pro peptide bond. The occurrence of diproline segments in turns and other secondary structures has been studied and compared to that of Xaa-Pro-Yaa segments in proteins which gives us a better understanding on the restriction imposed on other residues by the diproline segment and the single proline residue. The study indicates that P(II) -P(II) and P(II) -α are the most favorable conformational states for the diproline segment. The analysis on Xaa-Pro-Yaa sequences reveals that the Xaa-Pro peptide bond exists preferably as the trans conformer rather than the cis conformer. The present study may lead to a better understanding of the behavior of proline occurring in diproline segments which can facilitate various designed diproline-based synthetic templates for biological and structural studies.

Use of secondary structure element information in drug design: polypharmacology and conserved motifs in protein–ligand binding and protein–protein interfaces

The structure-based design of small-molecule inhibitors of protein–ligand and protein–protein interfaces is a key component of drug discovery. The underlying protein interactions can be regarded based on structural similarity of the secondary structure elements: similarities around the binding site (‘ligand-sensing cores’) or in the protein interface (‘interface-sensing surfaces’) in otherwise unrelated proteins can be useful in predicting polypharmacology and identifying new lead structures. Even small conserved motifs can provide similar interaction patterns in proteins with a completely different fold and function. The identification of these structural similarities can help in the design of new drugs by guiding further optimization. Here, the concepts and ideas based on secondary structure element similarities and their successful applications in drug design are reviewed and discussed.

Crystal structure of human FAF1 UBX domain reveals a novel FcisP touch-turn motif in p97/VCP-binding region

UBX domain is a general p97/VCP-binding module found in an increasing number of proteins including FAF1, p47, SAKS1 and UBXD7. FAF1, a multi-functional tumor suppressor protein, binds to the N domain of p97/VCP through its C-terminal UBX domain and thereby inhibits the proteasomal protein degradation in which p97/VCP acts as a co-chaperone. Here we report the crystal structure of human FAF1 UBX domain at 2.9Å resolution. It reveals that the conserved FP sequence in the p97/VCP-binding region adopts a rarely observed cis-Pro touch-turn structure. We call it an FcisP touch-turn motif and suggest that it is the conserved structural element of the UBX domain. Four FAF1 UBX molecules in an asymmetric unit of the crystal show two different conformations of the FcisP touch-turn motif. The phenyl ring of F(619) in the motif stacks partly over cis-Pro(620) in one conformation, whereas it is swung out from cis-P(620), in the other conformation, and forms hydrophobic contacts with the residues of the neighboring molecule. In addition, the entire FcisP touch-turn motif is pulled out in the second conformation by about 2Å in comparison to the first conformation. Those conformational differences observed in the p97/VCP-binding motif caused by the interaction with neighboring molecules presumably represent the conformational change of the UBX domain on its binding to the N domain of p97/VCP.

Phage display and peptide mapping of an immunoglobulin light chain fibril-related conformational epitope

Amyloid fibrils and partially unfolded intermediates can be distinguished serologically from native amyloidogenic precursor proteins or peptides. In this regard, we previously had reported that mAb 11-1F4, generated by immunizing mice with a thermally denatured variable domain (VL) fragment of the human kappa4 Bence Jones protein Len, bound to a non-native conformational epitope located within the N-terminal 18 residues of fibrillar, as well as partially denatured, Ig light chains (O'Nuallain, B., et al. (2006) Biochemistry 46, 1240-1247). To define further the antibody binding site, we used random peptide phage display and epitope mapping of VL Len using wild-type and alanine-mutated Len peptides where it was shown that the antibody epitope was reliant on up to 10 of the first 15 residues of protein Len. Comparison of Vkappa and Vlambda N-terminal germline consensus sequences with protein Len and 11-1F4-binding phages indicated that this antibody's cross-reactivity with light chains was related to an invariant proline at position(s) 7 and/or 8, bulky hydrophobic residues at positions 11 and 13, and additionally, to the ability to accommodate amino acid diversity at positions 1-4. Sequence alignments of the phage peptides revealed a central proline, often flanked by aromatic residues. Taken together, these results have provided evidence for the structural basis of the specificity of 11-1F4 for both kappa and lambda light chain fibrils. We posit that the associated binding site involves a rare type VI beta-turn or touch-turn that is anchored by a cis-proline residue. The identification of an 11-1F4-related mimotope should facilitate development of pan-light chain fibril-reactive antibodies that could be used in the diagnosis and treatment of patients with AL amyloidosis.

MolProbity: all-atom contacts and structure validation for proteins and nucleic acids

MolProbity is a general-purpose web server offering quality validation for 3D structures of proteins, nucleic acids and complexes. It provides detailed all-atom contact analysis of any steric problems within the molecules as well as updated dihedral-angle diagnostics, and it can calculate and display the H-bond and van der Waals contacts in the interfaces between components. An integral step in the process is the addition and full optimization of all hydrogen atoms, both polar and nonpolar. New analysis functions have been added for RNA, for interfaces, and for NMR ensembles. Additionally, both the web site and major component programs have been rewritten to improve speed, convenience, clarity and integration with other resources. MolProbity results are reported in multiple forms: as overall numeric scores, as lists or charts of local problems, as downloadable PDB and graphics files, and most notably as informative, manipulable 3D kinemage graphics shown online in the KiNG viewer. This service is available free to all users at http://molprobity.biochem.duke.edu.

Localization of a conformational epitope common to non-native and fibrillar immunoglobulin light chains

Amyloid fibrils and partially unfolded intermediates may be distinguished serologically from native amyloidogenic precursor proteins or peptides. In this regard, we had previously reported that the IgG1 mAb 11-1F4, generated by immunizing mice with a thermally denatured variable region fragment of the human Igkappa4 Bence Jones protein Len, reacted specifically with light chain (LC) fibrils, irrespective of kappa or lambda isotype but, notably, did not with native molecules (Hrncic, R. et al. (2000) Am. J. Pathol. 157, 1239-1246). To elucidate the molecular basis of this specificity, we have used a europium-linked fluorescent immunoassay, where it was demonstrated through epitope mapping that mAb 11-1F4 recognizes a conformational determinant contained within the first (N-terminal) 18 amino acids of misfolded LCs. The nature of this epitope was evidenced in competition studies where the peptide Len (1-18), but not the intact protein or other LCs, inhibited the binding of the antibody to fibrils. This unique reactivity was dependent on the structural integrity of this portion of the molecule, particularly the presence of a highly conserved prolyl residue at position 8. On the basis of our experimental data, we posit that the mAb 11-1F4 binding site found on partially denatured and fibrillar LCs involves an inducible N-terminal main chain reversal that results in the formation of a proline anchored beta-turn. Our delineation of this LC fibril-associated epitope provides the rationale for the design of novel amyloid-reactive antibodies with diagnostic and therapeutic potential for patients with LC-associated and other forms of amyloidosis.

A test of enhancing model accuracy in high-throughput crystallography

The high throughput of structure determination pipelines relies on increased automation and, consequently, a reduction of time spent on interactive quality control. In order to meet and exceed current standards in model accuracy, new approaches are needed for the facile identification and correction of model errors during refinement. One such approach is provided by the validation and structure-improvement tools of the MOL: PROBITY: web service. To test their effectiveness in high-throughput mode, a large subset of the crystal structures from the SouthEast Collaboratory for Structural Genomics (SECSG) has used protocols based on the MOL: PROBITY: tools. Comparison of 29 working-set and 19 control-set SECSG structures shows that working-set outlier scores for updated Ramachandran-plot, sidechain rotamer, and all-atom steric criteria have been improved by factors of 5- to 10-fold (relative to the control set or to a Protein Data Bank sample), while quality of covalent geometry, R(work), R(free), electron density and difference density are maintained or improved. Some parts of this correction process are already fully automated; other parts involve manual rebuilding of conformations flagged by the tests as trapped in the wrong local minimum, often altering features of functional significance. The ease and effectiveness of this technique shows that macromolecular crystal structures from either traditional or high-throughput determinations can feasibly reach a new level of excellence in conformational accuracy and reliability.

The advantages and limitations of protein crystal structures

Crystal structure analysis using X-ray diffraction is, in many cases, the most advanced method available for obtaining high-resolution structural information about biological macromolecules. The ways in which X-ray diffraction data are collected and refined have a strong impact on the final quality of the structural models and the type and magnitude of their associated errors. It is becoming increasingly necessary for both structural and non-structural biologists to judge the reliability and accuracy of these models, which are being used in many aspects of research, including structure-based drug design, and to address detailed functional biological questions. In this article, we discuss how errors in these models arise and how they can be evaluated, and we argue for even more stringent validation checks and documentation of structures before deposition with the Protein Data Bank.