Rigid and flexible region in lysozyme and the invariant features in its hydration shell

Water-mediated transformations provide a useful handle for exploring the flexibility in protein molecules and the invariant features in their hydration shells. Low-humidity monoclinic hen egg white lysozyme, resulting from such a transformation, has perhaps the lowest solvent content observed in any protein crystal so far and has a well-ordered structure. A detailed comparison involving this structure, low-humidity tetragonal lysozyme, and the other available refined crystal structures of the enzyme permits the delineation of the relatively rigid, moderately flexible and highly flexible regions of the molecule. The relatively rigid region forms a contiguous structural unit close to the molecular centroid and encompasses parts of of the main beta-structure and three alpha-helices. The hydration shell of the protein contains 30 invariant water molecules. Many of them are involved in holding different parts of the molecule together or in stabilizing local structure. Five of the six invariant water molecules attached to the substrate-binding region form part of a water cluster contiguous with the side-chains of the catalytic residues Glu-35 and Asp-52.

Visualisation of a 2′-5′ parallel stranded double helix at atomic resolution: crystal structure of cytidylyl-2′,5′-adenosine

X-ray crystallographlc studies on 3′–5′ ollgomers have provided a great deal of information on the stereochemistry and conformational flexibility of nucleic acids and polynucleotides. In contrast, there is very little Information available on 2′–5′ polynucleotides. We have now obtained the crystal structure of Cytidylyl-2′,5′-Adenoslne (C2′p5′A) at atomic resolution to establish the conformational differences between these two classes of polymers. The dlnucleoside phosphate crystallises in the monocllnlc space group C2, with a = 33.912(4)Å, b =16.824(4)Å, c = 12.898(2)Å and 0 = 112.35(1) with two molecules in the asymmetric unit. Spectacularly, the two independent C2′p5′A molecules in the asymmetric unit form right handed miniature parallel stranded double helices with their respective crystallographic two fold (b axis) symmetry mates. Remarkably, the two mini duplexes are almost indistinguishable. The cytosines and adenines form self-pairs with three and two hydrogen bonds respectively. The conformation of the C and A residues about the glycosyl bond is anti same as in the 3′–5′ analog but contrasts the anti and syn geometry of C and A residues in A2′p5′C. The furanose ring conformation is C3′endo, C2′endo mixed puckering as in the C3′p5′A-proflavine complex. A comparison of the backbone torsion angles with other 2′–5′ dinucleoside structures reveals that the major deviations occur in the torsion angles about the C3′–C2′ and C4′-C3′ bonds. A right-handed 2′–5′ parallel stranded double helix having eight base pairs per turn and 45° turn angle between them has been constructed using this dinucleoside phosphate as repeat unit. A discussion on 2′–5′ parallel stranded double helix and its relevance to biological systems is presented.

Protein Hydration and Water Structure: X-ray Analysis of a Closely Packed Protein Crystal with Very Low Solvent Content

Low-humidity monoclinic lysozyme, resulting from a water-mediated transformation, has one of the lowest solvent contents (22% by volume) observed in a protein crystal. Its structure has been solved by the molecular replacement method and refined to an R value of 0.175 for 7684 observed reflections in the 10–1.75 Å resolution shell. 90% of the solvent in the well ordered crystals could be located. Favourable sites of hydration on the protein surface include side chains with multiple hydrogen-bonding centres, and regions between short hydrophilic side chains and the main-chain CO or NH groups of the same or nearby residues. Major secondary structural features are not disrupted by hydration. However, the free CO groups at the C terminii and, to a lesser extent, the NH groups at the N terminii of helices provide favourable sites for water interactions, as do reverse turns and regions which connect β-structure and helices. The hydration shell consists of discontinuous networks of water molecules, the maximum number of molecules in a network being ten. The substrate-binding cleft is heavily hydrated, as is the main loop region which is stabilized by water interactions. The protein molecules are close packed in the crystals with a molecular coordination number of 14. Arginyl residues are extensively involved in intermolecular hydrogen bonds and water bridges. The water molecules in the crystal are organized into discrete clusters. A distinctive feature of the clusters is the frequent occurrence of three-membered rings. The protein molecules undergo substantial rearrangement during the transformation from the native to the low-humidity form. The main-chain conformations in the two forms are nearly the same, but differences exist in the side-chain conformation. The differences are particularly pronounced in relation to Trp 62 and Trp 63. The shift in Trp 62 is especially interesting as it is also known to move during inhibitor binding.

