Estimation of the number of alpha-helical and beta-strand segments in proteins using circular dichroism spectroscopy

A simple approach to estimate the number of alpha-helical and beta-strand segments from protein circular dichroism spectra is described. The alpha-helix and beta-sheet conformations in globular protein structures, assigned by DSSP and STRIDE algorithms, were divided into regular and distorted fractions by considering a certain number of terminal residues in a given alpha-helix or beta-strand segment to be distorted. The resulting secondary structure fractions for 29 reference proteins were used in the analyses of circular dichroism spectra by the SELCON method. From the performance indices of the analyses, we determined that, on an average, four residues per alpha-helix and two residues per beta-strand may be considered distorted in proteins. The number of alpha-helical and beta-strand segments and their average length in a given protein were estimated from the fraction of distorted alpha-helix and beta-strand conformations determined from the analysis of circular dichroism spectra. The statistical test for the reference protein set shows the high reliability of such a classification of protein secondary structure. The method was used to analyze the circular dichroism spectra of four additional proteins and the predicted structural characteristics agree with the crystal structure data.     

Design and characterization of peptides with amphiphilic beta-strand structures

To extend our studies on peptides and proteins with amphiphilic secondary structures, a series of peptides designed to form amphiphilic beta-strand structures was designed, synthesized, and characterized by circular dichroism and infrared spectroscopy. Amphiphilic beta-strand conformations may be likely to appear in a variety of surface-active proteins, including apolipoprotein B and fibronectin. In a beta-strand conformation, the synthetic peptides will possess a hydrophobic face composed of valine side chains and a hydrophilic face composed of alternating acidic (glutamic acid) and basic (ornithine or lysine) residues. The peptides studied had a variety of chain lengths (5, 9, and 13 residues), and had the amino groups either free or protected with the trifluoroacetyl group. While the peptides did not possess a high potential for beta-sheet formation based on the Chou Fasman parameters, they possessed significant beta-sheet content, with up to 90% beta-sheet calculated for the 13-residue protected peptide. The driving force for beta-sheet formation is the potential amphiphilicity of this conformation. The beta-strand conformation of the 13-residue deprotected peptide was stable in 50% trifluoroethanol, 6 M guanidine hydrochloride, and octanol. The peptides are strongly self-associating in water, which would reduce the unfavorable contacts of the hydrophobic residues with water. It is clear that small peptides can be designed to form stable beta-strand conformations.     

Site-directed spin labeling of a bacterial chemoreceptor reveals a dynamic, loosely packed transmembrane domain

We used site-directed spin labeling and electron paramagnetic resonance spectroscopy to investigate dynamics and helical packing in the four-helix transmembrane domain of the homodimeric bacterial chemoreceptor Trg. We focused on the first transmembrane helix, TM1, particularly on the nine-residue sequence nearest the periplasm, because patterns of disulfide formation between introduced cysteines had identified that segment as the region of closest approach among neighboring transmembrane helices. Along this sequence, mobility and accessibility of the introduced spin label were characteristic of loosely packed or solvent-exposed side chains. This was also the case for eight additional positions around the circumference and along the length of TM1. For the continuous nine-residue sequence near the periplasm, mobility and accessibility varied only modestly as a function of position. We conclude that side chains of TM1 that face the interior of the four-helix domain interact with neighboring helices but dynamic movement results in loose packing. Compared to transmembrane segments of other membrane proteins reconstituted into lipid bilayers and characterized by site-directed spin labeling, TM1 of chemoreceptor Trg is the most dynamic and loosely packed. A dynamic, loosely packed chemoreceptor domain can account for many experimental observations about the transmembrane domains of chemoreceptors.     

A functional protein pore with a "retro" transmembrane domain.

