Secondary Structure Preferences Research Articles

Interface of Apoptotic Protein Complexes Has Distinct Properties

Apoptosis is a programmed mechanism of cell death that is a normal component of development and health of multi-cellular organisms. In this study, we ask if interface properties of apoptotic protein complexes are different from protein complexes in general. We find that although in apoptotic protein complexes the overall distribution of interface size, surface complementarity, hydrogen bonding, hydrophobicity are similar to general interface properties, apoptotic complexes tend to have more fragmented interfaces and different secondary structural preferences. The statistics on the number of interfaces where specific amino acid(s) occur with significantly enhanced frequency suggest that Arg, Met and Asp are most important functional residues. The role of Met is believed to be unique, as evidenced from the existing data on hot spot potential of residues. These findings together provide insight into the possible role of various physico-chemical attributes at the protein interface in regulation of the apoptosis process.

A reexamination of the propensities of amino acids towards a particular secondary structure: classification of amino acids based on their chemical structure

The correlation between the primary and secondary structures of proteins was analysed using a large data set from the Protein Data Bank. Clear preferences of amino acids towards certain secondary structures classify amino acids into four groups: alpha-helix preferrers, strand preferrers, turn and bend preferrers, and His and Cys (the latter two amino acids show no clear preference for any secondary structure). Amino acids in the same group have similar structural characteristics at their Cbeta and Cgamma atoms that predicts their preference for a particular secondary structure. All alpha-helix preferrers have neither polar heteroatoms on Cbeta and Cgamma atoms, nor branching or aromatic group on the Cbeta atom. All strand preferrers have aromatic groups or branching groups on the Cbeta atom. All turn and bend preferrers have a polar heteroatom on the Cbeta or Cgamma atoms or do not have a Cbeta atom at all. These new rules could be helpful in making predictions about non-natural amino acids.

Secondary structural investigations into homo-oligomers of δ-2,4- cis oxetane amino acids

Investigations into the secondary structural preferences of homo-oligomers of a new class of sugar amino acids, δ-2,4- cis-oxetane amino acids, are reported. The oligomers were seen to adopt a well-defined repeating β-turn structure in solution.

Oxetane amino acids: synthesis of tetrameric and hexameric carbopeptoids derived from l-ribo 4-(aminomethyl)-oxetan-2-carboxylic acid

The synthesis of methyl 2,4-anhydro-5-azido-3- O-benzyl-5-deoxy- l-ribonate, a δ-2,4- cis-oxetane-azido ester scaffold derived from l-arabinose, is reported. Iterative coupling methods were utilised to form homo-oligomers up to the hexamer in order to investigate the secondary structural preferences of these systems.

Self-contacts in Asx and Glx residues of high-resolution protein structures: Role of local environment and tertiary interactions

In protein structures, side-chains of asparagine and aspartic acid (Asx) and glutamine and glutamic acid (Glx) can approach their own backbone nitrogen or carbonyl group. We have systematically analyzed intra-residue contacts in Asx and Glx residues and their secondary structure preferences in two different datasets consisting of 500 and 1506 high-resolution structures. Intra-residue contact in an Asx/Glx residue between the heavy atoms of side-chain and main-chain functional groups of the same residue was investigated irrespective of whether such contacts are due to hydrogen bonding or not. Our search yielded 563 and 1462 cases of self-contacting Asx and Glx residues from the two datasets. Two important observations have been made in this analysis. First, self-contacts involving side-chain oxygen and backbone nitrogen atoms in majority of Asx residues are not due to hydrogen bonds. In the second instance, surprisingly, side-chain and backbone carbonyl oxygens of a significant number of Asx and Glx residues approach each other. For a wide-range of accessible surface areas, self-contacting residues are surrounded by less number of polar groups compared to all other Asx/Glx residues. In buried and partially buried regions, side-chain and main-chain functional groups of these residues together participate in simultaneous interactions with the available polar groups or water molecules. Asx/Glx residues with self-contacts are rarely observed in the middle of an α-helix or a β-strand. Asx/Glx side-chain having contact with its own backbone nitrogen shows different capping preferences compared to those having contact with its backbone oxygen. Examples of proteins with multiple self-contacting Asx/Glx residues are found. We speculate that mutation of a self-contacting residue in the buried or partially buried region of a protein will destabilize the structure. The results of this analysis will help in engineering protein structures and site-directed mutagenesis experiments.

