Fast Structure-Based Assignment of 15N HSQC Spectra of Selectively 15N-Labeled Paramagnetic Proteins

A novel strategy for fast NMR resonance assignment of N-15 HSQC spectra of proteins is presented. It requires the structure coordinates of the protein, a paramagnetic center, and one or more residue-selectively N-15-labeled samples. Comparison of sensitive undecoupled N-15 HSQC spectra recorded of paramagnetic and diamagnetic samples yields data for every cross-peak on pseudocontact shift, paramagnetic relaxation enhancement, cross-correlation between Curie-spin and dipole-dipole relaxation, and residual dipolar coupling. Comparison of these four different paramagnetic quantities with predictions from the three-dimensional structure simultaneously yields the resonance assignment and the anisotropy of the susceptibility tensor of the paramagnetic center. The method is demonstrated with the 30 kDa complex between the N-terminal domain of the epsilon subunit and the theta subunit of Escherichia Coll DNA polymerase III. The program PLATYPUS was developed to perform the assignment, provide a measure of reliability of the assignment, and determine the susceptibility tensor anisotropy.

Structure determination of the three disulfide bond isomers of alpha-conotoxin GI: A model for the role of disulfide bonds in structural stability

The three possible disulfide bonded isomers of alpha-conotoxin GI have been selectively synthesised and their structures determined by H-1 NMR spectroscopy. alpha-Conotoxin GI derives from the venom of Conus geographus and is a useful neuropharmacological tool as it selectively binds to the nicotinic acetylcholine receptor (nAChR), a ligand-gated ion channel involved in nerve signal transmission. The peptide has the sequence ECCNPACGRHYSC-NH2, and the three disulfide bonded isomers are referred to as GI(2-7;3-13), GI(2-13;3-7) and GI(2-3;7-13). The NMR structure for the native isomer GI(2-7;3-13) is of excellent quality, with a backbone pairwise RMSD of 0.16 Angstrom for a family of 35 structures, and comprises primarily a distorted 3(10),, helix between residues 5 to 11. The two non-native isomers exhibit multiple conformers in solution, with the major populated forms being different in structure both from each other and from the native form. Structure-activity relationships for the native GI(2-7;3-13) as well as the role of the disulfide bonds on folding and stability of the three isomers are examined. It is concluded that the disulfide bonds in alpha-conotoxin GI play a crucial part in determining both the structure and stability of the peptide. A trend for increased conformational heterogeneity was observed in the order of GI(2-7;3-13) < GI(2-13;3-7) < GI(2-3;7-13). It was found that the peptide bond joining Cys2 to Cys3 in GI(2-3;7-13) is predominantly trans, rather than cis as theoretically predicted. These structural data are used to interpret the varying nAChR binding of the non-native forms. A model for the binding of native GI(2-7;3-13) to the mammalian nAChR is proposed, with an alpha-subunit binding face made up of Cys2, Asn4, Pro5, Ala6 and Cys7 and a selectivity face, comprised of Arg9 and His10. These two faces orient the molecule between the alpha and delta subunits of the receptor. The structure of the CCNPAC sequence of the native GI(2-7;3-13) is compared to the structure of the identical sequence from the toxic domain of heat-stable enterotoxins, which forms part of the receptor binding region of the enterotoxins, but which has a different disulfide connectivity. (C) 1998 Academic Press Limited.

Flush air data system calibration using numerical simulation

