Noncontact dipole effects on channel permeation. II. Trp conformations and dipole potentials in gramicidin A

The four Trp dipoles in the gramicidin A (gA) channel modulate channel conductance, and their side chain conformations should therefore be important, but the energies of different conformations are unknown. A conformational search for the right-handed helix based on molecular mechanics in vacuo yielded 46 conformations within 20 kcal/mol of the lowest energy conformation. The two lowest energy conformations correspond to the solid-state and solution-state NMR conformations, suggesting that interactions within the peptide determine the conformation. For representative conformations, the electrostatic potential of the Trp side chains on the channel axis was computed. A novel application of the image-series method of. Biophys. J. 9:1160-1170) was introduced to simulate the polarization of bulk water by the Trp side chains. For the experimentally observed structures, the CHARm toph19 potential energy (PE) of a cation in the channel center is -1.65 kcal/mol without images. With images, the PE is -1.9 kcal/mol, demonstrating that the images further enhance the direct dipole effect. Nonstandard conformations yielded less favorable PEs by 0.4-1.1 kcal/mol.     

Structure-fluorescence correlations in a single tryptophan mutant of carp parvalbumin: solution structure, backbone and side-chain dynamics

Heterogeneous fluorescence intensity decays of tryptophan in proteins are often rationalized using a model which proposes that different rotameric states of the indole alanyl side-chain are responsible for the observed fluorescence lifetime heterogeneity. We present here the study of a mutant of carp parvalbumin bearing a single tryptophan residue at position 102 (F102W) whose fluorescence intensity decay is heterogeneous and assess the applicability of a rotamer model to describe the fluorescence decay data. We have determined the solution structure of F102W in the calcium ligated state using multi-dimensional nuclear magnetic resonance (NMR) and have used the minimum perturbation mapping technique to explore the possible existence of multiple conformations of the indole moiety of Trp102 of F102W and, for comparison, Trp48 of holo-azurin. The maps for parvalbumin suggest two potential conformations of the indole side-chain. The high energy barrier for rotational isomerization between these conformers implies that interwell rotation would occur on time-scales of milliseconds or greater and suggests a rotamer basis for the heterogeneous fluorescence. However, the absence of alternate Trp102 conformers in the NMR data (to within 3 % of the dominant species) suggests that the heterogeneous fluorescence of Trp102 may arise from mechanisms independent of rotameric states of the Trp side-chain. The map for holo-azurin has only one conformation, and suggests a rotamer model may not be required to explain its heterogeneous fluorescence intensity decay. The backbone and Trp102 side-chain dynamics at 30 degrees C of F102W has been characterized based on an analysis of (15)N NMR relaxation data which we have interpreted using the Lipari-Szabo formalism. High order parameter (S(2)) values were obtained for both the helical and loop regions. Additionally, the S(2) values imply that the calcium binding CD and EF loops are not strictly equivalent. The S(2) value for the indole side-chain of Trp102 obtained from the fluorescence, NMR relaxation and minimum perturbation data are consistent with a Trp moiety whose motion is restricted.     

β‐Hairpin folding and stability: molecular dynamics simulations of designed peptides in aqueous solution

The structural properties of a 10‐residue and a 15‐residue peptide in aqueous solution were investigated by molecular dynamics simulation. The two designed peptides, SYINSDGTWT and SESYINSDGTWTVTE, had been studied previously by NMR at 278 K and the resulting model structures were classified as 3:5 β‐hairpins with a type I + G1 β‐bulge turn. In simulations at 278 K, starting from the NMR model structure, the 3:5 β‐hairpin conformers proved to be stable over the time period evaluated (30 ns). Starting from an extended conformation, simulations of the decapeptide at 278 K, 323 K and 353 K were also performed to study folding. Over the relatively short time scales explored (30 ns at 278 K and 323 K, 56 ns at 353 K), folding to the 3:5 β‐hairpin could only be observed at 353 K. At this temperature, the collapse to β‐hairpin‐like conformations is very fast. The conformational space accessible to the peptide is entirely dominated by loop structures with different degrees of β‐hairpin character. The transitions between different types of ordered loops and β‐hairpins occur through two unstructured loop conformations stabilized by a single side‐chain interaction between Tyr2 and Trp9, which facilitates the changes of the hydrogen‐bond register. In agreement with previous experimental results, β‐hairpin formation is initially driven by the bending propensity of the turn segment. Nevertheless, the fine organization of the turn region appears to be a late event in the folding process. Copyright © 2004 European Peptide Society and John Wiley & Sons, Ltd.     

