Solution structures of cyclic and dicyclic analogues of growth hormone releasing factor as determined by two-dimensional NMR and CD spectroscopies and constrained molecular dynamics.

Solution structures were determined for a linear analogue of growth hormone releasing factor (GRF), and cyclic and dicyclic analogues in which the side chains of aspartyl and lysyl residues spaced at positions i-(i + 4) were joined to form a lactam. The four analogues were [Ala15]-GRF-(1-29)-NH2 and its cyclo8-12, cyclo21-25, and dicyclo8-12;21-25 derivatives. The peptides were studied in two solvent systems: 75% methanol/25% water at pH 6.0; and 100% water at pH 3.0. CD spectroscopy was used to assess the overall alpha-helical content. Nuclear magnetic resonance spectroscopy was used to determine the structures in more detail. Nearly complete proton resonance assignments were made for each of the peptides, in both solvents. Nuclear Overhauser effects were converted into distance constraints and applied in the molecular dynamics program CHARMM to evaluate the range of low-energy structures that satisfied the nmr data. In 75% methanol, all of the peptides are comprised of a single alpha-helical segment with fraying of one to three residues at each end. The linear analogue has a tendency to kink. In water, the analogues have two helical segments with flexible regions between them and at the termini of the peptides. The linear analogue is helical at residues 7-14 and 21-28. In the cyclo8-12 analogue, the N-terminal helical region extends to include residues 7-19, while the other helical region is slightly shortened. In the cyclo21-25 analogue, the C-terminal helical region is extended to include residues 19-28, while the N-terminal helical region is destabilized. The dicyclic analogue has the largest N-terminal helix, spanning residues 7-20, but its helical segment at residues 21-28 is not well ordered. All of the analogues exhibit substantial biological activity. The cyclic and dicyclic analogues show dramatically increased resistance to degradation during incubation with human plasma. The i-(i + 4) lactam, therefore, appears to be a synthetic means of stabilizing a local alpha-helical conformation, which may be of general use in the design of active, stable peptides.     

Solution structure of alpha-conotoxin ImI determined by two-dimensional NMR spectroscopy.

The three-dimensional structure of alpha-conotoxin ImI, a potent antagonist targeting the neuronal alpha7 subtype of nicotinic acetylcholine receptor (nAChR), has been investigated by NMR spectroscopy. On the basis of 181 experimental constraints, a total of 25 converged structures were obtained. The average pairwise atomic root mean square difference is 0.40+/-0.11 A for the backbone atoms. The resulting structure indicates the presence of two successive type I beta-turns and a 310 helix for residues Cys2-Cys8 and Ala9-Arg11, respectively, and shows a significant structural similarity to that of alpha-conotoxin PnIA, which is also selective for the neuronal nAChR.     

Solution structure of duplex DNA containing an extrahelical abasic site analog determined by NMR spectroscopy and molecular dynamics.

Translesional DNA synthesis past abasic sites proceeds with the preferential incorporation of dAMP opposite the lesion and, depending on the sequence context, one or two base deletions. High-resolution NMR spectroscopy and molecular dynamics simulations were used to determine the three-dimensional structure of a DNA heteroduplex containing a synthetic abasic site (tetrahydrofuran) residue positioned in a sequence that promotes one base deletions. Analysis of NMR spectra indicates that the stem region of the duplex adopts a right-handed helical structure and the glycosidic torsion angle is in anti orientation for all residues. NOE interactions establish Watson-Crick alignments for all canonical base pairs of the duplex. Measurement of distance interactions at the lesion site shows the abasic residue excluded from the helix. Restrained molecular dynamics simulations generated three-dimensional models in excellent agreement with the spectroscopic data. These structures show a regular duplex region and a slight bend at the lesion site. The tetrahydrofuran residue extrudes from the helix and is highly flexible. The model reported here, in conjunction with a previous study performed on abasic sites, explains the structural bias of one-base deletion mutations.     

Solution structure of the Lewis x oligosaccharide determined by NMR spectroscopy and molecular dynamics simulations.

