Identification of disulfide bridges in a cardiotoxic peptide by electrospray ionization

A new method was developed for subjecting a peptide to a specified number of Edman degradation cycles on an automated polypeptide sequencer and desorbing the residual peptide for further investigations. The procedure was applied in combination with electrospray ionization mass spectrometry to identify the four disulfide bridges present in a small tightly bound peptide. The task was accomplished using only a few nanomoles of the intact peptide.

[D-TRP32]neuropeptide Y: a competitive antagonist of NPY in rat hypothalamus

Neuropeptide Y (NPY) is a potent orexigenic peptide. Structure-activity studies have revealed that nearly the entire sequence of NPY is required to elicit feeding responses. Therefore, in order to develop antagonistic peptides for NPY-induced feeding, we synthesized full-length analogs of NPY, substituting D-Trp in the C-terminal receptor binding region, and screened their activity in rat hypothalamus. Although [D-Trp36]NPY and [D-Trp34]NPY inhibited isoproterenol-stimulated hypothalamic membrane adenylate cyclase activity, [D-Trp32]NPY exhibited no intrinsic activity. Furthermore, [D-Trp32]NPY inhibited [125I]NPY binding to rat hypothalamic membranes with a potency comparable to that of NPY. The presence of 30 and 300 nM concentrations of [D-Trp32]NPY shifted the inhibitory dose-response curve of NPY on isoproterenol-stimulated hypothalamic membrane adenylate cyclase activity parallel to the right with comparable KB values. Moreover, in vivo experiments in rats revealed that [D-Trp32]NPY (10 micrograms) significantly attenuated the 1-h feeding response induced by NPY (1 microgram). Several other substitutions at position 32 including 2-D-Nal resulted in agonist activity, suggesting that there are strict structural requirements to induce the antagonistic property in NPY. These findings confirm that [D-Trp32]NPY is a competitive antagonist of NPY in both in vitro and in vivo systems. Analogs based on [D-Trp32]NPY may have potential clinical application, since NPY has been implicated in the pathophysiology of a number of feeding disorders including obesity, anorexia, and bulimia.

Molecular dynamics of sickle and normal hemoglobins

Molecular dynamics (MD) simulations have been carried out for 62.5 ps on crystal structures of deoxy sickle cell hemoglobin (HbS) and normal deoxy hemoglobin (HbA) using the CHARMM MD algorithm, with a time step of 0.001 ps. In the trajectory analysis of the 12.5-62.5 (50 ps) simulation, oscillations of the radius of gyration and solvent-accessible surface area were calculated. HbS exhibited a general contraction during the simulation, while HbA exhibited a nearly constant size. The average deviations of simulated structures from the starting structures were found to be 1.8 A for HbA and 2.3 A for HbS. The average rms amplitudes of atomic motions (atomic flexibility) were about 0.7 A HbA and about 1.0 A for HbS. The amplitudes of backbone motion correlate well with temperature factors derived from x-ray crystallography. A comparison of flexibility between the alpha- and beta-chains in both HbA and HbS indicates that the beta-chains generally exhibited greater flexibility than the alpha-chains, and that the HbS beta-chains exhibit greater flexibility in the N-terminal and D- and F-helix regions than do those of HbA. The average amplitude of backbone torsional oscillations was about 9 degrees, a value comparable with that of other simulations, with enhanced torsional oscillation occurring primarily at the ends of helices or in loop regions between helices. Comparison of atomic flexibility and torsional oscillation results suggests that the increased beta-chain flexibility results from relatively concerted motions of secondary structure elements. The increased flexibility may play an important role in HbS polymerization.(ABSTRACT TRUNCATED AT 250 WORDS)

Location of potential binding sites on deoxy hemoglobin for the design of antigelling agents

The binding sites of indole-based gelation inhibitors with sickle cell hemoglobin were investigated by two parallel theoretical approaches. A geometric approach originated by Kuntz and co-workers uses a spatial buildup scheme to locate potential binding regions, while a hybrid grid/geometric search method searches for specific indole ring binding pockets over the hemoglobin surface. The binding sites derived from these calculations were tested for their ability to accommodate indole rings by means of accessibility calculations with probes of various radii. These sites were further scanned for van der Waals' overlap and electrostatic interactions. A full 5BrTrp residue was built in each indole ring binding site, and its conformational energy of association with sickle hemoglobin was calculated at that site. Our theoretical results predict a total of 14 potential binding regions, including all of the sites observed from X-ray crystallography, and sites that are consistent with solution nuclear magnetic resonance studies.

Conformational space of alpha 1-4 glycosidic linkage: a molecular dynamics study

