Determination of productive conformations of acetylcholinesterase substrates using theoretical conformational analysis

All equilibrium conformations of twenty-three acetylcholinesterase effectors were calculated by the molecular mechanics method, nonbonded interactions, torsion energy and energy of bond angles deformation being taken into account. In a series of conformationally flexible derivatives of acetylcholine the correlation was found between hydrolysis rate and population of the completely extended tt-conformation. In a series of cyclic analogues of acetylcholine the high hydrolysis rate occurs only for substrates sterically corresponding to tt-conformation of acetylcholine with regard to disposition of ammonium group, carbonyl oxygen and carbonyl carbon. The hydrolysis rate of acetylcholine derivatives with elongated chain between acetyl and cationic groups is directly proportional to the population of the conformations similar to tt-conformation of acetylcholine. It is concluded that tt-conformation of acetylcholine is productive for acetylcholinesterase hydrolysis.     

Vibrational circular dichroism and IR spectral analysis as a test of theoretical conformational modeling for a cyclic hexapeptide

A model cyclohexapeptide, cyclo-(Phe-(D)Pro-Gly-Arg-Gly-Asp) was synthesized and its IR and VCD spectra were used as a test of density functional theory (DFT) level predictions of spectral intensities for a peptide with a nonrepeating but partially constricted conformation. Peptide structure and flexibility was estimated by molecular dynamics (MD) simulations and the spectra were simulated using full quantum mechanical (QM) approaches for the complete peptide and for simplified models with truncated side chains. After simulated annealing, the backbone conformation of the ring structure is relatively stable, consisting of a normal beta-turn and a tight loop (no H-bond) which does not vary over short trajectories. Only in quite long MD runs at high temperatures do other conformations appear. MD simulations were carried out for the cyclic peptide in water and in TFE, which match experimental solvents, as well as with and without protonation of the Asp carboxyl group. DFT spectral simulations were made using the annealed structure and were extended to include basis set variation, to determine an optimal computational approach, and solvent simulation with a polarized continuum model (PCM). Stepwise full DFT simulation of spectra was done for various sequences with the same backbone geometry but based on (1) solely Gly residues, (2) Ala substitution except Gly and Pro, and (3) complete sequences with side chains. Additionally, a selection of structures was used to compute IR and VCD spectra with the optimal method to determine structural variation effects. The side chains, especially the Asp-COOH and Arg-NH(2) transitions, had an impact on the computed amide frequencies, IR intensities and VCD pattern. Since experimentally these groups would have little chirality, due to conformational variation, they do not impact the observed VCD spectra. Correcting for frequency shifts, the Ala model for the cyclopeptide gives the clearest representation of the amide VCD. The experimental sign pattern for the amide I' band in D(2)O and also the sharper, more intense amide I VCD band in TFE was seen to some degree in one conformer with Type II' turns, but the data favor a mix of structures.     

Crystal structure and conformational analysis of ampullosporin A.

Ampullosporin A is a 15-mer peptaibol type polypeptide that induces pigment formation by the fungus Phoma destructiva, forms voltage-dependent ion channels in membranes and exhibits hypothermic effects in mice. The structure of ampullosporin A has been determined by x-ray crystallography. This is the first three-dimensional (3D) structure of the peptaibol subfamily SF6. From the N-terminus to residue 13 the molecule adopts an approximate right-handed alpha-helical geometry, whereas a less regular structure pattern with beta-turn characteristics is found in the C-terminus. Even though ampullosporin A does not contain a single proline or hydroxyproline it is significantly bent. It belongs to both the shortest and the most strongly bent peptaibol 3D structures. The straight structure part encompasses residues Ac-Trp(1)-Aib(10) and is thus less extended than the alpha-helical subunit. The 3D structure of ampullosporin A is discussed in relation to other experimentally determined peptaibol structures and in the context of its channel-forming properties. As a part of this comparison a novel bending analysis based on a 3D curvilinear axis describing the global structural characteristics has been proposed and applied to all 3D peptaibol structures. A sampling of 2500 conformations using different molecular dynamics protocols yields, for the complete ampullosporin A structure, an alpha-helix as the preferred conformation in vacuo with almost no bend. This indicates that solvent or crystal effects may be important for the experimentally observed peptide backbone bending characteristics of ampullosporin A.     