Extended retro (reversed) peptide sequences have not previously been accommodated within functional proteins. Here, we show that the entire transmembrane portion of the beta-barrel of the pore-forming protein alpha-hemolysin can be formed by retrosequences comprising a total of 175 amino acid residues, 25 contributed by the central sequence of each subunit of the heptameric pore. The properties of wild-type and retro heptamers in planar bilayers are similar. The single-channel conductance of the retro pore is 15% less than that of the wild-type heptamer and its current-voltage relationship denotes close to ohmic behavior, while the wild-type pore is weakly rectifying. Both wild-type and retro pores are very weakly anion selective. These results and the examination of molecular models suggest that beta-barrels may be especially accepting of retro sequences compared to other protein folds. Indeed, the ability to form a retro domain could be diagnostic of a beta-barrel, explaining, for example, the activity of the retro forms of many membrane-permeabilizing peptides. By contrast with the wild-type subunits, monomeric retro subunits undergo premature assembly in the absence of membranes, most likely because the altered central sequence fails to interact with the remainder of the subunit, thereby initiating assembly. Despite this difficulty, a technique was devised for obtaining heteromeric pores containing both wild-type and retro subunits. Most probably as a consequence of unfavorable interstrand side-chain interactions, the heteromeric pores are less stable than either the wild-type or retro homoheptamers, as judged by the presence of subconductance states in single-channel recordings. Knowledge about the extraordinary plasticity of the transmembrane beta-barrel of alpha-hemolysin will be very useful in the de novo design of functional membrane proteins based on the beta-barrel motif.     

Evaluation of transmembrane helix predictions in 2014

Experimental structure determination continues to be challenging for membrane proteins. Computational prediction methods are therefore needed and widely used to supplement experimental data. Here, we re‐examined the state of the art in transmembrane helix prediction based on a nonredundant dataset with 190 high‐resolution structures. Analyzing 12 widely‐used and well‐known methods using a stringent performance measure, we largely confirmed the expected high level of performance. On the other hand, all methods performed worse for proteins that could not have been used for development. A few results stood out: First, all methods predicted proteins in eukaryotes better than those in bacteria. Second, methods worked less well for proteins with many transmembrane helices. Third, most methods correctly discriminated between soluble and transmembrane proteins. However, several older methods often mistook signal peptides for transmembrane helices. Some newer methods have overcome this shortcoming. In our hands, PolyPhobius and MEMSAT‐SVM outperformed other methods. Proteins 2015; 83:473–484. © 2014 Wiley Periodicals, Inc.     

TMSEG: Novel prediction of transmembrane helices

Transmembrane proteins (TMPs) are important drug targets because they are essential for signaling, regulation, and transport. Despite important breakthroughs, experimental structure determination remains challenging for TMPs. Various methods have bridged the gap by predicting transmembrane helices (TMHs), but room for improvement remains. Here, we present TMSEG, a novel method identifying TMPs and accurately predicting their TMHs and their topology. The method combines machine learning with empirical filters. Testing it on a non‐redundant dataset of 41 TMPs and 285 soluble proteins, and applying strict performance measures, TMSEG outperformed the state‐of‐the‐art in our hands. TMSEG correctly distinguished helical TMPs from other proteins with a sensitivity of 98 ± 2% and a false positive rate as low as 3 ± 1%. Individual TMHs were predicted with a precision of 87 ± 3% and recall of 84 ± 3%. Furthermore, in 63 ± 6% of helical TMPs the placement of all TMHs and their inside/outside topology was correctly predicted. There are two main features that distinguish TMSEG from other methods. First, the errors in finding all helical TMPs in an organism are significantly reduced. For example, in human this leads to 200 and 1600 fewer misclassifications compared to the second and third best method available, and 4400 fewer mistakes than by a simple hydrophobicity‐based method. Second, TMSEG provides an add‐on improvement for any existing method to benefit from. Proteins 2016; 84:1706–1716. © 2016 Wiley Periodicals, Inc.     

The turn sequence directs beta-strand alignment in designed beta-hairpins

