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.     

Molecular mechanics investigation of the flexibility of some cyclic peptides

Conformation space near the crystal conformations of proline-containing cyclic octapeptides and cyclic hexapeptides of C2 sequence symmetry, e.g. cyclo-(Gly-Pro-D-Phe)2 and cyclo-(D-Ala-Gly-Pro-D-Phe)2, was explored using molecular mechanics. Conformations found in crystals were energy minimized, distortions were introduced by systematically fixing backbone dihedral angles at individual residues, and nearby energy-minimized conformations were then located. Interatomic distances and dihedral angles were examined in the conformations within a few kilocalories of the most stable conformation. A common form of flexibility was found to involve libration of amide planes. Among the peptides examined, the cyclic hexapeptides were found to have greater freedom than the cyclic octapeptides, and cyclo-(D-Ala-Gly-Pro-D-Phe)2 was found to be more rigid than cyclo-(D-Ala-Gly-Pro-Phe)2.     

The beta-turn scaffold of tripeptide containing an azaphenylalanine residue

The conformational preferences of azaphenylalanine-containing peptide were investigated using a model compound, Ac-azaPhe-NHMe with ab initio method at the HF/3-21G and HF/6-31G(*) levels, and the seven minimum energy conformations with trans orientation of acetyl group and the 4 minimum energy conformations with cis orientation of acetyl group were found at the HF/6-31G(*) level if their mirror images were not considered. An average backbone dihedral angle of the 11 minimum energy conformations is phi=+/-91 degrees +/-24 degrees , psi =+/-18 degrees +/-10 degrees (or +/-169 degrees +/-8 degrees ), corresponding to the i+2 position of beta-turn (delta(R)) or polyproline II (beta(P)) structure, respectively. The chi(1) angle in the aromatic side chain of azaPhe residue adopts preferentially between +/-60 degrees and +/-130 degrees, which reflect a steric hindrance between the N-terminal carbonyl group or the C-terminal amide group and the aromatic side chain with respect to the configuration of the acetyl group. These conformational preferences of Ac-azaPhe-NHMe predicted theoretically were compared with those of For-Phe-NHMe to characterize the structural role of azaPhe residue. Four tripeptides containing azaPhe residue, Boc-Xaa-azaPhe-Ala-OMe [Xaa=Gly(1), Ala(2), Phe(3), Asn(4)] were designed and synthesized to verify whether the backbone torsion angles of azaPhe reside are still the same as compared with theoretical conformations and how the preceding amino acids of azaPhe residue perturb the beta-turn skeleton in solution. The solution conformations of these tripeptide models containing azaPhe residue were determined in CDCl(3) and DMSO solvents using NMR and molecular modeling techniques. The characteristic NOE patterns, the temperature coefficients of amide protons and small solvent accessibility for the azapeptides 1-4 reveal to adopt the beta-turn structure. The structures of azapeptides containing azaPhe residue from a restrained molecular dynamics simulation indicated that average dihedral angles [(phi(1), psi(1)), (phi(2), psi(2))] of Xaa-azaPhe fragment in azapeptide, Boc-Xaa-azaPhe-Ala-OMe were [(-68 degrees, 135 degrees ), (116 degrees, -1 degrees )], and this implies that the intercalation of an azaPhe residue in tripeptide induces the betaII-turn conformation, and the volume change of a preceding amino acid of azaPhe residue in tripeptides would not perturb seriously the backbone dihedral angle of beta-turn conformation. We believe such information could be critical in designing useful molecules containing azaPhe residue for drug discovery and peptide engineering.     

Free energy landscapes of peptides by enhanced conformational sampling

