Aqueous channels within apolar peptide aggregates. Solvated helix of Boc-Aib-Ala-Leu-Aib-Ala-Leu-Aib- -Ala-Leu-Aib-OMe. H2O.CH3OH

Although the peptide Boc-Aibl-Ala2-Leu3- Aib4-Alas Leu'-Aib7-Ala8-Leu9-Aib'0-OMe [with a t-butoxycarbonyl(Boc) blocking group at the amino terminus, a methyl ester (OMe) at the carboxyl terminus, and four a-aminoisobutyric (Aib) residues] has a 3-fold repeat of residues, the helix formed by the peptide backbone is irregular. The carboxyl-terminal half assumes an at-helical form with torsion angles ) and r of approximately -60° and -45°, respectively, whereas the amino-terminal half is distorted by an insertion of a water molecule between the amide nitrogen of Ala5 [N(5)] and the carbonyl oxygen of Ala2 [0(2)]. The water molecule W(1) acts as a bridge by forming hydrogen bonds N(5).W(1) (2.93 A) and W(1)---0(2) (2.86 A). The distortion of the helix exposes the carbonyl oxygens of Aib' and Aib4 to the outside environment, with the consequence that the helix assumes an amphiphilic character despite having all apolar residues. Neighboring helices in the crystal run in antiparallel directions. On one side of a helix there are only hydrophobic contacts with efficient interdigitation of leucine side chains with those from the neighboring helix. On the other side of the helix there are hydrogen bonds between protruding carbonyl oxygens and four water molecules that separate two neighboring helices. Along the helix axis the helices bind head-to-tail with a direct hydrogen bond N(2)-0(9) (3.00 A). Crystals grown from methanol/water solution are in space group P2, with a = 15.778 ± 0.004 A, b = 11.228 ± 0.002 A, c = 18.415 ± 0.003 A, = 102.10 ± 0.02ur and two formula units per cell for C49HON1003 2H2OCH3OH. The overall agreement factorR is 7.5% for 3394 reflections observed with intensities >3a(F), and the resolution is 0.90 A.

Differences in hydration and association of helical Boc-Val-Ala-Leu-Aib-Val-Ala-Leu-(Val-Ala-Leu-Aib)2-OMe. xH2O in two crystalline polymorphs

The 15-residue apolar peptide, Boc-Val-Ala-Leu-Aib-Val-Ala-Leu-(Val-Ala-Leu-Aib)h2a-sO Mebeen crystallized from 2-propanol-water (form I). The crystal parameters for I are as follows:C74H133N15018*2H20s,p ace group P21, a = 9.185 (6) A, b = 47.410 (3) A, c = 10.325 (9) A, @ = 91.47(2)O, 2 = 2, R = 6.3% for 4532 reflections observed >3aQ, resolution 0.94 A. The structure isalmost completely a-helical with eleven 5-1 hydrogen bonds and one 441 hydrogen bond nearthe N-terminus. The structure has been compared with a polymorph (form 11) obtained frommethanol-water (Karle, I. L.; Flippen-Anderson, J. L.; Uma, K.; Sukumar, M.; Balaram, P., J. An.Chem. SOC19. 90,112,9350-9356). The two forms differ in the extent of hydration; form I contains two water molecules in the head-to-tail region of helical columns, while form I1 is more extensively solvated, with the equivalent of 7.5 water molecules. The three-dimensional packing of helices is completely parallel in I and antiparallel in 11.

b-Turn conformations in crystal structures of model peptides containing a,a-di-n-propylglycine (Dpg) and a,a-di-n-butyl-glycine (Dbg).

