Three conserved C-terminal residues of influenza fusion peptide alter its behavior at the membrane interface

The N-terminal fragment of the viral hemagglutinin HA2 subunit is termed a fusion peptide (HAfp). The 23-amino acid peptide (HAfp1-23) contains three C-terminal W21-Y22-G23 residues which are highly conserved among serotypes of influenza A and has been shown to form a tight helical hairpin very distinct from the boomerang structure of HAfp1-20. We studied the effect of peptide length on fusion properties, structural dynamics, and binding to the membrane interface. We developed a novel fusion visualization assay based on FLIM microscopy on giant unilamellar vesicles (GUV). By means of molecular dynamics simulations and spectroscopic measurements, we show that the presence of the three C-terminal W21-Y22-G23 residues promotes the hairpin formation, which orients perpendicularly to the membrane plane and induces more disorder in the surrounding lipids than the less structured HAfp1-20. Moreover, we report cholesterol-enriched domain formation induced exclusively by the longer fusion peptide.     

Exploring the free energy landscape of a model β‐hairpin peptide and its isoform

Secondary structural transitions from α‐helix to β‐sheet conformations are observed in several misfolding diseases including Alzheimer's and Parkinson's. Determining factors contributing favorably to the formation of each of these secondary structures is therefore essential to better understand these disease states. β‐hairpin peptides form basic components of anti‐parallel β‐sheets and are suitable model systems for characterizing the fundamental forces stabilizing β‐sheets in fibrillar structures. In this study, we explore the free energy landscape of the model β‐hairpin peptide GB1 and its E2 isoform that preferentially adopts α‐helical conformations at ambient conditions. Umbrella sampling simulations using all‐atom models and explicit solvent are performed over a large range of end‐to‐end distances. Our results show the strong preference of GB1 and the E2 isoform for β‐hairpin and α‐helical conformations, respectively, consistent with previous studies. We show that the unfolded states of GB1 are largely populated by misfolded β‐hairpin structures which differ from each other in the position of the β‐turn. We discuss the energetic factors contributing favorably to the formation of α‐helix and β‐hairpin conformations in these peptides and highlight the energetic role of hydrogen bonds and non‐bonded interactions. Proteins 2014; 82:2394–2402. © 2014 Wiley Periodicals, Inc.     

Physicochemical properties of thiol proteinase inhibitor isolated from goat pancreas

Thiol proteinase inhibitors are crucial to proper functioning of all living tissues consequent to their cathepsin regulatory and myriad important biologic properties. Equilibrium denaturation of dimeric goat pancreas thiol proteinase inhibitor (PTPI), a cystatin superfamily variant has been studied by monitoring changes in the protein's spectroscopic and functional characteristics. Denaturation of PTPI in guanidine hydrochloride and urea resulted in altered intrinsic fluorescence emission spectrum, diminished negative circular dichroism, and loss of its papain inhibitory potential. Native like spectroscopic properties and inhibitory activity are only partially restored when denaturant is diluted from guanidine hydrochloride unfolded samples demonstrating that process is partially reversible. Coincidence of transition curves and dependence of transition midpoint (3.2M) on protein concentration in guanidine hydrochloride‐induced denaturation are consistent with a two‐state model involving a native like dimer and denatured monomer. On the contrary, urea‐induced unfolding of PTPI is a multiphasic process with indiscernible intermediates. The studies demonstrate that functional conformation and stability are governed by both ionic and hydrophobic interactions. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 708–717, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Unusual bipartite mode of interaction between the nonsense‐mediated decay factors,UPF1 and UPF2

Nonsense‐mediated decay (NMD) is a eukaryotic quality control mechanism that degrades mRNAs carrying premature stop codons. In mammalian cells, NMD is triggered when UPF2 bound to UPF3 on a downstream exon junction complex interacts with UPF1 bound to a stalled ribosome. We report structural studies on the interaction between the C‐terminal region of UPF2 and intact UPF1. Crystal structures, confirmed by EM and SAXS, show that the UPF1 CH‐domain is docked onto its helicase domain in a fixed configuration. The C‐terminal region of UPF2 is natively unfolded but binds through separated α‐helical and β‐hairpin elements to the UPF1 CH‐domain. The α‐helical region binds sixfold more weakly than the β‐hairpin, whereas the combined elements bind 80‐fold more tightly. Cellular assays show that NMD is severely affected by mutations disrupting the beta‐hairpin binding, but not by those only affecting alpha‐helix binding. We propose that the bipartite mode of UPF2 binding to UPF1 brings the ribosome and the EJC in close proximity by forming a tight complex after an initial weak encounter with either element.     

