De novo determination of protein backbone structure from residual dipolar couplings using Rosetta

As genome-sequencing projects rapidly increase the database of protein sequences, the gap between known sequences and known structures continues to grow exponentially, increasing the demand to accelerate structure determination methods. Residual dipolar couplings (RDCs) are an attractive source of experimental restraints for NMR structure determination, particularly rapid, high-throughput methods, because they yield both local and long-range orientational information and can be easily measured and assigned once the backbone resonances of a protein have been assigned. While very extensive RDC data sets have been used to determine the structure of ubiquitin, it is unclear to what extent such methods will generalize to larger proteins with less complete data sets. Here we incorporate experimental RDC restraints into Rosetta, an ab initio structure prediction method, and demonstrate that the combined algorithm provides a general method for de novo determination of a variety of protein folds from RDC data. Backbone structures for multiple proteins up to approximately 125 residues in length and spanning a range of topological complexities are rapidly and reproducibly generated using data sets that are insufficient in isolation to uniquely determine the protein fold de novo, although ambiguities and errors are observed for proteins with symmetry about an axis of the alignment tensor. The models generated are not high-resolution structures completely defined by experimental data but are sufficiently accurate to accelerate traditional high-resolution NMR structure determination and provide structure-based functional insights.     

Computational design of the sequence and structure of a protein-binding peptide

The de novo design of protein-binding peptides is challenging because it requires the identification of both a sequence and a backbone conformation favorable for binding. We used a computational strategy that iterates between structure and sequence optimization to redesign the C-terminal portion of the RGS14 GoLoco motif peptide so that it adopts a new conformation when bound to Gα(i1). An X-ray crystal structure of the redesigned complex closely matches the computational model, with a backbone root-mean-square deviation of 1.1 ?.     

Isolated monohydrates of a model peptide chain: effect of a first water molecule on the secondary structure of a capped phenylalanine

The formation of monohydrates of capped phenylalanine model peptides, CH(3)-CO-Phe-NH(2) and CH(3)-CO-Phe-NH-CH(3), in a supersonic expansion has been investigated using laser spectroscopy and quantum chemistry methods. Conformational distributions of the monohydrates have been revealed by IR/UV double-resonance spectroscopy and their structures assigned by comparison with DFT-D calculations. A careful analysis of the final hydrate distribution together with a detailed theoretical investigation of the potential energy surface of the monohydrates demonstrates that solvation occurs from the conformational distribution of the isolated peptide monomers. The distribution of the monohydrates appears to be strongly dependent on both the initial monomer conformation (extended or folded backbone) and the solvation site initially occupied by the water molecule. The solvation processes taking place during the cooling can be categorized as follows: (a) solvation without significant structural changes of the peptide, (b) solvation inducing significant distortions of the backbone but retaining the secondary structure, and (c) solvation triggering backbone isomerizations, leading to a modification of the peptide secondary structure. It is observed that solvation by a single water molecule can fold a β-strand into a γ-turn structure (type c) or induce a significant opening of a γ-turn characterized by an elongated C(7) hydrogen bond (type b). These structural changes can be considered as a first step toward the polyproline II condensed-phase structure, illustrating the role played by the very first water molecule in the solvation process.     

Stability of cyclic beta-hairpins: asymmetric contributions from side chains of a hydrogen-bonded cross-strand residue pair

