NMR assignments for the cis and trans forms of the hemolysin II C-terminal domain

Pathogenic bacteria secrete pore-forming toxins (PFTs) to selectively defend against immune cells and to break through cellular barriers in the host. Understanding how PFTs attack cell membranes is not only essential for therapeutic intervention but for designing agents to deliver drugs to specific human cell subtypes, for example in anti-cancer or anti-viral therapies. Many toxins contain accessory domains that help recognize specific molecular epitopes on the membranes of target cells, including proteins, carbohydrates, and lipids. Here we report NMR assignments for the 94-residue 10 kDa C-terminal accessory domain of Bacillus cereus hemolysin II, HlyIIC, that has no known structural or functional homologues. The HlyIIC domain exists in a dynamic equilibrium due to cis/trans isomerization of its Gly86–Pro87 peptide bond. The cis and trans forms are about equally populated and are in slow exchange on the NMR timescale, giving rise to separate signals for approximately half of the residues in the domain. Assignments for the cis and trans forms were achieved with the aid of a P87M mutant that stabilizes the trans form, and separate sequential walks for the two forms in 3D NMR spectra of the wild-type HlyIIC. Based on backbone chemical shifts, the domain has a α1–α2–β1–β2–β3–β4–α3–β5 order of secondary structure elements. The last strand in the trans form and in the P87M mutant is shortened near Pro87 compared to the cis form. Both cis/trans isomerization and the P87M mutation cause large chemical shift changes throughout HlyIIC, suggesting that the proline is important in stabilizing the structure of the domain. The NMR assignments pave the way for solving the structures of the multiple conformational forms of HlyIIC and establishing their mechanism of interconversion.     

Effect of proline residues on protein folding

Conformational energy calculations have been used to study the role of the proline residues in the folding of bovine pancreatic trypsin inhibitor. In the calculation, each of the four proline residues of this small protein is forced from the trans to cis peptide isomer while still part of the native folded structure. The cis proline residue can always be accommodated by small changes of the native conformation (< 1 Å root-mean-square deviation). For three of the four proline residues, Pro2, Pro9 and Pro 13, being in the cis form is calculated to destabilize the folded conformation by less than 11 kcal/mol, suggesting that rapid folding to a stable native-like conformation can occur with either isomeric form. For one of these three, Pro13, the destabilization is only 1 kcal/mol, suggesting the existence of an alternative folded native conformation with Pro13 cis. The fourth proline residue, Pro8, is calculated to destabilize the native conformation by so much (33 kcal/mol) that it will block folding in the manner proposed by Brandts et al. (1975).     

The effect of phosphorylation on the conformation of oligo-peptides with Ser–Pro motif: a molecular dynamics simulation

Model tetrapeptide system was designed to investigate the cis/trans isomerization of peptidyl-prolyl imide bond of Ser–Pro motif. To establish the side-chain O-phosphorylation effect in regulating the peptides conformations, molecular dynamics (MD) simulations where carried out on the designed tetrapeptides and their corresponding phosphorylated forms by MD Insight II Discovery3 approach. The most stable configurations and the statistic cis/trans concentration distribution demonstrated that the phosphorylation evidently influences the peptidyl-prolyl imide bond isomerization and works as a key effect in regulating the peptide conformations. The charge state and the site provided for the charge of the phosphate moiety might be an important key. The results also demonstrated that phosphorylation changes the cis conformation ratio of the peptide and the maximum cis value is obtained when the phosphate group has no negative charge.     

Nmr studies of the rates of proline cis–trans isomerization in oligopeptides

¹H-Nmr was used to measure the rate of cis–trans interconversion of X-Pro bonds in linear and cyclic oligopeptides. k(cis → trans) = 2.5 × 10^?3 s^?1 at 25°C was found for the zwitterionic form of H-Ala-Pro-OH, in good agreement with earlier measurements. Replacement of Ala by Phe, Tyr, or Trp resulted in a 10-fold slower interconversion rate, whereas after substitution of Ala by His or Glu, the rate decreased only slightly. Independent of the residues X, the interconversion rate was increased by a factor of ca. 20 when the peptide chain was elongated by addition of Ala to the C-terminal Pro. An additional increase by a factor of 6 was observed when going from the protected linear peptide CF₃CO-Gly-Gly-Pro-Ala-OCH₃ to the closely related cyclic compound c[-Gly-Gly-Pro-Gly-Ala-]. These data are evaluated with regard to their possible use in future studies on the role of X-Pro cis–trans isomerization in the kinetics of protein folding.     

