Lincomycin and clindamycin conformations. A fragment shared by macrolides, ketolides and lincosamides determined from TRNOE ribosome-bound conformations

Two important lincosamide antibiotics, lincomycin and clindamycin were studied in the complex state with the bacterial ribosome after a conformational analysis by 1H and 13C NMR spectroscopy and molecular modelling of the unbound molecules. Lincosamide-ribosome interactions were investigated using two-dimensional transferred nuclear Overhauser effect spectroscopy (TRNOESY), resulting in a bound structure compatible with the experimental NMR data. The results compared with the conformational analysis of the substrates in solution indicate that specific conformations are preferred in the bound state. Clindamycin, the more bioactive antibiotic studied, displayed a stronger NMR response than lincomycin showing that in lincosamide-ribosome interactions, a low affinity binding level is associated to the tight binding one and is related to biological activity. This study shows that conformation plays an essential role for the low affinity binding site. Superimposition of lincosamide, macrolide and ketolide bound structures exhibited conformational similarities in a particular fragment which is in agreement with a hypothesis of partial overlapping lincosamide and macrolide binding sites.     

Molecular mechanics (MM3) study of the conformations of ethyl esters of diastereoisomeric 3-substituted 4,4,4-trichloro-2-cyano-butanoic acids

Molecular mechanics calculation (MM3 force field) were used to study the conformations of diastereoisomeric pairs of ethyl esters of 3-substituted 4,4,4-trichloro-2-cyano-butanoic acids in order to test an assumption from ¹H NMR spectra about their preferred conformations. On the basis of the most probable preferred conformation in each case, an assignment of the relative configurations of these newly synthesized compounds was made in the earlier study. The present computed most favorable conformations are in accord with the earlier qualitative considerations in all cases, thus validating the assignment of the configurations of the diastereoisomers made previously.     

Conformation of the O-specific polysaccharide of Shigella dysenteriae type 1: molecular modeling shows a helical structure with efficient exposure of the antigenic determinant alpha-L-Rhap-(1-->2)-alpha-D-Galp.

