Expanded turn conformations: characterization and sequence-structure correspondence in alpha-turns with implications in helix folding

Like the beta-turns, which are characterized by a limiting distance between residues two positions apart (i, i+3), a distance criterion (involving residues at positions i and i+4) is used here to identify alpha-turns from a database of known protein structures. At least 15 classes of alpha-turns have been enumerated based on the location in the phi,psi space of the three central residues (i+1 to i+3)-one of the major being the class AAA, where the residues occupy the conventional helical backbone torsion angles. However, moving towards the C-terminal end of the turn, there is a shift in the phi,psi angles towards more negative phi, such that the electrostatic repulsion between two consecutive carbonyl oxygen atoms is reduced. Except for the last position (i+4), there is not much similarity in residue composition at different positions of hydrogen and non-hydrogen bonded AAA turns. The presence or absence of Pro at i+1 position of alpha- and beta-turns has a bearing on whether the turn is hydrogen-bonded or without a hydrogen bond. In the tertiary structure, alpha-turns are more likely to be found in beta-hairpin loops. The residue composition at the beginning of the hydrogen bonded AAA alpha-turn has similarity with type I beta-turn and N-terminal positions of helices, but the last position matches with the C-terminal capping position of helices, suggesting that the existence of a "helix cap signal" at i+4 position prevents alpha-turns from growing into helices. Our results also provide new insights into alpha-helix nucleation and folding.     

Side-chain entropy effects on protein secondary structure formation

Loss of conformational entropy is one of the primary factors opposing protein folding. Both the backbone and side-chain of each residue in a protein will have their freedom of motion restricted in the final folded structure. The type of secondary structure of which a residue is part will have a significant impact on how much side-chain entropy is lost. Side-chain conformational entropies have previously been determined for folded proteins, simple models of unfolded proteins, alpha-helices, and a dipeptide model for beta-strands, but not for polyproline II (PII) helices. In this work, we present side-chain conformational estimates for the three regular secondary structure types: alpha-helices, beta-strands, and PII helices. Entropies are estimated from Monte Carlo computer simulations. Beta-strands are modeled as two structures, parallel and antiparallel beta-strands. Our data indicate that restraining a residue to the PII helix or antiparallel beta-strand conformations results in side-chain entropies equal to or higher than those obtained by restraining residues to the parallel beta-strand conformation. Side-chains in the alpha-helix conformation have the lowest side-chain entropies. The observation that extended structures retain the most side-chain entropy suggests that such structures would be entropically favored in unfolded proteins under folding conditions. Our data indicate that the PII helix conformation would be somewhat favored over beta-strand conformations, with antiparallel beta-strand favored over parallel. Notably, our data imply that, under some circumstances, residues may gain side-chain entropy upon folding. Implications of our findings for protein folding and unfolded states are discussed.     

Dynamic alpha-helices: conformations that do not conform

Structural transitions are important for the stability and function of proteins, but these phenomena are poorly understood. An extensive analysis of Protein Data Bank entries reveals 103 regions in proteins with a tendency to transform from helical to nonhelical conformation and vice versa. We find that these dynamic helices, unlike other helices, are depleted in hydrophobic residues. Furthermore, the dynamic helices have higher surface accessibility and conformational mobility (P-value = 3.35e-07) than the rigid helices. Contact analyses show that these transitions result from protein-ligand, protein-nucleic acid, and crystal-contacts. The immediate structural environment differs quantitatively (P-value = 0.003) as well as qualitatively in the two alternate conformations. Often, dynamic helix experiences more contacts in its helical conformation than in the nonhelical counterpart (P-value = 0.001). There is differential preference for the type of short contacts observed in two conformational states. We also demonstrate that the regions in protein that can undergo such large conformational transitions can be predicted with a reasonable accuracy using logistic regression model of supervised learning. Our findings have implications in understanding the molecular basis of structural transitions that are coupled with binding and are important for the function and stability of the protein. Based on our observations, we propose that several functionally relevant regions on the protein surface can switch over their conformation from coil to helix and vice-versa, to regulate the recognition and binding of their partner and hence these may work as "molecular switches" in the proteins to regulate certain biological process. Our results supports the idea that protein structure-function paradigm should transform from static to a highly dynamic one.     

A detailed unfolding pathway of a (beta/alpha)8-barrel protein as studied by molecular dynamics simulations