Computational study on the mechanisms of action of the potential anticancer drug trans-isopropylaminedimethylaminedichloroplatinum (trans-IPADMADP) and its cis isomer with DNA purine bases

The monofunctional and bifunctional bindings of the potential anticancer drug trans-isopropylaminedimethylaminedichloroplatinum (trans-IPADMADP) and its cis isomer to purine base in DNA are explored by using density functional theory and IEF-PCM solvation models. The computed lowest free energy barrier in the aqueous solution is 14.0/11.6 kcal/mol (from trans-Pt-chloroaqua complex to trans-/cis-monoadduct) for guanine(G), and 11.7/13.3 kcal/mol (from trans-Pt-chloroaqua complex to trans-/cis-monoadduct) for adenine(A). Our calculations demonstrate that the trans reactant complexes (or isolated reactants) can generate trans- or cis-monoadducts via similar trigonal bipyramidal transition state structures, suggesting that the monoadducts can subsequently close to form the bifunctional intrastrand Pt-DNA adducts and simultaneously distort DNA in the similar way as cisplatin. Our calculations show that Pt(isopropylamine)(dimethylamine)G₂²⁺ head-to-head path has the lowest free energy of activation at 17.6 kcal/mol, closely followed by the Pt(isopropylamine)(dimethylamine)GA²⁺ head-to-head path at 19.6 kcal/mol when the monofunctional cis-Pt-G complex serves as the reactant; while the Pt(isopropylamine)(dimethylamine)G₂²⁺ head-to-tail adduct has the lowest barrier of 20.5 kcal/mol, closely followed by the Pt(isopropylamine)(dimethylamine)GA²⁺ head-to-tail adduct at 23.0 kcal/mol if the monofunctional trans-Pt-G complex is the reactant.The calculated relatively lower activation energy barrier than that of cisplatin theoretically confirm that trans-[PtCl₂(isopropylamine)(dimethylamine)] is a potential anticancer drug as described by experiment.     

Molecular Determinants of a Native-State Prolyl Isomerization

Prolyl cis/trans isomerizations determine the rates of many protein-folding reactions, and they can serve as molecular switches and timers. The energy required to shift the prolyl cis/trans equilibrium during these processes originates from conformational reactions that are linked structurally and energetically with prolyl isomerization. We used the N2 domain of the gene-3-protein of phage fd to elucidate how such an energetic linkage develops in the course of folding. The Asp160-Pro161 bond at the tip of a β hairpin of N2 is cis in the crystal structure, but in fact, it exists as a mixture of conformers in folded N2. During refolding, about 10 kJ mol^− 1 of conformational energy becomes available for a 75-fold shift of the cis/trans equilibrium constant at Pro161, from 7/93 in the unfolded to 90/10 in the folded form. We combined single- and double-mixing kinetic experiments with a mutational analysis to identify the structural origin of this proline shift energy and to elucidate the molecular path for the transfer of this energy to Pro161. It originates largely, if not entirely, from the two-stranded β sheet at the base of the Pro161 hairpin. The two strands improve their stabilizing interactions when Pro161 is cis, and this stabilization is propagated to Pro161, because the connector peptides between the β strands and Pro161 are native-like folded when Pro161 is cis. In the presence of a trans-Pro161, the connector peptides are locally unfolded, and thus, Pro161 is structurally and energetically uncoupled from the β sheet. Such interrelations between local folding and prolyl isomerization and the potential modulation by prolyl isomerases might also be used to break and reestablish slow communication pathways in proteins.     

Conformational energy studies of linear dipeptides H-X-L-Pro-OH