Amino acid structural propensities measured in "host-guest" model studies are often used in protein structure prediction or to choose appropriate residues in de novo protein design. While this concept has proven useful for helical structures, it is more difficult to apply successfully to beta-sheets. We have developed a cyclic beta-hairpin scaffold as a host for measurement of individual residue contributions to hairpin structural stability. Previously, we have characterized substitutions in non-backbone-hydrogen-bonded strand sites; relative stability differences measured in the cyclic host are highly predictive of changes in folding free energy for linear beta-hairpin peptides. Here, we examine the hydrogen-bonded strand positions of our host. Surprisingly, we find a large favorable contribution to stability from a valine (or isoleucine) substitution immediately preceding the C-terminal cysteine of the host peptide, but not at the cross-strand position of the host or in either strand of a folded linear beta-hairpin (trpzip peptide). Further substitutions in the peptides and NMR structural analysis indicate that the stabilizing effect of valine is general for CX(8)C cyclic hairpins and cannot be explained by particular side-chain-side-chain interactions. Instead, a localized decrease in twist of the peptide backbone on the N-terminal side of the cysteine allows the valine side chain to adopt a unique conformation that decreases the solvent accessibility of the peptide backbone. The conformation differs from the highly twisted (coiled) conformation of the trpzip hairpins and is more typical of conformations present in multistranded beta-sheets. This unexpected structural fine-tuning may explain why cyclic hairpins selected from phage-displayed libraries often have valine in the same position, preceding the C-terminal cysteine. It also emphasizes the diversity of structures accessible to beta-strands and the importance of considering not only "beta-propensity", but also hydrogen-bonding pattern and strand twist, when designing beta structures. Finally, we observe correlated, cooperative stabilization from side-chain substitutions on opposite faces of the hairpin. This suggests that cooperative folding in beta-hairpins and other small beta-structures is driven by cooperative strand-strand association.     

Conformational landscape of the HIV-V3 hairpin loop from all-atom free-energy simulations

Small beta hairpins have many distinct biological functions, including their involvement in chemokine and viral receptor recognition. The relevance of structural similarities between different hairpin loops with near homologous sequences is not yet understood, calling for the development of methods for de novo hairpin structure prediction and simulation. De novo folding of beta strands is more difficult than that of helical proteins because of nonlocal hydrogen bonding patterns that connect amino acids that are distant in the amino acid sequence and there is a large variety of possible hydrogen bond patterns. Here we use a greedy version of the basin hopping technique with our free-energy forcefield PFF02 to reproducibly and predictively fold the hairpin structure of a HIV-V3 loop. We performed 20 independent basin hopping runs for 500 cycles corresponding to 7.4 x 10(7) energy evaluations each. The lowest energy structure found in the simulation has a backbone root mean square deviation (bRMSD) of only 2.04 A to the native conformation. The lowest 9 out of the 20 simulations converged to conformations deviating less than 2.5 A bRMSD from native.     

Side chain orientation from methyl 1H-1H residual dipolar couplings measured in highly deuterated proteins

High-level deuteration is a prerequisite for the study of high molecular weight systems using liquid-state NMR. Here, we present new experiments for the measurement of proton-proton dipolar couplings in CH(2)D methyl groups of (13)C labeled, highly deuterated (70-80%) proteins. (1)H-(1)H residual dipolar couplings (RDCs) have been measured in two alignment media for 57 out of 70 possible methyl containing residues in the 167-residue flavodoxin-like domain of the E. coli sulfite reductase. These data yield information on the orientation of the methyl symmetry axis with respect to the molecular alignment frame. The alignment tensor characteristics were obtained very accurately from a set of backbone RDCs measured on the same protein sample. To demonstrate that accurate structural information is obtained from these data, the measured methyl RDCs for Valine residues are analyzed in terms of chi(1) torsion angles and stereospecific assignment of the prochiral methyl groups. On the basis of the previously determined backbone solution structure of this protein, the methyl RDC data proved sufficient to determine the chi(1) torsion angles in seven out of nine valines, assuming a single-rotamer model. Methyl RDCs are complementary to other NMR data, for example, methyl-methyl NOE, to determine side chain conformation in high molecular weight systems.     

Protein backbone structure determination using RDC: An inverse kinematics approach with fast and exact solutions

Residual dipolar couplings (RDC) of proteins dissolved in anisotropic media promise to speed up the determination of protein structures. We consider the backbone as a robotic mechanism and formulate inverse kinematics problems using RDC restraints from two media. The φ, ψ of each secondary structure element (SSE) are computed from oriented vectors in consecutive peptide planes. We search for the optimum conformation joining the solutions of two independent backbone halves. The matrix transforming the vector Z of a global frame from one SSE into the other determines their orientation. Three distance constraints between two oriented SSE determine their relative position by solving nine polynomial equations. The benefit of this method is that complete and accurate solutions are obtained overcoming the local minima problems of heuristic procedures. The algorithm is implemented on MAPLE using the least number of experimental data; the runtimes take an order of seconds on a common PC. © 2013 Wiley Periodicals, Inc.     