The (beta/alpha)(8)-barrel is the most common protein fold. Similar structural properties for folding intermediates of (beta/alpha)(8)-barrel proteins involved in tryptophan biosynthesis have been reported in a number of experimental studies; these intermediates have the last two beta-strands and three alpha-helices partially unfolded, with other regions of the polypeptide chain native-like in conformation. To investigate the detailed folding/unfolding pathways of these (beta/alpha)(8)-barrel proteins, temperature-induced unfolding simulations of N-(5'-phosphoribosyl)anthranilate isomerase from Escherichia coli were carried out using a special-purpose parallel computer system. Unfolding simulations at five different temperatures showed a sequential unfolding pathway comprised of several events. Early events in unfolding involved disruption of the last two strands and three helices, producing an intermediate ensemble similar to those detected in experimental studies. Then, denaturation of the first two betaalpha units and separation of the sixth strand from the fifth took place independently. The remaining central betaalphabetaalphabeta module persisted the longest during all simulations, suggesting an important role for this module as the incipient folding scaffold. Our simulations also predicted the presence of a nucleation site, onto which several hydrophobic residues condensed forming the foundation for the central betaalphabetaalphabeta module.     

Noncharged amino acid residues at the solvent-exposed positions in the middle and at the C terminus of the alpha-helix have the same helical propensity

It was established previously that helical propensities of different amino acid residues in the middle of α‐helix in peptides and in proteins are very similar. The statistical analysis of the protein helices from the known three‐dimensional structures shows no difference in the frequency of noncharged residues in the middle and at the C terminus. Yet, experimental studies show distinctive differences for the helical propensities of noncharged residues in the middle and in the C terminus in model peptides. Is this a general effect, and is it applicable to protein helices or is it specific to the model alanine‐based peptides? To answer this question, the effects of substitutions at positions 28 (middle residue) and 32 (C2 position at the C terminus) of the α‐helix of ubiquitin on the stability of this protein are measured by using differential scanning calorimetry. The two data sets produce similar values for intrinsic helix propensity, leading to a conclusion that noncharged amino acid residues at the solvent‐exposed positions in the middle and at the C terminus of the α‐helix have the same helical propensity. This conclusion is further supported with an excellent correlation between the helix propensity scale obtained for the two positions in ubiquitin with the experimental helix propensity scale established previously and with the statistical distribution of the residues in protein helices.     

Exploring the protein G helix free-energy surface by solute tempering metadynamics

The free-energy landscape of the alpha-helix of protein G is studied by means of metadynamics coupled with a solute tempering algorithm. Metadynamics allows to overcome large energy barriers, whereas solute tempering improves the sampling with an affordable computational effort. From the sampled free-energy surface we are able to reproduce a number of experimental observations, such as the fact that the lowest minimum corresponds to a globular conformation displaying some degree of beta-structure, that the helical state is metastable and involves only 65% of the chain. The calculations also show that the system populates consistently a pi-helix state and that the hydrophobic staple motif is present only in the free-energy minimum associated with the helices, and contributes to their stabilization. The use of metadynamics coupled with solute tempering results then particularly suitable to provide the thermodynamics of a short peptide, and its computational efficiency is promising to deal with larger proteins.     

Designed molecules that fold to mimic protein secondary structures

Molecules that fold to mimic protein secondary structures have emerged as important targets of bioorganic chemistry. Recently, a variety of compounds that mimic helices, turns, and sheets have been developed, with notable advances in the design of beta-peptides that mimic each of these structures. These compounds hold promise as a step toward synthetic molecules with protein-like properties and as drugs that block protein-protein interactions.     

Cooperative effects in hydrogen-bonding of protein secondary structure elements: a systematic analysis of crystal data using Secbase

A systematic analysis of the hydrogen-bonding geometry in helices and beta sheets has been performed. The distances and angles between the backbone carbonyl O and amide N atoms were correlated considering more than 1500 protein chains in crystal structures determined to a resolution better than 1.5 A. They reveal statistically significant trends in the H-bond geometry across the different secondary structural elements. The analysis has been performed using Secbase, a modular extension of Relibase (Receptor Ligand Database) which integrates information about secondary structural elements assigned to individual protein structures with the various search facilities implemented into Relibase. A comparison of the mean hydrogen-bond distances in alpha helices and 3(10) helices of increasing length shows opposing trends. Whereas in alpha helices the mean H-bond distance shrinks with increasing helix length and turn number, the corresponding mean dimension in 3(10) helices expands in a comparable series. Comparing similarly the hydrogen-bond lengths in beta sheets there is no difference to be found between the mean H-bond length in antiparallel and parallel beta sheets along the strand direction. In contrast, an interesting systematic trend appears to be given for the hydrogen bonds perpendicular to the strands bridging across an extended sheet. With increasing number of accumulated strands, which results in a growing number of back-to-back piling hydrogen bonds across the strands, a slight decrease of the mean H-bond distance is apparent in parallel beta sheets whereas such trends are obviously not given in antiparallel beta sheets. This observation suggests that cooperative effects mutually polarizing spatially well-aligned hydrogen bonds are present either in alpha helices and parallel beta sheets whereas such influences seem to be lacking in 3(10) helices and antiparallel beta sheets.     

Novel approach for alpha-helical topology prediction in globular proteins: generation of interhelical restraints

The protein folding problem represents one of the most challenging problems in computational biology. Distance constraints and topology predictions can be highly useful for the folding problem in reducing the conformational space that must be searched by deterministic algorithms to find a protein structure of minimum conformational energy. We present a novel optimization framework for predicting topological contacts and generating interhelical distance restraints between hydrophobic residues in alpha-helical globular proteins. It should be emphasized that since the model does not make assumptions about the form of the helices, it is applicable to all alpha-helical proteins, including helices with kinks and irregular helices. This model aims at enhancing the ASTRO-FOLD protein folding approach of Klepeis and Floudas (Journal of Computational Chemistry 2003;24:191-208), which finds the structure of global minimum conformational energy via a constrained nonlinear optimization problem. The proposed topology prediction model was evaluated on 26 alpha-helical proteins ranging from 2 to 8 helices and 35 to 159 residues, and the best identified average interhelical distances corresponding to the predicted contacts fell below 11 A in all 26 of these systems. Given the positive results of applying the model to several protein systems, the importance of interhelical hydrophobic-to-hydrophobic contacts in determining the folding of alpha-helical globular proteins is highlighted.     

The structure of the ends of α-helices in globular proteins: effect of additional hydrogen bonds and implications for helix formation

We prepared a set of about 2000 α-helices from a relational database of high-resolution three-dimensional structures of globular proteins, and identified additional main chain i ← i+3 hydrogen bonds at the ends of the helices (i.e., where the hydrogen bonding potential is not fulfilled by canonical i ← i+4 hydrogen bonds). About one-third of α-helices have such additional hydrogen bonds at the N-terminus, and more than half do so at the C-terminus. Although many of these additional hydrogen bonds at the C-terminus are associated with Schellman loops, the majority are not. We compared the dihedral angles at the termini of α-helices having or lacking the additional hydrogen bonds. Significant differences were found, especially at the C-terminus, where the dihedral angles at positions C2 and C1 in the absence of additional hydrogen bonds deviate substantially from those occurring within the α-helix. Using a novel approach we show how the structure of the C-terminus of the α-helix can emerge from that of constituent overlapping α-turns and β-turns, which individually show a variation in dihedral angles at different positions. We have also considered the direction of propagation of the α-helix using this approach. If one assumes that helices start as a single α-turn and grow by successive addition of further α-turns, the paths for growth in the N → C and C → N directions differ in a way that suggests that extension in the C → N direction is favored.     

Variants of 3(10)-helices in proteins

An analysis of the shortest 3(10)-helices, containing three helical residues and two flanking capping residues that participate in two consecutive i + 3 --> i hydrogen bonds, shows that not all helices belong to the classic 3(10)-helix, where the three central residues adopt the right-handed helical conformation (alpha(R)). Three variants identified are: 3L10-helix with all residues in the left-handed helical region (alpha(L)), 3EL10-helix where the first residue is in the extended region followed by two residues in the alpha(L) conformation, and its mirror-image, the 3E'R10-helix. In the context of these helices, as well as the equivalent variants of alpha-helices, the length dependence of the handedness of secondary structures in protein structure is discussed. There are considerable differences in the amino acid preferences at different positions in the various types of 3(10)-helices. Each type of 3(10)-helix can be thought to be made up of an extension of a particular type of beta-turn (made up of residues i to i + 3) such that the (i + 3)th residue assumes the same conformation as the preceding residue. Distinct residue preferences at i and i + 3 positions seem to decide whether a particular stretch of four residues will be a beta-turn or a 3(10)-helix in the folded structure.     

Crystal-state 3D-structural characterization of novel 3(10)-helical peptides.

The crystal-state conformations of two octapeptides, pBrBz-(D-Iva)8-OtBu (8I) and Ac-[L-(alphaMe)Val]8-OH (8II), the heptapeptide Z-[L-(alphaMe)Val]7-OH (7), the hexapeptide Z-[L-(alphaMe)Leu]6-OtBu (6) and the tetrapeptide alkylamide Z-(Aib)2-L-Glu(OMe)-L-Ala-L-Lol (5) were assessed by x-ray diffraction analyses. Two independent molecules are observed in the asymmetric unit of each L-(alphaMe)Val homo-peptide. All four homo-peptides are folded in a regular 3(10)-helical structure (only the C-terminal H-bonded conformation of the D-Iva octapeptide is distorted to a type-I beta-turn). The hydroxyl groups of the C-terminal carboxyl moieties of the two L-(alphaMe)Val homo-peptides participate in an oxy-analogue of the type-III beta-turn conformation. While the two L-(alphaMe)Val 3(10)-helices are right-handed, the D-Iva and L-(alphaMe)Leu helices are left-handed. The tetrapeptide alkylamide is 3(10)-helical at the N-terminus, but it is mixed 3(10)/alpha-helical at the C-terminus.     

Steroidal constituents of rice (Rryza sativa) hulls with algicidal and herbicidal activity against blue-green algae and duckweed

Two new compounds, 14-methyl stigmast-9(11)-en-3alpha-ol-3beta-D-glucopyranoside (1) and cholest-11-en-3beta, 6beta, 7alpha, 22beta-tetraol-24-one-3beta-palmitoleate (2), along with the known compound beta-sitosteryl-3beta-D-glucopyranosyl-6'-linoleiate (3), were isolated from the methanolic extract of rice (Oryza sativa) hulls. The structures of the two new compounds were elucidated using one- and two-dimensional NMR in combination with IR, EI/MS, FAB/MS, HR-EI/MS and HR-FAB/MS. In bioassays with blue-green algae, Microcystis aeruginosa UTEX 2388 and duckweed, Lemna paucicostata Hegelm 381, the efficacy of bioactivity of the two new compounds linearly increased as the concentration increased from 0.3 to 300 IgM. Compared with momilactone A, compounds 1 and 2 showed similar and higher inhibitory activities against the growth of M. aeruginosa at a concentration of 300 microM. However, compound 2 was similar to momilactone A in inhibiting L. paucicostata growth at a concentration of 300 microM. As a result, compound 2 appears to have a strong potential for the environmentally friendly control of weed and algae that are harmful to water-logged rice.     

aV integrins and TGF-β-induced EMT: a circle of regulation

Transforming growth factor-β (TGF-β) has roles in embryonic development, the prevention of inappropriate inflammation and tumour suppression. However, TGF-β signalling also regulates pathological epithelial-to-mesenchymal transition (EMT), inducing or progressing a number of diseases ranging from inflammatory disorders, to fibrosis and cancer. However, TGF-β signalling does not proceed linearly but rather induces a complex network of cascades that mutually influence each other and cross-talk with other pathways to successfully induce EMT. Particularly, there is substantial evidence for cross-talk between αV integrins and TGF-β during EMT, and anti-integrin therapeutics are under development as treatments for TGF-β-related disorders. However, TGF-β's complex signalling network makes the development of therapeutics to block TGF-β-mediated pathology challenging. Moreover, despite our current understanding of integrins and TGF-β function during EMT, the precise mechanism of their role during physiological versus pathological EMT is not fully understood. This review focuses on the circle of regulation between αV integrin and TGF-β signalling during TGF-β induced EMT, which pose as a significant driver to many known TGF-β-mediated disorders.     

Local structure formation in simulations of two small proteins

Massively parallel all-atom, explicit solvent molecular dynamics simulations were used to explore the formation and existence of local structure in two small alpha-helical proteins, the villin headpiece and the helical fragment B of protein A. We report on the existence of transient helices and combinations of helices in the unfolded ensemble, and on the order of formation of helices, which appears to largely agree with previous experimental results. Transient local structure is observed even in the absence of overall native structure. We also calculate sets of residue-residue pairs that are statistically predictive of the formation of given local structures in our simulations.     

Hydrogen bonding is the prime determinant of carboxyl pKa values at the N-termini of alpha-helices

Experimentally determined mean pK(a) values of carboxyl residues located at the N-termini of alpha-helices are lower than their overall mean values. Here, we perform three types of analyses to account for this phenomenon. We estimate the magnitude of the helix macrodipole to determine its potential role in lowering carboxyl pK(a) values at the N-termini. No correlation between the magnitude of the macrodipole and the pK(a) values is observed. Using the pK(a) program propKa we compare the molecular surroundings of 18 N-termini carboxyl residues versus 233 protein carboxyl groups from a previously studied database. Although pK(a) lowering interactions at the N-termini are similar in nature to those encountered in other protein regions, pK(a) lowering backbone and side-chain hydrogen bonds appear in greater number at the N-termini. For both Asp and Glu, there are about 0.5 more hydrogen bonds per residue at the N-termini than in other protein regions, which can be used to explain their lower than average pK(a) values. Using a QM-based pK(a) prediction model, we investigate the chemical environment of the two lowest Asp and the two lowest Glu pK(a) values at the N-termini so as to quantify the effect of various pK(a) determinants. We show that local interactions suffice to account for the acidity of carboxyl residues at the N-termini. The effect of the helix dipole on carboxyl pK(a) values, if any, is marginal. Backbone amide hydrogen bonds constitute the single biggest contributor to the lowest carboxyl pK(a) values at the N-termini. Their estimated pK(a) lowering effects range from about 1.0 to 1.9 pK(a) units.     

Glycogen synthase kinase-3beta regulates DeltaNp63 gene transcription through the beta-catenin signaling pathway

Overexpression of DeltaNp63 has been observed in a number of human cancers, suggesting a role for DeltaNp63 in carcinogenesis. In the present study, we show that inhibition of glycogen synthase kinase-3beta (GSK-3beta) by lithium chloride (LiCl) elicited a stimulatory effect on DeltaNp63 promoter activity in HEK 293T cells. Exposure to LiCl induced DeltaNp63 promoter activation as well as DeltaNp63 protein expression in the cells. The effect of GSK-3beta on DeltaNp63 expression was further confirmed by the use of two highly specific GSK-3beta inhibitors, SB216763 and SB415286. Further study showed the presence of a putative beta-catenin responsive element (beta-catenin-RE) in the DeltaNp63 promoter region, and the stimulation of DeltaNp63 promoter activity by GSK-3beta inhibitor is markedly abolished by mutation or deletion of the putative beta-catenin-RE. Data are also presented to show that beta-catenin acts together with Lef-1 to influence DeltaNp63 promoter activity and protein expression.     

Structure and hydration of the amylopectin trisaccharide building blocks--Synthesis, NMR, and molecular dynamics

To gain insight into the molecular details and hydration of amylopectin, the five constituting trisaccharides have been chemically synthesized as their methyl alpha-glycosides. All five trisaccharides were subjected to 950 MHz NMR spectroscopy for complete assignment and nanosecond molecular dynamics trajectories were calculated to study the structure and dynamics of the trisaccharides in aqueous solution. Systematic analysis of the simulation data revealed several examples of bridging water molecules playing an important role in the stabilization of specific amylopectin conformations, which was also supported by the experimental NMR data such as interresidue NOE's and heteronuclear scalar couplings between nuclei from neighboring residues. Although alpha-maltotriose, alpha-iso-maltotriose, alpha-panose and alpha-isopanose are relatively well characterized structures, the study also includes one less characterized trisaccharide with the structure alphaGlcp(1-->4)alphaGlcp(1-->6)alphaGlcp. This trisaccharide, tentatively labelled alpha-forkose, is located at the branch point of amylopectin, forking the amylopectin into two strands that align into double-helical segments. The results show that the conformation of alpha-forkose takes a natural bend form which fits well into the structure of the double-helical segment of amylopectin. As the only trisaccharide in this study the structure of alpha-forkose is not significantly influenced by the hydration. In contrast, alpha-isopanose takes a restricted, but rather extended form due to an exceptionally strong localized water density. The two homo-linkage oligomers, alpha-maltotriose and alpha-iso-maltotriose, showed to be the most extended and the most flexible trimers, respectively, providing regular structure for crystalline domains and maximum linker flexibility for amorphous domains.