How stable is a collagen triple helix? An ab initio study on various collagen and β‐sheet forming sequences

Collagen forms the well characterized triple helical secondary structure, stabilized by interchain H-bonds. Here we have investigated the stability of fully optimized collagen triple helices and beta-pleated sheets by using first principles (ab initio and DFT) calculations so as to determine the secondary structure preference depending on the amino acid composition. Models composed of a total of 18 amino acid residues were studied at six different amino acid compositions: (i) L-alanine only, (ii) glycine only, (iii) L-alanines and glycine, (iv) L-alanines and D-alanine, (v) L-prolines with glycine, (vi) L-proline, L-hydroxyproline, and glycine. The last two, v and vi, were designed to mimic the core part of collagen. Furthermore, ii, iii, and iv model the binding and/or recognition sites of collagen. Finally, i models the G-->A replacement, rare in collagen. All calculated structures show great resemblance to those determined by X-ray crystallography. Calculated triple helix formation affinities correlate well with experimentally determined stabilities derived from melting point (T(m)) data of different collagen models. The stabilization energy of a collagen triple helical structure over that of a beta-pleated sheet is 2.1 kcal mol(-1) per triplet for the [(-Pro-Hyp-Gly-)(2)](3) collagen peptide. This changes to 4.8 kcal mol(-1) per triplet of destabilization energy for the [(-Ala-Ala-Gly-)(2)](3) sequence, known to be disfavored in collagen. The present study proves that by using first principles methods for calculating stabilities of supramolecular complexes, such as collagen and beta-pleated sheets, one can obtain stability data in full agreement with experimental observations, which envisage the applicability of QM in molecular design.

Endomorphin-2 with a β-Turn Backbone Constraint Retains the Potent μ-Opioid Receptor Agonist Properties

The constitutional similarity with different secondary structure preference between the Aba-Gly and the spiro-Aba-Gly scaffolds were exploited to design the novel endomorphin-2 analogs Tyr-spiro-( R/ S)-Aba-Gly-Phe-NH(2) ( 1 and 2) and Tyr-( R/ S)-Aba-Gly-Phe-NH(2) ( 3 and 4). The ( R)-spiro analog 1 was found to be a potent and selective micro-opioid agonist/partial agonist ( K (imicro) = 29.3 nM, IC(50) = 50 nM, K(e) = 0.57). NMR experiments and molecular modeling indicated that its backbone adopts mainly a beta-turn in aqueous solution.

Local interactions in protein folding determined through an inverse folding model

We adopt a model of inverse folding in which folding stability results from the combination of the hydrophobic effect with local interactions responsible for secondary structure preferences. Site-specific amino acid distributions can be calculated analytically for this model. We determine optimal parameters for the local interactions by fitting the complete inverse folding model to the site-specific amino acid distributions found in the Protein Data Bank. This procedure reduces drastically the influence on the derived parameters of the preference of different secondary structures for buriedness, which affects local interaction parameters determined through the standard approach based on amino acid propensities. The quality of the fit is evaluated through the likelihood of the observed amino acid distributions given the model and the Bayesian Information Criterion, which indicate that the model with optimal local interaction parameters is strongly preferable to the model where local interaction parameters are determined through propensities. The optimal model yields a mean correlation coefficient r = 0.96 between observed and predicted amino acid distributions. The local interaction parameters are then tested in threading experiments, in combination with contact interactions, for their capacity to recognize the native structure and structures similar to the native against unrelated ones. In a challenging test, proteins structurally aligned with the Mammoth algorithm are scored with the effective free energy function. The native structure gets the highest stability score in 100% of the cases, a high recognition rate comparable to that achieved against easier decoys generated by gapless threading. We then examine proteins for which at least one highly similar template exists. In 61% of the cases, the structure with the highest stability score excluding the native belongs to the native fold, compared to 60% if we use local interaction parameters derived from the usual amino acid propensities and 52% if we use only contact interactions. A highly similar structure is present within the five best stability scores in 82%, 81%, and 76% of the cases, for local interactions determined through inverse folding, through propensity, and set to zero, respectively. These results indicate that local interactions improve substantially the performances of contact free energy functions in fold recognition, and that similar structures tend to get high stability scores, although they are often not high enough to discriminate them from unrelated structures. This work highlights the importance to apply more challenging tests, as the recognition of homologous structures, for testing stability scores for protein folding.