Conformationally restricted mimics of vitamin D rotamers

Drug developments in the vitamin D field have continued to focus on structure-function studies of analogs produced by chemically modifying the structure of 1alpha,25-dihydroxyvitamin D(3) (1,25-D3) and its metabolites. Direct structural information gleaned from X-ray crystallographic or NMR studies regarding the ligand-receptor complex and other guest-host systems, which are likely involved in initiating biologic responses, also offers potential insight into drug design. Evidence has accrued suggesting that topologically different conformers of 1,25-D3 may bind to proteins in different ways, including the induction of different conformations of protein. This paper concerns our progress on the chemical synthesis of analogs (e.g. ansa-steroids, suprasterols, vinylallenes and other analogs) conformationally locked or at least rotationally restricted to mimic higher energy conformers of 1,25-D3.     

Conformational effects of C(alpha,alpha)-dipropargylglycine as a constrained residue

A useful synthon to approach artificial phenylalanyl peptides in a [2 + 2 + 2] cycloaddition reaction, C(alpha,alpha)-dipropargylglycine (Dprg) is examined for its conformational preferences as a constrained residue. Crystal structure analysis and preliminary NMR results establish possible preference of the residue for folded (alpha) rather than extended (beta) region of the straight phi,psi conformational space. Boc-Dprg-L-Leu-OMe (1) displays two molecular conformations within the same crystallographic asymmetric unit, with Dprg in the alpha(R) or alpha(L) conformation, participating in a type I beta-turn or an alpha(L)-alpha(R)-type fold, in which Leu(2) assumes the alpha(R) conformation stereochemically favored for an L-chiral residue. Boc-Dprg-D-Val-L-Leu-OMe (2) displays a type I' beta-turn conformation in crystal, with both Dprg(1) and D-Val(2) assuming the alpha(L) conformation stereochemically favored for a D-chiral residue, with 4 --> 1 type hydrogen bond linking L-Leu(3) NH with Boc CO. NMR analysis using temperature variation, solvent titration, and a spin probe study suggests a fully solvent-exposed nature of Dprg NH, ruling out a fully extended C(5)-type conformation for this residue, and solvent sequestered nature of L-Leu(3) NH, suggesting possibility of a beta-turn due to Dprg assuming a folded conformation.     

Refinement of the spatial structure of the gramicidin A ion channel]

The spatial structure of the gramicidin A (GA) transmembrane ion-channel was refined on the base of cross-peak volumes measured in NOESY spectra (mixing time tau m = 100 and 200 ms). The refinement methods included the comparison of experimental cross-peak volumes with those calculated for low-energy GA conformations, dynamic averaging of the low-energy conformation set and restrained energy minimization. Accuracy of the spatial structure determination was estimated by the penalty function Fr defined as a root mean square deviation of interproton distances corresponding to the calculated and experimental cross-peak volumes. As the initial conformation we used the right-handed pi 6,3 LD pi 6,3 LD helix established on the base of NMR data regardless of the cross-peak volumes. The conformation is in a good agreement with NOE cross-peak volumes (Fr 0.2 to 0.5 A depending on NOESY spectrum). For a number of NOEs formed by the side chain protons, distances errors were found as much as 0.5-2.0 A. Restrained energy minimization procedure had little further success. However some of these errors were eliminated by the change in torsional angle chi 2 of D-Leu12 and dynamic averaging of the Val7 side chain conformations. Apparently, majority of deviations of the calculated and experimental cross-peak volumes are due to the intramolecular mobility of GA and cannot be eliminated within the framework of rigid globule model. In summary the spatial structure of GA ion-channel can be thought as a set of low-energy conformations, differing by the side chain torsion angles chi 1 Val7 and chi 2 D-Leu4 and D-Leu10 and the orientation of the C-terminal ethanolamine group. Root mean square differences between the atomic coordinates of conformations are in the range of 0.3-0.8 A.     

Empirical solvation models in the context of conformational energy searches: application to bovine pancreatic trypsin inhibitor.

