The influence of helix morphology on co-operative polyamide backbone conformational flexibility in peptide nucleic acid complexes.

A systematic analysis of peptide nucleic acid (PNA) complexes deposited in the Protein Data Bank has been carried out using a set of contiguous atom torsion angle definitions. The analysis is complemented by molecular mechanics adiabatic potential energy calculations on hybrid PNA-nucleic acid model systems. Hitherto unobserved correlations in the values of the (alpha and epsilon) dihedral angles flanking the backbone secondary amide bond are found. This dihedral coupling forms the basis of a PNA backbone conformation classification scheme. Six conformations are thus characterised in experimental structures. Helix morphology is found to exert a significant influence on backbone conformation and flexibility: Watson-Crick PNA strands in complexes with DNA and RNA, that possess A-like base-pair stacking, adopt backbone conformations distinct from those in PNA.DNA-PNA triplex and PNA-PNA duplex P-helix forms. Solvation effects on Watson-Crick PNA backbone conformation in heterotriplexes are discussed and the possible involvement of inter-conformational transitions and dihedral angle uncoupling in asymmetric heteroduplex base-pair breathing is suggested.     

On the stability of peptide nucleic acid duplexes in the presence of organic solvents

Nucleic acid double helices are stabilized by hydrogen bonding and stacking forces (a combination of hydrophobic, dispersive and electrostatic forces) of the base pairs in the helix. One would predict the hydrogen bonding contributions to increase and the stacking contributions to decrease as the water activity in the medium decreases. Study of nucleobase paired duplexes in the absence of water and ultimately in pure aprotic, non-polar organic solvents is not possible with natural phosphodiester nucleic acids due to the ionic phosphate groups and the associated cations, but could be possible with non-ionic nucleic acid analogues or mimics such as peptide nucleic acids. We now report that peptide nucleic acid (PNA) (in contrast to DNA) duplexes show almost unaffected stability in up to 70% dimethylformamide (DMF) or dioxane, and extrapolation of the data to conditions of 100% organic solvents indicates only minor (or no) destabilization of the PNA duplexes. Our data indicate that stacking forces contribute little if at all to the duplex stability under these conditions. The differences in behaviour between the PNA and the DNA duplexes are attributed to the differences in hydration and counter ion release rather than to the differences in nucleobase interaction. These results support the possibility of having stable nucleobase paired double helices in organic solvents.     

Changes in interhelical hydrogen bonding upon rhodopsin activation

Hydrogen bonding interactions between transmembrane helices stabilize the visual pigment rhodopsin in an inactive conformation in the dark. The crystal structure of rhodopsin has previously revealed that Glu122 and Trp126 on transmembrane helix H3 form a complex hydrogen bonding network with Tyr206 and His211 on H5, while the indole nitrogen of Trp265 on H6 forms a water-mediated hydrogen bond with Asn302 on H7. Here, we use solid-state magic angle spinning NMR spectroscopy to probe the changes in hydrogen bonding upon rhodopsin activation. The NMR chemical shifts of 15N-labeled tryptophan are consistent with the indole nitrogens of Trp126 and Trp265 becoming more weakly hydrogen bonded between rhodopsin and metarhodopsin II. The NMR chemical shifts of 15N-labeled histidine show that His211 is neutral; the unprotonated imidazole nitrogen is not coordinated to zinc in rhodopsin and becomes more strongly hydrogen bonded in metarhodopsin II. Moreover, measurements of rhodopsin containing 13C-labeled histidine show that a strong hydrogen bond between the side-chain of Glu122 and the backbone carbonyl of His211 is disrupted in metarhodopsin II. The implications of these observations for the activation mechanism of rhodopsin are discussed.     

Probing multiple effects on 15N, 13C alpha, 13C beta,and 13C' chemical shifts in peptides using density functional theory