Investigations on unconventional hydrogen bonds in RNA binding proteins: The role of CH⋯O [dbnd]C interactions

We have investigated the roles played by C H⋯O C interactions in RNA binding proteins. There was an average of 78 C H⋯O C interactions per protein and also there was an average of one significant C H⋯O C interaction for every 6 residues in the 59 RNA binding proteins studied. Main chain–Main chain (MM) C H⋯O C interactions are the predominant type of interactions in RNA binding proteins. The donor atom contribution to C H⋯O C interactions was mainly from aliphatic residues. The acceptor atom contribution for MM C H⋯O C interactions was mainly from Val, Phe, Leu, Ile, Arg and Ala. The secondary structure preference analysis of C H⋯O C interacting residues showed that, Arg, Gln, Glu and Tyr preferred to be in helix, while Ala, Asp, Cys, Gly, Ile, Leu, Lys, Met, Phe, Trp and Val preferred to be in strand conformation. Most of the C H⋯O C interacting polar amino acid residues were solvent exposed while, majority of the C H⋯O C interacting non polar residues were excluded from the solvent. Long and medium-range C H⋯O C interactions are the predominant type of interactions in RNA binding proteins. More than 50% of C H⋯O C interacting residues had a higher conservation score. Significant percentage of C H⋯O C interacting residues had one or more stabilization centers. Sixty-six percent of the theoretically predicted stabilizing residues were also involved in C H⋯O C interactions and hence these residues may also contribute additional stability to RNA binding proteins.

Characterization of interfacial solvent in protein complexes and contribution of wet spots to the interface description