The Lewis x (Lex) determinant is a trisaccharide fragment that has been implicated as a specific differentiation antigen, as a tumor antigen, and as a key component of the ligand for the endothelial leukocyte adhesion molecule, ELAM-1. High-resolution nuclear magnetic resonance spectroscopy shows it to have a relatively rigid structure. Only a small range of glycosidic dihedral angles in the trisaccharide produce simulated nuclear Overhauser effect spectra agreeing with data measured for the human milk pentasaccharide, lacto-N-fucopentaose-3, which contains the Lex determinant. Independently, the same average structure for the Lex determinant arises from in vacuo molecular dynamics simulations. The proposed conformation of the Lex trisaccharide is very similar to that recently determined for the closely related Lea trisaccharide. In agreement with the recent finding that both sialylated Lea and Lex react with ELAM-1, the results presented here show that the Lea and Lex determinants contain very similar carbohydrate domains.     

Conformational studies of a benzodiazepine-like peptide in SDS micelles by circular dichroism, 1H NMR and molecular dynamics simulation

Summary The conformation of a benzodiazepine-like decapeptide corresponding to the YLGYLEQLLR fragment of a casein has been examined in a sodium dodecyl sulfate micellar medium using circular dichroism, two-dimensional¹H NMR spectroscopy and restrained molecular dynamics simulation. The decapeptide adopts an amphipathic 3₁₀-helicoid structure in which the E⁶...R¹⁰ ionic bridge stabilizes the C-terminus.     

Solution structure of a cathelicidin-derived antimicrobial peptide, CRAMP as determined by NMR spectroscopy.

CRAMP was identified from a cDNA clone derived from mouse femoral marrow cells as a member of cathelicidin-derived antimicrobial peptides. This peptide shows potent antimicrobial activity against gram-positive and gram-negative bacteria but no hemolytic activity against human erythrocytes. CRAMP was known to cause rapid permeabilization of the inner membrane of Escherichia coli. In this study, the structure of CRAMP in TFE/H2O (1 : 1, v/v) solution was determined by CD and NMR spectroscopy. CD spectra showed that CRAMP adopts a mainly alpha-helical conformation in TFE/H2O solution, DPC micelles, SDS micelles and liposomes, whereas it has a random structure in aqueous solution. The tertiary structure of CRAMP in TFE/H2O (1 : 1, v/v), as determined by NMR spectroscopy, consists of two amphipathic alpha-helices from Leu4 to Lys10 and from Gly16 to Leu33. These two helices are connected by a flexible region from Gly11 to Gly16. Previous analysis of series of fragments composed of various portion of CRAMP revealed that an 18-residue fragment with the sequence from Gly16 to Leu33 was found to retain antibacterial activity. Therefore, the amphipathic alpha-helical region from Gly16 to Leu33 of CRAMP plays important roles in spanning the lipid bilayers as well as its antibiotic activity. Based on this structure, novel antibiotic peptides having strong antibiotic activity, with no hemolytic effect will be developed.     

Equilibrium structure and folding of a helix-forming peptide: circular dichroism measurements and replica-exchange molecular dynamics simulations

We have performed experimental measurements and computer simulations of the equilibrium structure and folding of a 21-residue alpha-helical heteropeptide. Far ultraviolet circular dichroism spectroscopy is used to identify the presence of helical structure and to measure the thermal unfolding curve. The observed melting temperature is 296 K, with a folding enthalpy of -11.6 kcal/mol and entropy of -39.6 cal/(mol K). Our simulations involve 45 ns of replica-exchange molecular dynamics of the peptide, using eight replicas at temperatures between 280 and 450 K, and the program CHARMM with a continuum solvent model. In a 30-ns simulation started from a helical structure, conformational equilibrium at all temperatures was reached after 15 ns. This simulation was used to calculate the peptide melting curve, predicting a folding transition with a melting temperature in the 330-350 K range, enthalpy change of -10 kcal/mol, and entropy change of -30 cal/(mol K). The simulation results were also used to analyze the peptide structural fluctuations and the free-energy surface of helix unfolding. In a separate 15-ns replica-exchange molecular dynamics simulation started from the extended structure, the helical conformation was first attained after approximately 2.8 ns, and equilibrium was reached after 10 ns of simulation. These results showed a sequential folding process with a systematic increase in the number of hydrogen bonds until the helical state is reached, and confirmed that the alpha-helical state is the global free-energy minimum for the peptide at low temperatures.     

