Theoretical modeling of peptide α-helical circular dichroism in aqueous solution

Reliable modeling of protein and peptide circular dichroism (CD) spectra in the far UV presents a challenge for current theoretical approaches. In this study, the time-dependent density functional theory (TDDFT), configuration interaction with single excitation (CIS), and transition dipole coupling (TDC) were used to assess the most important factors contributing to the CD spectra of the α-helical secondary structure. The dependence on the peptide chain length and also the role of the flexibility and solvent environment were investigated with a model oligopeptide Ac-(Ala)(N)-NH-Me, (N = 1, ..., 18). Both the TDDFT and TDC-like methods suggest that the CD curve typical for the α-helix arises gradually, but its basic characteristic is discernible already for peptides with 4-5 amino acid residues. The calculated dependence was in a qualitative agreement with experimental spectra of short α-helices stabilized by the histidine-metal binding. The TDDFT computations of the CD were found to be unusually sensitive to the basis set and solvent model. Explicit hydration and temperature fluctuations of the peptide geometry, simulated with the aid of molecular dynamics (MD), significantly influenced the CD and absorption spectral shapes. An extensive averaging over MD configurations is thus required to obtain a converged spectral profile in cluster simulations. On the other hand, both the TDDFT and TDC models indicate only a minor influence of the alanine side chains. The CIS and TDC calculations also point toward a relatively small effect of the helix-helix interaction on the CD spectral profiles. For a model system of two helices, the CIS method predicted larger changes in the spectra than TDC. This suggests other than interactions between peptide chains, such as mutual polarization, can have a minor, but measurable, effect on the CD spectrum.     

A new “hydrophobic template” method detects segments forming transmembrane α-helical bundles in ion channels

We present a “hydrophobic template” method enabling recognition of α-helix bundles in membrane channels from sequence analysis. Inspection of hydrophobic properties of pore-forming helices in proteins with known structure (A-B₅ toxins) permits delineation of a common polarity motif: two hydrophobic surface stretches separated by polar areas. The bundles are stabilized by nonpolar interhelical contacts. A number of transmembrane segments were checked for presence of this motif, and it was detected for pore-forming helices of several ion transporters (segments M2 of acetylcholine and GABA_A receptors, α5 peptide of δ-endotoxin), which reveal five α-helix bundle architecture. Applications of the method to modeling of membrane channels are discussed. Received: 24 April 1998 / Accepted: 3 September 1998 / Published online: 7 December 1998     

Unravelling the interaction between α-cyclodextrin with the thaumatin protein and a peptide mimic

G-protein-coupled receptors (GPCRs) are responsible for signal transduction; through these transmembrane proteins, our senses are evoked: sight, smell and taste. Thaumatin is a natural sweet-tasting protein that is 100,000 times sweeter than sucrose but its use in food products has been hampered due to a liquorice aftertaste. Thaumatin has been shown to bind to a class C GPCR and the active binding site of the thaumatin protein is known. Here, we report on the binding of a well-known food grade host: α-cyclodextrin to thaumatin. We show through a combination of one- and two-dimensional NMR experiments that α-cyclodextrin binds to aromatic residues on thaumatin with K_a = 8.5 ± 2.4 M ^? 1. We also synthesise a heptapeptide KTGDRGF that mimics the active binding site of thaumatin and show that α-cyclodextrin binds to the C-terminal solvent accessible phenylalanine residue of this peptide with K_a = 8.8 ± 3.1 M ^? 1. This indicates that α-cyclodextrin may interact with the active binding site on thaumatin, suggesting that α-cyclodextrin could be used to modify the interaction of thaumatin with GPCRs and hence its sweet-taste profile.     

The Preparation of α-Phosphonovinylzirconocenes and their Application in the Stereospecific Synthesis of α-Halo-l-alkenylphosphonates

Hydrozirconation of 1-alkylnylphosphonates gives the organozirconium(Ⅳ) complexes 2 in syn-addition way.Complexs 2 was trapped with NCS,NBS or I2 to afford stereodifined α-halo-1-alkenylphosphonates in moderate to high yields.     

A direct link between the Passerini reaction and α-lactams