Continuum solvation models that estimate free energies of solvation as a function of solvent accessible surface area are computationally simple enough to be useful for predicting protein conformation. The behavior of three such solvation models has been examined by applying them to the minimization of the conformational energy of bovine pancreatic trypsin inhibitor. The models differ only with regard to how the constants of proportionality between free energy and surface area were derived. Each model was derived by fitting to experimentally measured equilibrium solution properties. For two models, the solution property was free energy of hydration. For the third, the property was NMR coupling constants. The purpose of this study is to determine the effect of applying these solvation models to the nonequilibrium conformations of a protein arising in the course of global searches for conformational energy minima. Two approaches were used: (1) local energy minimization of an ensemble of conformations similar to the equilibrium conformation and (2) global search trajectories using Monte Carlo plus minimization starting from a single conformation similar to the equilibrium conformation. For the two models derived from free energy measurements, it was found that both the global searches and local minimizations yielded conformations more similar to the X-ray crystallographic structures than did searches or local minimizations carried out in the absence of a solvation component of the conformational energy. The model derived from NMR coupling constants behaved similarly to the other models in the context of a global search trajectory. For one of the models derived from measured free energies of hydration, it was found that minimization of an ensemble of near-equilibrium conformations yielded a new ensemble in which the conformation most similar to the X-ray determined structure PTI4 had the lowest total free energy. Despite the simplicity of the continuum solvation models, the final conformation generated in the trajectories for each of the models exhibited some of the characteristics that have been reported for conformations obtained from molecular dynamics simulations in the presence of a bath of explicit water molecules. They have smaller root mean square (rms) deviations from the experimentally determined conformation, fewer incorrect hydrogen bonds, and slightly larger radii of gyration than do conformations derived from search trajectories carried out in the absence of solvent.     

Protein side chain conformation predictions with an MMGBSA energy function

The prediction of protein side chain conformations from backbone coordinates is an important task in structural biology, with applications in structure prediction and protein design. It is a difficult problem due to its combinatorial nature. We study the performance of an “MMGBSA” energy function, implemented in our protein design program Proteus, which combines molecular mechanics terms, a Generalized Born and Surface Area (GBSA) solvent model, with approximations that make the model pairwise additive. Proteus is not a competitor to specialized side chain prediction programs due to its cost, but it allows protein design applications, where side chain prediction is an important step and MMGBSA an effective energy model. We predict the side chain conformations for 18 proteins. The side chains are first predicted individually, with the rest of the protein in its crystallographic conformation. Next, all side chains are predicted together. The contributions of individual energy terms are evaluated and various parameterizations are compared. We find that the GB and SA terms, with an appropriate choice of the dielectric constant and surface energy coefficients, are beneficial for single side chain predictions. For the prediction of all side chains, however, errors due to the pairwise additive approximation overcome the improvement brought by these terms. We also show the crucial contribution of side chain minimization to alleviate the rigid rotamer approximation. Even without GB and SA terms, we obtain accuracies comparable to SCWRL4, a specialized side chain prediction program. In particular, we obtain a better RMSD than SCWRL4 for core residues (at a higher cost), despite our simpler rotamer library. Proteins 2016; 84:803–819. © 2016 Wiley Periodicals, Inc.     

Backbone dependency further improves side chain prediction efficiency in the Energy‐based Conformer Library (bEBL)

Side chain optimization is an integral component of many protein modeling applications. In these applications, the conformational freedom of the side chains is often explored using libraries of discrete, frequently occurring conformations. Because side chain optimization can pose a computationally intensive combinatorial problem, the nature of these conformer libraries is important for ensuring efficiency and accuracy in side chain prediction. We have previously developed an innovative method to create a conformer library with enhanced performance. The Energy‐based Library (EBL) was obtained by analyzing the energetic interactions between conformers and a large number of natural protein environments from crystal structures. This process guided the selection of conformers with the highest propensity to fit into spaces that should accommodate a side chain. Because the method requires a large crystallographic data‐set, the EBL was created in a backbone‐independent fashion. However, it is well established that side chain conformation is strongly dependent on the local backbone geometry, and that backbone‐dependent libraries are more efficient in side chain optimization. Here we present the backbone‐dependent EBL (bEBL), whose conformers are independently sorted for each populated region of Ramachandran space. The resulting library closely mirrors the local backbone‐dependent distribution of side chain conformation. Compared to the EBL, we demonstrate that the bEBL uses fewer conformers to produce similar side chain prediction outcomes, thus further improving performance with respect to the already efficient backbone‐independent version of the library. Proteins 2014; 82:3177–3187. © 2014 Wiley Periodicals, Inc.     