Water networks in protein interfaces can complement direct interactions contributing significantly to molecular recognition, function, and stability of protein association. Thus, water can be seen as an extension or addition of protein structural features, which may add plenty of information to protein interfacial definition. However, solvent is frequently neglected in protein interaction studies. Analysis of the interfacial information contained in the PDB is essential to achieve more accurate descriptions of protein interfaces. With this aim, we have used the SCOWLP database (http://www.scowlp.org) and applied computational geometry methods to extract and analyze interfacial information of a high-resolution nonredundant dataset of 176 protein complexes containing obligate and transient interfaces. We have identified all interfacial residues and characterized them in terms of temperature factors, secondary structure, residue composition, and pairing preferences to understand their contribution to the interface description. We have paid special attention to water-bridged residues; focusing on those that interact only mediated by a water molecule called wet spots. Our results show that 40.1% of the interfacial residues are interacting through water and that wet spots represent a 14.5% of the total, emphasizing the importance of the inclusion of solvent in protein interaction studies, and the contribution of wet spots to interfacial description. Wet spots present similar characteristics to residues binding buried water molecules in the core or cavities of proteins; being preferably located in nonregular secondary structures and establishing hydrogen bonds by their main-chains. We observe that obligate and transient interfaces present a comparable amount of solvent. Moreover, the role of solvent in both complex types differs according to the different nature of their interfaces. The information obtained in our studies will assist in the process of accomplishing more accurate descriptions of protein interfaces and may be helpful to improve comparison of protein family interfaces, to facilitate rational ligand design, and to guide protein docking.

Investigations on C–H⋯π interactions in RNA binding proteins

We have investigated the roles played by C–H⋯π interactions in RNA binding proteins. There was an average of 55 C–H⋯π interactions per protein and also there was an average of one significant C–H⋯π interaction for every nine residues in the 59 RNA binding proteins studied. Main-chain to side-chain C–H⋯π interactions is the predominant type of interactions in RNA binding proteins. The donor atom contribution to C–H⋯π interactions was mainly from Phe, Tyr, Trp, Pro, Gly, Lys, His and Ala residues. The acceptor atom contribution to main-chain to side-chain C–H⋯π and side-chain to side-chain C–H⋯π interactions was mainly from Phe and Tyr residues. On the contrary, the acceptor atoms of Trp residues contributed to all the four types of C–H⋯π interactions. Also, Trp contributed both donor and acceptor atoms in main-chain to side-chain, main-chain to side-chain five-member aromatic ring and side-chain to side-chain C–H⋯π interactions. The secondary structure preference analysis of C–H⋯π interacting residues showed that, Arg, Gln, Glu, His, Ile, Leu, Lys, Met, Phe and Tyr preferred to be in helix, while Ala, Asp, Cys, Gly, Trp and Val preferred to be in strand conformation. Long-range C–H⋯π interactions are the predominant type of interactions in RNA binding proteins. More than 50% of C–H⋯π interacting residues had a higher conservation score. Significant percentage of C–H⋯π interacting residues had one or more stabilization centers. Seven percent of the theoretically predicted stabilizing residues were also involved in C–H⋯π interactions and hence these residues may also contribute additional stability to RNA binding proteins.

The preferences of orientations between the pairs of amino acids

In this work, we make an investigation on the preferences of orientations between amino acids using the orientation defined based on the local geometry of the amino acids concerned. It is found that there are common preferences of orientations (70°, 30°, 140°) and (110°, 340°, 100°) for various pairs of amino acids. Different side chains may strengthen or weaken the common preferences, which is related to the effect of packing. Some amino acids having specific local flexibility may possess some preferences of orientations besides the common ones, such as (10°, 280°, 210°). Another analysis on the pairs of the amino acids with different secondary-structure preferences shows that the directional interaction may affect the distribution of orientation more effectively than the packing or local flexibility. All these results provide us some insight of the organization of amino acids in protein, and their relation with some related interactions.

Tryptophan rich peptides: Influence of indole rings on backbone conformation

Synthetic peptides with defined secondary structure scaffolds, namely hairpins and helices, containing tryptophan residues, have been investigated in this study to probe the influence of a large number of aromatic amino acids on backbone conformations. Solution NMR investigations of Boc-W-L-W-(D)P-G-W-L-W-OMe (peptide 1), designed to form a well-folded hairpin, clearly indicates the influence of flanking aromatic residues at the (D)Pro-Gly region on both turn nucleation and strand propagation. Indole-pyrrolidine interactions in this peptide lead to the formation of the less-frequent type I' turn at the (D)Pro-Gly segment and frayed strand regions, with the strand residues adopting local helical conformations. An analog of peptide 1 with an Aib-Gly turn-nucleated hairpin (Boc-W-L-W-U-G-W-L-W-OMe (peptide 2)) shows a preference for helical structures in solution, in both chloroform and methanol. Peptides with either one (Boc-W-L-W-U-W-L-W-OMe (peptide 3)) or two (Boc-U-W-L-W-U-W-L-W-OMe (peptide 4)) helix-nucleating Aib residues give rise to the well-folded helical conformations in the chloroform solution. The results are indicative of a preference for helical folding in peptides containing a large number of Trp residues. Investigation of a tetrapeptide analog of peptide 2, Boc-W-U-G-W-OMe (peptide 5), in solution and in the crystal state (by X-ray diffraction), also indicates a preference for a helical fold. Additionally, peptide 5 is stabilized in crystals by both aromatic interactions and an array of weak interactions. Examination of Trp-rich sequences in protein structures, however, reveals no secondary structure preference, suggesting that other stabilizing interactions in a well-folded protein may offset the influence of indole rings on backbone conformations.

NMR characterizations of an amyloidogenic conformational ensemble of the PI3K SH3 domain

Amyloid formation is associated with structural changes of native polypeptides to monomeric intermediate states and their self-assembly into insoluble aggregates. Characterizations of the amyloidogenic intermediate state are, therefore, of great importance in understanding the early stage of amyloidogenesis. Here, we present NMR investigations of the structural and dynamic properties of the acid-unfolded amyloidogenic intermediate state of the phosphatidylinositol 3-kinase (PI3K) SH3 domain--a model peptide. The monomeric amyloidogenic state of the SH3 domain studied at pH 2.0 (35 degrees C) was shown to be substantially disordered with no secondary structural preferences. (15)N NMR relaxation experiments indicated that the unfolded polypeptide is highly flexible on a subnanosecond timescale when observed under the amyloidogenic condition (pH 2.0, 35 degrees C). However, more restricted motions were detected in residues located primarily in the beta-strands as well as in a loop in the native fold. In addition, nonnative long-range interactions were observed between the residues with the reduced flexibility by paramagnetic relaxation enhancement (PRE) experiments. These indicate that the acid-unfolded state of the SH3 domain adopts a partly folded conformation through nonnative long-range contacts between the dynamically restricted residues at the amyloid-forming condition.

Structural analysis and prediction of protein mutant stability using distance and torsion potentials: Role of secondary structure and solvent accessibility

Analyzing the factors behind protein stability is a key research topic in molecular biology, and has direct implications on protein structure prediction and protein-protein interactions. We have analyzed protein stability upon point mutations using a distance-dependant pair potential representing mainly through-space interactions, and torsion angle potential representing mainly neighboring effects as a basic statistical mechanical setup for the analysis. The synergetic effect of accessible surface area and secondary structure preferences was used as a classifier for the potentials. In addition, short-, medium-, and long-range interactions of the protein environment were also analyzed. Two datasets of point mutations were taken for the comparison of theoretically predicted stabilizing energy values with experimental DeltaDeltaG and DeltaDeltaGH(2)O from thermal and chemical denaturation experiments. These include 1538 and 1603 mutations, respectively, and contain 101 proteins that share a wide range of sequence identity. The resulting force fields were carefully evaluated with different statistical tests. Results show a maximum correlation of 0.87 with a standard error of 0.71 kcal/mol between predicted and measured DeltaDeltaG values and a prediction accuracy of 85.3% (stabilizing or destabilizing) for all mutations together. A correlation of 0.77 (more than 80% prediction accuracy with a standard error of 0.95 kcal/mol) each for the test dataset of split-sample validation and fivefold crossvalidation was obtained and a correlation of 0.70 (77.4% prediction accuracy with a standard error of 1.17 kcal/mol) was shown by the jackknife test. The same model was implemented, and the results were analyzed for mutations with DeltaDeltaGH(2)O. A correlation of 0.78 (standard error 0.96 kcal/mol) was observed with a prediction efficiency of 84.65%. This model can be used for the future prediction of protein structural stability together with various experimental techniques.

Nearest‐neighbor effects and structural preferences in dipeptides are a function of the electronic properties of amino acid side‐chains

The electronic properties of amino acid side-chains are emerging as an important factor in the preference for secondary structure in proteins. These properties have not been fully characterized, nor has their role in the behavior of peptides been explored in any detail. The present studies sought to evaluate several possibilities: 1) that hydrophilicity can be expressed solely in electronic terms, 2) that substituent effects of side-chains extend across the peptide bond, and (3) nearest-neighbor effects in dipeptides correlate with secondary structural preferences. Quantum mechanics (QM) calculations were used to define the electronic properties of individual amino acids and dipeptides. It was found that the hydrophilicity of an amino acid side-chain can be accurately represented as a function of the electron densities of its component atoms. In addition, the nature of an amino acid in the second position of a dipeptide affects the electronic properties (Mulliken populations and electron densities) of the main-chain atoms of the first residue. Certain electronic features of the dipeptides strongly correlated with propensity for secondary structure. Specifically, Mulliken population data at the Calpha atom and N atom predicted preference for alpha-helices versus coil and strand conformations, respectively. Analysis of dipeptides arrayed in either helical or extended structures revealed lengthening of main-chain bonds in the alpha-helical conformations. A thorough characterization of the electronic properties of amino acids and short peptide segments may provide a better understanding of the forces that determine secondary structure in proteins.