The use of computational fluid dynamics simulations for calibrating a flush air data system is described, In particular, the flush air data system of the HYFLEX hypersonic vehicle is used as a case study. The HYFLEX air data system consists of nine pressure ports located flush with the vehicle nose surface, connected to onboard pressure transducers, After appropriate processing, surface pressure measurements can he converted into useful air data parameters. The processing algorithm requires an accurate pressure model, which relates air data parameters to the measured pressures. In the past, such pressure models have been calibrated using combinations of flight data, ground-based experimental results, and numerical simulation. We perform a calibration of the HYFLEX flush air data system using computational fluid dynamics simulations exclusively, The simulations are used to build an empirical pressure model that accurately describes the HYFLEX nose pressure distribution ol cr a range of flight conditions. We believe that computational fluid dynamics provides a quick and inexpensive way to calibrate the air data system and is applicable to a broad range of flight conditions, When tested with HYFLEX flight data, the calibrated system is found to work well. It predicts vehicle angle of attack and angle of sideslip to accuracy levels that generally satisfy flight control requirements. Dynamic pressure is predicted to within the resolution of the onboard inertial measurement unit. We find that wind-tunnel experiments and flight data are not necessary to accurately calibrate the HYFLEX flush air data system for hypersonic flight.

Determination of chemical shielding tensor of an indole carbon and application of tryptophan orientation of a membrane peptide

Using tryptophan C-13-enriched at the C-4 (C epsilon(3)) of the indole, the orientation of the C epsilon(3) chemical shift tensor relative to the C epsilon(3)-H dipolar axis was determined from the C-13 chemical shift/C-13-H-1 dipolar 2D NMR powder pattern. The principal values obtained were 208, 137 and 15 ppm with sigma(33) perpendicular to the indole plane, and sigma(11) (least shielded direction) 5 degrees off the C epsilon(3)-H bond toward C xi(3). The side off the C epsilon(3)-H bond was determined by comparing the reduced chemical shift anisotropies obtained by solid-state NMR and from molecular dynamics calculations of [4-C-13] tryptophans in gramicidin A aligned in phospholipid membranes. (C) 1999 Elsevier Science B.V. All rights reserved.

Solution structure of alpha-conotoxin ImI by H-1 nuclear magnetic resonance

alpha-Conotoxin ImI derives from the venom of Conus imperialis and is the first and only small-peptide ligand that selectively binds to the neuronal alpha(7) homopentameric subtype of the nicotinic acetylcholine receptor (nAChR). This receptor subtype is a possible drug target for several neurological disorders. The cysteines are connected in the pairs Cys2-Cys8 and Cys3-Cys12, To date it is the only alpha-conotoxin with a 4/3 residue spacing between the cysteines, The structure of ImI has been determined by H-1 NMR spectroscopy in aqueous solution, The NMR structure is of high quality, with a backbone pairwise rmsd of 0.34 Angstrom for a family of 19 structures, and comprises primarily a series of nested beta turns. Addition of organic solvent does not perturb the solution structure. The first eight residues of ImI are identical to the larger, but related, conotoxin EpI and adopt a similar structure, despite a truncated second loop. Residues important for binding of ImI to the alpha 7 nAChR are all clustered on one face of the molecule. Once further binding data for EPI and ImI are available, the ImI structure will allow for design of novel alpha(7) nAChR-specific agonists and antagonists with a wide range of potential pharmaceutical applications.

NMR studies of the low-density lipoprotein receptor-binding peptide of apolipoprotein E bound to dodecylphosphocholine micelles

Circular dichroism and NMR spectroscopy have been used to determine the structure of the low-density lipoprotein (LDL) receptor-binding peptide, comprising residues 130-152, of the human apolipoprotein E. This peptide has little persistent three-dimensional structure in solution, but when bound to micelles of dodecylphosphocholine (DPC) it adopts a predominantly alpha-helical structure. The three-dimensional structure of the DPC-bound peptide has been determined by using H-1-NMR spectroscopy: the structure derived from NOE-based distance constraints and restrained molecular dynamics is largely helical. The derived phi and psi angle order parameters show that the helical structure is well defined but with some flexibility that causes the structures not to be superimposable over the full peptide length. Deuterium exchange experiments suggest that many peptide amide groups are readily accessible to the solvent, but those associated with hydrophobic residues exchange more slowly, and this helix is thus likely to be positioned on the surface of the DPC micelles. In this conformation the peptide has one hydrophobic face and two that are rich in basic amino acid side chains. The solvent-exposed face of the peptide contains residues previously shown to be involved in binding to the LDL receptor.

