Three-dimensional structure of phoratoxin in solution: combined use of nuclear magnetic resonance, distance geometry, and restrained molecular dynamics

Deriving Structural Information from Experimentally Measured Data on Biomolecules

During the past half century, the number and accuracy of experimental techniques that can deliver values of observables for biomolecular systems have been steadily increasing. The conversion of a measured value Q^exp of an observable quantity Q into structural information is, however, a task beset with theoretical and practical problems: 1) insufficient or inaccurate values of Q^exp , 2) inaccuracies in the function Q(r→) used to relate the quantity Q to structure r→ , 3) how to account for the averaging inherent in the measurement of Q^exp , 4) how to handle the possible multiple-valuedness of the inverse r→(Q) of the function Q(r→) , to mention a few. These apply to a variety of observable quantities Q and measurement techniques such as X-ray and neutron diffraction, small-angle and wide-angle X-ray scattering, free-electron laser imaging, cryo-electron microscopy, nuclear magnetic resonance, electron paramagnetic resonance, infrared and Raman spectroscopy, circular dichroism, Förster resonance energy transfer, atomic force microscopy and ion-mobility mass spectrometry. The process of deriving structural information from measured data is reviewed with an eye to non-experts and newcomers in the field using examples from the literature of the effect of the various choices and approximations involved in the process. A list of choices to be avoided is provided.

Antimicrobial Peptides from Plants

Plant antimicrobial peptides (AMPs) have evolved differently from AMPs from other life forms. They are generally rich in cysteine residues which form multiple disulfides. In turn, the disulfides cross-braced plant AMPs as cystine-rich peptides to confer them with extraordinary high chemical, thermal and proteolytic stability. The cystine-rich or commonly known as cysteine-rich peptides (CRPs) of plant AMPs are classified into families based on their sequence similarity, cysteine motifs that determine their distinctive disulfide bond patterns and tertiary structure fold. Cystine-rich plant AMP families include thionins, defensins, hevein-like peptides, knottin-type peptides (linear and cyclic), lipid transfer proteins, α-hairpinin and snakins family. In addition, there are AMPs which are rich in other amino acids. The ability of plant AMPs to organize into specific families with conserved structural folds that enable sequence variation of non-Cys residues encased in the same scaffold within a particular family to play multiple functions. Furthermore, the ability of plant AMPs to tolerate hypervariable sequences using a conserved scaffold provides diversity to recognize different targets by varying the sequence of the non-cysteine residues. These properties bode well for developing plant AMPs as potential therapeutics and for protection of crops through transgenic methods. This review provides an overview of the major families of plant AMPs, including their structures, functions, and putative mechanisms.

Isolation and biochemical characterization of an antifungal peptide from Amaranthus hypochondriacus seeds

An antifungal peptide, Ay-AMP, was isolated from Amaranthus hypochondriacus seeds by acidic extraction and then purified by reverse-phase high-pressure liquid chromatography. The molecular mass of this peptide, as determined by mass spectrometry, is 3184 Da. The peptide belongs to the superfamily of chitin-binding proteins, containing a single cysteine/glycine-rich chitin-binding domain, and it was found that Ay-AMP degrades chitin. Ay-AMP inhibits the growth, at very low doses, of different pathogenic fungi, such as Candida albicans, Trichoderma sp., Fusarium solani, Penicillium chrysogenum, Geotrichum candidum, Aspergillus candidus, Aspergillus schraceus, and Alternaria alternata. Ay-AMP is very resistant to the effect of proteases and heating; however, it showed an antagonistic effect with CaCl2 and KCl.

Solution conformation of an immunodominant epitope in the hepatitis B virus preS2 surface antigen