The free energy landscapes of peptide conformations in water have been observed by the enhanced conformational sampling method, applying the selectively enhanced multicanonical molecular dynamics simulations. The conformations of the peptide dimers, -Gly-Gly-, -Gly-Ala-, -Gly-Ser-, -Ala-Gly-, -Asn-Gly-, -Pro-Gly-, -Pro-Ala-, and -Ala-Ala-, which were all blocked with N-terminal acetyl and C-terminal N-methyl groups, were individually sampled with the explicit TIP3P water molecules. From each simulation trajectory, we obtained the canonical ensemble at 300 K, from which the individual three-dimensional landscape was drawn by the potential of mean force using the three reaction coordinates: the backbone dihedral angle, psi, of the first amino acid, the backbone dihedral angle, phi, of the second amino acid, and the distance between the carbonyl oxygen of the N-terminal acetyl group and the C-terminal amide proton. The most stable state and several meta-stable states correspond to extended conformations and typical beta-turn conformations, and their free energy values were accounted for from the potentials of mean force at the states. In addition, the contributions from the intra-molecular energies of peptides and those from the hydration effects were analyzed. Consequently, the stable beta-turn conformations in the free energy landscape were consistent with the empirically preferred beta-turn types for each amino acid sequence. The thermodynamic values for the hydration effect were decomposed and they correlated well with the empirical values estimated from the solvent accessible surface area of each molecular conformation during the trajectories. The origin of the architecture of protein local fragments was analyzed from the viewpoint of the free energy and its decomposed factors.     

11-cis-retinal protonated Schiff base: influence of the protein environment on the geometry of the rhodopsin chromophore

Density functional theory (DFT) calculations based on the self-consistent-charge tight-binding approximation have been performed to study the influence of the protein pocket on the 3-dimensional structure of the 11-cis-retinal Schiff base (SB) chromophore. Starting with an effectively planar chromophore embedded in a protein pocket consisting of the 27 next-nearest amino acids, the relaxed chromophore geometry resulting from energy optimization and molecular dynamics (MD) simulations has yielded novel insights with respect to the following questions: (i) The conformation of the beta-ionone ring. The protein pocket tolerates both conformations, 6-s-cis and 6-s-trans, with a total energy difference of 0.7 kcal/mol in favor of the former. Of the two possible 6-s-cis conformations, the one with a negative twist angle (optimized value: -35 degrees ) is strongly favored, by 3.6 kcal/mol, relative to the one in which the dihedral is positive. (ii) Out-of-plane twist of the chromophore. The environment induces a nonplanar helical deformation of the chromophore, with the distortions concentrated in the central region of the chromophore, from C10 to C13. The dihedral angle between the planes formed by the bonds from C7 to C10 and from C13 to C15 is 42 degrees. (iii) The absolute configuration of the chromophore. The dihedral angle about the C12-C13 bond is +170 degrees from planar s-cis, which imparts a positive helicity on the chromophore, in agreement with earlier considerations based on theoretical and spectroscopic evidence.     

NMR investigations of proline and its derivatives. 4-Proton magnetic resonance parameters and structure of acetyl-proline amide.

The proton magnetic resonance (PMR) spectrum of acetyl-proline amide in D2O solution has been analysed by computer simulation. The spectra of the cis and the trans isomers have been separated and their PMR parameters (chemical shift and coupling constants) are given. Vicinal coupling constants of the pyrrolidine ring are interpreted by means of a Karplus zone relation. The chemical shift effect of the anisotropy of both peptide planes is considered. It follows that both isomers are puckered with Cgamma in an endo position, but the cis isomer is more rigid than the trans isomer, which moreover undergoes a small interconversion of the Cgamma and Cdelta atoms between two extreme spatial positions. The dihedral angle phi has different values in both isomers. Thus, the dihedral angle between the two peptide planes is smaller in the trans isomer than in the cis isomer.     

Proton nuclear magnetic resonance characterization of the aromatic residues in the variant-3 neurotoxin from Centruroides sculpturatus Ewing