We have used density functional calculations on model peptides to study conformational effects on (15)N, (13)C alpha, (13)C beta, and (13)C' chemical shifts, associated with hydrogen bonding, backbone conformation, and side-chain orientation. The results show a significant dependence on the backbone torsion angles of the nearest three residues. Contributions to (15)N chemical shifts from hydrogen bonding (up to 8 ppm), backbone conformation (up to 13 ppm), side-chain orientation and neighborhood residue effects (up to 22 ppm) are significant, and a unified theory will be required to account for their behavior in proteins. In contrast to this, the dependence on sequence and hydrogen bonding is much less for (13)C alpha and (13)C beta chemical shifts (<0.5 ppm), and moderate for carbonyl carbon shifts (<2 ppm). The effects of side-chain orientation are mainly limited to the residue itself for both nitrogen and carbon, but the chi(1) effect is also significant for the nitrogen shift of the following residue and for the (13)C' shift of the preceding residue. The calculated results are used, in conjunction with an additive model of chemical shift contributions, to create an algorithm for prediction of (15)N and (13)C shifts in proteins from their structure; this includes a model to extrapolate results to regions of torsion angle space that have not been explicitly studied by density functional theory (DFT) calculations. Crystal structures of 20 proteins with measured shifts have been used to test the prediction scheme. Root mean square deviations between calculated and experimental shifts 2.71, 1.22, 1.31, and 1.28 ppm for N, C alpha, C beta, and C', respectively. This prediction algorithm should be helpful in NMR assignment, crystal and solution structure comparison, and structure refinement.     

Sequence and structure patterns in proteins from an analysis of the shortest helices: implications for helix nucleation

The shortest helices (three-length 3(10) and four-length alpha), most abundant among helices of different lengths, have been analyzed from a database of protein structures. A characteristic feature of three-length 3(10)-helices is the shifted backbone conformation for the C-terminal residue (phi,psi angles: -95 degrees,0 degrees ), compared to the rest of the helix (-62 degrees,-24 degrees ). The deviation can be attributed to the release of electrostatic repulsion between the carbonyl oxygen atoms at the two C-terminal residues and further stabilization (due to a more linear geometry) of an intrahelical hydrogen bond. A consequence of this non-canonical C-terminal backbone conformation can be a potential origin of helix kinks when a 3(10)-helix is sequence-contiguous at the alpha-helix N-terminal. An analysis of hydrogen bonding, as well as hydrophobic interactions in the shortest helices shows that capping interactions, some of them not observed for longer helices, dominate at the N termini. Further, consideration of the distribution of amino acid residues indicates that the shortest helices resemble the N-terminal end of alpha-helices rather than the C terminus, implying that the folding of helices may be initiated at the N-terminal end, which does not get propagated in the case of the shortest helices. Finally, pairwise comparison of beta-turns and the shortest helices, based on correlation matrices of site-specific amino acid composition, and the relative abundance of these short secondary structural elements, leads to a helix nucleation scheme that considers the formation of an isolated beta-turn (and not an alpha-turn) as the helix nucleation step, with shortest 3(10)-helices as intermediates between the shortest alpha-helix and the beta-turn. Our results ascribe an important role played by shortest 3(10)-helices in proteins with important structural and folding implications.     

Multiple conformations of the acylenzyme formed in the hydrolysis of methicillin by Citrobacter freundii beta-lactamase: a time-resolved FTIR spectroscopic study