We have determined the solution conformation of the major B cell epitope (residues 123-145, adrl23 hereafter) in the preS2 region of hepatitis B virus known to be associated with infection neutralization. The adrl23 shows an "L" shaped helix-turn-helix topology with two beta-turns formed by residues Ala(130)-Asp(133) and Asp(133)-Val(136) intervening the N- and C-terminal helices. The N-terminal alpha-helix consists of residues Ser(124)-Gln(129) whereas the C-terminal 3(10) helix is formed by residues Val(136)-Tyr(140). The beta-turns overlap partially with the putative "conformational" epitope. The overall topology of adrl23 is primarily maintained by hydrophobic interactions involving Phe(127), Leu(131), Leu(132), Val(136), and Tyr(140) that are clustered on one side of the molecule. An additional hydrophobic stabilization comes from Phe(141) that is buried inside the concave side of the molecule. A network of hydrogen bonds formed among Thr(125), His(128), and Arg(137) further contribute to the "boomerang-shaped" architecture of adrl23. The N-terminus of adrl23 is immobile due to a hydrogen bond between the N-terminal amide proton of Asn(123) and the hydroxyl oxygen of Thr(126). The side chains of Asp(133), Arg(135), Val(136), Leu(139), and Tyr(140) that were shown to be important for binding to a monoclonal antibody H8 mAb are surface exposed.

Crystal structure of α-hordothionin at 1.9 Å resolution

Crystal structure of ubiquitous toxin from barley alpha-hordothionin (alpha-HT) has been determined at 1.9A resolution by X-ray crystallography. The primary sequence as well as the NMR solution structure of alpha-HT firmly established that alpha-HT belongs to a family of membrane active plant toxins-thionins. Since alpha-HT crystallized in a space group (P4(1)2(1)2) that is different from the space group (I422) of previously determined alpha(1)- and beta-purothionins, and visocotoxin A3, therefore, it provided independent information on protein-protein interactions that may be relevant to the toxin mechanism. The structure of alpha-HT not only confirms overall architectural features (crambin fold) but also provides an additional confirmation of the role for crucial solute molecules, that were postulated to be directly involved in the mechanism of toxicity for thionins.

Proposal for molecular mechanism of thionins deduced from physico-chemical studies of plant toxins

We propose a molecular model for phospholipid membrane lysis by the ubiquitous plant toxins called thionins. Membrane lysis constitutes the first major effect exerted by these toxins that initiates a cascade of cytoplasmic events leading to cell death. X-ray crystallography, solution nuclear magnetic resonance (NMR) studies, small angle X-ray scattering and fluorescence spectroscopy provide evidence for the mechanism of membrane lysis. In the crystal structures of two thionins in the family, alpha(1)- and beta-purothionins (MW: approximately 4.8 kDa), a phosphate ion and a glycerol molecule are modeled bound to the protein. (31)P NMR experiments on the desalted toxins confirm phosphate-ion binding in solution. Evidence also comes from phospholipid partition experiments with radiolabeled toxins and with fluorescent phospholipids. This data permit a model of the phospholipid-protein complex to be built. Further, NMR experiments, one-dimensional (1D)- and two-dimensional (2D)-total correlation spectroscopy (TOCSY), carried out on the model compounds glycerol-3-phosphate (G3P) and short chain phospholipids, supported the predicted mode of phospholipid binding. The toxins' high positive charge, which renders them extremely soluble (>300 mg/mL), and the phospholipid-binding specificity suggest the toxin-membrane interaction is mediated by binding to patches of negatively charged phospholipids [phosphatidic acid (PA) or phosphatidyl serine (PS)] and their subsequent withdrawal. The formation of proteolipid complexes causes solubilization of the membrane and its lysis. The model suggests that the oligomerization may play a role in toxin's activation process and provides insight into the structural principles of protein-membrane interactions.

Solution conformation of alphaA-conotoxin EIVA, a potent neuromuscular nicotinic acetylcholine receptor antagonist from Conus ermineus

