Molecular modelling of saccharide-lipid interactions

Certain saccharides, including trehalose, sucrose and glucose, stabilize lipid bilayers against dehydration. It has been suggested that these saccharides replace waters of hydration as the system is dried, thereby maintaining the headgroups at their hydrated spacing. The lipid acyl chains consequently have sufficient free volume to remain in the liquid crystallines state, and the processes that disrupt membrane integrity are inhibited. Initial molecular graphic investigations of a model trehalose/DMPC system supported this idea (Chandrasekhar, I. and Gaber, B.P. (1988) J. Biomol. Stereodyn, 5, 1163–1171). We have extended these studies to glucose and sucrose. A set of AMBER potential parameters has been established that reproduce simple saccharide conformations, including the anomeric effect. Extensive energy minimizations have been conducted on all three systems. The saccharide-lipid interaction energies become less stable in the order trehalose

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.

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

Conformational Polymorphism in G-tetraplex Structures: Strand Reversal by Base Flipover or Sugar Flipover

Guanine rich sequences adopt a variety of four stranded structures, which differ in strand orientation and conformation about the glycosidic bond even though they are all stabilised by Hoogsteen hydrogen bonded guanine tetrads. Detailed model building and molecular mechanics calculations have been carried out to investigate various possible conformations of guanines along a strand and different possible orientations of guanine strands in a G-tetraplex structure. It is found that for an oligo G stretch per se, a parallel four stranded structure with all guanines in anti conformation is favoured over other possible tetraplex structures. Hence an alternating syn-anti arrangement of guanines along a strand is likely to occur only in folded back tetraplex structures with antiparallel G strands. Our study provides a theoretical rationale for the observed alternation of glycosidic conformation and the inverted stacking arrangement arising from base flipover, in antiparallel G-tetraplex structures and also highlights the various structural features arising due to different types of strand orientations. The molecular mechanics calculations help in elucidating the various interactions which stabilize different G-tetraplex structures and indicate that screening of phosphate charge by counterions could have a dramatic effect on groove width in these four stranded structures.

Computer modeling studies on the binding of 2',5'-linked dinucleoside phosphates to ribonuclease T1-influence of subsite interactions on the substrate specificity

The modes of binding of Gp(2',5')A, Gp(2',5')C, Gp(2',5')G and Gp(2',5')U to RNase T1 have been determined by computer modelling studies. All these dinucleoside phosphates assume extended conformations in the active site leading to better interactions with the enzyme. The 5'-terminal guanine of all these ligands is placed in the primary base binding site of the enzyme in an orientation similar to that of 2'-GMP in the RNase T1-2'-GMP complex. The 2'-terminal purines are placed close to the hydrophobic pocket formed by the residues Gly71, Ser72, Pro73 and Gly74 which occur in a loop region. However, the orientation of the 2'-terminal pyrimidines is different from that of 2'-terminal purines. This perhaps explains the higher binding affinity of the 2',5'-linked guanine dinucleoside phosphates with 2'-terminal purines than those with 2'-terminal pyrimidines. A comparison of the binding of the guanine dinucleoside phosphates with 2',5'- and 3',5'-linkages suggests significant differences in the ribose pucker and hydrogen bonding interactions between the catalytic residues and the bound nucleoside phosphate implying that 2',5'-linked dinucleoside phosphates may not be the ideal ligands to probe the role of the catalytic amino acid residues. A change in the amino acid sequence in the surface loop region formed by the residues Gly71 to Gly74 drastically affects the conformation of the base binding subsite, and this may account for the inactivity of the enzyme with altered sequence i.e., with Pro, Gly and Ser at positions 71 to 73 respectively. These results thus suggest that in addition to recognition and catalytic sites, interactions at the loop regions which constitute the subsite for base binding are also crucial in determining the substrate specificity.

Conformational polymorphism in G-tetraplex structures: strand reversal by base flipover or sugar flipover