A Molecular Dynamics Simulation Study of Coaxial Stacking in RNA

Abstract

We report on unrestrained molecular dynamics simulations of an RNA tetramer binding to a tetra-nucleotide overhang at the 5′-end of an RNA hairpin (nicked structure) and of the corresponding continuous hairpin with Na⁺ as counterions. The simulations lead to stable structures and in this way a structural model for the coaxially stacked RNA hairpin is generated. The stacking interface in the coaxially stacked nicked hairpin structure is characterized by a reduced twist and shift and a slightly increased propeller twist as compared to the continuous system. This leads to an increased overlap between C22 and G23 in the stacking interface of the nicked structure. In the simulations the continuous RNA hairpin has an almost straight helical axis. On the other hand, the corresponding axis for the nicked structure exhibits a marked kink of 39°. The stacking interface exhibits no increased flexibility as compared to the corresponding base pair step in the continuous structure.     

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

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

Thermodynamics of the pseudo‐knot in helix 18 of 16S ribosomal RNA

A fragment of E. coli 16S rRNA formed by nucleotides 500 to 545 is termed helix 18. Nucleotides 505‐507 and 524‐526 form a pseudo‐knot and its distortion affects ribosome function. Helix 18 isolated from the ribosome context is thus an interesting fragment to investigate the structural properties and folding of RNA with pseudo‐knots. With all‐atom molecular dynamics simulations, spectroscopic and gel electrophoresis experiments, we investigated thermodynamics of helix 18, with a focus on its pseudo‐knot. In solution studies at ambient conditions we observed dimerization of helix 18. We proposed that the loop, containing nucleotides forming the pseudo‐knot, interacts with another monomer of helix 18. The native dimer is difficult to break but introducing mutations in the pseudo‐knot indeed assured a monomeric form of helix 18. Molecular dynamics simulations at 310 K confirmed the stability of the pseudo‐knot but at elevated temperatures this pseudo‐knot was the first part of helix 18 to lose the hydrogen bond pattern. To further determine helix 18 stability, we analyzed the interactions of helix 18 with short oligomers complementary to a nucleotide stretch containing the pseudo‐knot. The formation of higher‐order structures by helix 18 impacts hybridization efficiency of peptide nucleic acid and 2'‐O methyl RNA oligomers.     

Structural analysis of the N‐terminal fragment of the antiangiogenic protein endostatin: A molecular dynamics study

Endostatin is a potent antiangiogenic protein derived from the noncollagenous domain 1 (NC1) of collagen XVIII. The mechanism by which endostatin exerts its antiangiogenic effect is still incompletely understood. It has been shown that the 27 amino acid N‐terminal fragment of murine endostatin has antitumor, antimigration, and antipermeability activities comparable to the full soluble protein. To understand how this peptide can exert such elaborate function, we performed structural analysis using molecular dynamics to evaluate the behavior of this fragment in aqueous environment. Here, we show that the N‐terminal peptide of murine endostatin is able to assume a well‐defined structure, folding into a zinc‐dependent β‐hairpin conformation. Analyzing the folding mechanism, we were able to understand why the N‐terminal peptide of human endostatin with the same length failed to acquire a stable conformation. Conversely, we were able to predict the successful folding of the R4Q mutant and of a shorter form of the human peptide with 25 residues. Finally, we show that the β‐hairpin conformation assumed by the zinc‐bound peptide of murine endostatin has a high structural similarity with fragments of another family of angiogenesis inhibitors: the integrin‐binding portion of the NC1 domain of collagen IV. Indeed, our docking simulations show that arresten, canstatin, and the endostatin peptide bind to the same spot of αVβ3 integrin, suggesting similar interactions via a common binding site on this receptor. Proteins 2011;. © 2011 Wiley‐Liss, Inc.     

The influence of urea and guanidine chloride on the binding of the binding of the bacterial substrate and inhibitors to hen lysozyme at physiological temperature (40 degrees) (*).

The influence of urea and of guanidine chloride on the binding of the bacterial substrate and of inhibitors such as N-acetylglucosamine or chitotetraose to hen lysozyme were studied at 20 degrees and at 40 degrees C (physiological temperature). The action of urea did not prevent a certain degree of organization of the enzyme compatible with its usual behaviour in the presence of some inhibitors and with its crystallization ; guanidine chloride, already at low concentrations, seemed to have a more severe effect on lysozyme.     

