An Automatic Search for Similar Spatial Arrangements of α-Helices and β-Strands in Globular Proteins

Abstract

A fast search algorithm to reveal similar polypeptide backbone structural motifs in proteins is proposed. It is based on the vector representation of a polypeptide chain fold in which the elements of regular secondary structures are approximated by linear segments (Abagyan and Maiorov, J. Biomol. Struct. Dyn. 5, 1267–1279 (1988)). The algorithm permits insertions and deletions in the polypeptide chain fragments to be compared. The fast search algorithm implemented in FASEAR program is used for collecting βαβ supersecondary structure units in a number of α/β proteins of Brookhaven Data Bank. Variation of geometrical parameters specifying backbone chain fold is estimated. It appears that the conformation of the majority of the fragments, although almost all of them are right-handed, is quite different from that of standard βαβ units. Apart from searching for specific type of secondary structure motif, the algorithm allows automatically to identify new recurrent folding patterns in proteins. It may be of particular interest for the development of tertiary template approach for prediction of protein three-dimensional structure as well for constructing artificial polypeptides with goal-oriented conformation.     

The betaalphabetaalphabeta elementary supersecondary structure of the Rossmann fold from porcine lactate dehydrogenase exhibits characteristics of a molten globule

Protein classifications show that the Rossmann fold, which consists of two betaalphabetaalphabeta motifs (BABAB) related by a rough twofold axis, is the most populated alphabeta fold, and that the betaalphabeta submotif (BAB) is a widespread elementary structural arrangement. Herein, we report MD simulations, circular dichroism and NMR analyses on BAB and BABAB from porcine lactate dehydrogenase to evaluate their intrinsic stability. Our results demonstrate that BAB is not stable in solution and is not a folding nucleus. We also find that BABAB, despite its appearance of a functional and structural unit, is not an independent and thermodynamically stable folding unit. Rather, we show that BABAB retains most native secondary structure but very little tertiary structure, thus displaying characteristics of a molten globule.     

Fast in silico protein folding by introduction of alternating hydrogen bond potentials

We accelerate protein folding in all-atom molecular dynamics simulations by introducing alternating hydrogen bond potentials as a supplement to the force field. The alternating hydrogen bond potentials result in accelerated hydrogen bond reordering, which leads to rapid formation of secondary structure elements. The method does not require knowledge of the native state but generates the potentials based on the development of the tertiary structure in the simulation. In protein folding, the formation of secondary structure elements, especially α-helix and β-sheet, is very important, and we show that our method can fold both efficiently and with great speed.     

Molecular Dynamics Simulations of Peptides from the Central Domain of Smooth Muscle Caldesmon

Abstract

The central domain of smooth muscle caldesmon contains a highly charged region consisting of ten 13-residue repeats. Experimental evidence obtained from the intact protein and fragments thereof suggests that this entire region forms a single stretch of stable α-helix. We have carried out molecular dynamics simulations on peptides consisting of one, two and three repeats to examine the mechanism of α-helical stability of the central domain at the atomic level. All three peptides show high helical stability on the timescale of the MD simulations. Deviations from α-helical structure in all the simulations arise mainly from the formation of long stretches of π-helix. Interconversion between α-helical and π-helical conformations occurs through insertion of water molecules into α-helical hydrogen bonds and subsequent formation of reverse turns. The α-helical structure is stabilized by electrostatic interactions (salt bridges) between oppositely charged sidechains with i,i+4 spacings, while the π-helix is stabilized by i,i+5 salt bridge interactions. Possible i,i+3 salt bridges are of minor importance. There is a strong preference for salt bridges with a Glu residue N-terminal to a basic sidechain as compared to the opposite orientation. In the double and triple repeat peptides, strong i,i+4 salt bridges exist between the last Glu residue of one repeat and the first Lys residue of the next. This demonstrates a relationship between the repetitive nature of the central domain sequence and its ability to form very long stretches of α-helical structure.     

The effects of codon usage on the formation of secondary structures of nucleocapsid protein of peste des petits ruminants virus