The secondary structure of human amyloid deposits as determined by circular dichroism spectroscopy

The deposition of amyloid fibrils has been associated with a diversity of pathologies including plasma cell dyscrasias, chronic inflammatory diseases, and several types of neurological diseases including Alzheimer's disease. Using circular dichroism spectroscopy, the secondary structure of a human amyloid protein deposit from a patient with plasma cell (light chain)-associated amyloidosis (amyloidosis AL) has been determined. The protein contains 52% beta structure, which is consistent with these depositions arising from the aberrant catabolism of immunoglobulin light chains, which are rich in beta sheets. The protein was also found to contain 20% alpha-helix, suggesting that partial refolding may occur during amyloidogenesis.     

Solution conformations of cyclosporins and magnesium-cyclosporin complexes determined by vibrational circular dichroism

Vibrational circular dichroism (VCD) spectroscopy was used to investigate the solution conformations of cyclosporins A, C, D, G, and H in CDCl(3), in the amide I and NH/OH-stretching regions, and their corresponding magnesium complexes in CD(3)CN, in the amide I region. VCD spectra are sensitive to the chiral arrangement of Cdbond;O and NH bonds in this cyclic undecapeptide. Calculations of molecular geometries, as well as IR and VCD intensities of model cyclosporin fragments that include the intramolecular hydrogen bonds of the crystal conformations of cyclosporins A and H (CsA and CsH), were carried out at the density functional theory (DFT; BPW91 functional/6-31G* basis set) level. The good agreement between IR and VCD spectra from experiment and DFT calculations provides evidence that the crystal conformation of CsA is dominant in CDCl(3) solution; CsH, however, assumes both an intramolecularly hydrogen-bonded crystal conformation and more open forms in solution. Comparisons of the experimental and calculated VCD spectra in the NH/OH-stretching region of the noncomplexed cyclosporins indicate that conformers with both free and hydrogen-bonded NH and OH groups are present in solution. Differences between the IR and VCD spectra for the metal-free and magnesium-complexed cyclosporins are indicative of strong interactions between cyclosporins and magnesium ions.     

Structural determination of the vasoactive intestinal peptide by two-dimensional H-NMR spectroscopy

The structure of the vasoactive intestinal peptide 1-28 in 40% 2,2,2-trifluoroethanol was investigated by two-dimensional 1H-nmr spectroscopy. All 1H resonances, except the gamma, delta, and epsilon protons of the lysine residues, could be sequentially assigned. Numerous intraresidual as well as short-range interresidual nuclear Overhauser effect spectroscopy connectivities were observed. Using a variable-target function minimization, a molecular model consisting of two helical stretches involving residues 7-15 and 19-27 connected by a region of undefined structure was calculated. The existence of an undefined structure between residues 16 and 18 confers mobility to the peptide molecule.     

Solution structure of termite-derived antimicrobial peptide,spinigerin, as determined in SDS micelle by NMR spectroscopy

Spinigerin is a linear antibacterial peptide derived from a termite insect. It consists of 25 amino acids and is devoid of cysteines. Spinigerin displays good lytic activities against Gram-positive and Gram-negative bacteria, but has no hemolytic activities against human erythrocytes. In this study, we present a three-dimensional solution structure of spinigerin in SDS micelles. According to CD data spinigerin has an alpha-helical conformation in the presence of TFE, DPC micelles, and SDS micelles. The three-dimensional structure of spinigerin as determined by NMR spectroscopy contains a stable alpha-helix from Lys4 to Thr23. Spinigerin (4-21), an 18-residue fragment from Lys4 to Leu21, contains a similar content of alpha-helical structure compared to native spinigerin and was found to retain antibacterial activity, too. Therefore, this alpha-helical structure and the strong electrostatic attraction between four Lys and three Arg residues in spinigerin and the negatively charged polar head groups of the phospholipids on the membrane surface play important roles in disrupting membrane and subsequent cell death.     