Refinement of local and long-range structural order in theophylline-binding RNA using (13)C-(1)H residual dipolar couplings and restrained molecular dynamics.

13C-(1)H residual dipolar couplings (RDC) have been measured for the bases and sugars in the theophylline-binding RNA aptamer, dissolved in filamentous phage medium, and used to investigate the long-range structural and dynamic behavior of the molecule in the solution state. The orientation dependent RDC provide additional restraints to further refine the overall structure of the RNA-theophylline complex, whose long-range order was poorly defined in the NOE-based structural ensemble. Structure refinement using RDC normally assumes that molecular alignment can be characterized by a single tensor and that the molecule is essentially rigid. To address the validity of this assumption for the complex of interest, we have analyzed distinct domains of the RNA molecule separately, so that local structure and alignment tensors experienced by each region are independently determined. Alignment tensors for the stem regions of the molecule were allowed to float freely during a restrained molecular dynamics structure refinement protocol and found to converge to similar magnitudes. During the second stage of the calculation, a single alignment tensor was thus applied for the whole molecule and an average molecular conformation satisfying all experimental data was determined. Semirigid-body molecular dynamics calculations were used to reorient the refined helical regions to a relative orientation consistent with this alignment tensor, allowing determination of the global conformation of the molecule. Simultaneously, the local structure of the theophylline-binding core of the molecule was refined under the influence of this common tensor. The final ensemble has an average pairwise root mean square deviation of 1.50 +/- 0.19 A taken over all heavy atoms, compared to 3.5 +/- 1.1 A for the ensemble determined without residual dipolar coupling. This study illustrates the importance of considering both the local and long-range nature of RDC when applying these restraints to structure refinements of nucleic acids.     

Three‐dimensional structure of cyclic antibiotic teicoplanin aglycone using NMR distance and dihedral angle restraints in a DMSO solvation model

The three‐dimensional solution conformation of teicoplanin aglycone was determined using NMR spectroscopy. A combination of NOE and dihedral angle restraints in a DMSO solvation model was used to calculate an ensemble of structures having a root mean square deviation of 0.17 Å. The structures were generated using systematic searches of conformational space for optimal satisfaction of distance and dihedral angle restraints. Comparison of the NMR‐derived structure of teicoplanin aglycone with the X‐ray structure of a teicoplanin aglycone analog revealed a common backbone conformation with deviation of two aromatic side chain substituents. Experimentally determined backbone ¹³C chemical shifts showed good agreement with those computed at the density functional level of theory, providing a cross validation of the backbone conformation. The flexible portion of the molecule was consistent with the region that changes conformation to accommodate protein binding. The results showed that a hydrogen‐bonded DMSO molecule in combination with NMR‐derived restraints together enabled calculation of structures that satisfied experimental data. Copyright © 2015 John Wiley & Sons, Ltd.     

Template-directed assembly of a de novo designed protein

Many naturally occurring biomaterials are composed of laminated structures in which layers of beta-sheet proteins alternate with layers of inorganic mineral. These ordered laminates often have structural and mechanical properties that differ significantly from those of nonbiological materials. An important step in the construction of novel biomaterials is the creation of composites wherein a de novo designed protein assembles into an ordered structure. To achieve this goal, we layered a de novo protein onto the surface of highly ordered pyrolytic graphite (HOPG). The protein was derived from a combinatorial library of novel sequences designed to fold into amphiphilic beta-sheet structures. Atomic force microscopy reveals that the protein assembles on the HOPG surface into ordered fibers aligned in three orientations at 120 degrees to each other. The symmetry and extent of the ordered regions indicate that the hexagonal lattice underlying the graphite surface templates assembly of millions of protein molecules into a highly ordered structure.     

Magic angle spinning solid-state NMR spectroscopy for structural studies of protein interfaces. resonance assignments of differentially enriched Escherichia coli thioredoxin reassembled by fragment complementation