MiAMP1, a novel protein from Macadamia integrifolia adopts a Greek key beta-barrel fold unique amongst plant antimicrobial proteins

MiAMP1 is a recently discovered 76 amino acid residue, highly basic protein from the nut kernel of:Macadamia integrifolia which possesses no sequence homology to any known protein and inhibits the growth of several microbial plant pathogens in vitro while having no effect on mammalian or plant cells. It is considered to be a potentially useful tool for the genetic engineering of disease resistance in transgenic crop plants and for the design of new fungicides. The three-dimensional structure of MiAMP1 was determined through homonuclear and heteronuclear (N-15) 2D NMR spectroscopy and subsequent simulated annealing calculations with the ultimate aim of understanding the structure-activity relationships of the protein. MiAMP1 is made up of eight beta-strands which are arranged in two Greek key motifs. These Greek key motifs associate to form a Greek key beta-barrel. This structure is unique amongst plant antimicrobial proteins and forms a new class which we term the beta-barrelins. Interestingly, the structure of MiAMP1 bears remarkable similarity to a yeast killer toxin from Williopsis mrakii. This toxin acts by inhibiting beta-glucan synthesis and thereby cell wall construction in sensitive strains of yeast. The structural similarity of MiAMP1 and WmKT, which originate from plant and fungal phyla respectively, may reflect a similar mode of action. (C) 1999 Academic Press.

Computational chemistry on Fujitsu vector-parallel processors: Development and performance of applications software

In this and a preceding paper, we provide an introduction to the Fujitsu VPP range of vector-parallel supercomputers and to some of the computational chemistry software available for the VPP. Here, we consider the implementation and performance of seven popular chemistry application packages. The codes discussed range from classical molecular dynamics to semiempirical and ab initio quantum chemistry. All have evolved from sequential codes, and have typically been parallelised using a replicated data approach. As such they are well suited to the large-memory/fast-processor architecture of the VPP. For one code, CASTEP, a distributed-memory data-driven parallelisation scheme is presented. (C) 2000 Published by Elsevier Science B.V. All rights reserved.

Three-dimensional structure of RK-1: A novel alpha-defensin peptide

NMR spectroscopy and simulated annealing calculations have been used to determine the three-dimensional structure of RK-1, an antimicrobial peptide from rabbit kidney recently discovered from homology screening based on the distinctive physicochemical properties of the corticostatins/defensins. RK-1 consists of 32 residues, including six cysteines arranged into three disulfide bonds. It exhibits antimicrobial activity against Escherichia coli and activates Ca2+ channels in vitro. Through its physicochemical similarity, identical cysteine spacing, and linkage to the corticostatins/defensins, it was presumed to be a member of this family. However, RK-1 lacks both a large number of arginines in the primary sequence and a high overall positive charge, which are characteristic of this family of peptides. The three-dimensional solution structure, determined by NMR, consists of a triple-stranded antiparallel beta -sheet and a series of turns and is similar to the known structures of other alpha -defensins. This has enabled the definitive classification of RK-1 as a member of this family of antimicrobial peptides. Ultracentrifuge measurements confirmed that like rabbit neutrophil defensins, RK-1 is monomeric in solution, in contrast to human neutrophil defensins, which are dimeric.

Circular proteins in plants - Solution structure of a novel macrocyclic trypsin inhibitor from Momordica cochinchinensis