Guanlne rich sequences adopt a variety of four stranded structures, which differ in strand orientation and conformation about the glycosldic bond even though they are all stabilised by Hoogsteen hydrogen bonded guanlne tetrads. Detailed model building and molecular mechanics calculations have been carried out to investigate various possible conformations of guanlnes along a strand and different possible orientations of guanlne strands In a G-tetraplex structure. It is found that for an ollgo G stretch per se, a parallel four stranded structure with all guanines In anti conformation is favoured over other possible tetraplex structures. Hence an alternating syn-anti arrangement of guanlnes along a strand is likely to occur only in folded back tetraplex structures with antiparallel G strands. Our study provides a theoretical rationale for the observed alternation of glycosldic conformation and the inverted stacking arrangement arising from base filpover, In antlparallel G-tetraplex structures and also highlights the various structural features arising due to different types of strand orientations. The molecular mechanics calculations help in elucidating the various interactions which stabilize different G-tetraplex structures and indicate that screening of phosphate charge by counterions could have a dramatic effect on groove width in these four stranded structures.

Conformational Studies of Some Alpha-substituted Acetones CH3COCH2XCH3 Where X=CH2, O Or S And Phenacyl Sulfide

Both semiempirical and ab initio calculations are reported for conformational studies of a series of alpha-substituted acetones CH3COCH2XCH3 where X = CH2, O or S and of phenacyl sulfide PhCOCH2SCH3. For conformational studies in the lowest triplet state of these molecules, the MINDO/3 method was employed in the unrestricted Hartree-Fock frame. Results reveal that rotation around a bond alpha to the carbonyl group is more favourable than that around the beta bond. The preferred conformations in the lowest triplet state are nearly the same as in the ground state.

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.

Modelling multiple disulphide loop containing polypeptides by random conformation generation. The test cases of α-conotoxin GI and edothelin I

A general procedure for arriving at 3-D models of disulphiderich olypeptide systems based on the covalent cross-link constraints has been developed. The procedure, which has been coded as a computer program, RANMOD, assigns a large number of random, permitted backbone conformations to the polypeptide and identifies stereochemically acceptable structures as plausible models based on strainless disulphide bridge modelling. Disulphide bond modelling is performed using the procedure MODIP developed earlier, in connection with the choice of suitable sites where disulphide bonds could be engineered in proteins (Sowdhamini,R., Srinivasan,N., Shoichet,B., Santi,D.V., Ramakrishnan,C. and Balaram,P. (1989) Protein Engng, 3, 95-103). The method RANMOD has been tested on small disulphide loops and the structures compared against preferred backbone conformations derived from an analysis of putative disulphide subdatabase and model calculations. RANMOD has been applied to disulphiderich peptides and found to give rise to several stereochemically acceptable structures. The results obtained on the modelling of two test cases, a-conotoxin GI and endothelin I, are presented. Available NMR data suggest that such small systems exhibit conformational heterogeneity in solution. Hence, this approach for obtaining several distinct models is particularly attractive for the study of conformational excursions.

Molecular structure of Lantadene-B&C, triterpenoids ofantana camara, red variety: Lantadene-B, 22β-angeloyloxy-3-oxoolean-12-en-28-oic acid; Lantadene-C, 22 β-(S)-2′-methyl butanoyloxy-3-oxoolean-12-en-28-oic acid

(I)Lantadene-B: C35H52O5,M r =552.80, MonoclinicC2,a=25.65(1),b=6.819(9),c=18.75(1) Å,beta=100.61(9),V=3223(5) Å3,Z=4,D x =1.14 g cm–3 CuKagr (lambda=1.5418A),mgr=5.5 cm–1,F(000)=1208,R=0.118,wR=0.132 for 1527 observed reflections withF o ge2sgr(F o ). (II)Lantadene-C: C35H54O5·CH3OH,Mr=586.85, Monoclinic,P21,a=9.822(3),b=10.909(3),c=16.120(8)Å,beta=99.82(4),V=1702(1)Å3,Z=2,D x =1.145 g cm–3, MoKagr (lambda=0.7107Å), mgr=0.708 cm–1 F(000)=644,R=0.098, wR=0.094 for 1073 observed reflections. The rings A, B, C, D, and E aretrans, trans, trans, cis fused and are in chair, chair, sofa, half-chair, chair conformations, respectively, in both the structures. In the unit cell the molecules are stabilized by O-HctdotO hydrogen bonds in both the structures, however an additional C-HctdotO interaction is observed in the case of Lantadene-C.

Path integral description of polymers using fractional Brownian walks