The crystal state conformations of three peptides containing the a,a-dialkylated residues, a,adi n-propylglycine (Dpg) and a,@-di-n-butylglycine (Dbg), have been established by x-ray diffraction. Boc-Ala-Dpg-Ala-OMe ( I ) and Boc-Ala-Dbg-Ala-OMe (III) adopt distorted type II @-turn conformations with Ala ( I ) and Dpg/Dbg (2) as the corner residues. In both peptides the conformational angles at the Dxg residue (I: 4 = 66.23 J/ = 19.3'; III: 4 = 66S0, J. = 21 .la)deviate appreciablyfrom ideal values for the i + 2 residue in a type II @-turn. In both peptides the observed(N. 0) distances between the Boc CO andAla(3) NHgroups are far too long (I:3.44 k; III: 3.63 k) for an intramolecular 4 + 1 hydrogen bond. Boc-Ala-Dpg-Ala-NHMe (II)crystallizes with two independent molecules in the asymmetric unit. Both molecules IIA and IIB adopt consecutive @-turn (type III-III in IIA and type III-I in IIB) or incipient 3,,,-helical structures, stabilized by two intramolecular 4 --t I hydrogen bonds. In all four molecules the bond angle N-C"-C' ( T ) at the Dxg residues are 2 1109 The observation of conformational angles in the helical region of 4,J/ space at these residues is consistent with theoretical predictions

Diproline Templates as Folding Nuclei in Designed Peptides. Conformational Analysis of Synthetic Peptide Helices Containing Amino Terminal Pro-Pro Segments

The effect of N-terminal diproline segments in nucleating helical folding in designed peptides has been studied in two model sequences Piv-Pro-Pro-Aib-Leu-Aib-Phe-OMe (1) and Boc-Aib-Pro-Pro-Aib-Val-Ala-Phe-OMe (2). The structure of 1 in crystals, determined by X-ray diffraction, reveals a helical (RR) conformation for the segment residues 2 to 5, stabilized by one 4 -> 1 hydrogen bond and two 5 -> 1 interactions. The N-terminus residue, Pro(1) adopts a polyproline II (P-II) conformation. NMR studies in three different solvent systems support a conformation similar to that observed in crystals. In the apolar solvent CDCl3, NOE data favor the population of both completely helical and partially unfolded structures. In the former, the Pro-Pro segment adopts an alpha(R)-alpha(R) conformation, whereas in the latter, a P-II-alpha(R) structure is established. The conformational equilibrium shifts in favor of the P-II-alpha(R) structure in solvents like methanol and DMSO. A significant population of the Pro(1)- Pro(2) cis conformer is also observed. The NMR results are consistent with the population of at least three conformational states about Pro- Pro segment: trans alpha(R)-alpha(R), trans P-II-alpha(R) and cis P-II-alpha(R). Of these, the two trans conformers are in rapid dynamic exchange on the NMR time scale, whereas the interconversion between cis and trans form is slow. Similar results are obtained with peptide 2. Analysis of 462 diproline segments in protein crystal structures reveals 25 examples of the alpha(R)-alpha(R) conformation followed by a helix. Modeling and energy minimization studies suggest that both P-II-alpha(R) and alpha(R)-alpha(R) conformations have very similar energies in the model hexapeptide 1

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.

Steric and rotational constraints in the X-ray structure of the animalarial drug amodiaquine

In the crystal structure of the antimalarial drug amodiaquine, the bonds linking the quinoline and the phenyl groups show partial double-bond character. The partial double-bond character of the two exocyclic bonds, together with stereochemical constraints, reduce flexibility of the two ring systems of the molecule. The dihedral angle between the two ring planes is lowest compared to those in the antileukaemic drug amsacrine and its derivatives. CPK-modelling studies suggest the way amodiaquine can bind to DNA. Stacking interaction between the quinoline and phenyl groups of independent molecules and the hydrogen-bond network stabilize the crystal structure.

Conformational Variability in Modified Nucleosides. Crystal and Molecular Structure of 2',3'-O-Isopropylidene Inosine

The crystal structure of 2',3'-O-isopropylidene inosine shows a number of interesting features. The four independent molecules in the asymmetric unit exhibit significant conformational variations. Ribose puckers fall in the O(4')-exo region, unfavourable in unsubstituted nucleosides. Hypoxanthine bases show base-pairing (I.I) in a manner analogous to the guanine self pairs (G.G) in 2',3'-O-isopropylidene guanosine but with a C(2)-H…O(6) hydrogen bond instead of N(2)-H…O(6).

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.

Ab initio molecular orbital calculations on ion-molecule and ion pair-molecule complexes of the water-lithium cyanide system

Ab initio molecular orbital (MO) calculations with the 3-21G and 6-31G basis sets were performed on a series of ion-molecule and ion pair-molecule complexes for the H2O + LiCN system. Stabilisation energies (with counter-poise corrections), geometrical parameters, internal force constants and harmonic vibrational frequencies were evaluated for 16 structures of interest. Although the interaction energies are smaller, the geometries and relative stabilities of the monohydrated contact ion pair are reminiscent of those computed for the complexes of the individual ions. Thus, interaction of the oxygen lone pair with lithium leads to a highly stabilised C2v structure, while the coordination of water to the cyanide ion involves a slightly non-linear hydrogen bond. Symmetrical bifurcated structures are computed to be saddle points on the potential energy surface, and to have an imaginary frequency for the rocking mode of the water molecule. On optimisation the geometries of the solvent shared ion pair structures (e.g. Li+cdots, three dots, centered OH2cdots, three dots, centered CN−) revealed a proton transfer from the water molecule leading to hydrogen bonded forms such as Li-O-Hcdots, three dots, centered HCN. The variation in the force constants and harmonic frequencies in the various structures considered are discussed in terms of ion-molecular and ion pair-molecule interactions.

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.

Electronic excitations and alternation of conjugated polymers

The electronic excitations and fluorescence of conjugated polymers are related to large or small alternation ? of the transfer integrals t(1 ± ?) along the backbone. The fluorescence of polysilanes (PSs) and poly (para-phenylenevinylene (PPV) is linked to large ?, which places the one-photon gap Eg below the lowest two-photon gap Ea and reduces distortions due to electron-phonon (e-p) coupling. In contrast to small ? not, vert, similar 0.1 in ?-conjugated polymers, such as polyacetylene (PA), para-conjugated phenyls lead to an extended ?-system with increased alternation, to states localized on each ring and to charge-transfer excitations between them. Surprisingly good agreement is found between semiempirical parametric method 3 (PM3) bond lengths and exact Pariser-Parr-Pople (PPP) ?-bond orders for trans-stilbene, where the PPV bipolarons are confined to two phenyls. Stilbene spectra are consistent with increased alternation and small e-p distortions.

Exact eigenstates of the Pariser-Parr-Pople model for anthracene

The pi-electronic structure of anthracene is discussed by combining exact solutions of the Pariser-Parr-Pople (PPP) model and semiempirical PM3 calculations. Symmetry adaptation of the 2.8 million singlet valence-bond (VB) diagrams is explicitly demonstrated for D2h and electron-hole symmetry. Standard PPP parameters provide a comprehensive fit to one- and two-photon anthracene spectra and intensities up to the strong 1 B-1(3u)-absorption at 5.24 eV, the 10th excited state in the dense correlated spectrum, and indicate a reassignment of two-photon absorptions. The singlet-triplet gap and fine-structure constants also agree with experiment. Fully-relaxed PM3 geometries are obtained for the anthracene ground state and for singlet, triplet, and charged bipolarons. The PM3 bond lengths correlate well with PPP bond orders for the idealized structure. Single-determinantal PM3 excitation and relaxation energies for bipolarons are consistent with exact PPP results and contrast all-valence electron with pi-electron calculations. Several correlation effects are noted in the rich pi-spectra of anthracene in connection with improved PPP modeling of conjugated molecules and polymers.

Is tetramethyleneethane a ground state triplet?

We have examined a number of possible ways by which tetramethyleneethane (TME) can be a ground state triplet, as claimed by experimental studies, in violation of Ovchinnikov’s theorem for alternant hydrocarbons of equal bond lengths. Model exact π calculations of the low-lying states of TME, 3,4-dimethylenefuran and 3,4-dimethylenepyrrole were carried out using a diagrammatic valence bond approach. The calculations failed to yield a triplet ground state even after (a) tuning of electron correlation, (b) breaking alternancy symmetry, and (c) allowing for geometric distortions. In contrast to earlier studies of fine structure constants in other conjugated systems, the computedD andE values of all the low-lying triplet states of TME for various geometries are at least an order of magnitude different from the experimentally reported values. Incorporation of σ-π mixing by means of UHF MNDO calculations is found to favour a singlet ground state even further. A reinterpretation of the experimental results of TME is therefore suggested to resolve the conflict.

Lipid-amphotericin B complex structure in solution: A possible first step in the aggregation process in cell membranes

The interactions between the polyene antibiotic amphotericin B with dipalmitoylphosphatidylcholine were investigated in vesicles (using circular dichroism) and in chloroform solution (using circular dichroism and IH, I3C, and 31P nuclear magnetic resonance). The results show that amphotericin B readily aggregates in vesicles and that the extent of aggregation depends on the 1ipid:drug concentration ratio. Introduction of sterol molecules into the membrane hastens the process of aggregation of amphotericin B. In chloroform solutions amphotericin B strongly interacts with phospholipid molecules to form a stoichiometric complex. The results suggest that there are interactions between the conjugated heptene stretch of amphotericin B and the methylene groups of lipid acyl chains, while the sugar moiety interacts with the phosphate head group by the formation of a hydrogen bond. A model is proposed for the lipid-amphotericin B complex, in which amphotericin B interacts equally well with the two lipid acyl chains, forming a 1:l complex.