Much interest has been generated by recent reports on the discovery of circular (i.e. head-to-tail cyclized) proteins in plants. Here we report the three-dimensional structure of one of the newest such circular proteins, MCoTI-II, a novel trypsin inhibitor from Momordica cochinchinensis, a member of the Cucurbitaceae plant family. The structure consists of a small beta -sheet, several turns, and a cystine knot arrangement of the three disulfide bonds. Interestingly, the molecular topology is similar to that of the plant cyclotides (Craik, D. J., Daly, N. L., Bond, T., and Waine, C. (1999) J. Mol. Biol, 294, 1327-1336), which derive from the Rubiaceae and Violaceae plant families, have antimicrobial activities, and exemplify the cyclic cystine knot structural motif as part of their circular backbone. The sequence, biological activity, and plant family of MCoTI-II are all different from known cyclotides. However, given the structural similarity, cyclic backbone, and plant origin of MCoTI-II, we propose that MCoTI-II can be classified as a new member of the cyclotide class of proteins. The expansion of the cyclotides to include trypsin inhibitory activity and a new plant family highlights the importance and functional variability of circular proteins and the fact that they are more common than has previously been believed, Insights into the possible roles of backbone cyclization have been gained by a comparison of the structure of MCoTI-II with the homologous acyclic trypsin inhibitors CMTI-I and EETI-II from the Cucurbitaceae plant family.

Solid-State Experimental and Theoretical Investigation of the Ammonium Salt of Croconate Violet, a Pseudo-Oxocarbon Ion

The present work describes the crystal structure, vibrational spectra, and theoretical calculations of ammonium salts of 3,5-bis-(dicyanomethylene)cyclopentane-1,2,4-trionate, (NH(4))(2)(C(11)N(4)O(3)) [(NH(4))(2)CV], also known as ammonium croconate violet. This compound crystallizes in triclinic P (1) over bar and contains two water molecules per unit formula. The crystal packing is stabilized by hydrogen bonds involving water molecules and ammonium cations, giving rise to a 3D polymeric arrangement. In this structure, a pi-stacking interaction is not observed, as the smaller centroid-centroid distance is 4.35 angstrom. Ab initio electronic structure calculations under periodic boundary conditions were performed to predict vibrational and electronic properties. The vibrational analysis was used to assist the assignments of the Raman and infrared bands. The solid structure was optimized and characterized as a minimum in the potential-energy surface. The stabilizing intermolecular hydrogen bonds in the crystal Structure were characterized by difference charge-density analysis. The analysis of the density of states of (NH(4))(2)CV gives an energy gap of 1.4 eV with a significant contribution of carbon and nitrogen 2p states for valence and conduction bands.

Interfacial interactions and structure of organic-inorganic nanohybrids

Understanding the interfacial interactions and structure is important to better design and application of organic-inorganic nanohybrids. This paper presents our recent molecular dynamic studies on organoclays and polymer nanocomposites, including the layering behavior of organoclays, structural and dynamic properties of dioctadecyldimethyl ammoniums in organoclays, and interfacial interactions and structure of polyurethane nanocomposites. The results demonstrate that the layering behaviors of organoclays are closely related to the chain length of quaternary alkyl ammoniums and cation exchangeable capacity of clays. In addition to typical layered structures such as monolayer, bilayer and pseudo-trilayer, a pseudo-quadrilayer structure was also observed in organoclays modified with dioctadecyldimethyl ammoniums (DODDMA). In such a structure, alkyl chains do not lie flat within a single layer but interlace, and also jump to the next layer or even the next nearest layer. Moreover, the diffusion constants of nitrogen and methylene atoms increase with the temperature and methelene towards the tail groups. For polyurethane nanocomposite, the van der Waals interaction between apolar alkyl chains and soft segments of polyurethane predominates the interactions between organoclay and polyurethane. Different from most bulk polyurethane systems, there is no distinct phase-separated structure for the polyurethane.

Solid-state NMR structure determination of melittin in a lipid environment