A previous NMR investigation of model decapeptides with identical beta-strand sequences and different turn sequences demonstrated that, in these peptide systems, the turn residues played a more predominant role in defining the type of beta-hairpin adopted than cross-strand side-chain interactions. This result needed to be tested in longer beta-hairpin forming peptides, containing more potentially stabilizing cross-strand hydrogen bonds and side-chain interactions that might counterbalance the influence of the turn sequence. In that direction, we report here on the design and 1H NMR conformational study of three beta-hairpin forming pentadecapeptides. The design consists of adding two and three residues at the N- and C-termini, respectively, of the previously studied decapeptides. One of the designed pentadecapeptides includes a potentially stabilizing R-E salt bridge to investigate the influence of this interaction on beta-hairpin stability. We suggest that this peptide self-associates by forming intermolecular salt bridges. The other two pentadecapeptides behave as monomers. A conformational analysis of their 1H NMR spectra reveals that they adopt different types of beta-hairpin structure despite having identical strand sequences. Hence, the beta-turn sequence drives beta-hairpin formation in the investigated pentadecapeptides that adopt beta-hairpins that are longer than the average protein beta-hairpins. These results reinforce our previous suggestion concerning the key role played by the turn sequence in directing the kind of beta-hairpin formed by designed peptides.     

Reversible unfolding and refolding behavior of a monomeric aldolase from Staphylococcus aureus.

Thermal and GdmCl-induced unfolding transitions of aldolase from Staphylococcus aureus are reversible under a variety of solvent conditions. Analysis of the transitions reveals that no partially folded intermediates can be detected under equilibrium conditions. The stability of the enzyme is very low with a delta G0 value of -9 +/- 2 kJ/mol at 20 degrees C. The kinetics of unfolding and refolding of aldolase are complex and comprise at least one fast and two slow reactions. This complexity arises from prolyl isomerization reactions in the unfolded chain, which are kinetically coupled to the actual folding reaction. Comparison with model calculations shows that at least two prolyl peptide bonds give rise to the observed slow folding reactions of aldolase and that all of the involved bonds are presumably in the trans conformation in the native state. The rate constant of the actual folding reaction is fast with a relaxation time of about 15 s at the midpoint of the folding transition at 15 degrees C. The data presented on the folding and stability of aldolase are comparable to the properties of much smaller proteins. This might be connected with the simple and highly repetitive tertiary structure pattern of the enzyme, which belongs to the group of alpha/beta barrel proteins.     

Turn scanning by site-directed mutagenesis: application to the protein folding problem using the intestinal fatty acid binding protein

We have systematically mutated residues located in turns between beta-strands of the intestinal fatty acid binding protein (IFABP), and a glycine in a half turn, to valine and have examined the stability, refolding rate constants and ligand dissociation constants for each mutant protein. IFABP is an almost all beta-sheet protein exhibiting a topology comprised of two five-stranded sheets surrounding a large cavity into which the fatty acid ligand binds. A glycine residue is located in seven of the eight turns between the antiparallel beta-strands and another in a half turn of a strand connecting the front and back sheets. Mutations in any of the three turns connecting the last four C-terminal strands slow the folding and decrease stability with the mutation between the last two strands slowing folding dramatically. These data suggest that interactions between the last four C-terminal strands are highly cooperative, perhaps triggered by an initial hydrophobic collapse. We suggest that this trigger is collapse of the highly hydrophobic cluster of amino acids in the D and E strands, a region previously shown to also affect the last stage of the folding process (Kim et al., 1997). Changing the glycine in the strand between the front and back sheets also results in a unstable, slow folding protein perhaps disrupting the D-E strand interactions. For most of the other turn mutations there was no apparent correlation between stability and refolding rate constants. In some turns, the interaction between strands, rather than the turn type, appears to be critical for folding while in others, turn formation itself appears to be a rate limiting step. Although there is no simple correlation between turn formation and folding kinetics, we propose that turn scanning by mutagenesis will be a useful tool for issues related to protein folding.     

Discrimination of outer membrane proteins using machine learning algorithms