Synthesis and structural analyses of 3-acetamido-1,4-di-O-acetyl-2,3,5-trideoxy-5-C-(isopropylphosphinyl)-D-erythro-pentopyranoses

1H NMR spectroscopy of phosphorus containing hetero sugars (phospha sugars), revealed the alpha and beta configurations and chair conformations for 3-acetamido-1,4-di-O-acetyl-2,3,5-trideoxy-5-C-(isopropylphosphinyl)-alpha- and beta-D-erythro-pentopyranoses. The conformation of the title compounds was determined by 1H NMR as 1C4 in CDCl3 and the conformation was in accord with that in solid state determined by X-ray crystallographic analysis.     

NMR and molecular modeling characterization of RGD containing peptides.

The tripeptide sequence arginine-glycine-aspartic acid (RGD) has been shown to be the key recognition segment in numerous cell adhesion proteins. The solution conformation and dynamics in DMSO-d6 of the cyclic pentapeptides, [formula: see text], a potent fibrinogen receptor antagonist, and [formula: see text], a weak fibrinogen receptor antagonist, have been characterized by nuclear magnetic resonance (NMR) spectroscopy and molecular modeling. 1H-1H distance constraints derived from two-dimensional NOE spectroscopy and torsional angle constraints obtained from 3JNH-H alpha coupling constants, combined with computer-assisted modeling using conformational searching algorithms and energy minimization have allowed several low energy conformations of the peptides to be determined. Low temperature studies in combination with molecular dynamics simulations suggest that each peptide does not exist in a single, well-defined conformation, but as an equilibrating mixture of conformers in fast exchange on the NMR timescale. The experimental results can be fit by considering pairs of low energy conformers. Despite this inherent flexibility, distinct conformational preferences were found which may be related to the biological activity of the peptides.     

B3LYP/6-311++G** study of alpha- and beta-D-glucopyranose and 1,5-anhydro-D-glucitol: 4C1 and 1C4 chairs, (3,O)B and B(3,O) boats, and skew-boat conformations

Geometry optimization, at the B3LYP/6-311++G** level of theory, was carried out on 4C1 and 1C4 chairs, (3,O)B and B(3,O) boats, and skew-boat conformations of alpha- and beta-D-glucopyranose. Similar calculations on 1,5-anhydro-D-glucitol allowed examination of the effect of removal of the 1-hydroxy group on the energy preference of the hydroxymethyl rotamers. Stable minimum energy boat conformers of glucose were found, as were stable skew boats, all having energies ranging from approximately 4-15 kcal/mol above the global energy 4C1 chair conformation. The 1C4 chair electronic energies were approximately 5-10 kcal/mol higher than the 4C1 chair, with the 1C4 alpha-anomers being lower in energy than the beta-anomers. Zero-point energy, enthalpy, entropy, and relative Gibbs free energies are reported at the harmonic level of theory. The alpha-anomer 4C1 chair conformations were found to be approximately 1 kcal/mol lower in electronic energy than the beta-anomers. The hydroxymethyl gt conformation was of lowest electronic energy for both the alpha- and beta-anomers. The glucose alpha/beta anomer ratio calculated from the relative free energies is 63/37%. From a numerical Hessian calculation, the tg conformations were found to be approximately 0.4-0.7 kcal/mol higher in relative free energy than the gg or gt conformers. Transition-state barriers to rotation about the C-5-C-6 bond were calculated for each glucose anomer with resulting barriers to rotation of approximately 3.7-5.8 kcal/mol. No energy barrier was found for the path between the alpha-gt and alpha-gg B(3,O) boat forms and the equivalent 4C1 chair conformations. The alpha-tg conformation has an energy minimum in the 1S3 twist form. Other boat and skew-boat forms are described. The beta-anomer boats retained their starting conformations, with the exception of the beta-tg-(3,O)B boat that moved to a skew form upon optimization.     

Slow conformational dynamics of an endonuclease persist in its complex with its natural protein inhibitor