We report the solution three-dimensional structure of an alphaA-conotoxin EIVA determined by nuclear magnetic resonance spectroscopy and restrained molecular dynamics. The alphaA-conotoxin EIVA consists of 30 amino acids representing the largest peptide among the alpha/alphaA-family conotoxins discovered so far and targets the neuromuscular nicotinic acetylcholine receptor with high affinity. alphaA-Conotoxin EIVA consists of three distinct structural domains. The first domain is mainly composed of the Cys3-Cys11-disulfide loop and is structurally ill-defined with a large backbone root mean square deviation of 1.91 A. The second domain formed by residues His12-Hyp21 is extremely well defined with a backbone root mean square deviation of 0.52 A, thus forming a sturdy stem for the entire molecule. The third C-terminal domain formed by residues Hyp22-Gly29 shows an intermediate structural order having a backbone root mean square deviation of 1.04 A. A structurally ill-defined N-terminal first loop domain connected to a rigid central molecular stem seems to be the general structural feature of the alphaA-conotoxin subfamily. A detailed structural comparison between alphaA-conotoxin EIVA and alphaA-conotoxin PIVA suggests that the higher receptor affinity of alphaA-conotoxin EIVA than alphaA-conotoxin PIVA might originate from different steric disposition and charge distribution in the second loop "handle" motif.

The Xplor-NIH NMR molecular structure determination package

We announce the availability of the Xplor-NIH software package for NMR biomolecular structure determination. This package consists of the pre-existing XPLOR program, along with many NMR-specific extensions developed at the NIH. In addition to many features which have been developed over the last 20 years, the Xplor-NIH package contains an interface with a new programmatic framework written in C++. This interface currently supports the general purpose scripting languages Python and TCL, enabling rapid development of new tools, such as new potential energy terms and new optimization methods. Support for these scripting languages also facilitates interaction with existing external programs for structure analysis, structure manipulation, visualization, and spectral analysis.

Solution conformation of alpha-conotoxin EI, a neuromuscular toxin specific for the alpha 1/delta subunit interface of torpedo nicotinic acetylcholine receptor

A high resolution structure of alpha-conotoxin EI has been determined by (1)H NMR spectroscopy and molecular modeling. alpha-Conotoxin EI has the same disulfide framework as alpha 4/7 conotoxins targeting neuronal nicotinic acetylcholine receptors but antagonizes the neuromuscular receptor as do the alpha 3/5 and alpha A conotoxins. The unique binding preference of alpha-conotoxin EI to the alpha(1)/delta subunit interface of Torpedo neuromuscular receptor makes it a valuable structural template for superposition of various alpha-conotoxins possessing distinct receptor subtype specificities. Structural comparison of alpha-conotoxin EI with the gamma-subunit favoring alpha-conotoxin GI suggests that the Torpedo delta-subunit preference of the former originates from its second loop. Superposition of three-dimensional structures of seven alpha-conotoxins reveals that the estimated size of the toxin-binding pocket in nicotinic acetylcholine receptor is approximately 20 A (height) x 20 A (width) x 15 A (thickness).

Snake venom cardiotoxins-structure, dynamics, function and folding

Snake cardiotoxins are highly basic (pI > 10) small molecular weight (approximately 6.5 kDa), all beta-sheet proteins. They exhibit a broad spectrum of interesting biological activities. The secondary structural elements in these toxins include antiparallel double and triple stranded beta-sheets. The three dimensional structures of these toxins reveal an unique asymmetric distribution of the hydrophobic and hydrophilic amino acids. The 3D structures of closely related snake venom toxins such as neurotoxins and cardiotoxin-like basic proteins (CLBP) fail to show similar pattern(s) in the distribution of polar and nonpolar residues. Recently, many novel biological activities have been reported for cardiotoxins. However, to-date, there is no clear structure-function correlation(s) available for snake venom cardiotoxins. The aim of this comprehensive review is to summarize and critically evaluate the progress in research on the structure, dynamics, function and folding aspects of snake venom cardiotoxins.

CD and NMR structural characterization of ceratotoxins, natural peptides with antimicrobial activity

