Conformational analysis of the endogenous mu-opioid agonist endomorphin-1 using NMR spectroscopy and molecular modeling.

Endomorphin-1 (Tyr-Pro-Trp-Phe-NH2) is a highly selective and potent agonist of the mu-opioid receptor. To identify structural attributes unique to this opioid peptide and potential sites of recognition, a conformational analysis has been performed using multidimensional NMR and molecular modeling techniques. The spectroscopic results, derived from experiments in both DMSO and water, indicate that endomorphin-1 exists in the cis- and trans-configuration with respect to the Pro-omega bond in approximately 25% and 75% populations, respectively. In DMSO, the cis-configuration adopts a compact sandwich conformation in which the Tyr and Trp aromatic rings pack against the proline ring, whereas the trans-configuration adopts an extended conformation. Although non-random structure was not observed in water, condensed phase molecular dynamics calculations indicate that trans-isomers dominate the population in this higher dielectric medium. Structural comparison of the cis- and trans-configurations with morphine and selective mu-peptide ligands PL-017 and D-TIPP, as well as the delta-selective peptide ligands TIPP (delta-antagonist, mu-agonist) and DPDPE were also performed and suggest the trans-isomer is likely the bioactive form. A hypothesis is proposed to explain mu- and delta-selectivity based on the presence of spatially distinct selectivity pockets among these ligands.

Matrix metalloproteinase inhibitors containing a (carboxyalkyl)amino zinc ligand: modification of the P1 and P2' residues.

Systematic modification of the presumed P1 side chain in a series of (carboxyalkyl)amino-based inhibitors of matrix metalloproteinases enabled identification of the 2-(1,3-dihydro-1,3-dioxo-2H-benz[f]isoindol-2-yl)ethyl group as a preferred substituent imparting potent inhibition of the enzymes collagenase and gelatinase. It was subsequently found that the P2'-P3' residues in this series could be replaced by small non-peptide residues, while maintaining inhibitory potency. The imide group in this series of compounds can undergo autocatalytic hydrolysis under neutral conditions.

Isotope effects and alternative substrate reactivities for tryptophan 2,3-dioxygenase.

Tryptophan 2,3-dioxygenase (EC 1.13.1.12) is a hemoprotein which catalyzes the first step in the oxidative degradation of tryptophan. The reaction is believed to proceed by addition of O2 across the 2,3-bond of the indole ring, followed by decomposition of the resultant dioxetane to give N-formylkynurenine. A primary D2O isotope effect of 4.4 on Vmax/Km was observed at the pH optimum, pH 7.0. This implies that abstraction of the indole proton is at least partially rate-determining. An inverse secondary isotope effect of 0.96 was observed for L-[2-3H]tryptophan at this pH. The secondary isotope effect signals the formation of the C-O bond at C-2. As the rate of proton abstraction increased with increasing pH, the D2O isotope effect decreased to 1.2 at pH 8.5 and the secondary isotope effect increased to 0.92. The rate-determining steps therefore change with increasing pH, and bond formation at C-2 becomes more rate-limiting. The secondary isotope effect did not change significantly with varying O2 concentration so that substrate binding is primarily ordered with O2 binding first. The specificity of the enzyme towards substituted tryptophans shows that substitution of the phenyl ring of the indole is sterically unfavorable. Steric hindrance is highest at the 4- and 7-positions, while the 5- and 6-positions are less sensitive. 6-Fluoro-L-tryptophan was more reactive than tryptophan, and the increased reactivity can be explained by an electronic effect that enhances of the rate of C-O bond formation at C-2.

Solution conformation of human big endothelin-1.