The bacterial toxin colicin E9 is secreted by producing Escherichia coli cells with its 9.5 kDa inhibitor protein Im9 bound tightly to its 14.5 kDa C-terminal DNase domain. Double- and triple-resonance NMR spectra of the isolated DNase domain uniformly labeled with 13C/15N bound to unlabeled Im9 contain more signals than expected for a single DNase conformer, consistent with the bound DNase being present in more than one form. The presence of chemical exchange cross peaks in 750 MHz 15N-1H-15N HSQC-NOESY-HSQC spectra for backbone NH groups of Asp20, Lys21, Trp22, Leu23, Lys69, and Asn70 showed that the bound DNase was in dynamic exchange. The rate of exchange from the major to the minor form was determined to be 1.1 +/- 0.2 s(-1) at 298 K. Previous NMR studies have shown that the free DNase interchanges between two conformers with a forward rate constant of 1.61 +/- 0.11 s(-1) at 288 K, and that the bound Im9 is fixed in one conformation. The NMR studies of the bound DNase show that Im9 binds similarly to both conformers of the DNase and that the buried Trp22 is involved in the dynamic process. For the free DNase, all NH groups within a 9 A radius of any point of the Trp22 ring exhibit heterogeneity suggesting that a rearrangement of the position of this side chain is connected with the conformational interchange. The possible functional significance of this feature of the DNase is discussed.     

The effect of the amino-acid side chains on the energy profiles for ion transport in the gramicidin A channel

Computations on the energy profiles for Na+ in the gramicidin A (GA) channel have been extended by introducing the effect, previously neglected, of the amino acid side chains of GA, fixed in their most stable conformations. The calculations have been performed in two approximations: 1) with the ethanolamine tail fixed in its most stable conformation, 2) with the tail allowed to optimize its conformation upon the progression of the ion. In both approximations the overall shape of the energy profile is very similar to that obtained in the absence of the side chains. One observes, however, a general lowering of the profile upon the adjunction of the side chains. The analysis of the factors responsible for this energy lowering indicates that it is due essentially to the electrostatic and polarisation components of the interaction which interplay differently, however, in the different parts of the channel. A particular role is attributed in this respect to the tryptophan residues of GA. The role of the 4 tryptophans present, Trp 15, 13, 11 and 9, is individualized by stripping of one of them at a time. The strongest effect on the energy deepening is due to Trp 13 and is particularly prominent in the entrance zone at 14.5A from the center of the channel. The result indicates the possibility of investigating theoretically the effect on the energy profiles of the substitution of the "natural" side chain by others.     

Conformational solution studies of neuropeptide gamma using CD and NMR spectroscopy.

Neuropeptide gamma is one of the largest members of the tachykinin family of peptides, exhibiting strong agonistic activity towards the NK-2 tachykinin receptor. This peptide was synthesized by the solid-phase method using the Fmoc chemistry. Circular-dichroism spectroscopy (CD) investigations of this peptide were performed in phosphate buffer, in the presence of sodium dodecylsulphate (SDS) micelles and trifluoroethanol (TFE) solutions and in DMSO-d6 using the 2D NMR technique in conjunction with two different theoretical approaches. The first assumes multiconformational equilibrium of the peptide studied characterized by the values of statistical weights of low-energy conformations. These calculations were performed using three different force fields ECEPP/3, AMBER4.1 and CHARMM (implemented in the X-PLOR program). The second method incorporates interproton distance and dihedral angle constraints into the starting conformation using the Simulated Annealing algorithm (X-PLOR program). The CD experiments revealed that although the peptide studied is flexible in polar solvents, a tendency to adopt a helical structure was observed in the hydrophobic environment. The NMR data (NOE effects) indicate a helical or reverse structure in the Ile7-His12 fragment of the peptide studied in DMSO-d6 solution. The results obtained cannot be interpreted in terms of a single conformation. Most of the conformations obtained with the ECEPP/3 force field possess a high content of a helical structure. None of the conformers, obtained with the AMBER4.1 and CHARMM force fields, can be considered as the dominant one. In all conformations several beta-turns were detected and in some cases gamma-turns were also found. But in fact, it is rather difficult to select the position of the secondary element(s) present in the structure of NPgamma in solution. All conformers calculated with the X-PLOR program (with using NMR derived distance and torsion angle constraints) are stabilized by several beta-turns. Common structural motives are a type IV beta-turn in the Gln6-His12 fragment. All conformations obtained using two approaches adopt very similar turn shapes in the middle region of molecule and a random structure on the N- and C-terminal fragments.     