Antibacterial properties of the secretion from the female reproductive accessory glands of medfly Ceratitis capitata are mostly ascribed to the presence of two peptides, ceratotoxin A and B, which exhibit a strong activity against gram-positive and gram-negative bacterial strains, and show sequence and function homology with cecropins, melittin, and magainins. CD experiments performed in different solvents indicate the presence of a significant content of helical structures in organic solvent. Two-dimensional nmr results for ceratotoxin A in methanol show a helical behavior for the 8-25 region of the peptide. A ramachandran classification of each residue for the structures obtained from distance geometry calculations lead to the definition of four structural families in which the central segment 10-19 is always helical and differences refer to residues 8-9 and 19-23. A sequence analysis of the two ceratotoxins and a systematic search on the protein data bank revealed the occurrence of a KX-hydrophobic-hydrophobic-P motif that seems to be important for helix stabilization.

Synthesis and conformational analysis by 1H NMR and restrained molecular dynamics simulations of the cyclic decapeptide [Ser-Tyr-Ser-Met-Glu-His-Phe-Arg-Trp-Gly]

The design of enzyme mimics with therapeutic and industrial applications has interested both experimental and computational chemists for several decades. Recent advances in the computational methodology of restrained molecular dynamics, used in conjunction with data obtained from two-dimensional 1H NMR spectroscopy, make it a promising method to study peptide and protein structure and function. Several issues, however, need to be addressed in order to assess the validity of this method for its explanatory and predictive value. Among the issues addressed in this study are: the accuracy and generizability of the GROMOS peptide molecular mechanics force field; the effect of inclusion of solvent on the simulations; and the effect of different types of restraining algorithms on the computational results. The decapeptide Ser-Tyr-Ser-Met-Glu-His-Phe-Arg-Trp-Gly, which corresponds to the sequence of ACTH1-10, has been synthesized, cyclized, and studied by two-dimensional 1H NMR spectroscopy. Restrained molecular dynamics (RMD) and time-averaged restrained molecular dynamics (TARMD) simulations were carried out on four different distance-geometry starting structures in order to determine and contrast the behavior of cyclic ACTH1-10 in vacuum and in solution. For the RMD simulations, the structures did not fit the NOE data well, even at high values of the restraining potential. The TARMD simulation method, however, was able to give structures that fit the NOE data at high values of the restraining potential. In both cases, inclusion of explicit solvent molecules in the simulation had little effect on the quality of the fit, although it was found to dampen the motion of the cyclic peptide. For both simulation techniques, the number and size of the NOE violations increased as the restraining potential approached zero. This is due, presumably, to inadequacies in the force field. Additional TARMD vacuum-phase simulations, run with a larger memory length or with a larger sampling size (16 additional distance-geometry structures), yielded no significantly different results. The computed data were then analyzed to help explain the sparse NOE data and poor chymotryptic activity of the cyclic peptide. Cyclic ACTH1-10, which contains the functional moieties of the catalytic triad of chymotrypsin, was evaluated as a potential mimic of chymotrypsin by measurement of the rate of hydrolysis of esters of L- and D-phenylalanine. The poor rate of hydrolysis is attributed to the flexibility of the decapeptide, the motion of the side chains, which result in the absence of long-range NOEs, the small size of the macrocycle relative to that of the substrate, and the inappropriate orientation of the Gly, His, and Ser residues. The results demonstrate the utility of this method in computer-aided molecular design of cyclic peptides and suggest structural modifications for future work based on a larger and more rigid peptide framework.

Solution structure of the epsilon-aminohexanoic acid complex of human plasminogen kringle 1