The nucleocapsid (N) protein of peste des petits ruminants virus (PPRV) with a conserved amino acid usage pattern plays an important role in viral replication. The primary objective of this study was to estimate roles of synonymous codon usages of PPRV N gene and tRNA abundances of host in the formation of secondary structure of N protein. The potential effects of synonymous codon usages of N gene and tRNA abundances of host on shaping different folding units (α-helix, β-strand and the coil) in N protein were estimated, based on the information about the modeling secondary structure of PPRV N protein. The synonymous codon usage bias was found in different folding units in PPRV N protein. To better understand the role of translation speed caused by variant tRNA abundances in shaping the specific folding unit in N protein, we modeled the changing trends of tRNA abundance at the transition boundaries from one folding unit to another folding unit (β-strand → coil, coil → β-strand, α-helix → coil, coil → α-helix). The obvious fluctuations of tRNA abundance were identified at the two transition boundaries (β-strand → coil and coil → β-strand) in PPRV N protein. Our findings suggested that viral synonymous codon usage bias and cellular tRNA abundance variation might have potential effects on the formation of secondary structure of PPRV N protein.     

Controlled unfolding and refolding of a single sodium-proton antiporter using atomic force microscopy

Single-molecule force-spectroscopy was employed to unfold and refold single sodium-proton antiporters (NhaA) of Escherichia coli from membrane patches. Although transmembrane alpha-helices and extracellular polypeptide loops exhibited sufficient stability to individually establish potential barriers against unfolding, two helices predominantly unfolded pairwise, thereby acting as one structural unit. Many of the potential barriers were detected unfolding NhaA either from the C-terminal or the N-terminal end. It was found that some molecular interactions stabilizing secondary structural elements were directional, while others were not. Additionally, some interactions appeared to occur between the secondary structural elements. After unfolding ten of the 12 helices, the extracted polypeptide was allowed to refold back into the membrane. After five seconds, the refolded polypeptide established all secondary structure elements of the native protein. One helical pair showed a characteristic spring like "snap in" into its folded conformation, while the refolding process of other helices was not detected in particular. Additionally, individual helices required characteristic periods of time to fold. Correlating these results with the primary structure of NhaA allowed us to obtain the first insights into how potential barriers establish and determine the folding kinetics of the secondary structure elements.     

Conformational changes associated with L16P and T118M mutations in the membrane-embedded PMP22 protein,consequential in CMT-1A

Peripheral myelin protein 22 (PMP22) resides in the plasma membrane and is required for myelin formation in the peripheral nervous system. Excess PMP22 mutants accumulate in the endoplasmic reticulum (ER) resulting in the inherited neuropathies of Charcot–Marie–Tooth disease. However, there was no evidence of the structure of PMP22 or how mutations affect its folding. Therefore, in this study, we combined bioinformatics and homology modeling approaches to obtain three-dimensional native and mutated PMP22 models and its anchoring to a POPC membrane, submitted to .5-μs MD simulations, to determine how the L16P and T118M mutations affect the conformational behavior of PMP22. In addition, we investigated the ability of the native and mutated species to accumulate in the ER, via interaction with RER1, by combining protein–protein docking and MD simulations, taking the conformations that were most representative of the native and mutated PMP22 systems and RER1 conformations. Principal component analysis over MD simulations revealed that L16P and T118M mutations resulted in increased structural instability compared to the native form, which is consistent with previous experimental findings of increased structural fluctuations along a loop connecting transmembrane α-helix1 and α-helix2. Docking and MD simulations coupled with the MMGBSA approach allowed the identification that the binding interface for the PMP22-RER1 complex takes place through transmembrane α-helix1 and α-helix2, with higher effective binding free energy values between the mutated PMP22 systems and RER1 than for the native PMP22, mainly through van der Waals interactions.     

Asymmetric stability among the transmembrane helices of lactose permease