Molecular Dynamics simulations have been carried out for 100 ps on crystal structure of beta-cyclodextrin in vacuo and with explicit inclusion of solvent at constant pressure and constant temperature using the GROMOS MD algorithm, with a time step of 0.005 ps. The conformational space of the glycosidic linkage was studied by calculating two virtual dihedrals connecting the successive glucose units for the 2000 structures saved during the two simulations. Three preferred regions for alpha 1-4 glycosidic linkage were found in both the simulations. The use of these virtual dihedral angles in representing the glycosidic linkage is also brought out from these studies.

Nuclear magnetic resonance and molecular modeling studies on O-beta-D-galactopyranosyl-(1----4)-O-beta-D-xylopyranosyl-(1----0)-L-se rine, a carbohydrate-protein linkage region fragment from connective tissue proteoglycans

The solution conformation of O-beta-D-galactopyranosyl-(1----4)-O-beta-D-xylopyranosyl-(1----0)-L-ser ine (GXS), a carbohydrate-protein linkage region fragment from connective tissue proteoglycans, was investigated by two-dimensional NMR spectroscopy and molecular modeling calculations. Specifically, the 1H and 13C resonances were assigned by 2D-COSY and by 1H-13C heteronuclear correlation spectroscopy methods. 2D-NOESY was used to generate distance constraints between the galactose and xylose and between the xylose and serine residues. The 1H vicinal coupling constants for the sugars and the serine were also determined. A general molecular modeling methodology suitable for complex carbohydrates was developed. This methodology employed molecular dynamics and energy minimization procedures together with the application of inter-residue spatial constraints across the linkages derived from 2D-NOESY. The first step in this methodology is the generation of a wide variety of starting conformations that span the (phi, psi) space for each linkage. In the present study, nine such conformations were constructed for each linkage using the torsion angles phi and psi corresponding to the gauche+, gauche-, and trans configurations across each of the two bonds constituting the linkage. These conformations were subjected to a combined molecular dynamics/energy minimization refinement using the NOESY derived constraints as pseudoenergy functions. Families of conformations for the whole molecule were then constructed from the structures derived for each linkage. Characterization of GXS using this methodology identified a single family of conformations that are consistent with the solution phase NMR data on this molecule.

The distribution of physical, chemical and conformational properties in signal and nascent peptides

Signal peptides play a major role in an as-yet-undefined way in the translocation of proteins across membranes. The sequential arrangement of the chemical, physical and conformational properties of the signal and nascent amino acid sequences of the translocated proteins has been compiled and analysed in the present study. The sequence data of 126 signal peptides of length between 18 and 21 residues form the basis of this study. The statistical distribution of the following properties was studied hydrophobicity, Mr, bulkiness, chromatographic index and preference for adopting alpha-helical, beta-sheet and turn structures. The contribution of each property to the sequence arrangement was derived. A hydrophobic core sequence was found in all signal peptides investigated. The structural arrangement of the cleavage site was also clearly revealed by this study. Most of the physical properties of the individual sequences correlated (correlation coefficient approximately 0.4) very well with the average distribution. The preferred occupancy of amino acid residues in the signal and nascent sequences was also calculated and correlated with their property distribution. The periodic behaviour of the signal and nascent chains was revealed by calculating their hydrophobic moments for various repetitive conformations. A graphical analysis of average hydrophobic moments versus average hydrophobicity of peptides revealed the transmembrane characteristics of signal peptides and globular characteristics of the nascent peptides.

Secondary structure prediction for the spectrin 106-amino acid segment, and a proposed model for tertiary structure

A collective secondary structure prediction for the human erythrocyte spectrin 106-residue repeat segment is developed, based on the sequences of nine segments that have been reported in the literature, utilizing a consensus of several secondary structure prediction methods for locating turn regions. The analysis predicts a five-fold structure, with three alpha-helices and two beta-strand regions, and differs from previous models on the lengths of the helices and the existence of beta-strand structure. We also demonstrate that this structural motif can be folded into tertiary structures that satisfy the experimental spectrin data and several general principles of protein organization.

Tooth enamel protein, amelogenin, has a probable beta-spiral internal channel, Gln112-Leu138, within a single polypeptide chain: preliminary molecular mechanics and dynamics studies