Discriminating outer membrane proteins (OMPs) from other folding types of globular and membrane proteins is an important task both for identifying OMPs from genomic sequences and for the successful prediction of their secondary and tertiary structures. In this work, we have analyzed the performance of different methods, based on Bayes rules, logistic functions, neural networks, support vector machines, decision trees, etc. for discriminating OMPs. We found that most of the machine learning techniques discriminate OMPs with similar accuracy. The neural network-based method could discriminate the OMPs from other proteins [globular/transmembrane helical (TMH)] at the fivefold cross-validation accuracy of 91.0% in a dataset of 1,088 proteins. The accuracy of discriminating globular proteins is 88.8% and that of TMH proteins is 93.7%. Further, the neural network method is tested with globular proteins belonging to 30 different folding types and it could successfully exclude 95% of the considered proteins. The proteins with SAM domain such as knottins, rubredoxin, and thioredoxin folds are eliminated with 100% accuracy. These accuracy levels are comparable to or better than other methods in the literature. We suggest that this method could be effectively used to discriminate OMPs and for detecting OMPs in genomic sequences.     

Analysis of interactive packing of secondary structural elements in alpha/beta units in proteins

An alpha-helix and a beta-strand are said to be interactively packed if at least one residue in each of the secondary structural elements loses 10% of its solvent accessible contact area on association with the other secondary structural element. An analysis of all such 5,975 nonidentical alpha/beta units in protein structures, defined at < or = 2.5 A resolution, shows that the interaxial distance between the alpha-helix and the beta-strand is linearly correlated with the residue-dependent function, log[(V/nda)/n-int], where V is the volume of amino acid residues in the packing interface, nda is the normalized difference in solvent accessible contact area of the residues in packed and unpacked secondary structural elements, and n-int is the number of residues in the packing interface. The beta-sheet unit (beta u), defined as a pair of adjacent parallel or antiparallel hydrogen-bonded beta-strands, packing with an alpha-helix shows a better correlation between the interaxial distance and log(V/nda) for the residues in the packing interface. This packing relationship is shown to be useful in the prediction of interaxial distances in alpha/beta units using the interacting residue information of equivalent alpha/beta units of homologous proteins. It is, therefore, of value in comparative modeling of protein structures.     

Stoichiometry of lipid interactions with transmembrane proteins--Deduced from the 3D structures

The stoichiometry of the first shell of lipids interacting with a transmembrane protein is defined operationally by the population of spin-labeled lipid chains whose motion is restricted directly by the protein. Interaction stoichiometries have been determined experimentally for a wide range of alpha-helical integral membrane proteins by using spin-label ESR spectroscopy. Here, we determine the spatially defined number of first-shell lipids at the hydrophobic perimeter of integral membrane proteins whose 3D structure has been determined by X-ray crystallography and lipid-protein interactions characterized by spin-labeling. Molecular modeling is used to build a single shell of lipids surrounding transmembrane structures derived from the PDB. Constrained energy optimization of the protein-lipid assemblies is performed by molecular mechanics. For relatively small proteins (up to 7-12 transmembrane helices), the geometrical first shell corresponds to that defined experimentally by perturbation of the lipid-chain dynamics. For larger, multi-subunit alpha-helical proteins, the lipids perturbed directly by the protein may either exceed or be less in number than those that can be accommodated at the intramembranous perimeter. In these latter cases, the motionally restricted spin-labeled lipids can be augmented by intercalation, or can correspond to a specific subpopulation at the protein interface, respectively. For monomeric beta-barrel proteins, the geometrical lipid stoichiometry corresponds to that determined from lipid mobility for a 22-stranded barrel, but fewer lipids are motionally restricted than can be accommodated around an eight-stranded barrel. Deviations from the geometrical first shell, in the beta-barrel case, are for the smaller protein with a highly curved barrel.     

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.     

Motional restrictions of membrane proteins: a site-directed spin labeling study