Time-resolved infrared difference spectroscopy has been used to show that the carbonyl group of the acylenzyme reaction intermediate in the Citrobacter freundii beta-lactamase-catalyzed hydrolysis of methicillin can assume at least four conformations. A single-turnover experiment shows that all four conformations decline during deacylation with essentially the same rate constant. The conformers are thus in exchange on the reaction time scale, assuming that deacylation takes place only from the conformation which is most strongly hydrogen bonded or from a more minor species not visible in these experiments. All conformers have the same (10 cm-1) narrow bandwidth compared with a model ethyl ester in deuterium oxide (37 cm-1) which shows that all conformers are well ordered relative to free solution. The polarity of the carbonyl group environment in the conformers varies from 'ether-like' to strongly hydrogen bonding (20 kJ/mol), presumably in the oxyanion hole of the enzyme. From the absorption intensities, it is estimated that the conformers are populated approximately proportional to the hydrogen bonding strength at the carbonyl oxygen. A change in the difference spectrum at 1628 cm-1 consistent with a perturbation (relaxation) of protein beta-sheet occurs slightly faster than deacylation. Consideration of chemical model reactions strongly suggests that neither enamine nor imine formation in the acyl group is a plausible explanation of the change seen at 1628 cm-1. A turnover reaction supports the above conclusions and shows that the conformational relaxation occurs as the substrate is exhausted and the acylenzymes decline. The observation of multiple conformers is discussed in relation to the poor specificity of methicillin as a substrate of this beta-lactamase and in terms of X-ray crystallographic structures of acylenzymes where multiple forms are not apparently observed (or modeled). Infrared spectroscopy has shown itself to be a useful method for assessment of the uniqueness of enzyme-substrate interactions in physiological turnover conditions as well as for determination of ordering, hydrogen bonding, and protein perturbation.     

Molecular structures for microbial polysaccharides. X-ray diffraction from the capsular polysaccharide of Klebsiella K-type 5

Fibre Type X-ray diffraction patterns have been obtained from oriented, semicrystalline films prepared from the sodium salt form of the capsular polysaccharide of K5. The molecule has a linear trisaccharide repeating sequence containing a 1,3 linked β-d-glucuronic acid and a 1,4 linked β-d-glucose residue, resulting in a backbone linkage geometry of Man(1 eq-4e1)-GIcUA (1-eg-4eg) Gle (1 eq-3 eq) Man. It also contains an O-acetyl group and two charged groups, namely a uronic acid and a pyruvate. Analysis of the diffraction results gives rise to an extended two-fold helical conformation with an axially projected advance of 1.35 nm which correlates directly with the covalent repeating sequence. X-ray diffraction patterns from preparations of deacetylated K5 polysaccharide showed similar conformations for the individual helices but the interchain packing arrangements are different. In each case, isolated helices have been computer generated using molecular model building procedures and the most favourable conformations for both preparations were those which contained three stabilizing interresidue hydrogen bonds, one across each of the glycosidic linkages.     

Information transfer from DNA to peptide nucleic acids by template-directed syntheses.

Peptide nucleic acids (PNAs) are analogs of nucleic acids in which the ribose-phosphate backbone is replaced by a backbone held together by amide bonds. PNAs are interesting as models of alternative genetic systems because they form potentially informational base paired helical structures. Oligocytidylates have been shown to act as templates for formation of longer oligomers of G from PNA G2 dimers. In this paper we show that information can be transferred from DNA to PNA. DNA C4T2C4 is an efficient template for synthesis of PNA G4A2G4 using G2 and A2 units as substrates. The corresponding synthesis of PNA G4C2G4 on DNA C4G2C4 is less efficient. Incorporation of PNA T2 into PNA products on DNA C4A2C4 is the least efficient of the three reactions. These results, obtained using PNA dimers as substrates, parallel those obtained using monomeric activated nucleotides.     

Staggered molecular packing in crystals of a collagen-like peptide with a single charged pair

The crystal structure of the triple-helical peptide, (Pro-Hyp-Gly)(4)-Glu-Lys-Gly-(Pro-Hyp-Gly)(5) has been determined to 1.75 A resolution. This peptide was designed to examine the effect of a pair of adjacent, oppositely charged residues on collagen triple-helical conformation and intermolecular interactions. The molecular conformation (a 7(5) triple helix) and hydrogen bonding schemes are similar to those previously reported for collagen triple helices and provides a second instance of water mediated N--H . . . O==C interchain hydrogen bonds for the amide group of the residue following Gly. Although stereochemically capable of forming intramolecular or intermolecular ion pairs, the lysine and glutamic acid side-chains instead display direct interactions with carbonyl groups and hydroxyproline hydroxyl groups or interactions mediated by water molecules. Solution studies on the EKG peptide indicate stabilization at neutral pH values, where both Glu and Lys are ionized, but suggest that this occurs because of the effects of ionization on the individual residues, rather than ion pair formation. The EKG structure suggests a molecular mechanism for such stabilization through indirect hydrogen bonding. The molecular packing in the crystal includes an axial stagger between molecules, reminiscent of that observed in D-periodic collagen fibrils. The presence of a Glu-Lys-Gly triplet in the middle of the sequence appears to mediate this staggered molecular packing through its indirect water-mediated interactions with backbone C==O groups and side chains.     