α-Lactams (aziridinones) can function to replace two of the three reactants, the oxo-compound and the isonitrile, in the Passerini reaction. Four α-lactams (5a-d) were reacted with mono- and dicarboxylic acids of positive pK_a values to give 2-acyloxycarboxamides (4) and bis-2-acyloxycarboxamide products 12 and 13, respectively. The same compounds were also prepared via the Passerini reaction. Acids with a negative pK_a decarbonylate α-lactams to give immonium salts. The main path of the reaction depends on the pK_a of the acid component, the reactivity of the α-lactam, and the reaction conditions.     

The agonistic binding site at the histamine H2 receptor. II. Theoretical investigations of histamine binding to receptor models of the seven α-helical transmembrane domain

Summary In the first part (pp. 461–478 in this issue) of this study regarding the histamine H₂ receptor agonistic binding site, the best possible interactions of histamine with an -helical oligopeptide, mimicking a part of the fifth transmembrane -helical domain (TM5) of the histamine H₂ receptor, were considered. It was established that histamine can only bind via two H-bonds with a pure -helical TM5, when the binding site consists of Tyr¹⁸²/Asp¹⁸⁶ and not of the Asp¹⁸⁶/Thr¹⁹⁰ couple. In this second part, two particular three-dimensional models of G-protein-coupled receptors previously reported in the literature are compared in relation to agonist binding at the histamine H₂ receptor. The differences between these two receptor models are discussed in relation to the general benefits and limitations of such receptor models. Also the pros and cons of simplifying receptor models to a relatively easy-to-deal-with oligopeptide for mimicking agonistic binding to an agonistic binding site are addressed. Within complete receptor models, the simultaneous interaction of histamine with both TM3 and TM5 can be analysed. The earlier suggested three-point interaction of histamine with the histamine H₂ receptor can be explored. Our results demonstrate that a three-point interaction cannot be established for the Asp⁹⁸/Asp¹⁸⁶/Thr¹⁹⁰ binding site in either of the investigated receptor models, whereas histamine can form three H-bonds in case the agonistic binding site is constituted by the Asp⁹⁸/Tyr¹⁸²/Asp¹⁸⁶ triplet. Furthermore this latter triplet is seen to be able to accommodate a series of substituted histamine analogues with known histamine H₂ agonistic activity as well.     

The Preparation of α-Phosphonovinylzirconocenes and their Application in the Stereospecific Synthesis of α-Halo-1 -alkenylphosphonates

Hydrozirconation of 1-alkylnylphosphonates gives the organozirconium(Ⅳ) complexes 2 in syn-addition way. Complexs 2 was trapped with NCS, NBS or I2 to afford stereodifined α-halo-1-alkenylphosphonates in moderate to high yields.     

The Application of the Trichloroacetimidate Method to the Synthesis of α-D-Gluco- and α-D-Galactopyranosides

Abstract

The trichloroacetimidate method has been applied to the construction of α-d-galacto- and α-d-glucopyranosides. The readily available β-trichloroacetimidates of 2,3,4,6-tetra-O-benzyl-d-galacto- and glucopyranose (1-β and 3-β, respectively) have been employed in glycosidations with several monosaccharides (either A, B, C or D) under varying experimental conditions. With the galactose derivative 1-β as a donor and each of the monosaccharides A-D as acceptors, the corresponding disaccharides 1A-1D, were obtained in high yield and with good α-stereoselectivity when employing diethyl ether as solvent and either trimethylsilyl- or tert-butyldimethylsilyl trifluoro-methane sulphonate as catalyst. Glycosidations with the glucose derivative 3-β, as donor, and with the monosaccharide acceptors A, B or D, gave the corresponding disaccharides 3A, 3B and 3D, in high yield but with somewhat lower α-diastereoselectivity than observed with the galactose derivative 1-β. The stereochemical outcome of the reactions is rationalised in terms of possible reaction mechanisms.     

The fast and the slow: folding and trapping of λ6-85