Disulfide connectivity prediction using secondary structure information and diresidue frequencies

We describe a stand-alone algorithm to predict disulfide bond partners in a protein given only the amino acid sequence, using a novel neural network architecture (the diresidue neural network), and given input of symmetric flanking regions of N-terminus and C-terminus half-cystines augmented with residue secondary structure (helix, coil, sheet) as well as evolutionary information. The approach is motivated by the observation of a bias in the secondary structure preferences of free cysteines and half-cystines, and by promising preliminary results we obtained using diresidue position-specific scoring matrices. As calibrated by receiver operating characteristic curves from 4-fold cross-validation, our conditioning on secondary structure allows our novel diresidue neural network to perform as well as, and in some cases better than, the current state-of-the-art method. A slight drop in performance is seen when secondary structure is predicted rather than being derived from three-dimensional protein structures.

The structure of antimicrobial pexiganan peptide in solution probed by Fourier transform infrared absorption, vibrational circular dichroism, and electronic circular dichroism spectroscopy

Pexiganan (Gly-Ile-Gly-Lys-Phe-Leu-Lys-Lys-Ala-Lys-Lys-Phe-Gly-Lys-Ala-Phe-Val-Lys-Ile-Leu-Lys-Lys), a 22 amino acid peptide, is an analogue of the magainin family of antimicrobial peptides present in the skin of the African clawed frog. Conformational analysis of pexiganan was carried out in different solvent environments for the first time. Organic solvents, trifluoroethanol (TFE) and methanol, were used to study the secondary structural preferences of this peptide in the membrane-mimicking environments. In addition, aqueous (D2O) and dimethyl sulfoxide (DMSO) solutions were also investigated to study the role of hydrogen bonding involved in the secondary structure formation. Fourier transform infrared absorption, vibrational circular dichroism (VCD), and electronic circular dichroism (ECD) measurements were carried out under the same conditions to ascertain the conformational assignments in different solvents. All these spectroscopic measurements suggest that the pexiganan peptide has the tendency to adopt different structures in different environments. Pexiganan appears to adopt an alpha-helical conformation in TFE, a sheet-stabilized beta-turn structure in methanol, a random coil with beta-turn structure in D2O, and a solvated beta-turn structure in DMSO.