Empirical conformational energy calculations with the use of ECEPP energy functions have been carried out for linear dipeptides H-X-L -Pro-OH, with X = Gly, L -Ala, D -Ala, L -Leu, D -Leu, L -Phe, and D -Phe, in different states of protonation of the end groups. The results of these calculations are compared with the previously reported experimental equilibrium populations for the cis and trans isomers of the X-Pro bond in the different species. For all the protonation states of the seven dipeptides, the calculated nonbonded interactions and the conformational entropy term lead to a preference of the trans forms over the cis isomers by at least 1 kcal/mol. The electrostatic interactions stabilize the cis conformations in all species except the cationic forms of the D ,L -peptides, and it could further be shown that only the carbonyl group of X and the two end groups contribute significantly to the total electrostatic energy. One of the principal results of the experimental studies, i.e., the occurrence of 5–15% cis-proline in all the peptides with an uncharged C-terminus, was corroborated by our investigation of the cationic species. A detailed assessment of the electrostatic contribution to the total energy of the different conformations of H-Gly-L -Pro-OH indicates that the standard ECEPP parameters tend to overestimate the electrostatic interactions in aqueous solutions of the X-Pro dipeptides.     

Electronic structures of collagen model polymers: (Gly-Pro)n, (Gly-Hyp)n, (Ala-Pro)n, (Ala-Hyp)n, (Gly-Pro-Gly)n, (Gly-Hyp-Gly)n, (Gly-Pro-Pro)n, and (Gly-Pro-Hyp)n

In order to investigate the relationship existing between the electronic structures of collagen and its biochemical functions in vivo, the semiempirical CNDO/2 SCF MO calculations were carried out on several model polymers of collagen, (Gly-Pro)_n, (Gly-Hyp)_n, (Ala-Pro)_n, (Ala-Hyp)_n, (Gly-Pro-Gly)_n, (Gly-Hyp-Gly)_n, (Gly-Pro-Pro)_n and (Gly-Pro-Hyp)_n. Geometries of the skeleton of these polymers were assumed to be the same as those of poly(l-proline) I (cis) and II (trans) and the calculations were performed only on infinite polymers in a single chain. The results show that the cis form is always more stable than the trans form for all the polymers treated. This energy difference between the cis and trans forms depends, for example, on the kind of amino acid residue, Gly or Ala, but this could not be seen in the Pro or Hyp residue. The flexibility or mobility of the collagen structure was explained using the energy difference between the cis and trans forms of the polymers, i.e. the cis-trans conversion of the collagen was discussed in connection with the energy difference. The reason why the collagen has the constitution of (Gly-Pro-Hyp)_n is briefly discussed.     

Cis–trans peptide-bond equilibria in cyclic hexapeptides

Data are presented on the position of the equilibria of cyclo(Xxx-Pro-Yyy)₂ backbones between forms with two cis Xxx-Pro peptide bonds and forms with only trans peptide bonds. These data are interpreted in terms of two factors: a solvent-independent steric interaction between the Xxx and Pro side chains, and the ability of solvent to influence the transannular electrostatic interaction between N? H and C?O groups of the Xxx units in the all-trans form.     

Structure of H/ACA RNP protein Nhp2p reveals cis/trans isomerization of a conserved proline at the RNA and Nop10 binding interface