The amino acid sequence for the variant-3 (CsE-v3) toxin from the venom of the scorpion Centruroides sculpturatus Ewing contains eight aromatic residues. By use of 2D NMR spectroscopic methods, the resonances from the individual protons (NH, C alpha H, C beta H',H", and the ring) for each of the individual aromatic residues have been completely assigned. The spatial arrangement of the aromatic ring systems with respect to each other has been qualitatively analyzed by 2D-NOESY techniques. The results show that Trp-47, Tyr-4, and Tyr-42 are in close spatial proximity to each other. The NOESY contacts and the ring current induced shifts in the resonances of the individual protons of Tyr-4 and Trp-47 suggest that the aromatic ring planes of these residues are in an orthogonal arrangement. In addition, the spatial proximity of the rings in the pairs Tyr-4, Tyr-58; Tyr-42, Tyr-40; and Tyr-40, Tyr-38 has also been established. A comparison with the published crystal structure suggests that there is a minor rearrangement of the aromatic rings in the solution phase. No 2D-NOESY contacts involving Phe-44 and Tyr-14 to any other aromatic ring protons have been observed. The pH dependence of the aromatic ring proton chemical shifts has also been studied. These results suggest that the Tyr-58 phenolic group is experiencing a hydrogen-bonding interaction with a positively charged group, while Tyr-4, -14, -38, and -40 are experiencing through-space interactions with proximal negatively charged groups. The Trp-47 indole NH is interacting with the carboxylate groups of two proximal acidic residues. These studies define the microenvironment of the aromatic residues in the variant-3 neurotoxin in aqueous solution.     

Spatial relationship between the copper and carbonyl cofactors in the active site of pig plasma amine oxidase

The ortho-, meta-, and para-trifluoromethylphenylhydrazine inhibitors of porcine plasma amine oxidase were synthesized. Titrations of plasma amine oxidase with these inhibitors demonstrated that 1 mol of trifluoromethylphenylhydrazine completely and irreversibly modifies 1 mol of enzyme by covalently binding to the active carbonyl cofactor. NMR relaxation measurements on the fluorine nuclei were obtained at 188.22 and 74.84 MHz for each inhibitor-enzyme adduct. These measurements were used to calculate the exact distance and orientation between the inhibitor-binding site and the copper cofactor. The copper lies in the plane of the aromatic ring of the inhibitor at distances of 10.9, 14.3, and 15.5 A from the fluorines in the ortho-, meta-, and para- positions of the ring, respectively. Since the inhibitors react with the active carbonyl cofactor, this defines the relationship between the copper and the active carbonyl cofactor in the enzyme, and provides a basis for choosing between mechanisms for the transfer of electrons from the amine substrate to oxygen.     

The Relationship Between Structure and Activity Among Opioid Peptides

Summary Since the discovery and isolation of the endogenous opioid peptides Leu- and Met-enkephalin, structural studies have been focused on deducing the bioactive conformation of the peptide ligands. Theoretically, linear peptides can have many different backbone conformations, yet early, X-ray studies on enkephalin and its analogues showed only two different backbone conformations: extended and single β-bend. More recent reports include a third conformation for Leu-enkephalin and constrained opioid peptides from two ‘new’ classes (i.e. cyclic and ‘allaromatic’ peptides). In this report the relationship between solid-state X-ray structure and opioid peptide activity is examined. The N-terminal amine nitrogen and the two aromatic rings have previously been identified as structural features important to the biological activity of opioid peptides. From X-ray studies we find that the distances between the centroids of the aromatic rings, and between the N-terminal amine nitrogen and the centroid of the phenylalanine ring, vary over a large range. There is a discernible relationship, however, between the separation of the two rings and their orientation that correlates with activity.     

Hydrogen bonding monitored by deuterium isotope effects on carbonyl 13C chemical shift in BPTI: intra-residue hydrogen bonds in antiparallel beta-sheet.

Deuterium isotope effects on carbonyl 13C magnetic shielding were measured for the backbone carbonyl groups in BPTI (basic pancreatic trypsin inhibitor), and interpreted as a measure of hydrogen bond energies. The effects originate from peptide amide proton deuterium substitution and were observed on carbonyl carbons separated by two or three covalent bonds from the amide H/D. Two-bond isotope effects depend on the energy of the hydrogen bond donated by NH/D. Calibration of the effect with model compound data leads to hydrogen bond enthalpies less than 4.7 kcal/mol. Isotope effects over three bonds from the amide H/D to the carbonyl carbon of the same amino acid residue are observed for seven carbonyl resonances in BPTI. The three-bond isotope effects are highly related to the various backbone conformations. The largest effects are observed for residues with an approximate syn- periplanar conformation of the H-N-C alpha-C = O atoms, as realized for many residues in the BPTI antiparallel beta-sheet. The residues that show measurable three-bond effects have unusually short distances between H and O. The size of this effect decreases rapidly with increased O..H distance in the open five-membered ring. This observation suggests appreciable interactions in these rings.     

Conformation and absolute configuration of 11,12-cyclopropyl retinal analogs – an RHF/DFT/CIS study

Ab initio RHF and DFT/B3LYP calculations at the 6-31G** level have been performed to study possible conformations of the cyclopropyl retinal Schiff base analog 3 of known absolute configuration. In both the free base and the protonated form, the geometries are determined on the diene side by optimum conjugative interaction with the three-membered ring, on the triene side by repulsive interaction with the 9-methyl group. There are three low energy conformations, in which the seven-membered ring is either in a chair or in a twist-chair conformation. To decide between these alternatives, chiroptical parameters were calculated employing the GAUSSIAN/CIS routines and compared with the CD spectrum obtained by Nakanishi et al. Of the energy-minimized geometries only two fit the experimental data. In both, the dihedral angle C12-C13, which is indicative of the relative orientation of the two chromophores, is positive.     

Computational analysis of the first biheterocyclization site of the antibiotic microcin B17

Microcin B 17 (MccB17) undergoes an enzyme catalyzed posttranslational modification to form four oxazole and four thiazole rings. Four of these rings form 4,2 - connected biheterocyclic functionalities. In this study, the hexapeptide sequence surrounding the first biheterocyclization site of microcin B17 was examined using computational calculations and database analysis to see if it was preorganized for cyclization in a manner similar to that found in the autocatalytic posttranslational cyclization of Green Fluorescent Protein (GFP). Attention was focused on the intermolecular distances between the sulfur and oxygen atoms of the cysteine and serine residues and the carbonyl carbons which they attack in the ring formation. Conformational searches located some low energy conformations that contained relatively short oxygen to carbonyl carbon distances, which indicated that the oxazole forming fragment in microcin B17 is preorganized for cyclization. However, the lack of any clear patterns for the sulfur to carbon distances show that the side-chain of cysteine does not adopt any low energy conformations that are geometrically preorganized for cyclization. The MccB17 synthetase enzyme complex which catalyzes the cyclization process therefore has both steric and electronic functions. The data obtained in this investigation is in agreement with empirical data which shows that biheterocyclization will only occur if the thiazole forms before the oxazole.     

Ab initio quantum mechanical calculations of the variation of the 1H and 13C nuclear magnetic shielding constants in proline as a function of the angle psi

The variation of the nuclear magnetic shielding constant of the different protons and carbons of trans HCO-L-Pro-NH2 with the value of the angle psi is calculated by a non-empirical method for three conformations of the proline ring. The results concerning the CH protons show that the chemical shift of the alpha, beta and gamma endo hydrogens can vary by more than 1 ppm when psi goes from -30 degrees to 180 degrees. The theoretical variation of the chemical shift difference between alpha and gamma or beta and gamma carbons is found to be sensitive to the puckering of the proline ring. For the second of these differences the theoretical results are in agreement with Siemion's relation only for a limited range of molecular conformations. Additional calculations show that the variations of the proton shifts with the value of psi are due to the magnetic anisotropy of the proline carbonyl group and to the polarization of the CH bonds by the multipolar charge distribution carried by this carbonyl. The results are discussed in relation to experiment and the possibility of using 1H and 13C chemical shifts for the determination of the value of the torsion angle about the C alpha C' bond.     

New Insights into the Interdependence between Amino Acid Stereochemistry and Protein Structure

To successfully design new proteins and understand the effects of mutations in natural proteins, we must understand the geometric and physicochemical principles underlying protein structure. The side chains of amino acids in peptides and proteins adopt specific dihedral angle combinations; however, we still do not have a fundamental quantitative understanding of why some side-chain dihedral angle combinations are highly populated and others are not. Here we employ a hard-sphere plus stereochemical constraint model of dipeptide mimetics to enumerate the side-chain dihedral angles of leucine (Leu) and isoleucine (Ile), and identify those conformations that are sterically allowed versus those that are not as a function of the backbone dihedral angles ? and ψ. We compare our results with the observed distributions of side-chain dihedral angles in proteins of known structure. With the hard-sphere plus stereochemical constraint model, we obtain agreement between the model predictions and the observed side-chain dihedral angle distributions for Leu and Ile. These results quantify the extent to which local, geometrical constraints determine protein side-chain conformations.     

The structure of 2-methoxy-5-(2',3',4'-trimethoxyphenyl)tropone, an effective analog of colchicine

-Methoxy-5-(2',3',4'-trimethoxyphenyl) tropone is an active analog of colchicine, a mitotic spindle inhibitor, which is missing the middle "B" ring. This compound crystallizes in the triclinic system, space group P1, with Z = 2; a = 10.135(2), b = 10.166 (4), and c = 7.863(2) A; alpha = 82.15(3), beta = 103.49(3), and gamma = 107.16(2); degrees and V = 750.7(4) A. The structure was solved by direct methods and refined by full-matrix least-squares to a final R = 0.063, using 2503 observed reflections and 271 parameters. Despite the absence of the middle ring, the conformation of the molecule is similar to that of colchicine, isocolchicine , and their derivatives. The troponoid ring is dissimilar to the phenyl ring in that it is not aromatic and does have alternating short and long bond lengths. The dihedral angle between the least-squares planes of the two rings is -57.4 degrees. Van der Waals surface representations of the analog and colchicine are presented to demonstrate the similarity and differences of these two molecules . The structural information of the analog is consistent with the interpretation of thermodynamic parameters which govern the interactions between brain tubulin and the analog.     

The interaction with DNA of unfused aromatic systems containing terminal piperazino substituents. Intercalation and groove-binding

A number of unfused tricyclic aromatic intercalators have shown excellent activity as amplifiers of the anticancer activity of the bleomycins and the 4',6-diphenylpyrimidines, 2a and 2b, with terminal basic functions (4-methylpiperazino groups) have been synthesized to test the structural requirements for amplifier-DNA interactions. The terminal piperazine rings are bulky, have limited flexibility, and are twisted out of the phenyl ring plane in both 2a and 2b. With 2a the pyrimidine is unsubstituted at position 5 and the conformation predicted by molecular mechanics calculations has a 25-30 degrees twist between the phenyl and pyrimidine ring planes. With 2b the 5-position is substituted with a methyl group and this causes a larger twist angle (50-60 degrees) between the phenyl and pyrimidine planes. These conformational variations lead to markedly different DNA interactions for 2a and 2b. Absorption, CD and NMR spectral, viscometric, flow dichroism and kinetics results indicate that 2a binds strongly to DNA by intercalation while 2b binds more weakly in a groove complex. The general structure and conformation of 2a, a slightly twisted, unfused-aromatic system with terminal piperazino groups is more similar to groove-binding agents such as Hoechst 33258 than to intercalators. The fact that 2a forms a strong intercalation complex with DNA is unusual but in agreement with studies on other amplifiers of anticancer drug action. Molecular modeling studies provide a second unusual feature of the 2a intercalation complex. While most well-characterized intercalators bind with their bulky and/or cationic substitutents in the DNA minor groove, the cationic piperazino groups of 2a are too large to bind in the minor groove in an intercalation complex but can form strong interactions with DNA in the major groove. The tricyclic aromatic ring system of 2a stacks well with adjacent base-pairs in the major-groove complex and the piperazino groups have good electrostatic and van der Waals interactions with the DNA backbone.     

The normal modes of the gramicidin-A dimer channel.

The dynamics of the gramicidin-A dimer channel is studied in the harmonic approximation by a vibrational analysis of the atomic motions relative to their equilibrium positions. The system is represented by an empirical potential energy function, and all degrees of freedom (bonds lengths, bond angles, and torsional angles) are allowed to vary. The thermal fluctuations in the backbone dihedral angles phi and psi, atomic root mean square displacements, and the correlations between the different amide planes are computed. It is found that only adjacent dihedral psi i and phi i+1 are strongly correlated, while different hydrogen-bonded amide planes are only weakly correlated. Modes with relatively low vibrational frequencies (75-175 cm-1) make the dominant contributions to the carbonyl librations. The general flexibility of the structure and the role of carbonyl librations in the ion transport mechanism are discussed.     

New Insights into the Interdependence between Amino Acid Stereochemistry and Protein Structure

To successfully design new proteins and understand the effects of mutations in natural proteins, we must understand the geometric and physicochemical principles underlying protein structure. The side chains of amino acids in peptides and proteins adopt specific dihedral angle combinations; however, we still do not have a fundamental quantitative understanding of why some side-chain dihedral angle combinations are highly populated and others are not. Here we employ a hard-sphere plus stereochemical constraint model of dipeptide mimetics to enumerate the side-chain dihedral angles of leucine (Leu) and isoleucine (Ile), and identify those conformations that are sterically allowed versus those that are not as a function of the backbone dihedral angles ϕ and ψ. We compare our results with the observed distributions of side-chain dihedral angles in proteins of known structure. With the hard-sphere plus stereochemical constraint model, we obtain agreement between the model predictions and the observed side-chain dihedral angle distributions for Leu and Ile. These results quantify the extent to which local, geometrical constraints determine protein side-chain conformations.     

Dihedral angles of tripeptides in solution directly determined by polarized Raman and FTIR spectroscopy

The amide I mode of the peptide linkage is highly delocalized in peptides and protein segments due to through-bond and through-space vibrationally coupling between adjacent peptide groups. J. Phys. Chem. B. 104:11316-11320) used coherent femtosecond infrared (IR) spectroscopy to determine the excitonic coupling energy and the orientational angle between the transition dipole moments of the interacting amide I modes of cationic tri-alanine in D(2)O. Recently, the same parameters were determined for all protonation states of tri-alanine by analyzing the amide I bands in the respective IR and isotropic Raman spectra (. J. Am. Chem. Soc. 119:1720-1726.). In both studies, the dihedral angles phi and psi were then obtained by utilizing the orientational dependence of the coupling energy obtained from ab initio calculations on tri-glycine in vacuo (. J. Raman Spectrosc. 29:81-86) to obtain an extended 3(1) helix-like structure for the tripeptide. In the present paper, a novel algorithm for the analysis of excitonic coupling between amide I modes is presented, which is based on the approach by Schweitzer-Stenner et al. but avoids the problematic use of results from ab initio calculations. Instead, the dihedral angles are directly determined from infrared and visible polarized Raman spectra. First, the interaction energy and the corresponding degree of wave-function mixing were obtained from the amide I profile in the isotropic Raman spectrum. Second, the depolarization ratios and the amide I profiles in the anisotropic Raman and IR-absorption spectra were used to determine the orientational angle between the peptide planes and the transition dipole moments, respectively. Finally, these two geometric parameters were utilized to determine the dihedral angles phi and psi between the interacting peptide groups. Stable extended conformations with dihedral angles in the beta-sheet region were obtained for all protonation states of tri-alanine, namely phi(+) = -126 degrees, psi(+) = 178 degrees; phi(+/-) = -110 degrees, psi(+/-) = 155 degrees; and phi(-) = -127 degrees, psi(-) = 165 degrees for the cationic, zwitterionic, and anionic state, respectively. These values reflect an extended beta-helix structure. Tri-glycine was found to be much more heterogeneous in that different extended conformers coexist in the cationic and zwitterionic state, which yield a noncoincidence between isotropic and anisotropic Raman scattering. Our study introduces vibrational spectroscopy as a suitable tool for the structure analysis of peptides in solution and tripeptides as suitable model systems for investigating the role of local interactions in determining the propensity of peptide segments for distinct secondary structure motifs.     

Localized solution structure refinement of an F45W variant of ubiquitin using stochastic boundary molecular dynamics and NMR distance restraints.

The local structure within an 8-A radius around residue 45 of a recombinant F45W variant of human ubiquitin has been determined using 67 interproton distance restraints measured by two-dimensional proton NMR. Proton chemical shift evidence indicates that structural perturbations due to the F45W mutation are minimal and limited to the immediate vicinity of the site of mutation. Simulated annealing implemented with stochastic boundary molecular dynamics was applied to refine the structure of Trp 45 and 10 neighboring residues. The stochastic boundary method allowed the entire protein to be reassembled from the refined coordinates and the outlying unrefined coordinates with little distortion at the boundary. Refinement began with four low-energy indole ring orientations of F45W-substituted wild-type (WT) ubiquitin crystal coordinates. Distance restraints were derived from mostly long-range NOE cross peaks with 51 restraints involving the Trp 45 indole ring. Tandem refinements of 64 structures were done using either (1) upper and lower bounds derived from qualitative inspection of NOE crosspeak intensities or (2) quantitative analysis of cross-peak heights using the program MARDIGRAS. Though similar to those based on qualitative restraint, structures obtained using quantitative NOE analysis were superior in terms of precision and accuracy as measured by back-calculated sixth-root R factors. The six-membered portion of the indole ring is nearly coincident with the phenyl ring of the WT and the indole NH is exposed to solvent. Accommodation of the larger ring is accompanied by small perturbations in the backbone and a 120 degrees rotation of the chi 2 dihedral angle of Leu 50.