Solid-state C-13 NMR spectroscopy was used to investigate the three-dimensional structure of melittin as lyophilized powder and in ditetradecylphosphatidylcholine (DTPC) membranes. The distance between specifically labeled carbons in analogs [1-C-13]Gly3-[2-C-13]Ala4, [1-C-13]Gly3-[2-C-13]Leu6, [1-C-13]Leu13-[2-C-13]Ala15, [2-C-13]Leu13-[1-C-13]Ala15, and [1-C-13]Leu13-[2-C-13]Leu16 was measured by rotational resonance. As expected, the internuclear distances measured in [1-C-13]Gly3-[2-C-13]Ala4 and [1-C-13]Gly3-[2-C-13]Leu6 were consistent with alpha -helical structure in the N-terminus irrespective of environment. The Internuclear distances measured in [1-C-13]Leu13-[2-C-13]Ala15, [2-C-13]Leu13-[1-C-13]Ala15, and [1-C-13]Leu13-[2-C-13]Leu16 revealed, via molecular modeling, some dependence upon environment for conformation in the region of the bend in helical structure induced by Pro14. A slightly larger interhelical angle between the N- and C-terminal helices was indicated for peptide in dry or hydrated gel state DTPC (139 degrees -145 degrees) than in lyophilized powder (121 degrees -139 degrees) or crystals (129 degrees). The angle, however, is not as great as deduced for melittin in aligned bilayers of DTPC in the liquid-crystalline state (similar to 160 degrees) (R. Smith, F. Separovic, T. J. Milne, A. Whittaker, F. M. Bennett, B. A. Cornell, and A. Makriyannis, 1994, J. Mol, Biol 241:456-466). The study illustrates the utility of rotational resonance in determining local structure within peptide-lipid complexes.

Ring flexibility within tricyclic antidepressant drugs

The internal flexibility of the central seven-membered ring of a series of tricyclic antidepressant drugs (TCAs), imipramine {l}, amitriptyline {2}, doxepin {3}, and dothiepin {4}, has been investigated by H-1 and C-13 nuclear magnetic (NMR) techniques. Two dynamic processes were examined: ring inversion and bridge flexing. H-1 NMR lineshape analysis was used to obtain ring inversion barriers for 2-4. These studies yielded energy barriers of 14.3, 16.7, and 15.7 +/- 0.6 kcal/mol for the hydrochloride salts of doxepin, dothiepin, and amitriptyline, respectively. The barriers for the corresponding free bases were lower by 0.6 kcal/mol on average. (CT1)-C-13 relaxation measurements were used to determine the degree of bridge flexing associated with the central seven-membered ring for all four compounds. By fitting the T-1 data to a two-state jump model, lifetimes and amplitudes of rapid bridge flexing motions were determined. The results show that imipramine has the fastest rate of bridge flexing, followed by amitriptyline, doxepin, and dothiepin. The pharmacological profiles of the TCAs are complex and they interact with many receptor sites, resulting in numerous side effects and a general lack of understanding of their precise mode of action in different anxiety-related disorders. They all have similar three-dimensional structures, which makes it difficult to rationalize their differing relative potency in different assays/clinical settings. However, the clear finding here that there are significantly different degrees of internal mobility suggests that molecular dynamics should be an additional factor considered when trying to understand the mode of action of this clinically important family of molecules. (C) 2001 Wiley-Liss, Inc. and the American Pharmaceutical Association J Pharm Sci 90:713-721, 2001.

Subspace wavepacket evolution with Newton polynomials

Time-dependent wavepacket evolution techniques demand the action of the propagator, exp(-iHt/(h)over-bar), on a suitable initial wavepacket. When a complex absorbing potential is added to the Hamiltonian for combating unwanted reflection effects, polynomial expansions of the propagator are selected on their ability to cope with non-Hermiticity. An efficient subspace implementation of the Newton polynomial expansion scheme that requires fewer dense matrix-vector multiplications than its grid-based counterpart has been devised. Performance improvements are illustrated with some benchmark one and two-dimensional examples. (C) 2001 Elsevier Science B.V. All rights reserved.