Combining structure determinations from nuclear magnetic resonance (NMR) data and molecular dynamics simulations (MD) under the same environmental conditions revealed a startling asymmetry in the intrinsic conformational stability of secondary structure in the transmembrane domain of lactose permease (LacY). Eleven fragments, corresponding to transmembrane segments (TMs) of LacY, were synthesized, and their secondary structure in solution was determined by NMR. Eight of the TMs contained significant regions of helical structure. MD simulations, both in DMSO and in a DMPC bilayer, showed sites of local stability of helical structure in these TMs, punctuated by regions of conformational instability, in substantial agreement with the NMR data. Mapping the stable regions onto the crystal structure of LacY reveals a marked asymmetry, contrasting with the pseudosymmetry in the static structure: the secondary structure in the C-terminal half is more stable than in the N-terminal half. The relative stability of secondary structure is likely exploited in the transport mechanism of LacY. Residues supporting proton conduction are in more stable regions of secondary structure, while residues key to substrate binding are found in considerably unstable regions of secondary structure.     

Dynamics and cooperativity of Trp-cage folding

In this study, multiple independent molecular dynamics (MD) simulations on Trp-cage folding were performed at 300, 325 and 375 K using generalized Born (GB) implicit solvent model. The orientational movement of the side-chain of Trp6 to form a hydrophobic core with 3₁₀-helix was observed. The breaking/formation of a salt bridge between Asp9 and Arg16 was proposed to be the prerequisite for Trp-cage folding/refolding. Our results demonstrate that the cooperation between the salt bridge and the Trp6 orientation leads to a stable tertiary structure of Trp-cage. Analyses on backbone concerted motions at different temperatures indicate that interactions between Trp6 and 3₁₀-helix & Pro18 and between Pro12 and Pro17 & Pro18 are weakened at 375 K but strengthened at lower temperatures, suggesting that they could be the potential driving force of hydrophobic collapse.     

Structure–function relationships in ShKT domain peptides: ShKT-Ts1 from the sea anemone Telmatactis stephensoni

Diverse structural scaffolds have been described in peptides from sea anemones, with the ShKT domain being a common scaffold first identified in ShK toxin from Stichodactyla helianthus. ShK is a potent blocker of voltage-gated potassium channels (K_V1.x), and an analog, ShK-186 (dalazatide), has completed Phase 1 clinical trials in plaque psoriasis. The ShKT domain has been found in numerous other species, but only a tiny fraction of ShKT domains has been characterized functionally. Despite adopting the canonical ShK fold, some ShKT peptides from sea anemones inhibit K_V1.x, while others do not. Mutagenesis studies have shown that a Lys–Tyr (KY) dyad plays a key role in K_V1.x blockade, although a cationic residue followed by a hydrophobic residue may also suffice. Nevertheless, ShKT peptides displaying an ShK-like fold and containing a KY dyad do not necessarily block potassium channels, so additional criteria are needed to determine whether new ShKT peptides might show activity against potassium channels. In this study, we used a combination of NMR and molecular dynamics (MD) simulations to assess the potential activity of a new ShKT peptide. We determined the structure of ShKT-Ts1, from the sea anemone Telmatactis stephensoni, examined its tissue localization, and investigated its activity against a range of ion channels. As ShKT-Ts1 showed no activity against K_V1.x channels, we used MD simulations to investigate whether solvent exposure of the dyad residues may be informative in rationalizing and potentially predicting the ability of ShKT peptides to block K_V1.x channels. We show that either a buried dyad that does not become exposed during MD simulations, or a partially exposed dyad that becomes buried during MD simulations, correlates with weak or absent activity against K_V1.x channels. Therefore, structure determination coupled with MD simulations, may be used to predict whether new sequences belonging to the ShKT family may act as potassium channel blockers.     

The side chain interaction index as a tool for predicting fast-folding elements and the structure and stability of engineered peptides

The side chain interaction index (SCII) is a method of calculating the propensity for short-range interactions among side chains within a peptide sequence. Here, it is shown that the SCII values of secondary structure elements that have been shown to fold early and independently cluster separately from those of structures that fold later and/or are dependent on long-range interactions. In addition, the SCII values of engineered peptides that spontaneously adopt a particular desired fold in solution are significantly different from those of engineered peptides that fail to exhibit a stable conformation. Thus, the SCII, as a measure of local structural stability, constitutes a useful tool in folding prediction and in protein/peptide engineering. A program that allows rapid calculation of SCII values is presented.     