Non‐equivalent roles of two periplasmic subunits in the function and assembly of triclosan pump TriABC from Pseudomonas aeruginosa

In Gram‐negative bacteria, multidrug efflux transporters function in complexes with periplasmic membrane fusion proteins (MFPs) that enable antibiotic efflux across the outer membrane. In this study, we analyzed the function, composition and assembly of the triclosan efflux transporter TriABC–OpmH from Pseudomonas aeruginosa. We report that this transporter possesses a surprising substrate specificity that encompasses not only triclosan but the detergent SDS, which are often used together in antibacterial soaps. These two compounds interact antagonistically in a TriABC‐dependent manner and negate antibacterial properties of each other. Unlike other efflux pumps that rely on a single MFP for their activities, two different MFPs, TriA and TriB, are required for triclosan/SDS resistance mediated by TriABC–OpmH. We found that analogous mutations in the α‐helical hairpin and membrane proximal domains of TriA and TriB differentially affect triclosan efflux and assembly of the complex. Furthermore, our results show that TriA and TriB function as a dimer, in which TriA is primarily responsible for stabilizing interactions with the outer membrane channel, whereas TriB is important for the stimulation of the transporter. We conclude that MFPs are engaged into complexes as asymmetric dimers, in which each protomer plays a specific role.     

Mistic's membrane association and its assistance in overexpression of a human GPCR are independent processes

The interaction of the Bacillus subtilis protein Mistic with the bacterial membrane and its role in promoting the overexpression of other membrane proteins are still matters of debate. In this study, we aimed to determine whether individual helical fragments of Mistic are sufficient for its interaction with membranes in vivo and in vitro. To this end, fragments encompassing each of Mistic's helical segments and combinations of them were produced as GFP‐fusions, and their cellular localization was studied in Escherichia coli. Furthermore, peptides corresponding to the four helical fragments were synthesized by solid‐phase peptide synthesis, and their ability to acquire secondary structure in a variety of lipids and detergents was studied by circular dichroism spectroscopy. Both types of experiments demonstrate that the third helical fragment of Mistic interacts only with LDAO micelles but does not partition into lipid bilayers. Interestingly, the other three helices interact with membranes in vivo and in vitro. Nevertheless, all of these short sequences can replace full‐length Mistic as N‐terminal fusions to achieve overexpression of a human G‐protein‐coupled receptor in E. coli, although with different effects on quantity and quality of the protein produced. A bioinformatic analysis of the Mistic family expanded the number of homologs from 4 to 20, including proteins outside the genus Bacillus. This information allowed us to discover a highly conserved Shine‐Dalgarno sequence in the operon mstX‐yugO that is important for downstream translation of the potassium ion channel yugO.     

Cooperative assembly of a nativelike ubiquitin structure through peptide fragment complexation: energetics of peptide association and folding

Peptide fragments corresponding to the N- and C-terminal portions of bovine ubiquitin, U(1-35) and U(36-76), are shown by NMR to associate in solution to form a complex of modest stability (Kassn approximately 1.4 x 10(5) M(-1) at pH 7.0), with NMR features characteristic of a nativelike structure. The complex undergoes cold denaturation, with temperature-dependent estimates of stability from NMR indicating a DeltaC(p) degrees for fragment complexation in good agreement with that determined for native ubiquitin, suggesting that fragment association results in the burial of a similar hydrophobic surface area. The stability of the complex shows appreciable pH dependence, suggesting that ionic interactions on the surface of the protein contribute significantly. However, denaturation studies of native ubiquitin in the presence of guanidine hydrochloride (Gdn.HCl) show little pH dependence, suggesting that ionic interactions may be "screened" by the denaturant, as recently suggested. Examination of the conformation of the isolated peptide fragments has shown evidence for a low population of nativelike structure in the N-terminal beta-hairpin (residues 1-17) and weak nascent helical propensity in the helical fragment (residues 21-35). In contrast, the C-terminal peptide (36-76) shows evidence in aqueous solution, from some Halpha chemical shifts, for nonnative phi and psi angles; nonnative alpha-helical structure is readily induced in the presence of organic cosolvents, indicating that tertiary interactions in both native ubiquitin and the folded fragment complex strongly dictate its structural preference. The data suggest that the N-terminal fragment (1-35), where interaction between the helix and hairpin requires the minimum loss of conformational entropy, may provide the nucleation site for fragment complexation.     

The Impact of Dyskeratosis Congenita Mutations on the Structure and Dynamics of the Human Telomerase RNA Pseudoknot Domain

Abstract