Molecular dynamics simulations combining many microsecond trajectories have recently predicted that a very fast folding protein like lambda repressor fragment λ(6-85) D14A could have a slow millisecond kinetic phase. We investigated this possibility by detecting temperature-jump relaxation to 5 ms. While λ(6-85) D14A has no significant slow phase, two even more stable mutants do. A slow phase of λ(6-85) D14A does appear in mild denaturant. The experimental data and computational modeling together suggest the following hypothesis: λ(6-85) takes only microseconds to reach its native state from an extensively unfolded state, while the latter takes milliseconds to reach compact β-rich traps. λ(6-85) is not only thermodynamically but also kinetically protected from reaching such "intramolecular amyloids" while folding.     

Confinement in nanopores can destabilize α-helix folding proteins and stabilize the β structures

Protein folding in confined media has attracted wide attention over the past decade due to its importance in both in vivo and in vitro applications. Currently, it is generally believed that protein stability increases by decreasing the size of the confining medium, if its interaction with the confining walls is repulsive, and that the maximum folding temperature in confinement occurs for a pore size only slightly larger than the smallest dimension of the folded state of a protein. Protein stability in pore sizes, very close to the size of the folded state, has not however received the attention that it deserves. Using detailed, 0.3-ms-long molecular dynamics simulations, we show that proteins with an α-helix native state can have an optimal folding temperature in pore sizes that do not affect the folded-state structure. In contradiction to the current theoretical explanations, we find that the maximum folding temperature occurs in larger pores for smaller α-helices. In highly confined pores the free energy surface becomes rough, and a new barrier for protein folding may appear close to the unfolded state. In addition, in small nanopores the protein states that contain the β structures are entropically stabilized, in contrast to the bulk. As a consequence, folding rates decrease notably and the free energy surface becomes rougher. The results shed light on many recent experimental observations that cannot be explained by the current theories, and demonstrate the importance of entropic effects on proteins' misfolded states in highly confined environments. They also support the concept of passive effect of chaperonin GroEL on protein folding by preventing it from aggregation in crowded environment of biological cells, and provide deeper clues to the α → β conformational transition, believed to contribute to Alzheimer's and Parkinson's diseases. The strategy of protein and enzyme stabilization in confined media may also have to be revisited in the case of tight confinement. For in silico studies of protein folding in confined media, use of non-Go potentials may be more appropriate.     

The A53T mutation is key in defining the differences in the aggregation kinetics of human and mouse α-synuclein

Despite a 95% sequence similarity, the aggregation of human and mouse α-synuclein is remarkably different, as the human form is slower than the mouse form in forming fibrils but is associated with Parkinson's disease in both humans and transgenic mice. Here, the amino acid code underlying these differences is investigated by comparing the lag times, growth rates, and secondary structure propensities of a systematic series of eight human-mouse chimeras. Fluorescence analysis of these variants shows that the A53T substitution dominates the growth kinetics, while the lag phase is affected by a combination of the A53T and S87N substitutions. The secondary structure propensities derived from an NMR chemical shift analysis of the monomeric forms of the human-mouse variants enable us to establish a link between the changes in the conformational properties in the region of position 53 upon mutation and the corresponding changes in growth rates. These results suggest that the presence of an alanine residue at position 53 may be an evolutionary adaptation to minimize Parkinson's disease in humans and indicates that effective drug development efforts may be directed to target this N-terminal region of the sequence.     

Dynamic structure of transmembrane α-helical fragments of bacteriorhodopsin in lipid bilayer characterized by 13C chemical shift tensor and hydrogen bond distance by REDOR NMR