Folding study of an Aib-rich peptide in DMSO by molecular dynamics simulations.

To evaluate the ability of molecular dynamics (MD) simulations using atomic force-fields to correctly predict stable folded conformations of a peptide in solution, we show results from MD simulations of the reversible folding of an octapeptide rich in alpha-aminoisobutyric acid (2-amino-2-methyl-propanoic acid, Aib) solvated in di-methyl-sulfoxide (DMSO). This solvent generally prevents the formation of secondary structure, whereas Aib-rich peptides show a high propensity to form secondary structural elements, in particular 3(10)- and alpha-helical structures. Aib is, moreover, achiral, so that Aib-rich peptides can form left- or right-handed helices depending on the overall composition of the peptide, the temperature, and the solvation conditions. This makes the system an interesting case to study the ensembles of peptide conformations as a function of temperature by MD simulation. Simulations involving the folding and unfolding of the peptide were performed starting from two initial structures, a right-handed alpha-helical structure and an extended structure, at three temperatures, 298 K, 340 K, and 380 K, and the results are compared with experimental nuclear magnetic resonance (NMR) data measured at 298 K and 340 K. The simulations generally reproduce the available experimental nuclear Overhauser effect (NOE) data, even when a wide range of conformations is sampled at each temperature. The importance of adequate statistical sampling in order to reliably interpret the experimental data is discussed.     

Alpha-hairpin stability and folding of transmembrane segments

Molecular Dynamics (MD) simulations at low dielectric constant have been carried out for peptides matching the double spanning segments of transmembrane proteins. Different folding dynamics have been observed. The peptides folded into the stable helix-turn-helix conformation-alpha-hairpin-with antiparallel-oriented strands or unstable alpha-hairpin conformation that unfolded later into the straight helical structure. The peptide having flexible residues in the TM helices often misfolded into a tangled structure that can be avoided by restricting the flexibility of these residues. General conclusions can be drawn from the observed folding dynamics. The stability and folding of some double spanning transmembrane fragments are self-assembling. The following and/or neighboring peptide chains of the protein may support the stability of the hairpin structure of other fragments. The stability of the TM helices containing flexible residues could be maintained due to contacts with neighboring TM segments.     

Molecular dynamics as a tool to help design mutants

A 200 ps MD trajectory of wild type PCI and a 120 ps one for the Pro36Gly putative mutant are studied and compared with the structure of PCI in its complex with carboxypeptidase A (CPA). It is first established that the structures of PCI from X-ray and from MD simulation are essentially equal. Thereafter, data from the PCI-CPA and average MD structures together with available biochemical information are used to identify possible structural factors that may determine the inhibitory power of PCI. These structural determinants are used to analyze the mutant structure. The fold of the mutant protein shows a large degree of stability. The N-terminal tail in PCIm differs from the X-ray structure as it does in PCIw, while the mutant's C-terminal tail (which is the primary binding site with CPA) and residues 13–17 present deviations. Differences in fluctuation patterns exist between PCIm and PCIw in residues 2–4 (the N-terminal tail), 13–17, 22–23, 28–81 (the secondary contact site with CPA), and 37–38 (the C-terminal tail); the latter region is rigidified in PCIm. Results show that the MD method is able to sense long-range as well as local perturbative effects produced by amino-acid substitutions in flexible regions of this protein. The simulations suggest that the conformation of the C-terminal tail is less favorable for interaction with the target protein in the mutant than it is in the wild type protein. The Pro-36-Gly mutant is predicted to be a less potent inhibitor.Abbreviations CPA carboxypeptidase A - MD molecular dynamics - NIS non-inertial solvent - PCI potato carboxypeptidase A inhibitor - PCIm mutated inhibitor - PCIw wild inhibitor     

Beta-hairpin folding mechanism of a nine-residue peptide revealed from molecular dynamics simulations in explicit water