HELANAL: A program to characterise helix geometry in proteins,

A detailed analysis of structural and position dependent characteristic features of helices will give a better understanding of the secondary structure formation in globular proteins. Here we describe an algorithm that quantifies the geometry of helices in proteins on the basis of their C-alpha atoms alone. The Fortran program HELANAL can extract the helices from the PDB files and then characterises the overall geometry of each helix as being linear, curved or kinked, in terms of its local structural features, viz. local helical twist and rise, virtual torsion angle, local helix origins and bending angles between successive local helix axes. Even helices with large radius of curvature are unambiguously identified as being linear or curved. The program can also be used to differentiate a kinked helix and other motifs, such as helix-loop-helix or a helix-turn-helix (with a single residue linker) with the help of local bending angles. In addition to these, the program can also be used to characterise the helix start and end as well as other types of secondary structures.

Analysis of sequence dependent variations in secondary and tertiary structure of tRNA molecules

The double helical regions of the five tRNA(Phe) and two tRNA(Asp) crystal structures have been analyzed using the local basepair step parameters. The sequence dependent effects in the mini double helices of tRNA are very similar to those observed in the crystal structures of oligonucleotides in the A-form, the purine-pyrimidine and purine-purine steps have small roll angles when compared to the fiber models of A-DNA as well as A-RNA, while the pyrimidine-purine doublet steps have large roll angles. The orientation of the basepairs in the D-stem is unusual and invariant i.e. they are different from the other three stems but are very similar in all the five tRNA(Phe) crystal structures, presumably due to tertiary interaction of the Watson-Crick basepairs with other bases, with all bases being highly conserved. The origin of the differences between the tertiary structures of tRNA(Phe) and tRNA(Asp) from yeast has also been investigated. It is found that even though the angle between the acceptor arm and the D-stem is very similar in the two structures, the angle subtended by the acceptor arm and the anticodon arm is smaller in the tRNA(Phe) structure (by more than 10 degrees). This is due to differences in the orientation of the two mini helices constituting the anticodon arm, which are inclined to each other by approximately 25 degrees in tRNA(Phe) and 16 degrees in tRNA(Asp). In addition, the acceptor arm, the D-stem and the anticodon stem are nearly coplanar in tRNA(Phe), while in tRNA(Asp) the anticodon stem projects out of the plane defined by the acceptor arm and the anticodon stem. These two features together lead to a larger separation between the acceptor and anticodon ends in tRNA(Asp) and indicate that the junction between the D-stem and the anticodon stem is quite variable, with features characteristic of a ball-and-socket type joint and determined for each tRNA molecule by the base sequence at the junction.

Double helix conformation, groove dimensions and ligand binding potential of a G/C stretch in B-DNA

The self-complementary DNA fragment CCGGCGCCGG crystallizes in the rhombohedral space group R3 with unit cell parameters a = 54.07 angstrom and c = 44.59 angstrom. The structure has been determined by X-ray diffraction methods at 2.2 angstrom resolution and refined to an R value of 16.7%. In the crystal, the decamer forms B-DNA double helices with characteristic groove dimensions: compared with B-DNA of random sequence, the minor groove is wide and deep and the major groove is rather shallow. Local base pair geometries and stacking patterns are within the range commonly observed in B-DNA crystal structures. The duplex bears no resemblance to A-form DNA as might have been expected for a sequence with only GC base pairs. The shallow major groove permits an unusual crystal packing pattern with several direct intermolecular hydrogen bonds between phosphate oxygens and cytosine amino groups. In addition, decameric duplexes form quasi-infinite double helices in the crystal by end-to-end stacking. The groove geometries and accessibilities of this molecule as observed in the crystal may be important for the mode of binding of both proteins and drug molecules to G/C stretches in DNA.

Definitions and nomenclature of nucleic acid structure parameters

At an EMBO Workshop on DNA Curvature and Bending, held at Churchill College, Cambridge, on 10-15 September 1988, two sessions were scheduled on definitions of parameters used to describe the geometry of nucleic acid chains and helices, and a common nomenclature for these parameters. The most widely used library of helix analysis programs, HELIB (Fratini et al., 1982; Dickerson, 1985) suffers from the fact that the translations and rotations as defined are not fully independent and depend to a certain extent upon the choice of overall helix axis. Several research groups have been engaged independently in developing alternative programs for the geometrical analysis of polynucleotide chains, but with different definitions of quantities calculated and with widely different nomenclature even when the same parameter was involved.