Boundaries of the butyrylcholinesterase anion site from data of a conformational analysis of substrates

By means of molecular mechanics, theoretical conformational analysis has been made of 19 substrates of butyrylcholinesterase - acetylcholine derivatives with different structure of the ammonium group. It was concluded that the anionic point is located in the cavity of the enzymic molecule. Dimensions and shape of this cavity were established which provide satisfactory correlation between its filling by substrate conformers and the rate of their enzymic hydrolysis. Some suggestions were made with respect to the mechanism of the effect of non-productive sorbtion of the substrates on the rate of their enzymic hydrolysis.     

Solution conformational study of Scyliorhinin I analogues with conformational constraints by two-dimensional NMR and theoretical conformational analysis.

Two analogues of Scyliorhinin I (Scyl), a tachykinin with N-MeLeu in position 8 and a 1,5-disubstituted tetrazole ring between positions 7 and 8, introduced in order to generate local conformational constraints, were synthesized using the solid-phase method. Conformational studies in water and DMSO-d6 were performed on these peptides using a combination of the two-dimensional NMR technique and theoretical conformational analysis. The algorithm of conformational search consisted of the following three stages: (i) extensive global conformational analysis in order to find all low-energy conformations; (ii) calculation of the NOE effects and vicinal coupling constants for each of the low energy conformations; (iii) determining the statistical weights of these conformations by means of a nonlinear least-squares procedure, in order to obtain the best fit of the averaged simulated spectrum to the experimental one. In both solvents the three-dimensional structure of the analogues studied can be interpreted only in terms of an ensemble of multiple conformations. For [MeLeu8]Scyl, the C-terminal 6-10 fragment adopts more rigid structure than the N-terminal one. In the case of the analogue with the tetrazole ring in DMSO-d6 the three-dimensional structure is characterized by two dominant conformers with similar geometry of their backbones. They superimpose especially well (RMSD = 0.28 A) in the 6-9 fragments. All conformers calculated in both solvents superimpose in their C-terminal fragments much better than those of the first analogue. The results obtained indicate that the introduction of the tetrazole ring into the Scyl molecule rigidifies its structure significantly more than that of MeLeu.     

An NMR conformational analysis of cyclic bradykinin mimics. Evidence for a beta-turn

A detailed NMR study is carried out in acetonitrile/water solutions on three novel cyclic bradykinin antagonist analogues, BKM-824, BKM-870, and BKM-872, to examine their solution structures, and to correlate the structures with bradykinin antagonist and anti-cancer activities. The solution structures of the cyclic peptides are correlated with the structural data for known linear bradykinin antagonists. The sequences are: BKM-824 c[Ava-Ig1-Ser-DF5F-Oic-Arg] where Ava is 5-aminovaleric acid, Ig1 is alpha-(2-indanyl)glycine, F5F is pentafluorophenylalanine, and Oic is (2S,3aS,7aS)-octahydroindole-2-carboxylic acid; BKM-870; c[DArg-Arg-Add-DF5F-Oic-Arg] where Add is 12-aminododecanoic acid; and BKM-872; c[DArg-Arg-Eac-Ser-DF5F-Oic-Arg] where Eac is 6-aminocaproic acid. BKM-824 was the only peptide within this series that possessed a discernable solution structure. The NMR data indicate the presence of a type I beta-turn between residues F5F4 and Ava1, a C-terminal-like end. Molecular dynamics calculations show that a type I beta-turn from DF5F4 to Ava1 does exist although the turn was somewhat distorted. This result differs from the structures seen in linear bradykinin antagonists, which usually possess a type II'beta-turn at the C-terminal end and the presence of a defined turn is correlated with bradykinin antagonist activity. There is no solution structure for BKM-870 and BKM-872 but a correlation between the primary sequence Arg(terminal)-DArg1-Arg2-long chain aliphatic amino acid and anti-cancer activity is evident.     