Two-dimensional 1H NMR of two chemically modified analogs of the basic pancreatic trypsin inhibitor. Sequence-specific resonance assignments and sequence location of conformation changes relative to the native protein

Two-dimensional nuclear magnetic resonance was used to obtain sequence specific assignments for the 1H NMR spectra of two chemically modified analogs of the basic pancreatic trypsin inhibitor. In one analog the disulfide bond 14-38 was cleaved, in the second derivative the N-terminus was transaminated. From measurements of the chemical shifts and determination of the sequence locations of slowly exchanging backbone amide protons it was found that conformational differences between the native inhibitor and the chemical modifications occur exclusively near the modification sites and that the internal hydrogen bonds are nearly fully preserved. Intriguing conformation differences with respect to the native protein are that for five residues in the transaminated inhibitor and for one residue in the reduced inhibitor multiple local conformers are indicated, and that the four internal water molecules observed in the crystal structure of the native inhibitor appear not to be preserved after reduction of the disulfide bond 14-38.     

On the dependence of molecular conformation on the type of solvent environment: a molecular dynamics study of cyclosporin A

The dependence of the conformation of cyclosporin A (CPA), a cyclic undecapeptide with potent immunosuppressive activity, on the type of solvent environment is examined using the computer simulation method of molecular dynamics (MD). Conformational and dynamic properties of CPA in aqueous solution are obtained from MD simulations of a CPA molecule dissolved in a box with water molecules. Corresponding properties of CPA in apolar solution are obtained from MD simulations of CPA in a box with carbontetrachloride. The results of these simulations in H2O and in CCl4 are compared to each other and to those of previous simulations of crystalline CPA and of an isolated CPA molecule. The conformation of the backbone of the cyclic polypeptide is basically independent of the type of solvent. In aqueous solution the beta-pleated sheet is slightly weaker and the gamma-turn is a bit less pronounced than in apolar solution. Side chains may adopt different conformations in different solvents. In apolar solution the hydrophobic side chain of the MeBmt residue is in an extended conformation with its hydroxyl group hydrogen bonded to the backbone carbonyl group. In aqueous solution this hydrophobic side chain folds over the core of the molecule and the mentioned hydrogen bond is broken in favor of hydrogen bonding to water molecules. The conformation obtained from the MD simulation in CCl4 nicely agrees with experimental atom-atom distance data as obtained from nmr experiments in chloroform. In aqueous solution the relaxation of atomic motion tends to be slower than in apolar solution.     

Molecular modelling of the xanthan chain conformations

Determination of the steric features displayed by the bacterial polysaccharide xanthan is reported. Molecular modelling was conducted assuming an isolated, single polymer chain, having a well defined repeating unit. Potential energy was computed by taking into account van der Waals, torsional, exo-anomeric, and hydrogen bond contributions. Several families of stable conformers are found for the polymeric skeleton. Their occurences depend upon the conformations at the branching point, which are coupled to the orientations of the primary acetate group of the 6-O-acetate mannopyranosyl residue. The calculated stable conformations of the repeating unit would generate either hollow or extended helices, which, with the exception of one have left-handed chirality. One of these helices appears to correspond to the ordered conformation of xanthan in the solid state. The side chains exhibit a wide range of accessible conformations, but do not show any preferred tendency to fold back onto the polymeric backbone.     

Sequence determinants of the capping box, a stabilizing motif at the N-termini of alpha-helices.