De novo site-specific backbone and side-chain resonance assignments are presented for U-15N(1-73)/U-13C,15N(74-108) reassembly of Escherichia coli thioredoxin by fragment complementation, determined using solid-state magic angle spinning NMR spectroscopy at 17.6 T. Backbone dihedral angles and secondary structure predicted from the statistical analysis of 13C and 15N chemical shifts are in general agreement with solution values for the intact full-length thioredoxin, confirming that the secondary structure is retained in the reassembled complex prepared as a poly(ethylene glycol) precipitate. The differential labeling of complementary thioredoxin fragments introduced in this work is expected to be beneficial for high-resolution structural studies of protein interfaces formed by protein assemblies by solid-state NMR spectroscopy.     

Prediction of Protein Structure Using Surface Accessibility Data

An approach to the de novo structure prediction of proteins is described that relies on surface accessibility data from NMR paramagnetic relaxation enhancements by a soluble paramagnetic compound (sPRE). This method exploits the distance‐to‐surface information encoded in the sPRE data in the chemical shift‐based CS‐Rosetta de novo structure prediction framework to generate reliable structural models. For several proteins, it is demonstrated that surface accessibility data is an excellent measure of the correct protein fold in the early stages of the computational folding algorithm and significantly improves accuracy and convergence of the standard Rosetta structure prediction approach.     

Parallel β-sheet secondary structure is stabilized and terminated by interstrand disulfide cross-linking

Disulfide bonds between Cys residues in adjacent strands of parallel β-sheets are rare among proteins, which suggests that parallel β-sheet structure is not stabilized by such disulfide cross-links. We report experimental results that show, surprisingly, that an interstrand disulfide bond can stabilize parallel β-sheets formed by an autonomously folding peptide in aqueous solution. NMR analysis reveals that parallel β-sheet structure is terminated beyond the disulfide bond, which causes deviation from the extended backbone conformation at one of the Cys residues.     

Peptide backbone fragmentation initiated by side‐chain loss at cysteine residue in matrix‐assisted laser desorption/ionization in‐source decay mass spectrometry

Matrix‐assisted laser desorption/ionization in‐source decay (MALDI‐ISD) is initiated by hydrogen transfer from matrix molecules to the carbonyl oxygen of peptide backbone with subsequent radical‐induced cleavage leading to c′/z? fragments pair. MALDI‐ISD is a very powerful method to obtain long sequence tags from proteins or to do de novo sequencing of peptides. Besides classical fragmentation, MALDI‐ISD also shows specific fragments for which the mechanism of formation enlightened the MALDI‐ISD process. In this study, the MALDI‐ISD mechanism is reviewed, and a specific mechanism is studied in details: the N‐terminal side of Cys residue (Xxx‐Cys) is described to promote the generation of c′ and w fragments in MALDI‐ISD. Our data suggest that for sequences containing Xxx‐Cys motifs, the N–C_α bond cleavage occurs following the hydrogen attachment to the thiol group of Cys side‐chain. The c?/w fragments pair is formed by side‐chain loss of the Cys residue with subsequent radical‐induced cleavage at the N–C_α bond located at the left side (N‐terminal direction) of the Cys residue. This fragmentation pathway preferentially occurs at free Cys residue and is suppressed when the cysteines are involved in disulfide bonds. Hydrogen attachment to alkylated Cys residues using iodoacetamide gives free Cys residue by the loss of ?CH₂CONH₂ radical. The presence of alkylated Cys residue also suppress the formation of c?/w fragments pair via the (C_β)‐centered radical, whereas w fragment is still observed as intense signal. In this case, the z? fragment formed by hydrogen attachment of carbonyl oxygen followed side‐chain loss at alkylated Cys leads to a w fragment. Hydrogen attachment on peptide backbone and side‐chain of Cys residue occurs therefore competitively during MALDI‐ISD process. Copyright © 2013 John Wiley & Sons, Ltd.     

Microcin J25 has a threaded sidechain-to-backbone ring structure and not a head-to-tail cyclized backbone