Solution structure of human calcitonin gene-related peptide by 1H NMR and distance geometry with restrained molecular dynamics

The structure of human calcitonin gene-related peptide 1 (hCGRP-1) has been determined by 1H NMR in a mixed-solvent system of 50% trifluoroethanol/50% H2O at pH 3.7 and 27 degrees C. Complete resonance assignment was achieved by using two-dimensional methods. Distance restraints for structure calculations were obtained by semiquantitative analysis of intra- and interresidue nuclear Overhauser effects; in addition, stereospecific or X1 rotamer assignments were obtained for certain side chains. Structures were generated from the distance restraints by distance geometry, followed by refinement using molecular dynamics, and were compared with experimental NH-C alpha H coupling constants and amide hydrogen exchange data. The structure of hCGRP-1 in this solvent comprises an amino-terminal disulfide-bonded loop (residues 2-7) leading into a well-defined alpha-helix between residues 8 and 18; thereafter, the structure is predominantly disordered, although there are indications of a preference for a turn-type conformation between residues 19 and 21. Comparison of spectra for the homologous hCGRP-2 with those of hCGRP-1 indicates that the conformations of these two forms are essentially identical.     

Conformational studies of a benzodiazepine-like peptide in SDS micelles by circular dichroism, 1H NMR and molecular dynamics simulation

The conformation of a benzodiazepine-like decapeptide corresponding tothe YLGYLEQLLR fragment of a casein has been examined in a sodium dodecylsulfate micellar medium using circular dichroism, two-dimensional¹H NMR spectroscopy and restrained molecular dynamicssimulation. The decapeptide adopts an amphipathic 3₁₀-helicoidstructure in which theE⁶···R¹⁰ ionic bridgestabilizes the C-terminus.     

The secondary structure of bacteriorhodopsin determined by Raman and circular dichroism spectroscopy

The secondary structure of bacterio-opsin (BO), the retinal free protein-component of bacteriorhodopsin (BR), has been determined by Raman spectroscopy. Additional circular dichroism (CD) measurements have revealed only negligible conformational differences between BO in apomembranes and BR in purple membranes. Therefore, the secondary structure of BR was derived from the Raman data of BO. The protein conformation was determined to consist of 72-82% helices, 2-11% beta-strands, and 11-17% beta-turns. Only about half of the helical structures correspond to alpha 1-helices, the other half possess non-alpha 1-helical structures. According to the analysis of the Raman data, the derived secondary structure of BR was obtained with high reliability for all structure classes which can be distinguished by this method within the given uncertainty range. This is a remarkable difference from recently published secondary structural data derived from CD measurements where the helix content was reported to be between 50 and 80%. The inherent experimental and methodological uncertainties of the CD-technique leading to such a range of variation are critically discussed in comparison to the method of Raman spectroscopy. The combined application of Raman and CD spectroscopy, as performed here, is demonstrated to be a substantial improvement in the secondary structure determination of retinal-containing membrane proteins. On the basis of our results, some of the recently proposed structural models of BR with a beta-strand content of more than 11% can be ruled out.     

Solution structure of neuronal bungarotoxin determined by two-dimensional NMR spectroscopy: calculation of tertiary structure using systematic homologous model building, dynamical simulated annealing, and restrained molecular dynamics.

Neuronal bungarotoxin has previously been shown, using two-dimensional 1H NMR spectroscopy, to have a triple-stranded antiparallel beta-sheet structure which dimerizes in solution [Oswald, R.E., Sutcliffe, M.J., Bamberger, M., Loring, R.H., Braswell, E., & Dobson, C.M. (1991) Biochemistry 30, 4901-4909]. In this paper, structural calculations are described which use the 582 experimentally measured NOE restraints in conjunction with 27 phi-angle restraints from J-value measurements. The positions of the N-terminal region and C-terminal region were poorly defined in the calculated structures with respect to the remainder of the structure. The region of the structure containing the triple-stranded beta-sheet was, however, well defined and similar to that found in the structure of homologous alpha-bungarotoxin (45% amino acid identity). The experimental restraints did not result in a well-defined dimer interface region because of the small number of NOEs which could be identified in this region. An approach was therefore adopted which produced model structures based to varying degrees on the alpha-bungarotoxin structure. Fourteen different structures were generated in this manner and subsequently used as starting points for refinement using dynamical simulated annealing followed by restrained molecular dynamics. This approach, which combines NMR data and homologous model building, has enabled a family of structures to be proposed for the dimeric molecule. In particular, Phe49 has been identified as possibly playing an important role in dimer formation, this residue in one chain interacting with the corresponding residue in the adjacent chain.     

