Structural Investigations on Poly(4-hydroxy-L-proline). 2. Physicochemical Studies

We report experimental studies which confirm our prediction, namely that the ordered structure of poly(hydroxypro1ine) in solution corresponds to a left-handed helical structure with intrachain hydrogen bonds. The CD studies show that the poly(hydroxypro1ine) molecule has essentially the same conformation in aqueous solution and in the film obtained subsequently by evaporation. X-ray diffraction patterns of the sample in this form (B form) have been recorded at different relative humidities. The patterns recorded at relative humidities over 66% can be interpreted in terms of a helical structure with intrachain hydrogen bonds. These results lead us to conclude that the ordered conformation of poly(hydroxypro1ine) in solution is form B and not form A. This offers a simple explanation for the greater stability of the poly(hydroxypro1ine) helix in solution as compared to the poly(pro1ine) form I1 helix and also for the absence of mutarotation for poly(hydroxypro1ine).

Structural Investigations on Poly(4-hydroxy-L-proline). 1. Theoretical Studies

Model building studies on poly(hydroxypro1ine) indicate that in addition to the well-known helical structure of form A, a left-handed helical structure with trans peptide units and with h = 2.86 A and n = 2.67 (i.e., 8 residues in 3 turns) is also possible. In this structure which is shown to be in agreement with X-ray data of the form B in the next paper, the y-hydroxyl group of an (i + 1)th Hyp residue is hydrogen bonded to the carbonyl oxygen of an (i - 1)th residue. The possibility of a structure with cis peptide units is ruled out. It is shown that both forms A and B are equally favorable from considerations of intramolecular energies. Since form B is further stabilized by intrachain hydrogen bonds, we believe that this is likely to be the ordered conformation for poly(hydroxypro1ine) in water.

Modes of binding of guanosine monophosphates to ribonuclease T1 - A computer modelling study

Computer-modelling studies on the modes of binding of the three guanosine monophosphate inhibitors 2'-GMP, 3'-GMP, and 5'-GMP to ribonuclease (RNase) T1 have been carried out by energy minimization in Cartesian-coordinate space. The inhibitory power was found to decrease in the order 2'-GMP > 3'-GMP > 5'-GMP in agreement with the experimental observations. The ribose moiety was found to form hydrogen bonds with the protein in all the enzyme-inhibitor complexes, indicating that it contributes to the binding energy and does not merely act as a spacer between the base and the phosphate moieties as suggested earlier. 2'-GMP and 5'-GMP bind to RNase T1 in either of the two ribose puckered forms (with C3'-endo more favoured over the C2'-endo) and 3'-GMP binds to RNase T1 predominantly in C3'-endo form. The catalytically important residue His-92 was found to form hydrogen bond with the phosphate moiety in all the enzyme-inhibitor complexes, indicating that this residue may serve as a general acid group during catalysis. Such an interaction was not found in either X-ray or two-dimensional NMR studies.

Structure of valinomycin by molecular dynamics studies

Valinomycin is an important ionophore which exhibits a high conformational flexibility. The study of various conformations adopted by this molecule together with the study of flexibility in a given conformation can throw light on the ion transport by the ionophore across the membrane. Molecular dynamics (MD) studies are ideal to characterize the flexibility in different parts of the molecule and can also give an idea of various conformations adopted by the molecule at a given temperature. Hence MD studies at 100K have been carried out on the minimized crystal structure of the molecule to scan the possible conformations in the neighbourhood of the well known 'bracelet' like structure of uncomplexed Valinomycin, Properties, like the flexibility, average values, r.m.s. fluctuations of the various intramolecular hydrogen bonds are discussed. Energy minimization has been carried out on selected MD simulated points to analyze the characteristics of the unique conformation adopted by this molecule at this temperature.

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.

Fourfold Helical Structures for Polypeptides

Fourfold helical structures for polypeptides and their association in regular lattices with interchain hydrogen bonds were investigated by model building studies. These studies revealed that stereochemically satisfactory fourfold helical sturctures are possible for polyglycine, polyproline, and copolymers of glycine and proline with two and four units in the monomer. In these structures the unit height h for the backbone has been found to be restricted from 2.7 to 3.1 k, with four peptide units per turn of the helix. Energetically both fourfold and threefold helical structures are equally favorable.

X-ray studies on crystalline complexes involving amino acids and peptides. XXI. Structure of a (1:1) complex between L-phenylalanine and D-valine

CsHllNO2.C9HilNO2, Mr = 282.3, P1, a = 5.245 (1), b = 5.424 (1), c = 14.414 (2) A, a = 97.86 (1), fl = 93-69 (2), y = 70-48 (2) °, V= 356 A 3, Z = 1, O m = 1-32 (2), Dx = 1.32 g cm-3, h(Mo Ka) = 0-7107 A, g = 5-9 cm-1, F(000) = 158, T= 298 K, R=0.035 for 1518 observed reflections with I>2tr(I). The molecules aggregate in double layers, one ayer made up of L-phenylalanine molecules and the other of D-valine molecules. Each double layer is stabilized by interactions involving main-chain atoms of both types of molecules. The interactions include hydrogen bonds which give rise to two head-to-tail sequences. The arrangement of molecules in the complex is almost the same as that in the structure of DL-valine (and DL-leucine and DL-isoleucine) except for the change in the side chain of L molecules. The molecules in crystals containing an equal number of L and O hydrophobic amino-acid molecules thus appear to aggregate in a similar fashion, irrespective of the precise details of the side chain.

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.