Microcin J25 is a 21 amino acid bacterial peptide that has potent antibacterial activity against Gram-negative bacteria, resulting from its interaction with RNA polymerase. The peptide was previously proposed to have a head-to-tail cyclized peptide backbone and a tight globular structure (Blond, A., Péduzzi, J., Goulard, C., Chiuchiolo, M. J., Barthélémy, M., Prigent, Y., Salomón, R. A., Farías, R. N., Moreno, F. & Rebuffat, S. Eur. J. Biochem. 1999, 259, 747-755). It exhibits remarkable thermal stability for a peptide of its size lacking disulfide bonds and in part this was previously proposed to derive from its macrocyclic structure. We show here that in fact the peptide does not have a head-to-tail cyclic structure but rather a side chain to backbone cyclization between Glu8 and the N-terminus. This creates an embedded ring that is threaded by the C-terminal tail of the molecule, forming a noose-like feature. The three-dimensional structure deduced from NMR data suggests that slippage of the noose is prevented by two aromatic residues flanking the embedded ring. Unthreading does not occur even when the molecule is enzymatically digested with thermolysin. The new structural interpretation fully accounts for previously reported NMR and biophysical data and is consistent with the remarkable stability of this potent antimicrobial peptide.     

Structural studies of model peptides containing beta-, gamma- and delta-amino acids

The crystal structures of five model peptides Piv-Pro-Gly-NHMe (1), Piv-Pro-betaGly-NHMe (2), Piv-Pro-betaGly-OMe (3), Piv-Pro-deltaAva-OMe (4) and Boc-Pro-gammaAbu-OH (5) are described (Piv: pivaloyl; NHMe: N-methylamide; betaGly: beta-glycine; OMe: O-methyl ester; deltaAva: delta-aminovaleric acid; gammaAbu: gamma-aminobutyric acid). A comparison of the structures of peptides 1 and 2 illustrates the dramatic consequences upon backbone homologation in short sequences. 1 adopts a type II beta-turn conformation in the solid state, while in 2, the molecule adopts an open conformation with the beta-residue being fully extended. Piv-Pro-betaGly-OMe (3), which differs from 2 by replacement of the C-terminal NH group by an O-atom, adopts an almost identical molecular conformation and packing arrangement in the solid state. In peptide 4, the observed conformation resembles that determined for 2 and 3, with the deltaAva residue being fully extended. In peptide 5, the molecule undergoes a chain reversal, revealing a beta-turn mimetic structure stabilized by a C-H...O hydrogen bond.     

Beta-Ala containing peptides: potentials in design and construction of bioactive peptides and protein secondary structure mimics