The statistical properties of fractional Brownian walks are used to construct a path integral representation of the conformations of polymers with different degrees of bond correlation. We specifically derive an expression for the distribution function of the chains’ end‐to‐end distance, and evaluate it by several independent methods, including direct evaluation of the discrete limit of the path integral, decomposition into normal modes, and solution of a partial differential equation. The distribution function is found to be Gaussian in the spatial coordinates of the monomer positions, as in the random walk description of the chain, but the contour variables, which specify the location of the monomer along the chain backbone, now depend on an index h, the degree of correlation of the fractional Brownian walk. The special case of h=1/2 corresponds to the random walk. In constructing the normal mode picture of the chain, we conjecture the existence of a theorem regarding the zeros of the Bessel function.

Stereochemistry of 2'-5' linked nucleic acids: crystal and molecular structure of ammonium adenylyl-2',5'-adenosine tetrahydrate: a core fragment of 2'-5' oligo A's produced by interferon induced adenylate synthetase

The preponderance of 3'-5' phosphodiester links in nucleic acids is well known. Albeit less prevalent, the 2'-5' links are specifically utilised in the formation of 'lariat' in group II introns and in the msDNA-RNA junction in myxobacterium. As a sequel to our earlier study on cytidylyl-2',5'-adenosine we have now obtained the crystal structure of adenylyl-2',5'-adenosine (A2'p5'A) at atomic resolution. This dinucleoside monophosphate crystallizes in the orthorhombic space group P2(1)2(1)2(1) with a = 7.956(3) A, b = 12.212(3) A and c = 36.654(3) A. CuK alpha intensity data were collected on a diffractometer. The structure was sloved by direct methods and refined by full matrix least squares methods to R = 10.8%. The 2' terminal adenine is in the commonly observed anti (chi 2 = 161 degrees) conformation and the 5' terminal base has a syn (chi 1 = 55 degrees) conformation more often seen in purine nucleotides. A noteworthy feature of A2'p5'A is the intranucleotide hydrogen bond between N3 and O5' atoms of the 5' adenine base. The two furanose rings in A2'p5'A show different conformations - C2' endo, C3' endo puckering for the 5' and 2' ends respectively. In this structure too there is a stacking of the purine base on the ribose O4' just as in other 2'-5' dinucleoside structures, a feature characteristically seen in the left handed Z DNA. In having syn, anti conformation about the glycosyl bonds, C2' endo, C3' endo mixed sugar puckering and N3-O5' intramolecular hydrogen bond A2'p5'A resembles its 3'-5' analogue and several other 2'-5' dinucleoside monophosphate structures solved so far. Striking similarities between the 2'-5' dinucleoside monophosphate structures suggest that the conformation of the 5'-end nucleoside dictates the conformation of the 2' end nucleoside. Also, the 2'-5' dimers do not favour formation of miniature classical double helical structures like the 3'-5' dimers. It is conceivable, 2-5(A) could be using the stereochemical features of A2'p5'A which accounts for its higher activity.

Structure of a steriod fungal metabolite

3a,7a, 14a-Trihydroxypregn-16-en-20-one, C21H3204, M r = 348.48, orthorhombic, P212121, a = 9.211 (1), b = 13.201 (1), c = 16.031 (1),~, V = 1949.28 (29)/~3, Z = 4, Dx = 1.187 g cm -3, A(Cu Ka), = 1.5418 ,~,/z = 6.07 cm-l, F(000) = 760, T= 293 K, R = 0.061 for 1337 observations. The A, B and C rings adopt normal chair conformations with the D ring in a 14a-envelope conformation. The molecules are held together by two hydrogen bonds [0(3)..'0(20) = 2.879 and O(7).--O(14)= 2.612 A,].

Structure of pregna-4,16-diene-7[alpha],14[alpha]-diol-3,20-dione

7a,14a-Dihydroxypregna-4,16-diene-3,20- dione, C21H2804, M r = 344.45, orthorhombic, P212121, a = 7.136 (1), b = 12.342 (1), c = 20.049 (3)/k, V= 1765.7 (3)/k 3, Z = 4, Dx = 1.295 g cm -3, A(Cu Kte) = 1.5418/k, /z = 6.7 cm- a, F(000) = 744, T = 293 K, R = 0.048 for 1345 observations. The A ring may be described as in a l a,2flhalf- chair conformation or a l a-sofa conformation. The B and C rings adopt normal chair conformations and the D ring has a 14a-envelope conforma tion. The molecules are held together by a hydrogen bond [0(3)...0(7)= 2.767 A].

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