The solution structure of the human plasminogen kringle 1 domain complexed to the antifibrinolytic drug 6-aminohexanoic acid (epsilon Ahx) was obtained on the basis of 1H-NMR spectroscopic data and dynamical simulated annealing calculations. Two sets of structures were derived starting from (a) random coil conformations and (b) the (mutated) crystallographic structure of the homologous prothrombin kringle 1. The two sets display essentially the same backbone folding (pairwise root-mean-square deviation, 0.15 nm) indicating that, regardless of the initial structure, the data is sufficient to locate a conformation corresponding to an essentially unique energy minimum. The conformations of residues connected to prolines were localized to energetically preferred regions of the Ramachandran map. The Pro30 peptide bond is proposed to be cis. The ligand-binding site of the kringle 1 is a shallow cavity composed of Pro33, Phe36, Trp62, Tyr64, Tyr72 and Tyr74. Doubly charged anionic and cationic centers configured by the side chains of Asp55 and Asp57, and Arg34 and Arg71, respectively, contribute to anchoring the zwitterionic epsilon Ahx molecule at the binding site. The ligand exhibits closer contacts with the kringle anionic centers (approximately 0.35 nm average O...H distance between the Asp55/Asp57 carboxylate and ligand amino groups) than with the cationic ones (approximately 0.52 nm closest O...H distances between the ligand carboxylate and the Arg34/Arg71 guanidino groups). The epsilon Ahx hydrocarbon chain rests flanked by Pro33, Tyr64, Tyr72 and Tyr74 on one side and Phe36 on the other. Dipolar (Overhauser) connectivities indicate that the ligand aliphatic moiety establishes close contacts with the Phe36 and Trp62 aromatic rings. The computed structure suggests that the epsilon Ahx molecule adopts a kinked conformation when complexed to kringle 1, effectively shortening its dipole length to approximately 0.65 nm.

Tertiary structure of two-electron reduced Megasphaera elsdenii flavodoxin and some implications, as determined by two-dimensional 1H-NMR and restrained molecular dynamics

The tertiary structure of the non-crystallizable two-electron-reduced Megasphaera elsdenii flavodoxin (15 kDa, 137 amino acid residues) has been determined using nuclear Overhauser enhancement restraints extracted from two-dimensional 1H-NMR spectra. A tertiary structure satisfying the experimental restraints very well (maximum NOE violation of 66 pm) was obtained with use of restrained molecular dynamics, using 509 distance restraints (including one non-NOE) on a starting structure modeled from the crystal structure of one-electron-reduced Clostridium MP flavodoxin. The protein consists of a central parallel beta-sheet surrounded on both sides by two alpha-helices. The flavin is positioned at the periphery of the molecule. The tertiary structure of the protein is highly defined with the exception of the flavin. The latter is expected to result from performing the restrained molecular dynamics simulation without water molecules and without proper charges on the flavin. The flavin, including the phosphate, the ribityl side chain and the isoalloxazine ring, is solvent accessible under the experimental conditions used and evidenced by a two-dimensional amide exchange experiment. This accessibility is expected to be important in the redox potential regulation of the semiquinone/hydroquinone couple of the protein. The amide exchange against deuterons and several typical line shapes in the two-dimensional NMR spectra are consistent with the structure generated. The structure is discussed in detail.

The sampling properties of some distance geometry algorithms applied to unconstrained polypeptide chains: a study of 1830 independently computed conformations

In this paper we study the statistical geometry of ensembles of poly (L-alanine) conformations computed by several different distance geometry algorithms. Since basic theory only permits us to predict the statistical properties of such ensembles a priori when the distance constraints have a very simple form, the only constraints used for these calculations are those necessary to obtain reasonable bond lengths and angles, together with a lack of short- and long-range atomic overlaps. The geometric properties studied include the squared end-to-end distance and radius of gyration of the computed conformations, in addition to the usual rms coordinate and phi/psi angle deviations among these conformations. The distance geometry algorithms evaluated include several variations of the well-known embed algorithm, together with optimizations of the torsion angles using the ellipsoid and variable target function algorithms. The conclusions may be summarized as follows: First, the distribution with which the trial distances are chosen in most implementations of the embed algorithm is not appropriate when no long-range upper bounds on the distances are present, because it leads to unjustifiably expanded conformations. Second, chosing the trial distances independently of one another leads to a lack of variation in the degree of expansion, which in turn produces a relatively low rms square coordinate difference among the members of the ensemble. Third, when short-range steric constraints are present, torsion angle optimizations that start from conformations obtained by choosing their phi/psi angles randomly with a uniform distribution between -180 degrees and +180 degrees do not converge to conformations whose angles are uniformly distributed over the sterically allowed regions of the phi/psi plane. Finally, in an appendix we show how the sampling obtained with the embed algorithm can be substantially improved upon by the proper application of existing methodology.