H/ACA small nucleolar and Cajal body ribonucleoproteins (RNPs) function in site-specific pseudouridylation of eukaryotic rRNA and snRNA, rRNA processing, and vertebrate telomerase biogenesis. Nhp2, one of four essential protein components of eukaryotic H/ACA RNPs, forms a core trimer with the pseudouridylase Cbf5 and Nop10 that binds to H/ACA RNAs specifically. Crystal structures of archaeal H/ACA RNPs have revealed how the protein components interact with each other and with the H/ACA RNA. However, in place of Nhp2p, archaeal H/ACA RNPs contain L7Ae, which binds specifically to an RNA K-loop motif absent from eukaryotic H/ACA RNPs, while Nhp2 binds a broader range of RNA structures. We report solution NMR studies of Saccharomyces cerevisiae Nhp2 (Nhp2p), which reveal that Nhp2p exhibits two major conformations in solution due to cis/trans isomerization of the evolutionarily conserved Pro83. The equivalent proline is in the cis conformation in all reported structures of L7Ae and other homologous proteins. Nhp2p has the expected α-β-α fold, but the solution structures of the major conformation of Nhp2p with trans Pro83 and of Nhp2p-S82W with cis Pro83 reveal that Pro83 cis/trans isomerization affects the positions of numerous residues at the Nop10 and RNA binding interface. An S82W substitution, which stabilizes the cis conformation, also stabilizes the association of Nhp2p with H/ACA snoRNPs expressed in vivo. We propose that Pro83 plays a key role in the assembly of the eukaryotic H/ACA RNP, with the cis conformation locking in a stable Cbf5-Nop10-Nhp2 ternary complex and positioning the protein backbone to interact with the H/ACA RNA.     

Nmr studies of the molecular conformations in the linear oligopeptides H-(L-Ala)n-L-Pro-OH

The molecular conformations of the linear oligopeptides H-(L -Ala)_n-L -Pro-OH, with n = 1,2 and 3, have been investigated. ¹³C nmr observation of the equilibrium between the cis and trans forms of the Ala-Pro peptide bond indicated the occurrence of nonrandom conformations in solutions of these flexible peptides. The formation of the nonrandom species containing the cis form of the Ala-Pro bond was found to depend on the deprotonation of the carboxylic acid group of proline, the solvent, and the ionic strength in aqueous solution. The influence of intramolecular hydrogen bonding on the relative conformational energies of the species containing the cis and trans Ala-Pro peptide bond was studied by comparison of the peptides H-(Ala)_n-Pro-OH with analogous molecules where hydrogen bond formation was excluded by the covalent structure. In earlier work a hydrogen bond between the protonated terminal carboxylic acid group and the carbonyl oxygen of the penultimate amino acid residue had been suggested to stabilize conformations including trans proline. For the systems described here this hypothesis can be ruled out, since the cis:trans ratio is identical for molecules with methyl ester protected and free protonated terminal carboxylic acid groups of proline. Direct evidence for hydrogen bond formation between the deprotonated terminal carboxylic acid group and the amide proton of the penultimate amino acid residue in the molecular species containing cis proline was obtained from ¹H nmr studies. However, the cis:trans ratio of the Ala-Pro bond was not affected by N-methylation of the penultimate amino acid residue, which prevents formation of this hydrogen bond. Overall the experimental observations lead to the conclusion that the relative energies of the peptide conformations including cis or trans proline are mainly determined by intramolecular electrostatic interactions, whereas in the molecules considered, intramolecular hydrogen bonding is a consequence of specific peptide backbone conformations rather than a cause for the occurrence of energetically favored species. Independent support for this conclusion was obtained from model consideration which indicated that electrostatic interactions between the terminal carboxylic acid group and the carbonyl oxygen of the penultimate amino acid residue could indeed account for the observed relative conformational energies of the species containing cis and trans proline, respectively.     

A theoretical estimate of the energy barriers between stable conformations of the proline dimer.

Semi-empirical energy calculations for an internal Pro-Pro dimer are presented that take into account the nature of the flexibility of the proline ring due to its puckering. Calculations show that three stable conformations are available for the dimer: the cis (ω = 0°, ψ = 160°); the trans (ω = 180°, ψ = 160°, also referred to as trans′); and the cis′ (ω = 180°, ψ = ?40°) conformations. The best conformational pathways between these stable conformations are determined. Calculations also show that the barrier for cis′–trans′ conversion is of the same order of magnitude as that for cis–trans conversion.     

Equilibrium of the Cis-Trans Isomerisation of the Peptide Bond with N-alkyl Amino Acids Measured by 2D NMR