Dynamic structure of [1-¹³C]Ala¹⁴, [¹⁵N]Ala¹⁸-labeled transmembrane fragment A (6–42) of bacteriorhodopsin in DMPC bilayer was examined by careful analysis of ¹³C chemical shift tensors of [1-¹³C]Ala¹⁴ residue and ¹⁵N–H⋯O¹³C hydrogen bond length, as determined by its side-band patterns and rotational echo double resonance method (REDOR), respectively. It was found that rigid-body motions of this helical segment, fluctuating with a time scale of 10⁻⁴ s and dominant at a temperature above 10°C [Biopolymers 58 (2001) 78], ceased at a temperature below 0°C, to give rise to a static spectral pattern similar to that of the solid sample. No evidence was obtained for the presence of the distorted α-helix form in the DMPC bilayer, in view of the similarity between the principal components of ¹³C chemical shift tensors of the lipid bilayer and the solid. Furthermore, accurate C⋯N interatomic distance was measured between ¹⁵N and ¹³C nuclei corresponding to the hydrogen bond length, ¹⁵N–H (Ala¹⁸)⋯O¹³C (Ala¹⁴), characteristic of the α-helix form by REDOR NMR. The C⋯N distance turned out to be 4.5 Å at −30°C, which is in good agreement with that in the solid state taking the ordinary α-helix as estimated from the ¹³C chemical shifts. This finding shows that there is no anomaly in the α-helix form in the lipid bilayer. In this connection, it is emphasized that the ¹³C C_β peak of Ala residue resonating at unusually lower field is a very unique probe to be able to detect the low frequency, local anharmonic motions of the order of 10⁻² s in the lipid bilayer.     

Competitive electron transfers from a tyrosyl side-chain and peptide bond in the photodegradation of N-tosyl alpha-aminomethylamides: an insight into photosynthesis and photodamage in the biological oxidation of water?

Photo-excited N-tosyl derivatives of phenylalanyl- and, more particularly, O-methyltyrosylmethylamides undergo electron transfer from aryl to tosyl groups whereas the photo-degradation of aliphatic analogues is initiated by electron transfer from the peptide bond, suggesting the latter as one possible reason for the rapid turnover of the D1 protein in biological water oxidation when the essential mediating role of tyrosine 116 in the PSII complex is inhibited.     

Comparison on effect of hydrophobicity on the antibacterial and antifungal activities of α-helical antimicrobial peptides

HPRP-A1, a 15-mer α-helical cationic peptide, was derived from N-terminus of ribosomal protein L1 (RpL1) of Helicobacter pylori. In this study, HPRP-A1 was used as a framework to obtain a series of peptide analogs with different hydrophobicity by single amino acid substitutions in the center of nonpolar face of the amphipathic helix in order to systematically study the effect of hydrophobicity on biological activities of -helical antimicrobial peptides. Hydrophobicity and net charge of peptides played key roles in the biological activities of these peptide analogs; HPRP-A1 and peptide analogs with relative higher hydrophobicity exerted broad spectrum antimicrobial activity against Gram-negative bacteria, Gram-positive bacteria and pathogenic fungi, but also showed stronger hemolytic activity; the change of hydrophobicity and net charge of peptides had similar effects with close trend and extent on antibacterial activities and antifungal activities. This indicated that there were certain correlations between the antibacterial mode of action and the antifungal mode of action of these peptides in this study. The peptides exhibited antimicrobial specificity for bacteria and fungi, which provided potentials to develop new antimicrobial drugs for clinical practices.     

The α-effect of oximate ions and solvent effect in the transfer of the toluenesulfonyl group in aqueous dimethylsulfoxide

The kinetic behaviour of oximate ions was studied in the transfer of the toluenesulfonyl group in aqueous dimethylsulfoxide (0–95 vol.% DMSO). The solvation effects of the solvent are a factor controlling the nucleophilicity and magnitude of the α-effect of the oximate ions. L. M. Litvinenko Institute of Physical Organic and Coal Chemistry, National Academy of Sciences of Ukraine, 70 R. Lyuksemburg ul., Donetsk 340114, Ukraine. Translated from Teoreticheskaya i éksperimental'naya Khimiya, Vol. 36, No. 1, pp. 30–35, January–February, 2000.     

A sequence-based computational model for the prediction of the solvent accessible surface area for α-helix and β-barrel transmembrane residues