The solution conformation of human big endothelin-1, a 38-residue peptide which serves as the putative precursor to the potent vasoconstrictor endothelin-1 has been examined by 1H NMR. NOEs were utilized as distance restraints in the distance geometry program DSPACE to generate initial structures. Further refinement of these structures was accomplished through molecular mechanics/molecular dynamics in an iterative process involving the incorporation of stereospecific assignments of prochiral centers and the use of back-calculation of NOESY spectra. A family of structures consisting of a type II beta-turn for residues 5-8 and an alpha-helix extending from residues 9-16 constitute a well-defined region, as reflected by the atomic root-mean-square (RMS) difference of 1.56 A about the mean coordinate positions of the backbone atoms (N, C, C alpha and O). This core region (residues 1-15) is very similar to the core residues of endothelin-1 (Donlan, M. et al. (1991) J. Cell. Biochemistry, S15G, 85). While the evidence from NOESY and coupling constant data suggests that the C-terminal region, residues 17-34, is not a mixture of randomly distributed chain conformations, it is also not consistent with a single chain conformation. Under the conditions studied, residues 17-38 in human big endothelin-1 in water at pH 3.0 between 20-30 degrees C appear to be represented by a series of conformers in dynamic equilibrium.

NMR restraint analysis of transforming growth factor alpha: a key component for NMR structure refinement.

Structures of the protein, transforming growth factor alpha (TGF-alpha), have been derived from NMR data using distance geometry and subsequent energy refinement. Analysis of the sequential NOE distance bounds using a template algorithm provides a check for consistency in the calculation of bounds, stereospecific assignment of prochiral centers, and secondary structure assignment. Application of the template algorithm to the long range NOEs found within the NMR data sets collected at pH 6.3 and pH 3.4 is used to assess the confidence levels for the accuracy of the structures obtained from modeling. The method also provides critical insight in differentiating regions of the structure that are well defined from those that are not. Use of the restraint analysis protocol is shown to be a powerful adjunct to currently used methods for the assignment of protein structures from NMR data.

Probing structural factors stabilizing antisense oligonucleotide duplexes: NMR studies of a DNA.DNA duplex containing a formacetal linkage.

The duplex formed by annealing the formacetal backbone modified dodecamer d-(CGCGTTOCH2OTTGCGC) to its complementary strand, d(GCGCAAAACGCG) (duplex I), has been studied by NMR techniques and analyzed with reference to its unmodified counterpart (duplex II). Comparison of parameters such as 2D cross-peak intensities, coupling constants, and spectral patterns indicates that structural perturbations caused by the incorporation of the formacetal linkage are minimal and localized to the central T4.A4 block. Duplex I adopts a B-type helical conformation with regular Watson-Crick base pairing and normal minor groove width. The methylene group is accommodated along the phosphate backbone in a conformation similar to that of the PO2 group found in the B-form DNA family. The central T6-T7 base pairs of duplex I melt simultaneously with the duplex, indicating a cooperative transition to single strands. Although the formacetal linkage affects global melting, as evidenced by a 3 degree C reduction in Tm for duplex I with respect to duplex II, the present study indicates that this is not the result of localized premelting at the formacetal site of duplex I but rather reflects the subtle interplay of several structural and energy factors which need to be further explored.

Do interhelical side chain-backbone hydrogen bonds participate in formation of leucine zipper coiled coils?

The leucine zipper proteins are a group of transcriptional regulators that dimerize to form a DNA binding domain. It has been proposed that this dimerization results from the hydrophobic association of the alpha-helices of two leucine zipper monomers into a coiled coil. We propose a model for a coiled coil based on a periodic hydrophobic-hydrophilic amino acid motif found in the leucine zipper regions of 11 transcriptional regulatory proteins. This model predicts the symmetrical formation of secondary hydrogen bonds between the polar side chains of one helix and the peptide carbonyls of the opposite chain, supplementing the interactions between hydrophobic side chains. Physical modeling (CPK) and in vacuo molecular mechanics calculations of the stability of the GCN4 leucine zipper coiled coil configured in accordance with this model demonstrate a greater stability for this conformer than for a conformer configured according to a current hydrophobic model. Molecular dynamics simulations show similar stability of the two models in vacuo but a higher stability of the hydrophobic model in water.

Solution structures of human transforming growth factor alpha derived from 1H NMR data.