Theoretical conformational analysis of methylamide N-acetyl-L-arginine]

The spatial structure of methylamide N-acetyl-L-argine was studied taking into account the non-valent and electrostati interactions, the torsion energy, and the distorsion of valency angles. Calculation of the favourable conformations of the molecule was carried out with the use of all the combinations of angles phi, psi, chi1 divided by chi4 as an intital approximation. These correspond to the low energy forms of the main chain and to the minima of the torsion potentials of the side chain. Conformational possibilities of arginine and lysine were compared. The calculated stable conformation of N-acetyl-L-arginine-methylamide are compared with the geometry of arginine residues in the proteins with known structure.     

Theoretical conformational analysis of doulble-stranded polynucleotides]

Calculations of intramolecular interaction energy of two-stranded helical homopolynucleotide in the function of nine conformational variables have been carried out by the method of atom-atom potential functions. Four of these variables determine mutual position of base pairs, other four--deoxiribose ring conformation and other one--orientation of this ring with respect to the base. For this purpose an algorythm connecting dependent variables with independent ones has been developed. The investigation of energy function has shown that in the space of conformational parameters there are two valleys, which correspond to A-and B-families of conformations. Experimentaly determined conformations of two-stranded helical polynucleotides are located along the bottoms of these valleys. Along the bottom of each valley the intramolecular interaction energy changes rather little when conformational parameters change within a wide range. The valleys are separated by an energetical barrier.     

A pharmacologic analysis of several adrenomimetic effects of serotonin]

Contractions of the isolated rat vas deferens in response to the addition of serotonin could be explained by the release of catecholamines from the nerve endings. As shown by the use of D-, M- and T-antagonists of serotonin (LSD-25, indocarb, typindole), symatholytic--bretilium, alpha-adrenolytic--droperidol and imipramine, this effect was not associated with the direct activation of serotonin- and adrenoreactive receptors.     

Theoretical conformational analysis of a μ-selective cyclic opioid peptide analog

The allowed conformations of the μ-receptor-selective cyclic opioid peptide analog were determined using a grid search through the entire conformational space. Energy minimization of the 13-membered ring structure lacking the exocyclic Tyr¹ residue and the Phe³ side chain using the molecular mechanics program Maximin resulted in only four low-energy conformations. These four ring structures served as templates for a further energy minimization study with the Tyr¹ residue and Phe³ side chain added to the molecule. The results indicated that the Tyr¹ and Phe³ side chains enjoy considerable orientational freedom, but nevertheless, only a limited number of low-energy side-chain configurations were found. The obtained low-energy conformers are discussed in relation to various proposed models of the bioactive conformation of enkephalins and morphiceptin.     

Theoretical conformational analysis of methylamide of N-acetyl-L-lysine]

The spatial structure of the methylamide of N-acetyl-L-lysine has been analysed taking into account non-bonded and electrostatic interactions, torsional energy, bond angles distortion and hydrogen bonding. Conformational capacities of the backbone and mutual dependence of spatial structures of the backbone and the side chain was described by conformational maps obtained by energy minimisation, the dihedral angles and the bond angles of the side chain being varied for every phi, psi point. Every possible combination for phi, psi, x1-x5-angles was used corresponding to the stable form of the backbone and to torsion potential minima of the initial approximations in the calculation of preferred conformations of the molecule. Comparisons are made between stable forms of the methylamide of N-acetyl-L-lysine and Lys residues in proteins with known structure.     