Refinement of the NMR structures for acyl carrier protein with scalar coupling data

Structure determination of small proteins using NMR data is most commonly pursued by combining NOE derived distance constraints with inherent constraints based on chemical bonding. Ideally, one would make use of a variety of experimental observations, not just distance constraints. Here, coupling constant constraints have been added to molecular mechanics and molecular dynamics protocols for structure determination in the form of a psuedoenergy function that is minimized in a search for an optimum molecular conformation. Application is made to refinement of a structure for a 77 amino acid protein involved in fatty acid synthesis, Escherichia coli acyl carrier protein (ACP). 54 3JHN alpha coupling constants, 12 coupling constants for stereospecifically assigned side chain protons, and 450 NOE distance constraints were used to calculate the 3-D structure of ACP. A three-step protocol for a molecular dynamics calculation is described, in analogy to the protocol previously used in molecular mechanics calculations. The structures calculated with the molecular mechanics approach and the molecular dynamics approach using a rigid model for the protein show similar molecular energies and similar agreement with experimental distance and coupling constant constraints. The molecular dynamics approach shows some advantage in overcoming local minimum problems, but only when a two-state averaging model for the protein was used, did molecular energies drop significantly.

A Monte Carlo simulation study of the influence of internal motions on the molecular conformation deduced from two-dimensional NMR experiments

Monte Carlo methods have been used to simulate internal motions of aromatic protons of an oligonucleotide at the nanosecond time scale. Each proton is allowed to fluctuate about its equilibrium position. The longitudinal cross-relaxation rates of such a system of spins have been determined by computing the appropriate correlation functions. Then the interproton distances have been deduced according to the procedure generally used in two-dimensional nmr techniques (nuclear Overhauser effect spectroscopy--NOESY) and compared to the true values. The influence of the amplitude A and of the internal rotational diffusion constant Dint characterizing the dynamics of the system has been checked for in-phase and for uncorrelated motions. It is shown that for the investigated models the distances deduced from NOESY experiments may be under- or overestimated, depending strongly on the values of A and Dint. Furthermore, the cross-relaxation rate of a couple of protons is very sensitive to the correlation level of the motions of both protons.

Determination of three-dimensional structures of proteins and nucleic acids in solution by nuclear magnetic resonance spectroscopy

Nuclear magnetic resonance (NMR) spectroscopy has evolved over the last decade into a powerful method for determining three-dimensional structures of biological macromolecules in solution. Key advances have been the introduction of two-dimensional experiments, high-field superconducting magnets, and computational procedures for converting the NMR-derived interproton distances and torsion angles into three-dimensional structures. This article outlines the methodology employed, describes the major NMR experiments necessary for the spectral analysis of macromolecules, and discusses the computational approaches employed to date. The present state of the art is illustrated using a variety of examples, and future developments are indicated.

Three-dimensional structure of acyl carrier protein in solution determined by nuclear magnetic resonance and the combined use of dynamical simulated annealing and distance geometry

The solution conformation of acyl carrier protein from Escherichia coli (77 residues) has been determined on the basis of 423 interproton-distance restraints and 32 hydrogen-bonding restraints derived from NMR measurements. A total of nine structures were computed using a hybrid approach combining metric matrix distance geometry and dynamic simulated annealing. The polypeptide fold is well defined with an average backbone atomic root-mean-square difference of 0.20 +/- 0.03 nm between the final nine converged structures and the mean structure obtained by averaging their coordinates. The principal structural motif is composed of three helices: 1 (residues 3-12), 2 (residues 37-47) and 4 (residues 65-75) which line a hydrophobic cavity. Helices 2 and 4 are approximately parallel to each other and anti-parallel at an angle of approximately equal to 150 degrees to helix 1. The smaller helix 3 (residues 56-63) is at an angle of approximately equal to 100 degrees to helix 4.