The 600-MHz 1H NMR spectrum of the des-Val-Val mutant of human transforming growth factor alpha (TGF-alpha) was reassigned at pH = 6.3. The conformation space of des-Val-Val TGF-alpha was explored by distance geometry embedding followed by restrained molecular dynamics refinement using NOE distance constraints and some torsion angle constraints derived from J-couplings. Over 80 long-range NOE constraints were found by completely assigning all resolved cross-peaks in the NOESY spectra. Low NOE constraint violations were observed in structures obtained with the following three different refinement procedures: interactive annealing in DSPACE, AMBER 3.0 restrained molecular dynamics, and dynamic simulated annealing in XPLOR. The segment from Phe15 to Asp47 was found to be conformationally well-defined. Back-calculations of NOESY spectra were used to evaluate the quality of the structures. Our calculated structures resemble the ribbon diagram presentations that were recently reported by other groups. Several side-chain conformations appear to be well-defined as does the relative orientation of the C loop to the N-terminal half of the protein.

Molecular dynamics simulations of "loop closing" in the enzyme triose phosphate isomerase.

We present molecular dynamics simulations on the active site region of dimeric triose phosphate isomerase (TIM) using the co-ordinates of native chicken muscle TIM as a starting point and performing simulations with no substrate, with dihydroxyacetone phosphate (DHAP), the natural substrate, and with dihydroxyacetone sulfate (DHAS), a substrate analog. Whereas most of the protein moves less than 1 A during the simulation, some residues in the active site loop move more than 8 A during the 10.5 picoseconds of dynamics for each of the simulations. Most interestingly, the nature of the loop motion depends on the substrate, with the largest motion found in the presence of DHAP, and only in the presence of DHAP does the loop move to "close off" the active site pocket. The final structure found for the DHAP-chicken TIM complex is qualitatively similar to that described by Alber et al. for DHAP-yeast TIM. Simulations on the monomeric protein gives insight into why the molecule is active only as a dimer.

Calculation of the relative change in binding free energy of a protein-inhibitor complex.

By means of a thermodynamic perturbation method implemented with molecular dynamics, the relative free energy of binding was calculated for the enzyme thermolysin complexed with a pair of phosphonamidate and phosphonate ester inhibitors. The calculated difference in free energy of binding was 4.21 +/- 0.54 kilocalories per mole. This compares well with the experimental value of 4.1 kilocalories per mole. The method is general and can be used to determine a change or "mutation" in any system that can be suitably represented. It is likely to prove useful for protein and drug design.

Theory and modeling of stereoselective organic reactions.

Theoretical investigations of the transition structures of additions and cycloadditions reveal details about the geometries of bond-forming processes that are not directly accessible by experiment. The conformational analysis of transition states has been developed from theoretical generalizations about the preferred angle of attack by reagents on multiple bonds and predictions of conformations with respect to partially formed bonds. Qualitative rules for the prediction of the stereochemistries of organic reactions have been devised, and semi-empirical computational models have also been developed to predict the stereoselectivities of reactions of large organic molecules, such as nucleophilic additions to carbonyls, electrophilic hydroborations and cycloadditions, and intramolecular radical additions and cycloadditions.

Stabilized venous distensibility of normotensive and hypertensive humans on high and low sodium intake.

Venous responses to stabilized orthostasis (45 degrees head-up tilt) were studied in seven normotensive subjects and eight hypertensive patients, when on high and low dietary sodium intake. Exchangeable sodium and blood volumes were determined to permit correlation with any significant changes in venous behavior. The intent of this study was to detect and analyze any diet-induced changes in responses of forearm veins to prolonged orthostasis. The pharmacological effects of sodium depletion by medication and diet on arteries and veins of hypertensives are discussed. The results of this study indicate that dietary sodium depletion did not have adverse effects on the ability to maintain stabilized venous tone during orthostasis. These results support recommendations that moderate dietary sodium restriction be included as part of antihypertensive regimens.