Predicting the solvent accessible surface area (ASA) of transmembrane (TM) residues is of great importance for experimental researchers to elucidate diverse physiological processes. TM residues fall into two major structural classes (α-helix membrane protein and β-barrel membrane protein). The reported solvent ASA prediction models were developed for these two types of TM residues respectively. However, this prevents the general use of these methods because one cannot determine which model is suitable for a given TM residue without information of its type. To conquer this limitation, we developed a new computational model that can be used for predicting the ASA of both TM α-helix and β-barrel residues. The model was developed from 78 α-helix membrane protein chains and 24 β-barrel membrane protein. Its prediction ability was evaluated by cross validation method and its prediction result on an independent test set of 20 membrane protein chains. The results show that our model performs well for both types of TM residues and outperforms other prediction model which was developed for the specific type of TM residues. The prediction results also proved that the random forest model incorporating conservation score is an effective sequence-based computational approach for predicting the solvent ASA of TM residues.     

On the orientation of a designed transmembrane peptide: toward the right tilt angle?

The orientation of the transmembrane peptide WALP23 under small hydrophobic mismatch has been assessed through long-time-scale molecular dynamics simulations of hundreds of nanoseconds. Each simulation gives systematically large tilt angles (>30 degrees). In addition, the peptide visits various azimuthal rotations that mostly depend on the initial conditions and converge very slowly. In contrast, small tilt angles as well as a well-defined azimuthal rotation were suggested by recent solid-state 2H NMR studies on the same system. To optimally compare our simulations with NMR data, we concatenated the different trajectories in order to increase the sampling. The agreement with 2H NMR quadrupolar splittings is spectacularly better when these latter are back-calculated from the concatenated trajectory than from any individual simulation. From these ensembled-average quadrupolar splittings, we then applied the GALA method as described by Strandberg et al. (Biophys J. 2004, 86, 3709-3721), which basically derives the peptide orientation (tilt and azimuth) from the splittings. We find small tilt angles (6.5 degrees), whereas the real observed tilt in the concatenated trajectory presents a higher value (33.5 degrees). We thus propose that the small tilt angles estimated by the GALA method are the result of averaging effects, provided that the peptide visits many states of different azimuthal rotations. We discuss how to improve the method and suggest some other experiments to confirm this hypothesis. This work also highlights the need to run several and rather long trajectories in order to predict the peptide orientation from computer simulations.     

Concurrent display of both α- and β-turns in a model peptide

This article describes a model peptide that concurrently displays both α- and β-turns, as demonstrated by structural investigations using single crystal X-ray crystallography and solution-state NMR studies. The motif reported herein has the potential for the design of novel conformationally ordered synthetic oligomers with structural architectures distinct from those classically observed.     

Structures of α-type ions formed in the 157 nm photodissociation of singly-charged peptide ions

One hundred fifty-seven nm photodissociation of singly-charged peptide ions induces the cleavage of alpha-carbon to carbonyl-carbon bonds along the backbone. a(n) + 1 radical ions are observed as the primary photolysis products of peptides with N-terminal arginines in a linear ion trap mass spectrometer. The radical elimination pathways undertaken by the a(n) + 1 radical ions to form more stable even-electron species are studied in hydrogen-deuterium (H/D) exchange experiments. Two types of a(n) ions along with d-type ions are observed as secondary elimination products. The relative abundance of each depends on the C-terminal residue of the radical fragment ion.     

The Molecular Basis of the Interaction of Cyclophilin A with α-Synuclein

Peptidylprolyl isomerases (PPIases) catalyze cis/trans isomerization of prolines. The PPIase CypA colocalizes with the Parkinson's disease (PD)-associated protein α-synuclein in cells and interacts with α-synuclein oligomers. Herein, we describe atomic insights into the molecular details of the α-synuclein/CypA interaction. NMR spectroscopy shows that CypA catalyzes isomerization of proline 128 in the C-terminal domain of α-synuclein. Strikingly, we reveal a second CypA-binding site formed by the hydrophobic sequence ⁴⁷GVVHGVATVA⁵⁶, termed PreNAC. The 1.38 Å crystal structure of the CypA/PreNAC complex displays a contact between alanine 53 of α-synuclein and glutamine 111 in the catalytic pocket of CypA. Mutation of alanine 53 to glutamate, as found in patients with early-onset PD, weakens the interaction of α-synuclein with CypA. Our study provides high-resolution insights into the structure of the PD-associated protein α-synuclein in complex with the most abundant cellular cyclophilin.