Three-dimensional solution structure of plastocyanin from the green alga Scenedesmus obliquus

The solution conformation of plastocyanin from the green alga Scenedesmus obliquus has been determined from distance and dihedral angle constraints derived by nuclear magnetic resonance (NMR) spectroscopy. Structures were generated with distance geometry and restrained molecular dynamics calculations. A novel molecular replacement method was also used with the same NMR constraints to generate solution structures of S. obliquus plastocyanin from the x-ray structure of the homologous poplar protein. Scenedesmus obliquus plastocyanin in solution adopts a beta-barrel structure. The backbone conformation is well defined and is similar overall to that of poplar plastocyanin in the crystalline state. The distinctive acidic region of the higher plant plastocyanins, which functions as a binding site for electron transfer proteins and inorganic complexes, differs in both shape and charge in S. obliquus plastocyanin.

Determination of three-dimensional structures of proteins from interproton distance data by hybrid distance geometry-dynamical simulated annealing calculations

A new hybrid distance space-real space method for determining three-dimensional structures of proteins on the basis of interproton distance restraints is presented. It involves the following steps: (i) the approximate polypeptide fold is obtained by generating a set of substructures comprising only a small subset of atoms by projection from multi-dimensional distance space into three-dimensional cartesian coordinate space using a procedure known as 'embedding'; (ii) all remaining atoms are then added by best fitting extended amino acids one residue at a time to the substructures; (iii) the resulting structures are used as the starting point for real space dynamical simulated annealing calculations. The latter involve heating the system to a high temperature followed by slow cooling in order to overcome potential barriers along the pathway towards the global minimum region. This is carried out by solving Newton's equations of motion. Unlike conventional restrained molecular dynamics, however, the non-bonded interactions are represented by a simple van der Waals repulsion term. The method is illustrated by calculations on crambin (46 residues) and the globular domain of histone H5 (79 residues). It is shown that the hybrid method is more efficient computationally and samples a larger region of conformational space consistent with the experimental data than full metric matrix distance geometry calculations alone, particularly for large systems.

Calculation of protein conformation by the build-up procedure. Application to bovine pancreatic trypsin inhibitor using limited simulated nuclear magnetic resonance data

Low-energy conformations of a set of tetrapeptides derived from the small protein bovine pancreatic trypsin inhibitor (BPTI) were generated by a build-up procedure from the low-energy conformations of single amino acid residues. At each stage, various-size fragments were built up from all combinations of smaller ones, the total energies were then minimized, and the low-energy conformations were retained for the next stage. The energies of the tetrapeptides were re-ordered by including the effects of hydration. No information other than the amino acid sequence was used to obtain the low-energy conformations of the hydrated tetrapeptides. The latter were then supplemented with a limited set of simulated NMR distance information, derived from the X-ray structure of BPTI, to provide a basis for building the rest of the whole protein molecule by the same procedure. A total of 189 upper bounds, plus 12 pairs of upper and lower bounds pertaining to the location of the three disulfide bonds in this molecule, were used. Four sets of conformations of the entire molecule were generated by utilizing different combinations of smaller fragments. It was possible to obtain low-energy conformations with small rms deviations, 1.1 to 1.4 A for the alpha-carbons, from the structure derived by X-ray diffraction. The average deviations of the backbone dihedral angles were also low, viz. 23 degrees to 26 degrees.

Variable-target-function and build-up procedures for the calculation of protein conformation. Application to bovine pancreatic trypsin inhibitor using limited simulated nuclear magnetic resonance data