Molecular dynamics simulation, with backbone constraints for 20 ps of equilibration and simulation, of a repeating polypeptide segment, Gln-Pro-His-Gln-Pro-Leu-Gln-Pro-His-Gln-Pro-Leu-Gln-Pro-Met-(Gln-Pro-Leu )4, constituting residues 112-138 of bovine amelolgenin, a 19.35 kD hydrophobic protein, are discussed. It is generally believed that the above polypeptide segment is important for the interaction of amelogenin with Ca++ ions, which occurs in the early phases of enamel mineralization. An energetically stable structure of the above polypeptide with recurrent beta-turns is observed and contains a pore of approximately 1 A radius along the helical that can accommodate an unhydrated Ca++ ion. The length of the polypeptide possesses correct dimensions to span a bilayer. The proposed structure is unique among known polypeptide and protein structures.

Identification of peptide hormones of the amphipathic helix class using the helical hydrophobic moment algorithm

Eisenberg's helical hydrophobic moment (less than mu H greater than) algorithm was applied to the analysis of the primary structure of amphipathic alpha-helical peptide hormones and an optimal method for identifying other peptides of this class determined. We quantitate and compare known amphipathic helical peptide hormones with a second group of peptides with proven nonamphipathic properties and determine the best method of distinguishing between them. The respective means of the maximum 11 residue less than mu H greater than for the amphipathic helical and control peptides were 0.46 (+/-/-0.07) and 0.33 (0.07) (P + 0.004). To better reflect the amphipathic potential of the entire peptide, the percent of 11 residue segments in each peptide above a particular less than mu H greater than was plotted vs less than mu H greater than. The resulting curves are referred to as HM-C. The mean HM-C (of the two groups) was highly significantly different such that the HM-C method was superior to others in its ability to distinguish amphipathic from nonamphipathic peptides. Several potential new members of this structural class were identified using this approach. Molecular modeling of a portion of one of these, prolactin inhibitory factor, reveals a strongly amphipathic alpha helix at residues 4-21. This computer-based method may enable rapid identification of peptides of the amphipathic alpha-helix class.

AUGUR: a program to predict, display and analyze the tertiary structure of B-DNA

AUGUR is a program to predict, display and analyze the three-dimensional structure of B-DNA. The user can choose one of six models to predict the helical parameters of a given sequence. These parameters are then used to generate the coordinates of the DNA model in three-dimensional space (trajectory). The trajectory can be displayed and rotated on a graphics terminal. The trajectory and helical parameters can also be searched for bends and structural homologues.

Phenylalanine transfer RNA: molecular dynamics simulation

Yeast phenylalanine transfer RNA was subjected to a 12-picosecond molecular dynamics simulation. The principal features of the x-ray crystallographic analysis are reproduced, and the amplitudes of atomic displacements appear to be determined by the degree of exposure of the atoms. An analysis of the hydrogen bonds shows a correlation between the average length of a bond and the fluctuation in that length and reveals a rocking motion of bases in Watson-Crick guanine X cytosine base pairs. The in-plane motions of the bases are generally of larger amplitude than the out-of-plane motions, and there are correlations in the motions of adjacent bases.

Spatial constraints and group behaviour in globular proteins

A comparative analysis has been attempted on the spatial placement of amino acid residues derived from radial, ellipsoidal and exposure arrangements. The group behaviour of residues and their restraining influence in protein folding have been brought out. A study is also made on the geometry of proteins, the exposure arrangement of residues and the spatial distribution of the physical properties of the residues in globular proteins. It has been shown that the group constraints along with the information on the shape of the globular proteins would be highly useful in assigning the spatial and exposure arrangements of residues in globular proteins.

Molecular dynamics of phenylalanine transfer RNA

The atomic motions of yeast phenylalanine transfer RNA have been simulated using the molecular dynamics algorithm. Two simulations were carried out for a period of 12 picoseconds, one with a normal Van der Waals potential and the other with a modified Van der Waals potential intended to mimic the effect of solvent. An analysis of large scale motions, surface exposure, root mean square displacements, helical oscillations and relaxation mechanisms reveals the maintenance of stability in the simulated structures and the general similarity of the various dynamic features of the two simulations. The regions of conformational flexibility and rigidity for tRNA(Phe) have been shown in a quantitative measure through this approach.

Geometrical parameters and shape anisotropy in globular proteins

The anisotropy in the shape of globular proteins is derived by a comparative analysis of three types of geometrical parameters. The role of secondary structures in the design of the shape of globular proteins is worked out.

Hydrophobic packing and spatial arrangement of amino acid residues in globular proteins

Amino acid residues acquire characteristic hydrophobic environments in globular proteins. Using the crystal data on 21 proteins, a new scale of hydrophobic indices for the residues is set up. This scale provides valuable information with regard to hydrophobic domains, nucleation sites, surface domains, loop sites and the spatial positions of residues in protein molecules.