The pseudoknot domain is a functionally crucial part of telomerase RNA and influences the activity and stability of the ribonucleoprotein complex. Autosomal dominant dyskeratosis congenita (DKC) is an inherited disease that is linked to mutations in telomerase RNA and impairs telomerase function. In this paper, we present a computational prediction of the influence of two base DKC mutations on the structure, dynamics, and stability of the pseudoknot domain. We use molecular dynamics simulations, MM-GBSA free energy calculations, static analysis, and melting simulations analysis. Our results show that the DKC mutations stabilize the hairpin form and destabilize the pseudoknot form of telomerase RNA. Moreover, the P3 region of the predicted DKC-mutated pseudoknot structure is unstable and fails to form as a defined helical stem. We directly compare our predictions with experimental observations by calculating the enthalpy of folding and melting profiles for each structure. The enthalpy values are in very good agreement with values determined by thermal denaturation experiments. The melting simulations and simulations at elevated temperatures show the existence of an intermediate structure, which involves the formation of two UU base pairs observed in the hairpin form of the pseudoknot domain.     

A molecular dynamics simulation study of coaxial stacking in RNA

We report on unrestrained molecular dynamics simulations of an RNA tetramer binding to a tetra-nucleotide overhang at the 5'-end of an RNA hairpin (nicked structure) and of the corresponding continuous hairpin with Na+ as counterions. The simulations lead to stable structures and in this way a structural model for the coaxially stacked RNA hairpin is generated. The stacking interface in the coaxially stacked nicked hairpin structure is characterized by a reduced twist and shift and a slightly increased propeller twist as compared to the continuous system. This leads to an increased overlap between C22 and G23 in the stacking interface of the nicked structure. In the simulations the continuous RNA hairpin has an almost straight helical axis. On the other hand, the corresponding axis for the nicked structure exhibits a marked kink of 39 degrees. The stacking interface exhibits no increased flexibility as compared to the corresponding base pair step in the continuous structure.     

Static and dynamic quenching of protein fluorescence. II. lysozyme.

The fluorescence parameters—lifetime, relative quantum yield, wavelength of maximum fluorescence intensity, half-width, and polarization—of 0.01% lysozyme were measured at 15°C in aqueous solution, in glycerol–water mixtures (0–90% v/v glycerol), in aqueous urea (0–8M) solutions, and in aqueous guanidine hydrochloride (0–6.4M) solutions. The changes in the static and dynamic quenching of lysozyme fluorescence, monitored by the quantum yield and lifetime measurements, were correlated with the other fluorescence parameters and compared with our earlier results with bovine serum albumin. The results were interpreted in terms of induced conformational changes. The various perturbants altered the fluorescence parameters of lysozyme and bovine serum albumin very differently. The differences were shown to be entirely consistent with our earlier conclusion that bovine serum albumin fluorophores are nonsurface residues and with the conclusion of others that lysozyme fluorophores are surface residues. Unlike their effects on bovine serum albumin, urea and guanidine hydrochloride affect lysozyme structure quite differently, both in nature and degree. We have suggested that the affect of urea on lysozyme fluorescence is an indirect result of reduction in the size of the cleft brought about by the structure-breaking action of urea on water in the cleft. 4M Urea is sufficient for this reaction. Large decreases in the polarization of the fluorescence of lysozyme in the 0.8–1.6M and 3.2–4.8M guanidine hydrochloride ranges demonstrated two guanidine hydrochloride-induced conformation changes. A red shift of the fluorescence maximum to 354 nm indicated that the second transition completely exposes all fluorescing tryptophan residues of lysozyme to mobile solvent water. However, even 6.4M guanidine hydrochloride did not completely unravel the lysozyme molecule at 15°C, as evidenced by its failure to cause any of the tyrosine residues to become fluorescent.     

Effect of glutathione on ribonuclease

GSH, but not GSSG, inhibits the reactivation by phosphate ion of ribonuclease activity inactivated by urea or guanidine. The effects of GSH are rather slow and pretreatment of ribonuclease with urea is a requisite for the inhibitory action of GSH on enzyme reactivation. GSH is more effective in urea than in guanidine and its action is greatly enhanced by EDTA. An optimum pH of about 9.0 was found for the inhibitory effect of GSH. Titration of the thiol groups formed after inactivation of ribonuclease by GSH strongly suggests that the reduction of only one disulphide linkage is involved. The reduction of this bond is sufficient to completely abolish the enzymic activity.     

Solution effects on the self-association of a water-soluble peptoid

