Helix stability of oligoglycine,oligoalanine, and oligo‐β‐alanine dodecamers reflected by hydrogen‐bond persistence

Helices are important structural/recognition elements in proteins and peptides. Stability and conformational differences between helices composed of α‐ and β‐amino acids as scaffolds for mimicry of helix recognition has become a theme in medicinal chemistry. Furthermore, helices formed by β‐amino acids are experimentally more stable than those formed by α‐amino acids. This is paradoxical because the larger sizes of the hydrogen‐bonding rings required by the extra methylene groups should lead to entropic destabilization. In this study, molecular dynamics simulations using the second‐generation force field, AMOEBA (Ponder, J.W., et al., Current status of the AMOEBA polarizable force field. J Phys Chem B, 2010. 114 (8): p. 2549–64.) explored the stability and hydrogen‐bonding patterns of capped oligo‐β‐alanine, oligoalanine, and oligoglycine dodecamers in water. The MD simulations showed that oligo‐β‐alanine has strong acceptor+2 hydrogen bonds, but surprisingly did not contain a large content of 3₁₂‐helical structures, possibly due to the sparse distribution of the 3₁₂‐helical structure and other structures with acceptor+2 hydrogen bonds. On the other hand, despite its backbone flexibility, the β‐alanine dodecamer had more stable and persistent <3.0 Å hydrogen bonds. Its structure was dominated more by multicentered hydrogen bonds than either oligoglycine or oligoalanine helices. The 3₁ (PII) helical structure, prevalent in oligoglycine and oligoalanine, does not appear to be stable in oligo‐β‐alanine indicating its competition with other structures (stacking structure as indicated by MD analyses). These differences are among the factors that shape helical structural preferences and the relative stabilities of these three oligopeptides. Proteins 2014; 82:3043–3061. © 2014 Wiley Periodicals, Inc.     

The structure of the aquaporin-1 water channel: a comparison between cryo-electron microscopy and X-ray crystallography

Three different medium-resolution structures of the human water channel aquaporin-1 (AQP1) have been solved by cryo-electron microscopy (cryo-EM) during the last two years. Recently, the structure of the strongly related bovine AQP1 was solved by X-ray crystallography at higher resolution, allowing a validation of the original medium-resolution structures, and providing a good indication for the strengths and limitations of state of the art cryo-EM methods. We present a detailed comparison between the different models, which shows that overall, the structures are highly similar, deviating less than 2.5 A from each other in the helical backbone regions. The two original cryo-EM structures, however, also show a number of significant deviations from the X-ray structure, both in the backbone positions of the transmembrane helices and in the location of the amino acid side-chains facing the pore. In contrast, the third cryo-EM structure that included information from the X-ray structure of the homologous bacterial glycerol facilitator GlpF and that was subsequently refined against cryo-EM AQP1 data, shows a root mean square deviation of 0.9A from the X-ray structure in the helical backbone regions.     

The elusive π-helix

Central to protein architecture is the local arrangement or secondary structure of the polypeptide backbone. Thirty to forty percent of protein domains are α-helices with 3.6 residues per turn. π-Helices, in which the peptide chain is more loosely coiled (4.4 residues per turn), have also been proposed. However, such structures necessitate an energetically unfavorable ～1 Å central helical hole. We show that rather than being composed of idealized π-helices, helical regions formed from putative π-helices actually consist of a series of concatenated wide turns with unique elliptical configurations. These structures have a larger helical radius akin to that of a π-helix, but without the loss of favorable cross-core van der Waals interactions. This not only obviates the helical void, but also endows proteins with important functionalities, including metal ion coordination, enhanced flexibility and specific enzyme-substrate binding interactions.     

A parallel stranded G‐quadruplex composed of threose nucleic acid (TNA)