Summary The conformationalcis-trans equilibrium around the peptide bond in model tripeptides has been determined by 2D NMR methods (HOHAHA, ROESY). The study was limited to three different N-substituted amino acids in position 2, namely Pro (proline), Tic (1,2,3,4-tetrahydroisoquinoline-3-carboxylic acid), and N-MePhe (N-methylphenylalanine). In all cases the amino acid in position 1 was tyrosine and in position 3, phenylalanine. The results of our studies show that thecis-trans ratio depends mostly on the configuration of the amino acids forming the peptide bond undergoing thecis-trans isomerisation. The amino acid following the sequence (in position 3) does not have much influence on thecis-trans isomerisation, indicating that there is no interaction of the side chains between these amino acids. The model peptides with the L-Tyr-L-AA-(L-or D-)Phe (where AA is N-substituted amino acid) chiralities give 80–100% more of thecis form in comparison to the corresponding peptides with the D-Tyr-L-AA-(L-or D-)Phe chiralities. These results indicate that the incorporation of N-substituted amino acids in small peptides with the same chirality as the precedent amino acid involved in the peptide bound undergoing thecis/trans isomerisation moves the equilibrium to a significant amount of thecis form.     

A method to detect nonproline cis peptide bonds in proteins

Based on the geometrical parameters around seventeen incorrectly assigned trans conformations of peptide bonds in protein structures and their correct cis counterparts, we have devised an algorithm that is capable of detecting these sites. The algorithm was optimized to reliably find all of the seventeen test cases. It can be used to quickly scan an atomic coordinate file or the complete Brookhaven Protein Data Base for more likely candidates for non‐Pro cis peptide bonds. Also, it can be of help to guide the crystallographer in intermediate stages of structure determination towards suspect areas. © 1999 John Wiley & Sons, Inc. Biopoly 50: 536–544, 1999     

Proline isomerization in the C-terminal region of HSP27

In mammals, small heat-shock proteins (sHSPs) typically assemble into interconverting, polydisperse oligomers. The dynamic exchange of sHSP oligomers is regulated, at least in part, by molecular interactions between the α-crystallin domain and the C-terminal region (CTR). Here we report solution-state nuclear magnetic resonance (NMR) spectroscopy investigations of the conformation and dynamics of the disordered and flexible CTR of human HSP27, a systemically expressed sHSP. We observed multiple NMR signals for residues in the vicinity of proline 194, and we determined that, while all observed forms are highly disordered, the extra resonances arise from cis-trans peptidyl-prolyl isomerization about the G193-P194 peptide bond. The cis-P194 state is populated to near 15% at physiological temperatures, and, although both cis- and trans-P194 forms of the CTR are flexible and dynamic, both states show a residual but differing tendency to adopt β-strand conformations. In NMR spectra of an isolated CTR peptide, we observed similar evidence for isomerization involving proline 182, found within the IPI/V motif. Collectively, these data indicate a potential role for cis-trans proline isomerization in regulating the oligomerization of sHSPs.     

Prediction ofcis/transisomerization in proteins using PSI-BLAST profiles and secondary structure information

Background  

The majority of peptide bonds in proteins are found to occur in thetransconformation. However, for proline residues, a considerable fraction of Prolyl peptide bonds adopt thecisform. Prolinecis/transisomerization is known to play a critical role in protein folding, splicing, cell signaling and transmembrane active transport. Accurate prediction of prolinecis/transisomerization in proteins would have many important applications towards the understanding of protein structure and function.     

Conformational energy calculations on glycosylated turns in glycoproteins