Interactions of tris buffer with nucleotides: the crystal structure of tris(hydroxymethyl)methylammonium adenosine 5'-diphosphate dihydrate

The crystal and molecular structures of the Tris salt of adenosine 5'-diphosphate were determined from X-ray diffraction data. The crystals are monoclinic, space P21, and Z = 2 with a=9.198 (2) A, b=6.894 (1) A, c=18.440 (4) A, and beta = 92.55 (2) degrees. Intensity data were collected on an automated diffractometer. The structure was solved by the heavy-atom technique and refined by least squares to R = 0.047. The ADP molecule adopts a folded conformation. The conformation about the glycosidic bond is anti. The conformation of the ribose ring is close to a perfect C(2')-endo-C-(3')-exo puckering. The conformation about C(4')-C(5') is gauche-gauche, similar to other nucleotide structures. The pyrophosphate chain displays a nearly eclipsed geometry when viewed down the P-P vector, unlike the staggered conformation observed in crystal structures of other pyrophosphates. The less favorable eclipsed conformation probably results from the observed association of Tris molecules with the polar diphosphate chain through electrostatic interactions and hydrogen bonds. Such interactions may play an important role in Tris-buffered aqueous solutions of nucleotides and metal ions.

Conformational studies on beta-bend containing a cis peptide unit

Conformational studies have been carried out on the X-cis-Pro tripeptide system (a system of three linked peptide units, in the trans-cis-trans configuration) using energy minimization techniques. For X, residues Gly, L-Ala, D-Ala and L-Pro have been used. The energy minima have been classified into different groups based upon the conformational similarity. There are 15, 20, 18 and 6 minima that are possible for the four cases respectively and these fall into 11 different groups. A study of these minima shows that, (i) some minima contain hydrogen bonds - either 4-->1 or 1-->2 type, (ii) the low energy minima qualify themselves as bend conformations, (iii) cis' and trans' conformations are possible for the prolyl residue as also the C(gamma)-endo and C(gamma)-exo puckerings, and (iv) for Pro-cis-Pro, cis' at the first prolyl residue is ruled out, due to the high energy. The available crystal structure data on proteins and peptides, containing cis-Pro segment have been examined with a view to find the minima that occur in solid state. The data from protein show that they fall under two groups. The conformation at X in X-cis-Pro is near extended when it is a non-glycyl residue. In both peptides and proteins there exists a preference for trans' conformation at prolyl residue over cis' when X is a non-glycyl residue. The minima obtained can be useful in modelling studies.

Molecular dynamics simulations on parallel and antiparallel C.G*G triplexes

Molecular dynamics (MD) studies have been carried out on the Hoogsteen hydrogen bonded parallel and the reverse Hoogsteen hydrogen banded antiparallel C.G*G triplexes. Earlier, the molecular mechanics studies had shown that the parallel structure was energetically more favourable than the antiparallel structure. To characterize the structural stability of the two triplexes and to investigate whether the antiparallel structure can transit to an energetically more favourable structure, due to the local fluctuations in the structure during the MD simulation, the two structures were subjected to 200ps of constant temperature vacuum MD simulations at 300K. Initially no constraints were applied to the structures and it was observed that for the antiparallel tripler, the structure showed a large root mean square deviation from the starting structure within the first 12ps and the N4-H41-O6 hydrogen bond in the WC duplex got distorted due to a high propeller twist and a moderate increase in the opening angle in the basepairs. Starting from an initial value of 30 degrees, helical twist of the average structure from this simulation had a value of 36 degrees, while the parallel structure stabilized at a twist of 33 degrees. In spite of the hydrogen bond distortions in the antiparallel tripler, it was energetically comparable to the parallel tripler. To examine the structural characteristics of an undistorted structure, another MD simulation was performed on the antiparallel tripler by constraining all the hydrogen bonds. This structure stabilized at an average twist of 33 degrees. In the course of the dynamics though the energy of the molecule - compared to the initial structure - improved, it did not become comparable to the parallel structure. Energy minimization studies performed in the presence of explicit water and counterions also showed the two structures to be equally favourable energetically Together these results indicate that the parallel C.G*G tripler with Hoogsteen hydrogen bonds also represents a stereochemically and energetically favourable structure for this class of triplexes.

Sequence Independent Recombination Triple Helices: A Molecular Dynamics Study

Recent experimental studies have shown that the Rec-A mediated homologous recombination reaction involves a triple helical intermediate, in which the third strand base forms hydrogen bonds with both the bases in the major groove of the Watson-Crick duplex. Such 'mixed' hydrogen bonds allow formation of sequence independent triplexes. DNA triple helices involving 'mixed' hydrogen bonds have been studied, using model building, molecular mechanics (MM) and molecular dynamics (MD). Models were built for a tripler comprising all four possible triplets viz., G.C*C, C.G*G, A.T*T and T.A*A. To check the stability of all the 'mixed' hydrogen bonds in such triplexes and the conformational preferences of such tripler structures, MD studies were carried out starting from two structures with 30 degrees and 36 degrees twist between the basepairs. It was observed that though the two triplexes converged towards a similar structure, the various hydrogen bonds between the WC duplex and the third strand showed differential stabilities. An MD simulation with restrained hydrogen bonds showed that the resulting structure was stable and remained close to the starting structure. These studies help us in defining stable hydrogen bond geometries involving the third strand and the WC duplex. It was observed that in the C.G*G triplets the N7 atom of the second strand is always involved in hydrogen bonding. In the G.C*C triplets, either N3 or O2 in the third strand cytosine can interchangeably act as a hydrogen bond acceptor.

Molecular Mechanics Studies of Homopolymeric and Mixed Sequence Py.Pu*Pu Triple Helices,

DNA triple helices containing two purine strands and one pyrimidine strand (C.G*G and T.A*A) have been studied, using model building followed by energy minimisation, for different orientations of the third strand resulting from variation in the hydrogen bonding between the Watson-Crick duplex and the third strand and the glycosidic torsion angle in the third strand. Our results show that in the C.G*G case the structure with a parallel orientation of the third strand, resulting from Hoogsteen hydrogen bonds between the third strand and the Watson-Crick duplex, is energetically the most favourable while in the T.A*A case the antiparallel orientation of the third strand, resulting from reverse Hoogsteen hydrogen bonds, is energetically the most favourable. These studies when extended to the mixed sequence triplexes, in which the second strand is a mixture of G and A, correspondingly the third strand is a mixture of G and APT, show that though the parallel orientation is still energetically more favourable, the antiparallel orientation becomes energetically comparable with an increasing number of thymines in the third strand. Structurally, for the mixed triplexes containing G and T in the third strand, it is seen that the basepair non-isomorphism between the C.G*G and the T.A*T triplets can be overcome with some changes in the base pair parameters without much distortion of either the backbone or the hydrogen bonds.

Structural Polymorphism in d(T).d(A)*d(T) Triple Helices

DNA triple helices containing two thymine strands and one adenine strand have been studied, using model building followed by energy minimisation, for different orientations of the third strand resulting from variation in the hydrogen bonding between the Watson-Crick duplex and the third strand and the glycosidic torsion angle in the third strand. Our results show that the structure with a parallel orientation of the third strand, in which the third strand base forms Hoogsteen hydrogen bonds with the adenine base in the Watson-Crick duplex, is energetically the most favourable. An antiparallel orientation of the third strand is also possible, in which the third strand base hydrogen bonds to both the bases in the Watson-Crick duplex. This structure is energetically comparable to the parallel structure. For the parallel triplex a 200ps molecular dynamics simulation starting from two different starting structures indicates that at 300K significant structural heterogeneity exists in this tripler structure. The results are compared with existing structural data on this class of triplexes derived from theoretical and NMR techniques.

Conformational Polymorphism in Telomeric Structures: Loop Orientation and Interloop Pairing in d(G4TnG4),

Sequence repeats constituting the telomeric regions of chromosomes are known to adopt a variety of unusual structures, consisting of a G tetraplex stem and short stretches of thymines or thymines and adenines forming loops over the stem. Detailed model building and molecular mechanics studies have been carried out for these telomeric sequences to elucidate different types of loop orientations and possible conformations of thymines in the loop. The model building studies indicate that a minimum of two thymines have to be interspersed between guanine stretches to form folded-back structures with loops across adjacent strands in a G tetraplex (both over the small as well as large groove), while the minimum number of thymines required to build a loop across the diagonal strands in a G tetraplex is three. For two repeat sequences, these hairpins, resulting from different types of folding, can dimerize in three distinct ways-i.e., with loops across adjacent strands and on same side, with loops across adjacent strands and on opposite sides, and with loops across diagonal strands and on opposite sides-to form hairpin dimer structures. Energy minimization studies indicate that all possible hairpin dimers have very similar total energy values, though different structures are stabilized by different types of interactions. When the two loops are on the same side, in the hairpin dimer structures of d(G(4)T(n)G(4)), the thymines form favorably stacked tetrads in the loop region and there is interloop hydrogen bonding involving two hydrogen bonds for each thymine-thymine pair. Our molecular mechanics calculations on various folded-back as well as parallel tetraplex structures of these telomeric sequences provide a theoretical rationale for the experimentally observed feature that the presence of intervening thymine stretches stabilizes folded-back structures, while isolated stretches of guanines adopt a parallel tetraplex structure