Conformational and structural analysis of the equilibrium between single- and double-strand beta-helix of a D,L-alternating oligonorleucine

Alternating sequences of D and L residues in peptides are directly related to the formation of several kinds of regular helical conformations usually called beta-helices. The major feature of these structures is that they can be associated with the transmembrane ion-conducting channel activity in some natural antibacterial peptides. The study of alternating D,L synthetic peptides is critical to understand how factors such as surrounding media, main chain length, type of side chain and terminal groups, among others, can determine the adoption of a specific kind of beta-helix. Early studies pointed out that the peptides Boc-(D-NLeu-L-NLeu)(6)-D-MeNLe-L-Nl-D-Nl-L-Nl-OMe (Boc: tert-butyloxycarbonyl) and Boc-L-Nle-(D-Nle-L-Nle)(5)-D-MeNle-L-Nle-D-Nle-L-Nle-OMe adopt in chloroform a unique detectable conformation single beta(4.4)- and double beta(5.6) upward arrow downward arrow -helix, respectively. The influence of terminal groups on the final stable conformation of N-formylated peptides has been studied in this work. The initial basic NMR data analysis of a synthetic alternating D,L-oligopeptide with ten norleucines, N-methylated on the residue 7 and having HCO- and -OMe as terminal groups clearly indicates the coexistence of two different conformations in equilibrium. NMR data and molecular dynamics calculations point to a dimeric antiparallel beta-helix structure beta(5.6) upward arrow downward arrow for the main conformation. On the other hand, NMR data suggest a single beta-helix structure beta(4.4) for the second conformation. Finally, a thermodynamic analysis of the equilibrium between both conformations has been carried out by one-dimensional NMR measurements at ten different temperatures. The temperature at which 50% of dimer conformation is dissociated is 319 K. In addition, the dimer-monomer equilibrium curve obtained shows a DeltaG>0 for the whole range of studied temperatures, and its behavior can be considered similar to the thermodynamic denaturation protein processes.     

New acylated triterpenoid saponins from Maesa lanceolata

Ten new acylated triterpenoid saponins were isolated from the leaves of Maesa lanceolata. For their structure elucidation extensive use was made of homo- and heteronuclear 2D NMR techniques such as COSY, NOESY, HSQC and HMBC. All saponins identified contained the same tetraglycosidic side chain, but the triterpenoid moiety showed a variable esterification pattern. Monoester, diester and triester derivatives were present. Maesasaponin I was a 21-monoester derivative, i.e. ?3 beta-O-[alpha-L-rhamnopyranosyl-(1-->2)-beta-D-galactopyranosyl- (1-->3)]-[beta-D-galactopyranosyl-(1-->2)]-beta-D-glucuronopyranosyl+ ++?-21 beta-angeloyloxy-13 beta, 28-oxidoolean-16 alpha, 22 alpha, 28 alpha-triol. Maesasaponins III, IV3, V3 and VI2 had an additional acetyl, propanoyl, n-butanoyl and angeloyl substituent, respectively, in position 22. Maesasaponins II, IV2, V2, VI3 and VII1 were characterised as the 16-acetyl derivatives of maesasaponins I, III, IV3, V3 and VI2, respectively. Structures of saponins previously reported in M. lanceolata had to be revised.     

Cloning,overexpression, and purification of aminoglycoside antibiotic 3-acetyltransferase-IIIb: conformational studies with bound substrates