Site-directed mutagenesis was used to produce 27 single cysteine mutants of bacteriophage M13 major coat protein spanning the whole primary sequence of the protein. Single-cysteine mutants were labeled with nitroxide spin labels and incorporated into phospholipid bilayers with increasing acyl chain length. The SDSL is combined with ESR and CD spectroscopy. CD spectroscopy provided information about the overall protein conformation in different mismatching lipids. The spin label ESR spectra were analyzed in terms of a new spectral simulation approach based on hybrid evolutionary optimization and solution condensation. This method gives the residue-level free rotational space (i.e., the effective space within which the spin label can wobble) and the diffusion constant of the spin label attached to the protein. The results suggest that the coat protein has a large structural flexibility, which facilitates a stable protein-to-membrane association in lipid bilayers with various degrees of hydrophobic mismatch.     

Asymmetric amino acid compositions of transmembrane beta-strands

In contrast to water-soluble proteins, membrane proteins reside in a heterogeneous environment, and their surfaces must interact with both polar and apolar membrane regions. As a consequence, the composition of membrane proteins' residues varies substantially between the membrane core and the interfacial regions. The amino acid compositions of helical membrane proteins are also known to be different on the cytoplasmic and extracellular sides of the membrane. Here we report that in the 16 transmembrane beta-barrel structures, the amino acid compositions of lipid-facing residues are different near the N and C termini of the individual strands. Polar amino acids are more prevalent near the C termini than near the N termini, and hydrophobic amino acids show the opposite trend. We suggest that this difference arises because it is easier for polar atoms to escape from the apolar regions of the bilayer at the C terminus of a beta-strand. This new characteristic of beta-barrel membrane proteins enhances our understanding of how a sequence encodes a membrane protein structure and should prove useful in identifying and predicting the structures of trans-membrane beta-barrels.     

Free energy barrier estimation of unfolding the alpha-helical surfactant-associated polypeptide C

Molecular dynamics simulations were conducted to estimate the free energy barrier of unfolding surfactant-associated polypeptide C (SP-C) from an alpha-helical conformation. Experimental studies indicate that while the helical fold of SP-C is thermodynamically stable in phospholipid micelles, it is metastable in a mixed organic solvent of CHCl3/CH3OH/0.1 M HCl at 32:64:5 (v/v/v), in which it undergoes an irreversible transformation to an insoluble aggregate that contains beta-sheet. On the basis of experimental observations, the free energy barrier was estimated to be approximately 100 kJ/mole by applying Eyring's transition state theory to the experimental rate of unfolding [Protein Sci 1998;7:2533-2540]. These studies prompted us to carry out simulations to investigate the unwinding process of two helical turns encompassing residues 25-32 in water and in methanol. The results give an upper bound estimation for the free energy barrier of unfolding of SP-C of approximately 20 kJ/mole. The results suggest a need to reconsider the applicability of a single-mode activated process theory to protein unfolding.     

Sequence specificity, statistical potentials, and three-dimensional structure prediction with self-correcting distance geometry calculations of beta-sheet formation in proteins

A statistical analysis of a representative data set of 169 known protein structures was used to analyze the specificity of residue interactions between spatial neighboring strands in beta-sheets. Pairwise potentials were derived from the frequency of residue pairs in nearest contact, second nearest and third nearest contacts across neighboring beta-strands compared to the expected frequency of residue pairs in a random model. A pseudo-energy function based on these statistical pairwise potentials recognized native beta-sheets among possible alternative pairings. The native pairing was found within the three lowest energies in 73% of the cases in the training data set and in 63% of beta-sheets in a test data set of 67 proteins, which were not part of the training set. The energy function was also used to detect tripeptides, which occur frequently in beta-sheets of native proteins. The majority of native partners of tripeptides were distributed in a low energy range. Self-correcting distance geometry (SECODG) calculations using distance constraints sets derived from possible low energy pairing of beta-strands uniquely identified the native pairing of the beta-sheet in pancreatic trypsin inhibitor (BPTI). These results will be useful for predicting the structure of proteins from their amino acid sequence as well as for the design of proteins containing beta-sheets.     