Solution structures of alpha-conotoxin G1 determined by two-dimensional NMR spectroscopy

Two-dimensional NMR data have been used to generate solution structures of alpha-conotoxin G1, a potent peptide antagonist of the acetylcholine receptor. Structural information was obtained in the form of proton-proton internuclear distance constraints, and initial structures were produced with a distance geometry algorithm. Energetically more favorable structures were generated by using the distance geometry structures as input for a constrained energy minimization program. The results of both of these calculations indicate that the overall backbone conformation of the molecule is well-defined by the NMR data whereas the side-chain conformations are generally less well-defined. The main structural features derived from the NMR data were the presence of tight turns centered on residues Pro5 and Arg9. The solution structures are compared with previous proposed models of conotoxin G1, and the NMR data are interpreted in conjunction with chemical modification studies and structural properties of other antagonists of the acetylcholine receptor to gain insight into structure-activity relationships in these peptide toxins.     

Folding-unfolding of goat alpha-lactalbumin studied by stopped-flow circular dichroism and molecular dynamics simulations

Folding reaction of goat alpha-lactalbumin has been studied by stopped-flow circular dichroism and molecular dynamics simulations. The effects of four single mutations and a double mutation on the stability of the protein under a native condition were studied. The mutations were introduced into residues located at a hydrophobic core in the alpha-domain of the molecule. Here we show that an amino acid substitution (T29I) increases the native-state stability of goat alpha-lactalbumin against the guanidine hydrochloride-induced unfolding by 3.5 kcal/mol. Kinetic refolding and unfolding of wild-type and mutant goat alpha-lactalbumin measured by stopped-flow circular dichroism showed that the local structure around the Thr29 side chain was not constructed in the transition state of the folding reaction. To characterize the local structural change around the Thr29 side chain to an atomic level of resolution, we performed high-temperature (at 400 K and 600 K) molecular dynamics simulations and studied the structural change at an initial stage of unfolding observed in the simulation trajectories. The Thr29 portion of the molecule experienced structural disruption accompanied with the loss of inter-residue contacts and with the water molecule penetration in the 400-K simulation as well as in four of the six 600-K simulations. Disruption of the N-terminal portion was also observed and was consistent with the results of kinetic refolding/unfolding experiments shown in our previous report.     

Solution structure of neuronal bungarotoxin determined by two-dimensional NMR spectroscopy: sequence-specific assignments, secondary structure, and dimer formation

The solution structure of neuronal bungarotoxin (nBgt) has been studied by using two-dimensional 1H NMR spectroscopy. Sequence-specific assignments for over 95% of the backbone resonances and 85% of the side-chain resonances have been made by using a series of two-dimensional spectra at four temperatures. From these assignments over 75% of the NOESY spectrum has been assigned, which has in turn provided 582 distance constraints. Twenty-seven coupling constants (NH-alpha CH) were determined from the COSY spectra, which have provided dihedral angle constraints. In addition, hydrogen exchange experiments have suggested the probable position of hydrogen bonds. The NOE constraints, dihedral angle constraints, and the rates of amide proton exchange suggest that a triple-stranded antiparallel beta sheet is the major component of secondary structure, which includes 25% of the amino acid residues. A number of NOE peaks were observed that were inconsistent with the antiparallel beta-sheet structure. Because we have confirmed by sedimentation equilibrium that nBgt exists as a dimer, we have reinterpreted these NOE constraints as intermolecular interactions. These constraints suggest that the dimer consists of a six-stranded antiparallel beta sheet (three from each monomer), with residues 55-59 forming the dimer interface.