Effect of constraints by threonine on proline containing αa-helix—A molecular dynamics approach

Proline plays an important role in the secondary structure of proteins. In the pursuit of understanding its structural role, Proline containing helices with constraints have been studied by employing molecular dynamics (MD) technique. In the present study, the constraint introduced is a threonine residue, whose sidechain has intramolecular hydrogen bond interaction with the backbone oxygen atom. The three systems that have been chosen for characterization are: (1) Ace-(Ala)12−Thr-Pro-(Ala)10−NHMe, (2) Ace-(Ala)13-Pro-Ala-Thr- (Ala)8-NHMe and (3) Ace-(Ala)13-Pro-(Ala)3-Thr-(Ala)6-NHMe. The equilibrium structures and structural transitions have been identified by monitoring the backbone dihedral angles, bend related parameters and the hydrogen bond interactions. The MD averages and root mean square (r.m.s.) fluctuations are compared and discussed. Energy minimization has been carried out on selected MD simulated points in order to analyze the characteristics of different conformations.

Maximum-entropy calculation of the electron density at 4 A resolution of Pf1 filamentous bacteriophage

A 4 A electron-density map of Pf1 filamentous bacterial virus has been calculated from x-ray fiber diffraction data by using the maximum-entropy method. This method produces a map that is free of features due to noise in the data and enables incomplete isomorphous-derivative phase information to be supplemented by information about the nature of the solution. The map shows gently curved (banana-shaped) rods of density about 70 A long, oriented roughly parallel to the virion axis but slewing by about 1/6th turn while running from a radius of 28 A to one of 13 A. Within these rods, there is a helical periodicity with a pitch of 5 to 6 A. We interpret these rods to be the helical subunits of the virion. The position of strongly diffracted intensity on the x-ray fiber pattern shows that the basic helix of the virion is right handed and that neighboring nearly parallel protein helices cross one another in an unusual negative sense.

Hybrid Polypeptides: Gabapentin as a Stereochemically Constrained gamma-Amino Acid Residue

The design of folded structures in peptides containing the higher homologues of alpha-amino acid residues requires the restriction of the range of local conformational choices In alpha-amino acids stereochemically constrained residues like alpha,alpha-dialkylated residue, aminoisobutyric acid (Aib), and D-Proline ((D)Pro) have proved extremely useful in the design of helices and hairpins in short peptides Extending this approach, backbone substitution and cyclization are anticipated to bc useful in generating conformationally constrained beta- and gamma-residues This brief review provides a survey of work on hybrid peptide sequences concerning the conformationally constrained gamma-amino acid residue 1-aminomethyl cyclohexane acetic acid, gabapentin (Gpn) This achiral, beta,beta-disubstituted, gamma-residue strongly favors gauche-gauche conformations about the C-alpha-C-beta (0(2)) and C-alpha-C-gamma (0(1)) bonds, facilitating local folding The Gpn residue can adopt both C-7 (NH1 -> CO1) and C-9 (CO1 (I)<- NH1+I) hydrogen bonds which are analogous to the C-5 and C7 (gamma-turn) conformations at alpha-residues In conjunction with adjacent residues, Gpn may be used in ay and gamma alpha segments to generate C-12 hydrogen bonded conformations which may be considered as expanded analogs of conventional beta-turns The structural characterization of C-12 helices, C-12/C-10 helices with mixed hydrogen bond directionalities and beta-hairpins incorporating Gpn residues at the turn segment is illustrated (C) 2010 Wiley Periodicals, Inc Biopolymers (Pept Sci) 94 733-741 2010

The Design and Construction of Synthetic Protein Mimics

A strategy for the modular construction of synthetic protein mimics based on the ability non-protein amino acids to act as stereochemical directors of polypeptide chain folding, is described. The use of alpha-aminoisobutyric acid (Aib) to construct stereochemically rigid helices has been exemplified by crystallographic and spectroscopic studies of several apolar peptides, ranging in length from seven to sixteen residues. The problem of linker design in elaborating alpha,alpha motifs has been considered. Analysis of protein crystal structure data provides a guide to choosing linking sequences. Attempts at constructing linked helical motifs using linking Gly-Pro segments have been described. The use of flexible linkers, like epsilon-aminocaproic acid has been examined and the crystallographic and solution state analysis of a linked helix motif has been presented. The use of bulky sidechain modifications on a helical scaffold, as a means of generating putative binding sites has been exemplified by a crystal structure of a peptide packed in a parallel zipper arrangement.