A desire to replicate the structural and functional complexity of proteins with structured, sequence-specific oligomers motivates study of the structural features of water-soluble peptoids (N-substituted glycine oligomers). Understanding the molecular-level details of peptoid self-assembly in water is essential to advance peptoids' application as novel materials. Peptoid 1 , an amphiphilic, putatively helical peptoid previously studied in our laboratory, shows evidence of self-association in aqueous solution. In this work, we evaluate how changes to aqueous solution conditions influence the self-association of 1 . We report that changes to pH influence the fluorescence and CD spectroscopic features as well as the peptoid's interaction with a solvatochromic fluorophore and its apparent size as estimated by size exclusion chromatography. Addition of guanidine hydrochloride and ammonium sulfate also modulate spectroscopic features of the peptoid, its interaction with a solvatochromic fluorophore, and its elution in size exclusion chromatography. These data suggest that the ordering of the self-assembly changes in response to pH and with solvent additives and is more ordered at higher pH and in the presence of guanidine hydrochloride. The deeper understanding of the self-association of 1 afforded by these studies informs the design of new stimuli-responsive peptoids with stable tertiary or quaternary structures.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin‐binding protein G from Streptococcus. IV. Implication for the mechanism of folding of the parent protein

A 34‐residue α/β peptide [IG(28–61)], derived from the C‐terminal part of the B3 domain of the immunoglobulin binding protein G from Streptoccocus, was studied using CD and NMR spectroscopy at various temperatures and by differential scanning calorimetry. It was found that the C‐terminal part (a 16‐residue‐long fragment) of this peptide, which corresponds to the sequence of the β‐hairpin in the native structure, forms structure similar to the β‐hairpin only at T = 313 K, and the structure is stabilized by non‐native long‐range hydrophobic interactions (Val47–Val59). On the other hand, the N‐terminal part of IG(28–61), which corresponds to the middle α‐helix in the native structure, is unstructured at low temperature (283 K) and forms an α‐helix‐like structure at 305 K, and only one helical turn is observed at 313 K. At all temperatures at which NMR experiments were performed (283, 305, and 313 K), we do not observe any long‐range connectivities which would have supported packing between the C‐terminal (β‐hairpin) and the N‐terminal (α‐helix) parts of the sequence. Such interactions are absent, in contrast to the folding pathway of the B domain of protein G, proposed recently by Kmiecik and Kolinski (Biophys J 2008, 94, 726–736), based on Monte‐Carlo dynamics studies. Alternative folding mechanisms are proposed and discussed. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 469–480, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

The impact of dyskeratosis congenita mutations on the structure and dynamics of the human telomerase RNA pseudoknot domain

The pseudoknot domain is a functionally crucial part of telomerase RNA and influences the activity and stability of the ribonucleoprotein complex. Autosomal dominant dyskeratosis congenita (DKC) is an inherited disease that is linked to mutations in telomerase RNA and impairs telomerase function. In this paper, we present a computational prediction of the influence of two base DKC mutations on the structure, dynamics, and stability of the pseudoknot domain. We use molecular dynamics simulations, MM-GBSA free energy calculations, static analysis, and melting simulations analysis. Our results show that the DKC mutations stabilize the hairpin form and destabilize the pseudoknot form of telomerase RNA. Moreover, the P3 region of the predicted DKC-mutated pseudoknot structure is unstable and fails to form as a defined helical stem. We directly compare our predictions with experimental observations by calculating the enthalpy of folding and melting profiles for each structure. The enthalpy values are in very good agreement with values determined by thermal denaturation experiments. The melting simulations and simulations at elevated temperatures show the existence of an intermediate structure, which involves the formation of two UU base pairs observed in the hairpin form of the pseudoknot domain.     

Dielectric screening effect of electronic polarization and intramolecular hydrogen bonding

Recent site‐resolved hydrogen exchange measurements have uncovered significant discrepancies between simulations and experimental data during protein folding, including the excessive intramolecular hydrogen bonds in simulations. This finding indicates a possibility that intramolecular charge–charge interactions have not included sufficient dielectric screening effect of the electronic polarization. Scaling down peptide atomic charges according to the optical dielectric constant is tested in this study. As a result, the number of intramolecular hydrogen bonds is lower than using unscaled atomic charges while reaching the same levels of helical contents or β‐hairpin backbone hydrogen bonds, because van der Waals interactions contribute substantially to peptide folding in water. Reducing intramolecular charge–charge interactions and hydrogen bonding increases conformational search efficiency. In particular, it reduces the equilibrium helical content in simulations using AMBER force field and the energy barrier in folding simulations using CHARMM force field.