Crystal structure and conformational analysis of angiotensinogen fragments

The tripeptide acetyl-L-prolyl-L-phenylalanyl-L-histidine crystallizes in the orthorhombic space group P2(1)2(1)2(1) with eight molecules in a unit cell of dimensions a = 9.028(2), b = 140.54(6) and c = 42.41(1)A. The structure has been solved by direct methods and refined to an R value of 0.056 for 2904 observed reflections. The molecule exists as a zwitterion with terminal (His)CO2- and (imidazole)H+ as charged groups. The two peptide molecules in the structure adopt a type I beta-turn with Pro and Phe as the corner residues. The main conformational difference between the two crystallographically independent molecules is seen to be in the histidine side-chain orientations. The molecules arrange themselves in sheets perpendicular to the c axis. All hydrophobic side chains lie on one side of the sheets thus generated, whereas the hydrophilic groups are located on the other side. An interesting feature of the crystal structure is the existence of a water layer between adjacent peptide sheets. The conformational study of the isolated Ac-His-Pro-Phe-His-MA using energy calculations gives a rather limited number of stable conformers. The most stable corresponds to a type I beta-turn stabilized through two hydrogen bonds, followed by a less stable type II beta-turn (delta E = 2.0 kcal) and a partly helical structure (delta E = 2.6 kcal).     

Synthesis, biological activity and conformational analysis of cyclic GRF analogs

A novel cyclic GRF analog, cyclo(Asp8-Lys12)-[Asp8,Ala15]-GRF(1-29)-NH2, i.e. cyclo8,12[Asp8,Ala15]-GRF(1-29)-NH2, was synthesized by the solid phase procedure and found to retain significant biological activity. Solid phase cyclization of Asp8 to Lys12 proceeded rapidly (approximately 2 h) using the BOP reagent. Substitution of Ala2 with D-Ala2 and/or NH2-terminal replacement (desNH2-Tyr1 or N-MeTyr1) in the cyclo8,12[Asp8,Ala15]-GRF(1-29)-NH2 system resulted in highly potent analogs that were also active in vivo. Conformational analysis (circular dichroism and molecular dynamics calculations based on NOE-derived distance constraints) demonstrated that cyclo8,12[Asp8,Ala15]-GRF(1-29)-NH2 contains a long alpha-helical segment even in aqueous solution. A series of cyclo8,12 stereoisomers containing D-Asp8 and/or D-Lys12 were prepared and also found to be highly potent and to retain significant alpha-helical conformation. The high biological activity of cyclo8,12[N-MeTyr1,D-Ala2,Asp8,Ala15]-GRF(1-29)- NH2 may be explained on the basis of retention of a preferred bioactive conformation.     

Nuclear magnetic resonance analysis and conformational characterization of a cyclic decapeptide antagonist of gonadotropin-releasing hormone

Two-dimensional proton nuclear magnetic resonance spectroscopy at 500 MHz has been carried out on the cyclic decapeptide antagonist of gonadotropin-releasing hormone: cyclo-(delta 3-Pro1-D-pClPhe2-D-Trp3-Ser4-Tyr5-D-Trp6-NMeLeu7-Arg8- Pro9-beta-Ala 10). The antagonist exists in two slowly interconverting conformations. All data are consistent with the conclusion that one form has all-trans peptide bonds and the other has a cis beta-Ala10-delta3-Pro1 bond. With the use of sequential assignment methods, chemical shift assignments were obtained for all backbone and side-chain protons of both conformational isomers except for the serine and tyrosine hydroxyl groups and the C gamma, C delta, and guanidinium group protons of the arginine. Temperature dependence of spectral parameters and magnitudes of observed nuclear Overhauser effects support the interpretation that both conformers of the antagonist consist of two beta-turns (type II', D-Trp6-NMeLeu7; type II, delta 3-Pro1-D-pClPhe2) connected by extended antiparallel beta-like strands.