Conformation of polymerizable diacetylenic lecithin, DC8,9PC, by 1H-NMR spectroscopy

Diacetylenic phospholipid, 1,2 bis-(tricosa-10,12-diynoyl)-sn-glycero-3-phosphocholine (DC8,9PC), forms helices and tubules in addition to liposomes. The diacetylenic moiety responsible for the transformation is probed by 2-D NMR correlated spectroscopy. Chemical shift assignments and the analysis of 2D-COSY measurements were done on the lipid in chloroform-d solution. Based on this analysis, a model for the lipid is proposed. The geometry of the headgroup, glycerol backbone and acyl chains up to three methylenes from glycerol backbone [-(CH2)(3)-] is similar to that of dipalmitoyl phosphatidylcholine. The estimated torsional angle for methylene groups adjacent to diacetylenic moieties suggested an overall tilt of the diacetylenic lipid molecule from the bilayer axis of 25-30 degrees. This tilt could be negative or positive depending on the handedness of the resultant microstructures.

Effect of the valine-threonine constraint on the dynamics of the proline helix — A molecular dynamics study

Proline residues in helices play an important role in the structure of proteins. The proline residue introduces a kink in the helix which varies from about 5-degrees to 50-degrees. The presence of other residues such as threonine or valine near the proline region can influence the flexibility exhibited by the kinked helix, which can have an important biological role. In the present paper, the constraint introduced by threonine and valine on a proline helix is investigated by molecular dynamics studies. The systems considered am (1) a poly-alanine helix with threonine-proline residues (TP) and (2) a poly-alanine helix with valine-threonine-proline residues (VTP), in the middle. Molecular dynamics simulations are carried out on these two systems for 500 ps. The results are analyzed in terms of structural transitions, bend-related parameters and sidechain orientations.

A Change in the 3(10)- to alpha-Helical Transition Point in the Heptapeptides Containing Sulfur and Selenium

Crystal structures of three heptapeptides Boc-Ala-Leu-Aib-XXX-Ala-Leu-Aib-OMe (where XXX = methionine in peptide A, selenomethionine in peptide B, and S-benzyl cysteine in peptide C) reveal mixed 3(10)-/alpha-helical conformations with R factors of 6.94, 5.79, and 5.98, respectively. All the structures were solved in the P2(1)2(1)2(1) space group. 3(10)- to a-helical transitions are observed in all of these peptides. The helices begin as a 3(10)-helical segment at the N-terminus and then transit for peptides A and C at residue Aib(3) carbonyl (O(3)), while for peptide B the transition occurs at residue Leu(2) carbonyl oxygen (O(2)). There are water molecules associated in the crystal of each of these peptides and they form different types of hydrogen bonding patterns in each crystal. The observations suggest that 3(10)- to alpha-helical transition is sequence dependent in these short heptapeptide sequences.

Free Energy Barriers for Escape of Water Molecules from Protein Hydration Layer

Free energy barriers separating interfacial water molecules from the hydration layer at the surface of a protein to the bulk are obtained by using the umbrella sampling method of free energy calculation. We consider hydration layer of chicken villin head piece (HP-36) which has been studied extensively by molecular dynamics simulations. The free energy calculations reveal a strong sensitivity to the secondary structure. In particular, we find a region near the junction of first and second helix that contains a cluster of water molecules which are slow in motion, characterized by long residence times (of the order of 100 ps or more) and separated by a large free energy barrier from the bulk water. However, these ``slow'' water molecules constitute only about 5-10% of the total number of hydration layer water molecules. Nevertheless, they play an important role in stabilizing the protein conformation. Water molecules near the third helix (which is the important helix for biological function) are enthalpically least stable and exhibit the fastest dynamics. Interestingly, barrier height distributions of interfacial water are quite broad for water surrounding all the three helices (and the three coils), with the smallest barriers found for those near the helix-3. For the quasi-bound water molecules near the first and second helices, we use well-known Kramers' theory to estimate the residence time from the free energy surface, by estimating the friction along the reaction coordinate from the diffusion coefficient by using Einstein relation. The agreement found is satisfactory. We discuss the possible biological function of these slow, quasi-bound (but transient) water molecules on the surface.