Aminoglycoside 3-acetyltransferase-IIIb (AAC3), which acetylates N3 amine of aminoglycoside antibiotics, was cloned from P. Aeruginosa and purified from overexpressing E. coli BL21 (DE3) cells. Bound conformations of kanamycin A and ribostamycin, in the active site of the enzyme that modifies the essential N3B of aminoglycoside antibiotics, were determined by NMR spectroscopy. Experimentally determined interproton distances were used in a simulated annealing protocol to determine enzyme-bound conformations of both antibiotics. Two conformations, consistent with the NOE restraints, were determined for ribostamycin. The only difference between the two conformers was the orientation of the A ring with respect to the rest of the molecule. The average glycosidic dihedral angles were Phi(1A) = -22 degrees +/- 3 and Psi(1A) = -42 degrees +/- 1 (conformer 1) and Phi(1A) = -67 degrees +/- 0.7 and Phi(1A) = -59 degrees +/- 0.8 (conformer 2). Three conformers were determined for the enzyme-bound kanamycin A. Two conformers of kanamycin A were matched well with the two conformers of ribostamycin when the A and the B rings of the antibiotics were superimposed. Conformations of kanamycin A and ribostamycin were compared to those of other aminoglycosides that are bound to different enzymes and RNA. The results lend further support to our earlier hypothesis that the A and B rings of aminoglycosides adopt a conformation that is recognized not only by the aminoglycoside-modifying enzymes but also by RNA (Serpersu, E. H., Cox, J. R., Digiammarino, E. L., Mohler, M. L., Akal, A., Ekman, D. R., and Owston, M. (2000) Cell Biochem. Biophys. 33, 309-321). These results may be useful in designing new antibiotics to combat the antibiotic resistance against infectious diseases.     

Solution conformation of a model hexapeptide containing RGD sequence.

The solution conformation of a model hexapeptide Asp-Arg-Gly-Asp-Ser-Gly (DRGDSG) containing the RGD sequence has been studied in DMSO-d6 as well as in aqueous solution (H2O:D2O/90:10%) by 1H NMR spectroscopy. The unambiguous identification of spin systems of various amino acid residues and sequence specific assignment of all proton resonances was achieved by a combination of two dimensional COSY and NOESY experiments. The temperature coefficient data of the amide proton chemical shifts in conjunction with the vicinal coupling constants, i.e. 3JNH-C alpha H, NOESY and ROESY results indicate that the peptide in both the solvents exists in a blend of conformers with beta-sheet like extended backbone structure and folded conformations. The folded conformers do not appear to be stabilised by intramolecular hydrogen bonding. Our results are consistent with the flexibility of RGD segment observed in the NMR studies on the protein echistatin containing the RGD motif (references 23-25).     

Comparative molecular dynamics simulations of the potent synthetic classical cannabinoid ligand AMG3 in solution and at binding site of the CB1 and CB2 receptors

The C-1'-dithiolane Delta(8)-tetrahydrocannabinol (Delta(8)-THC) amphiphilic analogue (-)-2-(6a,7,10,10a-tetrahydro-6,6,9-trimethylhydroxy-6H-dibenzo[b,d]pyranyl)-2-hexyl-1,3-dithiolane (AMG3) is considered as one of the most potent synthetic analgesic cannabinoid (CB) ligands. Its structure is characterized by rigid tricyclic and flexible alkyl chain segments. Its conformational properties have not been fully explored. Structure-activity relationship (SAR) studies on classical CBs showed that the alkyl side chain is the most critical structural part for the receptor activation. However, reported low energy conformers of classical CB analogues vary mainly in the conformation of their alkyl side chain segment. Therefore, comparative molecular dynamics (MD) simulations of low energy conformers of AMG3 were performed in order to investigate its structural and dynamical properties in two different systems. System-I includes ligand and amphoteric solvent DMSO, simulating the biological environment and system-II includes ligand at active site of the homology models of CB1 and CB2 receptors in the solvent. The trajectory analysis results are compared for the systems I and II. In system-I, the dihedral angle defined between aromatic ring and dithiolane ring of AMG3 shows more resistance to be transformed into another torsional angle and the dihedral angle adjacent to dithiolane ring belonging in the alkyl chain has flexibility to adopt gauche+/- and trans dihedral angles. The rest of the dihedral angles within the alkyl chain are all trans. These results point out that wrapped conformations are dynamically less favored in solution than linear conformations. Two possible plane angles defined between the rigid and flexible segments are found to be the most favored and adopting values of approximately 90 degrees and approximately 140 degrees. In system-II, these values are approximately 90 degrees and approximately 120 degrees. Conformers of AMG3 at the CB1 receptor favor to establish a cis conformation defined between aromatic and dithiolane ring and a trans conformation in the CB2 receptor. These different orientations of ligand inside the binding pocket of CB1 and CB2 receptors may explain its different binding affinity in the two receptors. The results of this study can be applied to other synthetic classical CB ligands to produce low energy conformations and can be of general use for the molecules possessing flexible alkyl chain(s). In addition, this study can be useful when restraint of the alkyl chain is sought for optimizing drug design.