The beta-hairpin fold mechanism of a nine-residue peptide, which is modified from the beta-hairpin of alpha-amylase inhibitor tendamistat (residues 15-23), is studied through direct folding simulations in explicit water at native folding conditions. Three 300-nanosecond self-guided molecular dynamics (SGMD) simulations have revealed a series of beta-hairpin folding events. During these simulations, the peptide folds repeatedly into a major cluster of beta-hairpin structures, which agree well with nuclear magnetic resonance experimental observations. This major cluster is found to have the minimum conformational free energy among all sampled conformations. This peptide also folds into many other beta-hairpin structures, which represent some local free energy minimum states. In the unfolded state, the N-terminal residues of the peptide, Tyr-1, Gln-2, and Asn-3, have a confined conformational distribution. This confinement makes beta-hairpin the only energetically favored structure to fold. The unfolded state of this peptide is populated with conformations with non-native intrapeptide interactions. This peptide goes through fully hydrated conformations to eliminate non-native interactions before folding into a beta-hairpin. The folding of a beta-hairpin starts with side-chain interactions, which bring two strands together to form interstrand hydrogen bonds. The unfolding of the beta-hairpin is not simply the reverse of the folding process. Comparing unfolding simulations using MD and SGMD methods demonstrate that SGMD simulations can qualitatively reproduce the kinetics of the peptide system.     

Analyzing the Interplay Between Secondary and Tertiary Structure Predictions in Folding Simulations with a Genetic Algorithm

Three different strategies to tackle mispredictions from incorrect secondary structure prediction are analysed using 21 small proteins (22-121 amino acids; 1-6 secondary structure elements) with known three dimensional structures: (1) Testing accuracy of different secondary structure predictions and improving them by combinations, (2) correcting mispredictions exploiting protein folding simulations with a genetic algorithm and (3) applying and combining experimental data to refine predictions both for secondary structure and tertiary fold. We demonstrate that predictions from secondary structure prediction programs can be efficiently combined to reduce prediction errors from missed secondary structure elements. Further, up to two secondary structure elements (helices, strands) missed by secondary structure prediction were corrected by the genetic algorithm simulation. Finally, we show how input from experimental data is exploited to refine the predictions obtained.Electronic Supplementary Material available.     

Conformational mapping of the N-terminal segment of surfactant protein B in lipid using 13C-enhanced Fourier transform infrared spectroscopy.

Synthetic peptides based on the N-terminal domain of human surfactant protein B (SP-B1-25; 25 amino acid residues; NH2-FPIPLPYCWLCRALIKRIQAMIPKG) retain important lung activities of the full-length, 79-residue protein. Here, we used physical techniques to examine the secondary conformation of SP-B1-25 in aqueous, lipid and structure-promoting environments. Circular dichroism and conventional, 12C-Fourier transform infrared (FTIR) spectroscopy each indicated a predominate alpha-helical conformation for SP-B1-25 in phosphate-buffered saline, liposomes of 1-palmitoyl-2-oleoyl phosphatidylglycerol and the structure-promoting solvent hexafluoroisopropanol; FTIR spectra also showed significant beta- and random conformations for peptide in these three environments. In further experiments designed to map secondary structure to specific residues, isotope-enhanced FTIR spectroscopy was performed with 1-palmitoyl-2-oleoyl phosphatidylglycerol liposomes and a suite of SP-B1-25 peptides labeled with 13C-carbonyl groups at either single or multiple sites. Combining these 13C-enhanced FTIR results with energy minimizations and molecular simulations indicated the following model for SP-B1-25 in 1-palmitoyl-2-oleoyl phosphatidylglycerol: beta-sheet (residues 1-6), alpha-helix (residues 8-22) and random (residues 23-25) conformations. Analogous structural motifs are observed in the corresponding homologous N-terminal regions of several proteins that also share the 'saposin-like' (i.e. 5-helix bundle) folding pattern of full-length, human SP-B. In future studies, 13C-enhanced FTIR spectroscopy and energy minimizations may be of general use in defining backbone conformations at amino acid resolution, particularly for peptides or proteins in membrane environments.