An implementation of the variable-target-function procedure, first introduced by Braun and Go [W. Braun and N. Go, J. Mol. Biol. 186, 611-626 (1985)], has been used to generate conformations of the small protein bovine pancreatic trypsin inhibitor (BPTI), given a limited set of simulated data that could be obtained by nuclear magnetic resonance (NMR) techniques. A hybrid strategy was also used to calculate conformations of BPTI, given the same information. In the hybrid strategy, low-energy structures of medium-size fragments (decapeptides) of BPTI were generated using the variable-target-function method, followed by restrained energy optimization. The low-energy conformations were used as a basis to build up the complete fifty-eight-residue BPTI molecule. By using the variable-target-function approach, in which energy considerations were not introduced until full conformations of the entire BPTI molecule had been generated, it was not possible to obtain calculated structures with rms deviations from the X-ray conformation of less than 1.6 A for the alpha-carbons. On the other hand, with the hybrid strategy, which involved the consideration of realistic energy terms in the early stages of the calculations, it was possible to calculate low-energy conformations of BPTI with rms deviations from the X-ray structure of 1.06 to 1.50 A for the alpha-carbons. When the rms deviations were computed along the amino acid sequence, it was found that there was a good correlation between deviations among the calculated structures and deviations from the X-ray structure.

Test of circular dichroism (CD) methods for crambin and CD-assisted secondary structure prediction of its homologous toxins

Methods that analyze protein circular dichroism (CD) spectra for fractions of secondary structure are evaluated for the plant protein crambin, which has a known high-resolution crystal structure. In addition, a two-step secondary structure prediction scheme is presented and used for the toxins homologous to crambin, shown by others to have secondary structures similar to crambin. The test of CD spectral analysis methods with the protein crambin employed two computer programs and several CD basis sets. Crambin's crystal structure, known to 0.945A resolution (Hendrickson, W.A., Teeter, M.M. Nature 290:107-113, 1981), allows accurate evaluation of results. Analysis with the protein spectra basis sets (Provencher, S.W., Glöckner, J. Biochemistry 20:33-37, 1981) as modified (Manavalan, P., Johnson, W.C., Jr. Anal. Biochem. 167:76-85, 1987) agreed most closely with crambin's crystal structure. This method was then applied to the CD spectra of the membrane-active toxins homologous to crambin (alpha 1- and beta-purothionin, phoratoxin A and B, and viscotoxin A3 and B). The new program SEQ (pronounced "seek") was developed to assign the secondary structure along the protein chain in a hierarchical fashion and applied to the plant toxins. The method constrained the secondary structure fractions to those from CD analysis and combined standard statistical methods with amphipathic helix location. Both CD-arrived secondary structure percentages and sequence assignment indicate that the viscotoxins are structurally most similar to crambin. Purothionin's secondary structure was predicted to be fundamentally similar to crambin's with a difference at the start of the first helix. This assignment agreed with Raman and NMR analyses of purothionin and lends validity to the method presented here. Differences from the NMR in the CD secondary structure fraction analysis for phoratoxin suggest interference in the CD from tryptophan residues.

Combined procedure of distance geometry and restrained molecular dynamics techniques for protein structure determination from nuclear magnetic resonance data: application to the DNA binding domain of lac repressor from Escherichia coli

The technique of two-dimensional nuclear magnetic resonance (2D-NMR) has recently assumed an active role in obtaining information on structures of polypeptides, small proteins, sugars, and DNA fragments in solution. In order to generate spatial structures from the atom-atom distance information obtained by the NMR method, different procedures have been developed. Here we introduce a combined procedure of distance geometry (DG) and molecular dynamics (MD) calculations for generating 3D structures that are consistent with the NMR data set and have reasonable internal energies. We report the application of the combined procedure on the lac repressor DNA binding domain (headpiece) using a set of 169 NOE and 17 "hydrogen bond" distance constraints. Eight of ten structures generated by the distance geometry algorithm were refined within 10 ps MD simulation time to structures with low internal energies that satisfied the distance constraints. Although the combination of DG and MD was designed to combine the good sampling properties of the DG algorithm with an efficient method of lowering the internal energy of the molecule, we found that the MD algorithm contributes significantly to the sampling as well.