In recent years, there has been increasing interest in de novo design and construction of novel synthetic peptides that mimic protein secondary structures, i.e., turns, helices and sheets. The unique structural influences exerted by unsubstituted, non-coded, non-chiral beta-amino acid, i.e., beta-alanine (beta-Ala; 3- or beta- aminopropionic acid) on peptide backbone, when inserted into peptide chain comprised alpha-amino acids, offer an excellent opportunity to design and construct diverse well-defined three-dimensional structures. Our current understanding of folding-unfolding behavior of the beta-Ala residues relies primarily from an examination of conformational preferences of a large number of short cyclic- as well as acyclic beta-Ala containing peptides investigated using single crystal X-ray diffraction analysis. In addition, theoretical conformational energy calculations and different spectroscopic techniques: 1H NMR, FT-IR and CD, have also been employed although, to a lesser extent. The obtainable results tend to reveal overwhelming preferences of the beta-Ala moiety for the folded gauche (mu approximately +/-65+/-10 degrees conformation in cyclic- and for an extended trans (mu approximately +/-165+/-10 degrees) as well as gauche (mu approximately +/-65+/-10 degrees) orientations in acyclic beta-Ala containing peptides. The results also indicate that in short linear beta-Ala containing peptides, the specific influence of selective neighboring side-chain substituents e.g. linear- or cyclic symmetrically C(alpha,alpha)-disubstituted glycines and other conformational constraints, may be significant in controlling the overall folded-unfolded topographical features across the two methylene units (-CbetaH2-CalphaH2-) of the beta-Ala residue. Taking into consideration the wide occurrence of beta-Ala moiety in animal and plant kingdoms and the remarkable structural versatility of the peptides incorporating beta-Ala residue(s), together with appreciable resistance towards enzymatic degradation, hold strong promise for biophysicists and biochemists not only to design molecules that fold to mimic protein secondary structures but also to develop potent peptide analogs and peptidomimetics displaying unique pharmaceutical properties.     

Examination of the folding of a short alanine-based helical peptide with salt bridges using molecular dynamics simulation

A molecular dynamics simulation of the folding of a short alanine-based helical peptide of 17 residues with three Glu...Lys (i, i + 4) salt bridge pairs, referred to as the AEK17 peptide, was carried out. The simulation gave an estimated simulation folding time of 2.5 ns, shorter than 12 ns for an alanine-based peptide of 16 residues with three Lys residues only, referred to as the AK16 peptide, simulated previously. After folded, the AEK17 peptide had a helical content of 77%, in excellent agreement with the experimentally determined value of 80%. An examination of the folding pathways of AEK17 indicated that the peptide proceeded via three-turn helix conformations more than the helix-turn-helix conformation in the folding pathways. An analysis of interactions indicated that the formation of hydrogen bonds between Lys residue side chains and backbone carbonyls is a major factor in the abundant conformation of the three-turn helix intermediate. The substitution of three Ala with Glu residues reduces the extent of hydrophobic interaction in alanine-based AK peptides with the result that the breaking of the interactions of Lys epsilon-NH3+(side chain)...C=O(backbone) is a major activation action for the AEK17 to achieve a complete fold, in contrast to the AK16 peptide, in which breaking non-native hydrophobic interaction is the rate-determining step.     

The gas-phase dipeptide analogue acetyl-phenylalanyl-amide: a model for the study of side chain/backbone interactions in proteins

The issue of the influence of the side chain/backbone interaction on the local conformational preferences of a phenylalanine residue in a peptide chain is addressed. A synergetic approach is used, which combines gas-phase UV spectroscopy as well as gas-phase IR/UV double-resonance experiments with DFT and post Hartree-Fock calculations. N-Acetyl-Phe-amide was chosen as a model system for which three different conformers were observed. The most stable conformer has been identified as an extended beta(L) conformation of the peptide backbone. It is stabilized by a weak but significant NH-pi interaction bridging the aromatic ring on the residue (i) with the NH group on residue (i+1), with the aromatic side chain being in an anti conformation. This stable conformation corresponds to the common NH(i+1)-aromatic(i) interaction encountered in proteins for the three aromatic residues (phenylalanine, tyrosine, and tryptophan), which illustrates the relevance of gas-phase investigations to structural biology issues. The two other less abundant conformers have been assigned to two gamma-folded backbone conformations that differ by the orientation of the side chain. In all cases, the IR data provided spectroscopic fingerprints of these interactions. Finally, the strong conformational dependence of the fluorescence yield found for N-acetyl-Phe-amide illustrates the role of the environment on the excited-state dynamics of these species, which is often exploited by biochemists to monitor protein structural changes from tryptophan lifetime measurements.     

Organic Ligand Binding by a Hydrophobic Cavity in a Designed Tetrameric Coiled‐Coil Protein

The design and characterization of a hydrophobic cavity in de novo designed proteins provides a wide range of information about the functions of de novo proteins. We designed a de novo tetrameric coiled‐coil protein with a hydrophobic pocketlike cavity. Tetrameric coiled coils with hydrophobic cavities have previously been reported. By replacing one Leu residue at the a position with Ala, hydrophobic cavities that did not flatten out due to loose peptide chains were reliably created. To perform a detailed examination of the ligand‐binding characteristics of the cavities, we originally designed two other coiled‐coil proteins: AM2, with eight Ala substitutions at the adjacent a and d positions at the center of a bundled structure, and AM2W, with one Trp and seven Ala substitutions at the same positions. To increase the association of the helical peptides, each helical peptide was connected with flexible linkers, which resulted in a single peptide chain. These proteins exhibited CD spectra corresponding to superhelical structures, despite weakened hydrophobic packing. AM2W exhibited binding affinity for size‐complementary organic compounds. The dissociation constants, K_d, of AM2W were 220 nM for adamantane, 81 μM for 1‐adamantanol, and 294 μM for 1‐adamantaneacetic acid, as measured by fluorescence titration analyses. Although it was contrary to expectations, AM2 did not exhibit any binding affinity, probably due to structural defects around the designed hydrophobic cavity. Interestingly, AM2W exhibited incremental structure stability through ligand binding. Plugging of structural defects with organic ligands would be expected to facilitate protein folding.