G‐rich sequences can adopt four‐stranded helical structures, called G‐quadruplexes, that self‐assemble around monovalent cations like sodium (Na⁺) and potassium (K⁺). Whether similar structures can be formed from xeno‐nucleic acid (XNA) polymers with a shorter backbone repeat unit is an unanswered question with significant implications on the fold space of functional XNA polymers. Here, we examine the potential for TNA (α‐l ‐threofuranosyl nucleic acid) to adopt a four‐stranded helical structure based on a planar G‐quartet motif. Using native polyacrylamide gel electrophoresis (PAGE), circular dichroism (CD) and solution‐state nuclear magnetic resonance (NMR) spectroscopy, we show that despite a backbone repeat unit that is one atom shorter than the backbone repeat unit found in DNA and RNA, TNA can self‐assemble into stable G‐quadruplex structures that are similar in thermal stability to equivalent DNA structures. However, unlike DNA, TNA does not appear to discriminate between Na⁺ and K⁺ ions, as G‐quadruplex structures form equally well in the presence of either ion. Together, these findings demonstrate that despite a shorter backbone repeat unit, TNA is capable of self‐assembling into stable G‐quadruplex structures.     

Geometry motivated alternative view on local protein backbone structures

We present an alternative to the classical Ramachandran plot (R-plot) to display local protein backbone structure. Instead of the (ϕ, ψ)-backbone angles relating to the chemical architecture of polypeptides generic helical parameters are used. These are the rotation or twist angle ϑ and the helical rise parameter d. Plots with these parameters provide a different view on the nature of local protein backbone structures. It allows to display the local structures in polar (d, ϑ)-coordinates, which is not possible for an R-plot, where structural regimes connected by periodicity appear disconnected. But there are other advantages, like a clear discrimination of the handedness of a local structure, a larger spread of the different local structure domains—the latter can yield a better separation of different local secondary structure motives—and many more. Compared to the R-plot we are not aware of any major disadvantage to classify local polypeptide structures with the (d, ϑ)-plot, except that it requires some elementary computations. To facilitate usage of the new (d, ϑ)-plot for protein structures we provide a web application (http://agknapp.chemie.fu-berlin.de/secsass), which shows the (d, ϑ)-plot side-by-side with the R-plot.     

Solution structure of peptides from HIV-1 Vpr protein that cause membrane permeabilization and growth arrest

Vpr, one of the accessory gene products encoded by HIV-1, is a 96-residue protein with a number of functions, including targeting of the viral pre-integration complex to the nucleus and inducing growth arrest of dividing cells. We have characterized by 2D NMR the solution conformations of bioactive synthetic peptide fragments of Vpr encompassing a pair of H(F/S)RIG sequence motifs (residues 71–75 and 78–82 of HIV-1 Vpr) that cause cell membrane permeabilization and death in yeast and mammalian cells. Due to limited solubility of the peptides in water, their structures were studied in aqueous trifluoroethanol. Peptide Vpr^59–86 (residues 59–86 of Vpr) formed an α-helix encompassing residues 60–77, with a kink in the vicinity of residue 62. The first of the repeated sequence motifs (HFRIG) participated in the well-defined α-helical domain whereas the second (HSRIG) lay outside the helical domain and formed a reverse turn followed by a less ordered region. On the other hand, peptides Vpr^71–82 and Vpr^71–96, in which the sequence motifs were located at the N-terminus, were largely unstructured under similar conditions, as judged by their C^αH chemical shifts. Thus, the HFRIG and HSRIG motifs adopt α-helical and turn structures, respectively, when preceded by a helical structure, but are largely unstructured in isolation. The implications of these findings for interpretation of the structure–function relationships of synthetic peptides containing these motifs are discussed. © 1998 European Peptide Society and John Wiley & Sons, Ltd.     

Intrinsic Bent DNA Sites in the Developmentally Amplified C3-22 Gene Promoter of Rhynchosciara americana (Diptera: Sciaridae)

The Rhynchosciara americana C3-22 gene is located in an amplified domain and is developmentally expressed. The aim of the present work was to identify intrinsically bent DNA sites in a segment containing the gene promoter and downstream sequence. The results indicated that this gene is flanked by intrinsically bent DNA sites. Three bent DNA sites (b^?3, b^?2, and b^?1) were localized in the promoter, and one was localized downstream of the gene (b⁺¹). These sites had helical parameters that confirmed the curved structure, as well as segments with left-handed superhelical writhe. In silico analysis of the promoters of four other insect genes, which encode secreted polypeptides, showed that they all had curved structures and similar helical parameters. Correlation with other results indicates that the detected intrinsically bent DNA sites that flank the C3-22 gene might be a consensus feature of the gene structure in the amplified domains.     

Effect of acetylation and permethylation on the conformation and candidacidal activity of salivary histatin-5

Summary Salivary histatins provide the non-immune defense against oral pathogens such as Candida albicans. The structural requirements of histatin-5 for anticandida activity were examined with respect to its ability to adopt helical structures, its electrostatic interactions and the hydrogen-bonding potency of its basic residues. For this purpose, the lysine and/or histidine residues of histatin-5 were chemically modified by acetylation and permethylation. Acetylated histatin-5 retained its ability to adopt helical structures in 2,2,2-trifluoroethanol, but completely lost its ability to kill yeast cells. In contrast, permethylated histatin-5 shows very little tendency to adopt a helical structure, but retained significant anticandida activity. The results suggest that the candidacidal activity can arise even when the histatin does not have the ability to adopt helical structures. The candidacidal activity of the derivatives is discussed in terms of electrostatic interactions and hydrogen-bonding potency.     

Synthèse et structure de la poly-O-carbobenzoxy L- tyrosine et de copolymères de O-carbobenzoxy L-tyrosine et de glutamate de benzyle ou d'aspartate de benzyle

The optical rotatory dispersion of copolymers of O-carbobenzoxy-L -tyrosine and benzyl L (or D )-glutamate as well as benzyl L -aspartate, dissolved in nonpolar solvent, has been studied. Moffitt's equation permits the determination of b₀ coefficients whose variation, with varying composition in amino acid residues, suggests that the molecules of poly-O-carbobenzoxy-L -tyrosine have a helical structure similar to that of poly-(benzyl L -glutamate). Results obtained from infrared spectroscopy and x-ray diffraction show that the copolymers possess a helical conformation in the solid state, even when they are very rich in carbobenzoxy-L -tyrosine residues. The value of the b₀, coefficient for poly-O-carbobenzoxy-L -tyrosine may be explained by a regular stacking of the chromophore groups around the helical backbone. The ordering of the molecules of this polymer in a purely helical structure seems favored by the insertion of a small number of foreign residues in the polypeptide chain.     

G219S mutagenesis as a means of stabilizing conformational flexibility in the bacterial sodium channel NaChBac

The NaChBac sodium channel from Bacillus halodurans is a homologue of eukaryotic voltage-gated sodium channels. It can be solubilized in a range of detergents and consists of four identical subunits assembled as a tetramer. Sodium channels are relatively flexible molecules, adopting different conformations in their closed, open and inactivated states. This study aimed to design and construct a mutant version of the NaChBac protein that would insert into membranes and retain its folded conformation, but which would have enhanced stability when subjected to thermal stress. Modelling studies suggested a G219S mutant would have decreased conformational flexibility due to the removal of the glycine hinge around the proposed gating region, thereby imparting increased resistance to unfolding. The mutant expressed in Escherichia coli and purified in the detergent dodecyl maltoside was compared to wildtype NaChBac prepared in a similar manner. The mutant was incorporated into the membrane fraction and had a nearly identical secondary structure to the wildtype protein. When the thermal unfolding of the G219S mutant was examined by circular dichroism spectroscopy, it was shown to not only have a Tm ～10°C higher than the wildtype, but also in its unfolded state it retained more ordered helical structure than did the wildtype protein. Hence the G219S mutant was shown to be, as designed, more thermally stable.     

Stabilization of Mikamycin B Lactonase

Myoglobin (Mb) purified from fast skeletal muscle of bluefin tuna Thunnus thynnus orientalis was subjected to thermal treatment, and the denaturation profiles were examined by thermodynamic analysis. Based on the ellipticity or helical content obtained by circular dichroism (CD) spectrometry, it was found that denaturation of tuna Mb consisted of three steps, and that slight structural changes of Mb started below 20 °C. However, major structural changes were observed at around 58 and 72 °C. Differential scanning calorimetry (DSC) analysis revealed a similar but somewhat different thermal denaturation profile of Mb. In comparison with the denaturing profiles of whale Mb under the same conditions, the thermal stability of tuna Mb was found to be much lower. In the modeled tertiary structures of these Mbs, they were roughly similar to each other, though minor conformational differences were recognized and the total energy was found to be lower for tuna Mb.     

Crystal structure of a tripeptide containing aminocyclododecane carboxylic acid: a supramolecular twisted parallel β‐sheet in crystals

The crystal structure of a tripeptide Boc‐Leu‐Val‐Ac₁₂c‐OMe ( 1 ) is determined, which incorporates a bulky 1‐aminocyclododecane‐1‐carboxylic acid (Ac₁₂c) side chain. The peptide adopts a semi‐extended backbone conformation for Leu and Val residues, while the backbone torsion angles of the C^α,α‐dialkylated residue Ac₁₂c are in the helical region of the Ramachandran map. The molecular packing of 1 revealed a unique supramolecular twisted parallel β‐sheet coiling into a helical architecture in crystals, with the bulky hydrophobic Ac₁₂c side chains projecting outward the helical column. This arrangement resembles the packing of peptide helices in crystal structures. Although short oligopeptides often assemble as parallel or anti‐parallel β‐sheet in crystals, twisted or helical β‐sheet formation has been observed in a few examples of dipeptide crystal structures. Peptide 1 presents the first example of a tripeptide showing twisted β‐sheet assembly in crystals. Copyright © 2016 European Peptide Society and John Wiley & Sons, Ltd.     

Backbone solution structures of proteins using residual dipolar couplings: Application to a novel structural genomics target

Structural genomics (or proteomics) activities are critically dependent on the availability of high-throughput structure determination methodology. Development of such methodology has been a particular challenge for NMR based structure determination because of the demands for isotopic labeling of proteins and the requirements for very long data acquisition times. We present here a methodology that gains efficiency from a focus on determination of backbone structures of proteins as opposed to full structures with all sidechains in place. This focus is appropriate given the presumption that many protein structures in the future will be built using computational methods that start from representative fold family structures and replace as many as 70% of the sidechains in the course of structure determination. The methodology we present is based primarily on residual dipolar couplings (RDCs), readily accessible NMR observables that constrain the orientation of backbone fragments irrespective of separation in space. A new software tool is described for the assembly of backbone fragments under RDC constraints and an application to a structural genomics target is presented. The target is an 8.7 kDa protein from Pyrococcus furiosus, PF1061, that was previously not well annotated, and had a nearest structurally characterized neighbor with only 33% sequence identity. The structure produced shows structural similarity to this sequence homologue, but also shows similarity to other proteins, which suggests a functional role in sulfur transfer. Given the backbone structure and a possible functional link this should be an ideal target for development of modeling methods. This revised version was published online in March 2005 with corrections to the references.     

Probing designability via a generalized model of helical bundle geometry

Because the space of folded protein structures is highly degenerate, with recurring secondary and tertiary motifs, methods for representing protein structure in terms of collective physically relevant coordinates are of great interest. By collapsing structural diversity to a handful of parameters, such methods can be used to delineate the space of designable structures (i.e., conformations that can be stabilized with a large number of sequences)—a crucial task for de novo protein design. We first demonstrate this on natural α-helical coiled coils using the Crick parameterization. We show that over 95% of known coiled-coil structures are within  1-Å C_α root mean square deviation of a Crick-ideal backbone. Derived parameters show that natural geometric space of coiled coils is highly restricted and can be represented by “allowed” conformations amidst a potential continuum of conformers. Allowed structures have (1) restricted axial offsets between helices, which differ starkly between parallel and anti-parallel structures; (2) preferred superhelical radii, which depend linearly on the oligomerization state; (3) pronounced radius-dependent a- and d-position amino acid propensities; and (4) discrete angles of rotation of helices about their axes, which are surprisingly independent of oligomerization state or orientation. In all, we estimate the space of designable coiled-coil structures to be reduced at least 160-fold relative to the space of geometrically feasible structures. To extend the benefits of structural parameterization to other systems, we developed a general mathematical framework for parameterizing arbitrary helical structures, which reduces to the Crick parameterization as a special case. The method is successfully validated on a set of non-coiled-coil helical bundles, frequent in channels and transporter proteins, which show significant helix bending but not supercoiling. Programs for coiled-coil parameter fitting and structure generation are provided via a web interface at http://www.gevorggrigoryan.com/cccp/, and code for generalized helical parameterization is available upon request.     

Interactions of molecules with nucleic acids. I. An algorithm to generate nucleic acid structures with an application to the B-DNA structure and a counterclockwise helix.

An algorithm is developed that enables the routine determination of backbone conformations of nucleic acids. All atomic positions including hydrogen are specified in accord with experimental bond lengths and angles but with theoretically determined conformational angles. For two Watson-Crick base pairs at a separation of 3.38 Å, and perpendicular to a common helical axis, minimum energy configurations are found for all 10 combinations at helical angles of α ～ 36°–38°, corresponding to the B-DNA structure with C(2′)-endo sugar puckers. Backbone configurations exist only within the range 35.5° ? α ? 42°, which suggests the origin of the 10-fold helix. Calculated stacking energies for the B-DNA structure increases for each of the clustered groups of base pairs: G·C with G·C, G·C with A·T, and A·T with A·T, and they are in approximate agreement with experimental observations. The counter-clockwise helix is examined, and physically meaningful structures are found only when the helical axes of successive base pairs are disjointed.     

Similarity of Domain Organization of Proteins in Phylogenetically Distant Organisms as a Basis of Meiosis Conservatism

The cytological mechanism of meiosis is very conservative in all eukaryotes. Some meiosis-specific structural proteins of yeast, nematode Caenorhabditis elegans, Drosophila, and mammals, which play identical roles in cells during meiosis, do not have homology of the primary structure, but their domain organization and conformation are similar. The enzymes of meiotic recombination in yeast and plants have similar epitopes. These facts suggest that the similarity of the higher level of organization of the meiosis-specific proteins allows these proteins to form similar subcellular structures and produce similar cytological picture of meiosis and similar functions of these subcellular structures. Finally, this leads to a conservative scheme of meiosis in evolutionally distant eukaryotes.     

A molecular dynamics simulation of the flavin mononucleotide–RNA aptamer complex

We report on an unrestrained molecular dynamics simulation of the flavin mononucleotide (FMN)–RNA aptamer. The simulated average structure maintains both cross‐strand and intermolecular FMN–RNA nuclear Overhauser effects from the nmr experiments and has all qualitative features of the nmr structure including the G10–U12–A25 base triple and the A13–G24, A8–G28, and G9–G27 mismatches. However, the relative orientation of the hairpin loop to the remaining part of the molecule differs from the nmr structure. The simulation predicts that the flexible phosphoglycerol part of FMN moves toward G27 and forms hydrogen bonds. There are structurally long‐lived water molecules in the FMN binding pocket forming hydrogen bonds within FMN and between FMN and RNA. In addition, long‐lived water is found bridging primarily RNA backbone atoms. A general feature of the environment of long‐lived “structural” water is at least two and in most cases three or four potential acceptor atoms. The 2′‐OH group of RNA usually acts as an acceptor in interactions with the solvent. There are almost no intrastrand O2′H(n)⋮O4′(n + 1) hydrogen bonds within the RNA backbone. In the standard case the preferred orientation of the 2′‐OH hydrogen atoms is approximately toward O3′ of the same nucleotide. However, a relatively large number of conformations with the backbone torsional angle γ in the trans orientation is found. A survey of all experimental RNA x‐ray structures shows that this backbone conformation occurs but is less frequent than found in the simulation. Experimental nmr RNA aptamer structures have a higher fraction of this conformation as compared to the x‐ray structures. The backbone conformation of nucleotide n + 1 with the torsional angle γ in the trans orientation leads to a relatively short distance between 2′‐OH(n) and O5′(n + 1), enabling hydrogen‐bond formation. In this case the preferred orientation of the 2′‐OH hydrogen atom is approximately toward O5′(n + 1). We find two relatively short and dynamically stable types of backbone–backbone next‐neighbor contacts, namely C2′(H)(n)⋮O4′(n + 1) and C5′(H)(n + 1)⋮O2′(n). These interactions may affect both backbone rigidity and thermodynamic stability of RNA helical structures. © 1999 John Wiley & Sons, Inc. Biopoly 50: 287–302, 1999     

Multiple helical conformations of the helix‐turn‐helix region revealed by NOE‐restrained MD simulations of tryptophan aporepressor,TrpR

The nature of flexibility in the helix‐turn‐helix region of E. coli trp aporepressor has been unexplained for many years. The original ensemble of nuclear magnetic resonance (NMR structures showed apparent disorder, but chemical shift and relaxation measurements indicated a helical region. Nuclear Overhauser effect (NOE) data for a temperature‐sensitive mutant showed more helical character in its helix‐turn‐helix region, but nevertheless also led to an apparently disordered ensemble. However, conventional NMR structure determination methods require all structures in the ensemble to be consistent with every NOE simultaneously. This work uses an alternative approach in which some structures of the ensemble are allowed to violate some NOEs to permit modeling of multiple conformational states that are in dynamic equilibrium. Newly measured NOE data for wild‐type aporepressor are used as time‐averaged distance restraints in molecular dynamics simulations to generate an ensemble of helical conformations that is more consistent with the observed NMR data than the apparent disorder in the previously reported NMR structures. The results indicate the presence of alternating helical conformations that provide a better explanation for the flexibility of the helix‐turn‐helix region of trp aporepressor. Structures representing these conformations have been deposited with PDB ID: 5TM0. Proteins 2017; 85:731–740. © 2016 Wiley Periodicals, Inc.     

Beyond the Twilight Zone: Automated prediction of structural properties of proteins by recursive neural networks and remote homology information

The prediction of 1D structural properties of proteins is an important step toward the prediction of protein structure and function, not only in the ab initio case but also when homology information to known structures is available. Despite this the vast majority of 1D predictors do not incorporate homology information into the prediction process. We develop a novel structural alignment method, SAMD, which we use to build alignments of putative remote homologues that we compress into templates of structural frequency profiles. We use these templates as additional input to ensembles of recursive neural networks, which we specialise for the prediction of query sequences that show only remote homology to any Protein Data Bank structure. We predict four 1D structural properties – secondary structure, relative solvent accessibility, backbone structural motifs, and contact density. Secondary structure prediction accuracy, tested by five‐fold cross‐validation on a large set of proteins allowing less than 25% sequence identity between training and test set and query sequences and templates, exceeds 82%, outperforming its ab initio counterpart, other state‐of‐the‐art secondary structure predictors (Jpred 3 and PSIPRED) and two other systems based on PSI‐BLAST and COMPASS templates. We show that structural information from homologues improves prediction accuracy well beyond the Twilight Zone of sequence similarity, even below 5% sequence identity, for all four structural properties. Significant improvement over the extraction of structural information directly from PDB templates suggests that the combination of sequence and template information is more informative than templates alone. Proteins 2009. © 2009 Wiley‐Liss, Inc.