A consensus algorithm to screen genomes for novel families of transmembrane beta barrel proteins

The ability to search sequence datasets for membrane spanning proteins is an important requirement for genome annotation. However, the development of algorithms to identify novel types of transmembrane beta-barrel (TMB) protein has proven substantially harder than for transmembrane helical proteins, owing to a shorter TM domain in which only alternate residues are hydrophobic. Although recent reports have described important improvements in the development of such algorithms, there is still concern over their ability to confidently screen genomes. Here we describe a new algorithm combining composition and hidden Markov model topology based classifiers (called TMB-Hunt2), which achieves a crossvalidation accuracy of >95%, with 96.7% precision and 94.2% recall. An overview is given of the algorithm design, with a thorough assessment of performance and application to a number of genomes. Of particular note is that TMB/extracellular protein discrimination is significantly more difficult than TMB/cytoplasmic protein discrimination, with the predictor correctly rejecting just 74% of extracellular proteins, in comparison to 98% of cytoplasmic proteins. Focus is given to directions for further improvements in TMB/non-TMB protein discrimination, with a call for the development of standardized tests and assessments of such algorithms. Tools and datasets are made available through a website called TMB-Web (http://www.bioinformatics.leeds.ac.uk/TMB-Web/TMB-Hunt2).     

Accommodation of insertion mutations on the surface and in the interior of staphylococcal nuclease.

Alignment of homologous amino acid sequences reveals that insertion mutations are fairly common in evolution. Hitherto, the structural consequences of insertion mutations on the surface and in the interior of proteins of known structures have received little attention. We report here the high-resolution X-ray crystal structures of 2 site-directed insertion mutants of staphylococcal nuclease. The structure of the first insertion mutant, in which 2 glycine residues were inserted on the protein surface in the amino-terminal beta-strand, has been solved to 1.70 A resolution and refined to a crystallographic R value of 0.182. The inserted residues are accommodated in a special 3-residue beta-bulge. A bridging water molecule in the newly created cavity satisfies the hydrogen bonding requirements of the beta-sheet by forming a bifurcated hydrogen bond to 1 beta-strand, and a single hydrogen bond to the other beta-strand. The second insertion mutant contains a single leucine residue inserted at the end of the third beta-strand. The structure was solved to 2.0 A resolution and refined to a final R value of 0.196. The insertion is accommodated in a register shift that changes the conformation of the flexible loop portion of the molecule, relaxing and widening the omega turn. This structural alteration results in changes in position and coordination of a bound calcium ion important for catalysis. These structures illustrate important differences in how amino acid insertions are accommodated: as localized bulges, and as extensive register shifts.     

Detection of drug‐induced conformational change of a transmembrane protein in lipid bilayers using site‐directed spin labeling

As a target of antiviral drugs, the influenza A M2 protein has been the focus of numerous structural studies and has been extensively explored as a model ion channel. In this study, we capitalize on the expanding body of high‐resolution structural data available for the M2 protein to design and interpret site‐directed spin‐labeling electron paramagnetic resonance spectroscopy experiments on drug‐induced conformational changes of the M2 protein embedded in lipid bilayers. We obtained data in the presence of adamantane drugs for two different M2 constructs (M2TM 22–46 and M2TMC 23–60). M2TM peptides were spin labeled at the N‐terminal end of the transmembrane domain. M2TMC peptides were spin labeled site specifically at cysteine residues substituted for amino acids within the transmembrane domain (L36, I39, I42, and L43) and the C‐terminal amphipathic helix (L46, F47, F48, C50, I51, Y52, R53, F54, F55, and E56). Addition of adamantane drugs brought about significant changes in measured electron paramagnetic resonance spectroscopy environmental parameters consistent with narrowing of the transmembrane channel pore and closer packing of the C‐terminal amphipathic helices.