AMBER-based hybrid force field for conformational sampling of polypeptides

AMBER-based hybrid force field was developed, with mixing the mainchain torsion energies of AMBER parm94 and parm96 force fields, with introducing a weighting factor to control the mixing ratio. First, dependence of an adiabatic Φ– Ψ potential–energy map of alanine dipeptide on the weighting factor was investigated, and the map was compared with an ab initio Φ– Ψ map. Second, we did enhanced conformational sampling of two six-residue peptides known to form a helical or turn conformation in solution. The hybrid force field reproduced the ab initio energy better than either parm94 or parm96, and exhibited the secondary-structure preference depending on the amino acid sequence.

Secondary-structure preferences of force fields for proteins evaluated by generalized-ensemble simulations

Secondary-structure forming tendencies are examined for six well-known protein force fields: AMBER94, AMBER96, AMBER99, CHARMM22, OPLS-AA/L, and GROMOS96. We performed generalized-ensemble molecular dynamics simulations of two peptides. One of these peptides is C-peptide of ribonuclease A, and the other is the C-terminal fragment from the B1 domain of streptococcal protein G. The former is known to form α-helix structure and the latter β-hairpin structure by experiments. The simulation results revealed significant differences of the secondary-structure forming tendencies among the force fields. Of the six force fields, the results of AMBER99 and CHARMM22 were in accord with experiments for C-peptide. For G-peptide, on the other hand, the results of OPLS-AA/L and GROMOS96 were most consistent with experiments. Therefore, further improvements on the force fields are necessary for studying the protein folding problem from the first principles, in which a single force field can be used for all cases.