Analysis of the amino acid sequence of glycoproteins has suggested the β-turn as a likely site of glycosylation in glycoproteins. According to this model, the peptide chain traverses the interior of a globular protein, reversing its direction at the protein surface, a likely point for the attachment of hydrophilic carbohydrate residues. In order to search for plausible conformations of glycosylated β-turns in asparagine-linked glycoproteins, we have adapted the conformational energy calculation method of Scheraga and coworkers for use in carbohydrates. The parameters for nonbonded and hydrogen-bonded interactions have been published, and electrostatic parameters are derived from a CNDO calculation on a model glycopeptide. Our results indicate that the orientation of the glycosyl amide bond having the amide proton nearly trans to the anomeric proton of the sugar has the lowest energy. Although CD and nmr experiments in our laboratory have consistently found this conformation, our calculations show the conformation having these two protons in a cis relationship to lie very close in energy. Calculations on the glycopeptide linkage model, α-N-acetyl, δ-N(2-acetamido-1,2-dideoxy-β-D -glucopyranosyl)-N′-methyl-L -asparaginyl amide show that several distinct geometries are allowed for glycosylated β-turns. For a type I β-turn, three conformations of the glycosylated side chain are found within 4 kcal of the minimum, while two conformations of the glycosylated side chain are allowed for a type II turn. The hydrogen-bonded C7 conformation is also allowed. Stereoviews of the low-energy conformations reveal no major hydrogen-bonding interaction between the peptide and sugar.     

Stereochemical aspects of the metabolism of the isomeric methylcyclohexanols and methylcyclohexanones

1. The seven isomeric optically inactive forms of methylcyclohexanol (i.e. 1-, and cis- and trans-2-, -3- and -4-) are excreted by rabbits mainly as glucuronides of the thermodynamically more stable forms of the alcohols. The eighth isomer, cyclohexylmethanol, however, undergoes aromatization in vivo, giving rise to benzoic acid and hippuric acid. The (±)-2-, (±-3- and 4-methylcyclohexanones are reduced in the rabbit and excreted mainly as the glucuronides of the thermodynamically more stable forms of the corresponding methylcyclohexanols. 2. Racemic cis- and trans-2-methylcyclohexanol and 2-methylcyclohexanone are all excreted as conjugated trans-2-methylcyclohexanol. However, when the (±)-cis-alcohol or the (±)-ketone is fed, (+)-trans-2-methylcyclohexanol is excreted, whereas when the (±)-trans-alcohol is fed it is excreted as the (±)-trans-alcohol. 3. Racemic cis- and trans-3-methylcyclohexanol and 3-methylcyclohexanone are all excreted as conjugated racemic cis-3-methylcyclohexanol. cis- and trans-4-Methylcyclohexanol and 4-methylcyclohexanone are all excreted as conjugated trans-4-methylcyclohexanol. 4. The metabolic differences between the various methylcyclohexanols are explicable in terms of their conformations and of Vennesland's (1958) hypothesis of the role of NADH in dehydrogenation reactions.     

Cis-trans isomerism of the peptide bond

The molecular orbital method PCILO is applied to eight. N-monsubstituted amides. Experimentally known geometric properties are reasonably predicted by minimization of total energy with respect to molecular geometry. The same procedure shows that molecular deformations during rotation around the peptide bond significantly lower calculated barriers. Experimental heats of activation and the free-energy changes associated with cis–trans isomerism are in good agreement with those calculated, which include qualitative estimates of configurational entropy contributions to the isomerism energies. Both the calculations and revised infrared data indicate that N-phenylurethane, which has been used as a model for the cis peptide bond, should be predominantly trans. However the variations in rotational barriers and cis–trans isomerism energies among the N-monosubstituted amides provide no reason to suppose that the cis peptide bond should be excluded from stable protein conformations.     

Conformational analysis of d-tubocurarine: Implications for minimal dimensions of its binding site within ion channels

All the minimum-energy conformations of d-tubocurarine were calculated by the method of molecular mechanics. The energy was minimized from 413 closed forms of the 18-member ring. The set of minimum-energy conformations includes 10 forms with energies less than 6 kcal/mol from the most stable one. Among the four lowest minimum-energy conformations, two forms correspond to those known from X-ray studies, whereas two conformations were not detected experimentally earlier. The flexibility of d-tubocurarine was estimated by calculating six paths of interconversion between the four lowest minimum-energy conformations. Using a molecular graphics technique, it was found that the most extended minimum-energy conformation of d-tubocurarine may fit in an ion channel of a rectangular profile of 8.7 × 11.2 Å, while one tetrahydroisoquinoline head may fit a profile as small as 6.9 × 11.0 Å. A possible model of d-tubocurarine location within the ion channel of the neuronal nicotinic acetylcholine receptor is suggested.