The O-specific polysaccharide of Shigella dysenteriae type 1, which has the repeating tetrasaccharide unit -->3)-alpha-L-Rhap-(1-->3)-alpha-L-Rhap-(1-->2)-alpha-D-Galp-(1-->3)-alpha-D-GlcNAcp-(1--> (A-B-C-D), is a major virulence factor, and it is believed that antibodies against this polysaccharide confer protection to the host. The conformational properties of fragments of this O-antigen were explored using systematic search with a modified HSEA method (GLYCAN) and with molecular mechanics MM3(96). The results show that the alpha-D-Gal-(1-->3)-alpha-D-GlcNAc linkage adopts two favored conformations, phi/psi approximately equal to -40 degrees /-30 degrees (I) and approximately 15 degrees /30 degrees (II), whereas the other glycosidic linkages only have a single favored phi/psi conformational range. MM3 indicates that the trisaccharide B-C-D and tetrasaccharides containing the B-C-D moiety exist as two different conformers, distinguished by the conformations I and II of the C-D linkage. For the pentasaccharide A-B-C-D-A' and longer fragments, the calculations show preference for the C-D conformation II. These results can explain previously reported nuclear magnetic resonance data. The pentasaccharide in its favored conformation II is sharply bent, with the galactose residue exposed at the vertex. This hairpin conformation of the pentasaccharide was successfully docked with the binding site of a monoclonal IgM antibody (E3707 E9) that had been homology modeled from known crystal structures. For fragments made of repetitive tetrasaccharide units, the hairpin conformation leads to a left-handed helical structure with the galactose residues protruding radially at the helix surface. This arrangement results in a pronounced exposure of the galactose and also the adjacent rhamnose in each repeating unit, which is consistent with the known role of the as alpha-L-Rhap-(1-->2)-alpha-D-Galp moiety as a major antigenic epitope of this O-specific polysaccharide.     

Addition reactions at the 16(17) double bond of 3-methoxy-13alpha-estra-1,3,5(10),16-tetraene

The epoxidation, the addition of hypobromous acid, and the hydroboration of 3-methoxy-13alpha-estra-1,3,5(10),16-tetraene 1 with diborane, catecholborane, and 9-BBN were investigated in order to determine the stereochemical outcome and to synthesize new 13alpha-estra-1,3,5(10)-trienes for biological and conformational investigations. It was shown that the sterically demanding reagent 9-BBN participated in a preferred beta attack (53% 16betaOH 10, 34% 17betaOH 8, 13% 16alphaOH 11). This stereochemical result is in agreement with that from another cis addition reaction, the recently described OsO4 dihydroxylation of 1 [Steroids 68 (2003) 113]. With smaller reagents such as B2H6, catecholborane, or magnesium monoperoxyphthalate, a diminished stereoselectivity was observed with only a slight excess of beta attack. The ionic trans addition of hypobromous acid gave two 17-bromo-16-alcohols with 16beta,17alpha (4, 76%) and 16alpha,17beta configuration (5, 24%) formed by trans cleavage of the 16,17alpha- and beta-bromonium ion at position 16. The same regioselective and stereoselective course was found for the cleavage of the 16alpha,17alpha- and 16beta,17beta-epoxides (3 and 2) with hydrazoic acid (3-->16betaN3,17alphaOH 7, 2-->16alphaN3,17betaOH 6). The stereochemistry of the addition reactions to 1 can be explained in terms of a twist-boat conformation involving the C ring of compound 1. From a synthetic viewpoint the synthesis of the beta-epoxide 2 from the bromohydrin 4, the cleavage of this epoxide to 16alpha-substituted-17beta-hydroxy compounds, such as 6, and hydroboration/oxidation with 9-BBN to the hitherto unknown 16beta-hydroxy compound 10 are useful procedures. The bromohydrin 5 is the first 13alpha-steroid with a 17beta-bromo substituent. X-ray analysis revealed twist-boat and 16beta-envelope conformations for rings C and D, respectively.     

Conformational studies of irreversible HIV-1 protease inhibitors containing cis-epoxide as an amide isostere.

We have carried out NMR and molecular modeling studies of peptidomimetic HIV-1 protease inhibitors, LB71116: Qc-Asn-Phepsi[(1R,2S)-cis-epoxide]Gly-NH-CH(isopropyl)2 where Qc stands for quinaldic acid and LB71148: Qc-(SMe)Pen(O)2-Phepsi[(1R,2S)-cis-epoxide]Gly-NH-CH(isoprop yl)2 where (SMe)Pen(O)2 stands for S-methyl-S-dioxo-penicillamine. Through conformational calculations and NMR data analysis, we have obtained preferred conformations of the two inhibitors in solution. To our knowledge, this work is one of the first extensive conformational studies of peptidomimetics containing cis-epoxide amide isostere. The resulting preferred conformations contain extended structures. In these conformations, the psi of Phe(cep) is maintained about 130 degrees and the phi angle of (cep)Gly prefers +/- 150 degrees [where Phe(cep) and (cep)Gly are the residues generated by the replacement of the Phe-Gly peptide bond with cis-epoxide]. Two conformations were commonly observed in the preferred conformations of each inhibitor. Through restrained molecular dynamics simulating the hydrogen bond formation between our inhibitor and a water molecule ('flap water'), one of the conformations is assumed as the conformation which can bind to the enzyme without large conformational changes. Recently, we had the opportunity to compare the selected preferred conformation with the binding conformation of LB71116 observed from the X-ray studies of the complex between LB71116 and HIV-1 protease. These two conformations are surprisingly similar to each other. Thus, we can explain high activity and selectivity of our inhibitors to the HIV-1 protease by the similarity between the preferred conformations in solution and the binding conformation.     

Conformation of the exopolysaccharide of Burkholderia cepacia predicted with molecular mechanics (MM3) using genetic algorithm search

We present a computational conformational analysis of the exopolysaccharide of Burkholderia cepacia, which is believed to play a role in colonization and persistence of B. cepacia in the lungs of cystic fibrosis patients. The repeating unit of the exopolysaccharide is a heptasaccharide with three branches, which cause significant steric restraints. Conformational searches using glygal, an in-house developed software using genetic algorithm search methods, were performed on fragments as well as on the complete repeating unit with wrap-over residues. The force field used for the calculations was MM3(96). The search showed four favored conformations for an isolated repeating unit. However, for a sequence of several repeating units, the calculations indicate a single, well-defined linear conformation.     

A nonhelical,multiple β-turn conformation in a glycine-rich heptapeptide fragment of trichogin A IV containing a single central α-aminoisobutyric acid residue

The conformational properties of the protected seven-residue C-terminal fragment the lipopeptaibol antibiotic Trichogin A IV (Boc-Gly-Gly-Leu-Aib-Gly-Ile-Leu-OMe) has been examined in CDCl₃ and (CD₃)₂SO by ¹H-nmr. Evidence for a multiple β-turn conformation [type I′ at Gly(1)-Gly(2), type II at Leu(3)-Aib(4), and a type I′ at Aib(4)-Gly(5)] suggests that Leu(3) has preferred an extended or semiextended conformation over a helical conformation in CDCl₃. This structure is thus in contrast to earlier observations of seven-residue peptides containing a single central Aib preferring helical conformations in both solution and crystalline slates. A structural transition to a frayed right-handed helix is absented in (CD₃)₂SO. These results suggest that nonhelical conformations may be important in Gly-rich peptides containing Aib. Further, the presence of amino acids with contradictory influences on backbone conformational freedom can lead to well-defined conformational transitions even in small peptides. © 1995 John Wiley & Sons, Inc.     

Pi-turns in proteins and peptides: Classification, conformation, occurrence, hydration and sequence.

The i + 5-->i hydrogen bonded turn conformation (pi-turn) with the fifth residue adopting alpha L conformation is frequently found at the C-terminus of helices in proteins and hence is speculated to be a "helix termination signal." An analysis of the occurrence of i + 5-->i hydrogen bonded turn conformation at any general position in proteins (not specifically at the helix C-terminus), using coordinates of 228 protein crystal structures determined by X-ray crystallography to better than 2.5 A resolution is reported in this paper. Of 486 detected pi-turn conformations, 367 have the (i + 4)th residue in alpha L conformation, generally occurring at the C-terminus of alpha-helices, consistent with previous observations. However, a significant number (111) of pi-turn conformations occur with (i + 4)th residue in alpha R conformation also, generally occurring in alpha-helices as distortions either at the terminii or at the middle, a novel finding. These two sets of pi-turn conformations are referred to by the names pi alpha L and pi alpha R-turns, respectively, depending upon whether the (i + 4)th residue adopts alpha L or alpha R conformations. Four pi-turns, named pi alpha L'-turns, were noticed to be mirror images of pi alpha L-turns, and four more pi-turns, which have the (i + 4)th residue in beta conformation and denoted as pi beta-turns, occur as a part of hairpin bend connecting twisted beta-strands. Consecutive pi-turns occur, but only with pi alpha R-turns. The preference for amino acid residues is different in pi alpha L and pi alpha R-turns. However, both show a preference for Pro after the C-termini. Hydrophilic residues are preferred at positions i + 1, i + 2, and i + 3 of pi alpha L-turns, whereas positions i and i + 5 prefer hydrophobic residues. Residue i + 4 in pi alpha L-turns is mainly Gly and less often Asn. Although pi alpha R-turns generally occur as distortions in helices, their amino acid preference is different from that of helices. Poor helix formers, such as His, Tyr, and Asn, also were found to be preferred for pi alpha R-turns, whereas good helix former Ala is not preferred. pi-Turns in peptides provide a picture of the pi-turn at atomic resolution. Only nine peptide-based pi-turns are reported so far, and all of them belong to pi alpha L-turn type with an achiral residue in position i + 4. The results are of importance for structure prediction, modeling, and de novo design of proteins.     

Theoretical studies of cis-Pt(II)-diammine binding to duplex DNA

The binding of cis-Pt(II) diammine (cis-DP) to double-stranded DNA was studied with several kinked conformations that can accommodate the formation of a square planar complex. Molecular mechanics (MM) calculations were performed to optimize the molecular fit. These results were combined with quantum mechanical (QM) calculations to ascertain the relative energetics of ligand binding through water vs direct binding of the phosphate to the ammine and platinum, and to guide the selection of DNA conformations to model complex formation. Based on QM and MM calculations, models are proposed that may be characterized by several general features. A structure involving hydrogen bonding between each ammine and distinct adjacent phosphate groups, referred to as closed conformation (CC), has already been reported. This is also found in the crystal structure of small dimers. We report alternative conformations that may be important in platination of duplex DNA. They are characterized by an intermediate conformation (IC), involving hydrogen bonding between one ammine and phosphate group, and an open conformation (OC), without ammine phosphate hydrogen bonding. The IC and OC can be stabilized by water bridges in the space between the ammine and the phosphate groups. Sugar puckers alternate from the type C(2')-endo or C(1')-exo (S), to the type C(3')-endo or C(2')-exo (N), with intermediate types near O(1')-endo (O). In general, the sugar puckers alternate from S to N to S through the platinated region (3'-TpG*pG*p-5'), with the complexed strand exhibiting, (3')-S*-N*-S-(5') alternation, while the complementary strand shows either (3')-S*-N*-S-(5') or (3')-S*-N*-O-(5') alternation. In both the OC and IC, a hydrogen bond is found between the ammine and O4(T) on thymine (T) at the (3') end, adjacent to the complex site. There is a continuous range of backbone conformations through the platinated region which relate the OC to the IC. The models presented suggest that the dynamics of the binding of the cis-Pt(II)-diammines to adjacent N7(G) in double-stranded DNA may encompass several conformational possibilities, and that water bridges may play a roll in supporting open and intermediate conformations. Proton-proton distances are reported to assist in the experimental determination of conformations.     

The utility of side-chain cyclization in determining the receptor-bound conformation of peptides: cyclic tripeptides and angiotensin II.

The effect of side-chain cyclization on accessible backbone conformations of tripeptides, X-Ala-Y (X and/or Y = Cys, Hcy (Hcy: homocysteine), cis 4-mercaptoproline (MPc), and trans 4-mercaptoproline (MPt)), was elucidated using two variants of systematic conformational search. In addition to cyclization through a disulfide bond, the thioether (-S-CH2-) and amide (-CO-NH-) side-chain analogues of Cys-Ala-Cys and Hcy-Ala-Hcy were evaluated. The number of valid backbone conformations and the allowed phi, psi space were evaluated for each compound, and the ability of the cyclic tripeptides to accommodate beta-turn conformations was examined in order to assess the value of cyclization in limiting conformational freedom. Based on the number of conformations, cyclization was highly effective in reducing the backbone degree of freedom: in order of decreasing number of conformations, Ala-Ala-Ala 1 > Hcy-Ala-Hcy 2 > Cys-Ala-Hcy 3 approximately equal to Hcy-Ala-Cys 4 > MPc-Ala-Hcy 5, 7 > Cys-Ala-Cys 6 > MPc-Ala-Cys 8 > Hcy-Ala-MPt 9 > Cys-Ala-MPt 10 approximately equal to MPc-Ala-MPt 11. Although Hcy-Ala-Hcy 2 had the greatest number of conformations of the cyclic peptides studied, it was still greatly constrained relative to its linear analogue 1. The bicyclic ring system introduced by MP was even more effective in constraining the cycle, having greater impact at position 3 than at position 1. Under the conditions of the study, cyclization of MP-containing analogues could be effected only with the cis isomer (MPc) at position 1 and/or the trans isomer (MPt) at position 3. Sterically allowed conformations of Ala2 for the cyclic tripeptides 2-4 were generally similar to those of the linear tripeptide 1, while those of Cys-Ala-Cys 6 and MPc-Ala-Hcy 7 were restricted to a smaller region of phi 2, psi 2 space: the right- and left-handed alpha-helical conformation and the beta-conformation. This trend was even more pronounced for Hcy-Ala-MPt 9, Cys-Ala-MPt 10, and MPc-Ala-MPt 11, in which Ala2 was severely restricted to a very small region of phi, psi space: the left-handed alpha-helical conformation for 9-11, plus the beta conformation for 9. This suggests that MP at the 3-position is incompatible with a right-handed alpha-helical conformation at position 2.(ABSTRACT TRUNCATED AT 400 WORDS)     

Synthesis and conformation of four 16,17-diols in the 3-methoxy-13alpha-estra-1,3,5(10)-triene series

All four diasteromeric 16,17-diols in the 3-methoxy-13alpha-estra-1,3,5(10)-triene series have been synthesized. The trans-diols 1 and 2 can be obtained by hydroborating the 17-enol acetate 6 (61%, ratio 27:73, preferred alpha attack). OsO(4) dihydroxylation of the olefin 7 yielded the cis-diols 3 and 4 (ratio 13:87). The dihydroxylation proceeds with preference for beta attack caused by a C-ring twist-boat form of 7. The conformations of the diols 2 and 4, the 17-benzyl-17-hydroxy compounds 9 and 10 (obtained by Grignard reaction), and the 16alpha-bromo-17beta-hydroxy compound 8 were determined by X-ray analysis and by 1H NMR spectroscopy in solution. Some compounds, in spite of a 17beta-hydroxy group, had a conformation with a ring C chair form (4, 8, 9) caused by intermolecular interaction in the solid state. The rest of the compounds studied here (2, 10) possessed a conformation with a ring C twist-boat form, which has been also found for all 17beta-substituted compounds in solution. The preferred conformation of the D-ring with 17beta-substituents seems to be the 16beta-envelope form or near this form, but the existence of the 16alpha-envelope form (inversion of the ring D) for some compounds showed great variance in the conformation of ring D, which is substituent dependent.     

Solution conformation of a pectin fragment disaccharide using molecular modelling and nuclear magnetic resonance

In the present study, the conformational behaviour of methylated pectic disaccharide 4-O-alpha-D-galactopyranurosyl 1-O-methyl-alpha-D-galactopyranuronic 6,6'-dimethyl diester 1 has been completely characterized through combined n.m.r. and molecular modelling studies. The 1H-1H n.O.e. across the glycosidic bond was measured by both steady-state and transient 1D and 2D experiments. In parallel, the complete conformational analysis of the disaccharide has been achieved with the MM3 molecular mechanics method. The conformation of the pyranose ring is confirmed by the excellent agreement between the experimental and calculated intracyclic scalar coupling constants. The iso-energy contours displayed on the 'relaxed' map indicate an important flexibility about the glycosidic linkage. There is no significant influence of the methoxyl group on the conformational behaviour of the disaccharide. The theoretical n.m.r. data were calculated taking into account all the accessible conformations and using the averaging methods appropriate for slow internal motions. 3JC-H coupling constants were calculated using an equation suitable for C-O-C-H segments. The agreement between experimental and theoretical data is excellent. Within the potential energy surface calculated for the disaccharide, several conformers can be identified. When these conformations are extrapolated to a regular polymer structure, they generate pectins with right- and left-handed chirality along with a two-fold helix. These different types of helical structure are the result of small changes in conformation, without any drastic variation of the fibre repeat.     

The use of a genetic algorithm search for molecular mechanics (MM3)-based conformational analysis of oligosaccharides

We have implemented a system called glygal that can perform conformational searches on oligosaccharides using several different genetic algorithm (GA) search methods. The searches are performed in the torsion angle conformational space, considering both the primary glycosidic linkages as well as the pendant groups (C-5-C-6 and hydroxyl groups) where energy calculations are performed using the MM3(96) force field. The system includes a graphical user interface for setting calculation parameters and incorporates a 3D molecular viewer. The system was tested using dozens of structures and we present two case studies for two previously investigated O-specific oligosaccharides of the Shigella dysenteriae type 2 and 4. The results obtained using glygal show a significant reduction in the number of structures that need to be sampled in order to find the best conformation, as compared to filtered systematic search.     

Conformational preferences of 1,4,7-trithiacyclononane: a molecular mechanics and density functional theory study

Conformational preferences of 1,4,7-trithiacyclononane were studied using a highly efficient sampling technique based on local nonstochastic deformations and the MM2(91) force field. The results show that conformers that the molecule adopts in the crystal state were found to be low-energy conformers (LECs) within 5 kcal mol(-1) of the global minimum. A conformation with C1 symmetry was the global minimum and the C3 and C2 conformations were calculated to be 0.03 and 1.78 kcal mol(-1) higher in energy, respectively. The structures were further minimized using Density Functional Theory (DFT) calculations with two different functionals. The C2 and the C1 conformations were found to be LECs with the C3 conformation more than 4.0 kcal mol(-1) above the global minimum. The relative energies and structural ordering obtained using the BP86 functional are in agreement with the previously reported relative energies calculated using second-order Moller-Plesset (MP2) ab initio calculations. With the energy ordering being dependent on the molecular mechanics force field used, the approach of MM-->DFT (searching exhaustively the available conformational space at the MM level followed by generating the energy ordering through DFT calculations) appears to be appropriate for thiacrown ethers.     

Structural Factors that Distinguish Dopamine D1 and D2 Agonists

To determine the structural features responsible for their selectivity as dopamine D1 agonists, a conformational analysis has been performed on an analog of nomifensine, dihydrexidine, a benzergoline, and an isochroman using the MM2-87 program. The preferred three dimensional structure of the hydroxylated phenyl ring of the nomifensine analog was found to differ from the other compounds with a substantial energy barrier to achieving the planar conformation of the other compounds which may explain its weak potency for D1 receptors. The preferred three dimensional structures of dihydrexidine and the benzergoline were found to differ significantly despite their molecular similarity. These conformational differences were also evident in crystal structures of the compounds or their analogs. The hypothesis that an equatorial ammonium hydrogen (or amine lone pair) is required for D1 agonist selectivity was tested by performing calculations on N-methyl equatorial and N-methyl axial analogs of the compounds. Calculations were also performed on nonselective dopamine agonists (apomorphine and 5,6-diOH- and 6,7-diOH aminotetralin) and dopamine D2-selective agonists ((+)-PHNO and an analog of quinpirole). The energy difference for the N-methyl axial conformations (or their equivalent) were found to be relatively small for the nonselective agonists and more substantial for the D2-selective agonists. This suggests that D2-selectivity may be associated with the relative unfavorability of the N-methyl axial conformation and provides an explanation for the decreased potency of tertiary amine analogs of the D1-selective agonists. In the benzergoline, where the energy difference is computed to be smaller, the addition of the N-methyl group appears to have a smaller deleterious effect on D1 activity. An N-methyl axial conformation has also been observed for the benzergoline in the crystal state suggesting that this conformation is energetically accessible.     

Stereochemical effects of non-ionic methylphosphonates on nucleotide conformations.

The preferred conformations of deoxyribo and ribonucleoside 3'-methylphosphonates are analysed by minimizing the conformational energy as a function of all the major parameters including the sugar ring for both the S- and R-isomers. The results show that neither the substitution nor the nature of the diastereomer affects significantly the preferred conformations compared to the naturally occurring nucleoside 3'-phosphates. The preferred range of C3'-O3' bond torsions or the phase angles of pseudorotation (P) of the sugar are unaffected. The chiral substitution on the phosphate always adopts a conformation distal to the secondary C3' carbon atom in the minimum energy conformational state. Further, it introduces certain restrictions on the preferred range of P-O3' torsions depending on the methylphosphonate configuration. Methylphosphonate, especially the S-isomer, renders the normal gauche- range of P-O3' bond torsions responsible for the stacked helical duplexes to be energetically unfavourable besides introducing a high energy barrier between trans and gauche conformations. Therefore it is suggested that duplexes with S-methylphosphonate may favour extended phosphodiester conformations. These factors explain the observed lower melting temperature as well as the downfield shifts in the 31P signals in duplexes containing the S-isomer.     

Prediction of the native conformation of a polypeptide by a statistical-mechanical procedure. III. Probable and average conformations of enkephalin

The program SMAPPS (Statistical-Mechanical Algorithm for Predicting Protein Structure) was originally designed to determine the probable and average backbone (?, ψ) conformations of a polypeptide by the application of equilibrium statistical mechanics in conjunction with an adaptive importance sampling Monte Carlo procedure. In the present paper, the algorithm has been extended to include the variation of all side-chain (χ) and peptide-bond (ω) dihedral angles of a polypeptide during the Monte Carlo search of the conformational space. To test the effectiveness of the generalized algorithm, SMAPPS was used to calculate the probable and average conformations of Met-enkephalin for which all dihedral angles of the pentapeptide were allowed to vary. The total conformational energy for each randomly generated structure of Met-enkephalin was obtained by summing over the interaction energies of all pairs of nonbonded atoms of the whole molecule. The interaction energies were computed by the program ECEPP /2 (Empirical Conformational Energy Program for Peptides). Solvent effects were not included in the computation. The results of the Monte Carlo calculation of the structure of Met-enkephalin indicate that the thermodynamically preferred conformation of the pentapeptide contains a γ-turn involving the three residues Gly²-Gly³-Phe⁴. The γ-turn conformation, however, does not correspond to the structure of lowest conformational energy. Rather, the global minimum-energy conformation, recently determined by a new optimization technique developed in this laboratory, contains a type II′ β-bend that is formed by the interaction of the four residues Gly²-Gly³-Phe⁴-Met⁵. A similar minimum-energy conformation is found by the SMAPPS procedure. The thermodynamically preferred γ-turn structure has a conformational energy of 4.93 kcal/mole higher than the β-bend structure of lowest energy but, because of the inclusion of entropy in the SMAPPS procedure, it is estimated to be ～ 9 kcal/mole lower in free energy. The calculation of the average conformation of Met-enkephalin was repeated until a total of ten independent average conformations were established. As far as the phenylalanine residue of the pentapeptide is concerned, the results of the ten independent average conformations were all found to lie in the region of conformational space corresponding to the γ-turn. These results further support the conclusion that the γturn conformation is thermodynamically favored.     

Computational studies of sialyllactones: methods and uses

N-Acetylneuraminic acid (1) is a common sugar in many biological recognition processes. Neuraminidase enzymes recognize and cleave terminal sialic acids from cell surfaces. Viral entry into host cells requires neuraminidase activity, thus inhibition of neuraminidase is a useful strategy for development of drugs for viral infections. A recent crystal structure for influenza viral neuraminidase with sialic acid bound shows that the sialic acid is in a boat conformation [Prot Struct Funct Genet 14: 327 (1992)]. Our studies seek to determine if structural pre-organization can be achieved through the use of sialyllactones. Determination of whether siallylactones are pre-organized in a binding conformation requires conformational analysis. Our inability to find a systematic study comparing the results obtained by various computational methods for carbohydrate modeling led us to compare two different conformational analysis techniques, four different force fields, and three different solvent models. The computational models were compared based on their ability to reproduce experimental coupling constants for sialic acid, sialyl-1,4-lactone, and sialyl-1,7-lactone derivatives. This study has shown that the MM3 forcefield using the implicit solvent model for water implemented in Macromodel best reproduces the experimental coupling constants. The low-energy conformations generated by this combination of computational methods are pre-organized toward conformations which fit well into the active site of neuraminidase. This revised version was published online in November 2006 with corrections to the Cover Date.     

Conformation-activity relationship of peptide T and new pseudocyclic hexapeptide analogs.

Peptide T (ASTTTNYT), a segment corresponding to residues 185-192 of gp120, the coat protein of HIV, has several important biological properties in vitro that have stimulated the search for simpler and possibly more active analogs. We have previously shown that pseudocyclic hexapeptide analogs containing the central residues of peptide T retain considerable chemotactic activity. We have now extended the design of this type of analogs to peptides containing different aromatic residues and/or Ser in lieu of Thr. The complex conformation-activity relationship of these analogs called for a reexamination of the basic conformational tendencies of peptide T itself. Here, we present an exhaustive NMR conformational study of peptide T in different media. Peptide T assumes a gamma-turn in aqueous mixtures of ethylene glycol, a type-IV beta-turn conformation in aqueous mixtures of DMF, and a type-II beta-turn conformation in aqueous mixtures of DMSO. The preferred conformations for the analogs were derived from modeling, starting from the preferred conformations of peptide T. The best models derived from the gamma-turn conformation of peptide T are those of peptides XII (DSNYSR), XIII (ETNYTK) and XVI (ESNYSR). The best models derived from the type-IV beta-turn conformation of peptide T are those of peptides XIV (KTTNYE) and XV (DSSNYR). No low-energy models could be derived starting from the type-II beta-turn conformation of peptide T. The analogs with the most favored conformations are also the most active in the chemotactic test.