The capping box, a recurrent hydrogen bonded motif at the N-termini of alpha-helices, caps 2 of the initial 4 backbone amide hydrogen donors of the helix (Harper ET, Rose GD, 1993, Biochemistry 32:7605-7609). In detail, the side chain of the first helical residue forms a hydrogen bond with the backbone of the fourth helical residue and, reciprocally, the side chain of the fourth residue forms a hydrogen bond with the backbone of the first residue. We now enlarge the earlier definition of this motif to include an accompanying hydrophobic interaction between residues that bracket the capping box sequence on either side. The expanded box motif--in which 2 hydrogen bonds and a hydrophobic interaction are localized within 6 consecutive residues--resembles a glycine-based capping motif found at helix C-termini (Aurora R, Srinivasan R, Rose GD, 1994, Science 264:1126-1130).     

Change in entropy of the free polypeptide chain during formation of hydrogen bonds]

Formation probabilities of different hydrogen bonds between carbonyl oxygen and amide hydrogen were determined by Monte Carlo simulations using a computer model in the space of sterically allowable conformations of alanine and glycine oligopeptides, and the corresponding entropy losses for the peptide backbone, T delta S, were calculated. The model was studied at different criteria of steric interactions. Comparison with the data of other authors showed the values of T delta S to be mainly determined by overall extent and type of the state space and to be only slightly dependent on its energy profile. Both short-range and long-range steric interactions were shown to prevent hydrogen bonding, especially in alanine peptides. In the model studied, the initiation of alpha(R)-helices is associated with T delta S = 8-10 kT, and prior formation of a 3/10-turn or one three-center H-bond does not appreciably decrease this entropy barrier. Elongation of the alpha(R)-helix by one residue leads to T delta S = 3.0-3.7 kT, the helices begin to stabilize after at least three sequential H-bonds are formed. The difference in the probability of insertion of Ala and Gly into the helix is lower than it follows from comparison of their mobility. The results could be explained assuming that factors different from helical H-bonds take part in the stabilization of the helices. One may suppose upon modeling of folding that even three sequential H-bonds are unable to fix the structure of a flexible peptide loop, while the elongation of alpha(R)-helices in the supersecondary helix-loop-helix structure is favorable as long as the loop conformation remains nearly optimal.     

The Ribbon of Hydrogen Bonds in the Three-Dimensional Structure of Globular Proteins

We report quantum mechanical computations and experimental evidence which suggest that the backbone conformation of globular proteins depends generally on the conservation of that part of the hydrogen bond network or ribbon which is joined, in general, directly to the backbone and is largely independent of the remainder of this whole network of hydrogen bonds. The familiar hydrogen bonds of the helix and the sheet form about one-half of this ribbon of hydrogen bonds. Both water molecules and hydrogen bonding side chain groups are involved in the formation of the ribbon.This view of the three-dimensional structure of globular proteins in terms of the `molecule' allows us to deal with the non-secondary structure as well as with the familiar secondary structure. It also suggests that the ribbon contains approximately the same number of hydrogen bonds within all three structures – the helix, the sheet and the coil – and that this is the reason for the ease of interconversion of these three structures.The quantum mechanical computations on hydrogen bonding suggest that delocalised water molecules which have substantial mobility are an essential part of the ribbon. This situation arises because the hydrogen bonding groups of the protein molecule are not free to move to optimise the hydrogen bonding geometries as are the oxygen atoms in the waters and ices. Such delocalised water molecules either have high B values or are invisible in the X-ray data and yet are able to form a structure which is as strong as a normal hydrogen bond.The experimental data on the point mutations of the THRI57 residue of the T4 phage lysome provides an initial test of this model. Both the local backbone conformation and the ribbon of hydrogen bonds are conserved throughout all the mutations of residue 157,providing that the delocalised water molecules are accepted as a genuine part of the structure. These mutations include the introduction of hydrocarbon side chains at position 157 when water molecules or other side chain groups take over the formation of the hydrogen bonds.We suggest that, provided steric effects are not important, many point mutations succeed because they leave the ribbon of hydrogen bonds (and so the backbone conformation) largely unchanged.