Three-dimensional structure of acyl carrier protein determined by NMR pseudoenergy and distance geometry calculations

Distance constraints from two-dimensional NMR cross-relaxation data are used to derive a three-dimensional structure for acyl carrier protein from Escherichia coli. Several approaches to structure determination are explored. The most successful proves to be an approach that combines the early stages of a distance geometry program with energy minimization in the presence of NMR constraints represented as pseudopotentials. Approximately 450 proton to proton distance constraints including 50 long-range constraints were included in these programs. Starting structures were generated at random by the distance geometry program and energies minimized by a molecular mechanics module to give final structures. Seven of the structures were deemed acceptable on the basis of agreement with experimentally determined distances. Root-mean-square deviations from the mean of these structures for backbone atoms range from 2 to 3 A. All structures show three roughly parallel helices with hydrophobic residues facing inward and hydrophilic residues facing outward. A hydrophobic cleft is recognizable and is identified as a likely site for acyl chain binding.     

Refinement of the thrombin-bound structure of a hirudin peptide by a restrained electrostatically driven Monte Carlo method.

Energy refinement of the structure of a linear peptide, hirudin56-65, bound to thrombin was carried out using a conformational search method in combination with restrained minimization. Five conformations originated from nmr data and distance geometry calculations having a similar global folding pattern but quite different backbone conformations were used as the starting structures. As a result of this approach, a series of low-energy conformations compatible with a set of upper and lower bounds of interproton distances determined from transferred nuclear Overhauser effects were found. A comparison among the lowest energy conformations of each run showed that the combination of energy refinement plus distance constraints led to a very well-defined structure for both the backbone and the side chains of the last 7 residues of the polypeptide. Furthermore, the low-energy conformations generated with this technique contain a segment of 3(10)-helix involving the last 5 residues at the COOH terminal end.     

Determining the three-dimensional fold of a protein from approximate constraints

We propose a new approach for calculating the three-dimensional (3D) structure of a protein from distance and dihedral angle constraints derived from experimental data. We suggest that such constraints can be obtained from experiments such as tritium planigraphy, chemical or enzymatic cleavage of the polypeptide chain, paramagnetic perturbation of nuclear magnetic resonance (NMR) spectra, measurement of hydrogen-exchange rates, mutational studies, mass spectrometry, and electron paramagnetic resonance. These can be supplemented with constraints from theoretical prediction of secondary structures and of buried/exposed residues. We report here distance geometry calculations to generate the structures of a test protein Staphylococcal nuclease (STN), and the HIV-1 rev protein (REV) of unknown structure. From the available 3D atomic coordinates of STN, we set up simulated data sets consisting of varying number and quality of constraints, and used our group's Self Correcting Distance Geometry (SECODG) program DIAMOD to generate structures. We could generate the correct tertiary fold from qualitative (approximate) as well as precise distance constraints. The root mean square deviations of backbone atoms from the native structure were in the range of 2.0 A to 8.3 A, depending on the number of constraints used. We could also generate the correct fold starting from a subset of atoms that are on the surface and those that are buried. When we used data sets containing a small fraction of incorrect distance constraints, the SECODG technique was able to detect and correct them. In the case of REV, we used a combination of constraints obtained from mutagenic data and structure predictions. DIAMOD generated helix-loop-helix models, which, after four self-correcting cycles, populated one family exclusively. The features of the energy-minimized model are consistent with the available data on REV-RNA interaction. Our method could thus be an attractive alternative for calculating protein 3D structures, especially in cases where the traditional methods of X-ray crystallography and multidimensional NMR spectroscopy have been unsuccessful.     

Structural comparison of acyl carrier protein in acylated and sulfhydryl forms by two-dimensional 1H NMR spectroscopy

Sequence-specific assignments of 1H NMR resonances are obtained for the backbone protons of Escherichia coli acyl carrier protein, acylated with an eight-carbon chain covalently attached to the prosthetic group thiol (octanoyl-ACP). Comparison of 1H-1H sequential connectivities in the NOESY spectra of octanoyl-ACP and the unacylated protein (ACPSH) indicates that secondary structure is largely conserved on acylation. Changes in resonance positions observed for certain groups of residues are interpreted in terms of a model that describes the spatial reorientation of secondary structural elements in the protein resulting from introduction of the acyl chain.     

Secondary structure of acyl carrier protein as derived from two-dimensional 1H NMR spectroscopy

Sequence-specific assignments of 1H NMR resonances were obtained for the backbone protons in acyl carrier protein (ACP) from Escherichia coli, a protein of 77 residues. The observations, in the NOESY spectra, of 1H-1H sequential and medium-range connectivities indicate the presence of three or four alpha-helical segments joined by short sequences of mixed conformations. The observations are used to refine a secondary structure model previously proposed on the basis of a Chou-Fasman algorithm [Rock, C. O., & Cronan, J. E., Jr. (1979) J. Biol. Chem. 254, 9778-9785].     

Biosynthesis of D-alanyl-lipoteichoic acid: the tertiary structure of apo-D-alanyl carrier protein.

The D-alanylation of lipoteichoic acid (LTA) allows the Gram-positive organism to modulate its surface charge, regulate ligand binding, and control the electromechanical properties of the cell wall. The incorporation of D-alanine into LTA requires the D-alanine:D-alanyl carrier protein ligase (AMP-forming) (Dcl) and the carrier protein (Dcp). The high-resolution solution structure of the 81-residue (8.9 kDa) Dcp has been determined by multidimensional heteronuclear NMR. An ensemble of 30 structures was calculated using the torsion angle dynamics approach of DYANA. These calculations utilized 3288 NOEs containing 1582 unique nontrivial NOE distance constraints. Superposition of residues 4-81 on the mean structure yields average atomic rmsd values of 0.43 +/- 0.08 and 0.86 +/- 0.09 A for backbone and non-hydrogen atoms, respectively. The solution structure is composed of three alpha-helices in a bundle with additional short 3(10)- and alpha-helices in intervening loops. Comparisons of the three-dimensional structure with the acyl carrier proteins involved in fatty acid, polyketide, and nonribosomal peptide syntheses support the conclusion that Dcp is a homologue in this family. While there is conservation of the three-helix bundle fold, Dcp has a higher enthalpy of unfolding and no apparent divalent metal binding site(s), features that distinguish it from the fatty acid synthase acyl carrier protein of Escherichia coli. This three-dimensional structure also provides insights into the D-alanine ligation site recognized by Dcl, as well as the site which may bind the poly(glycerophosphate) acceptor moiety of membrane-associated LTA.     

Calculation of protein backbone geometry from alpha-carbon coordinates based on peptide-group dipole alignment.

An algorithm is proposed for the conversion of a virtual-bond polypeptide chain (connected C alpha atoms) to an all-atom backbone, based on determining the most extensive hydrogen-bond network between the peptide groups of the backbone, while maintaining all of the backbone atoms in energetically feasible conformations. Hydrogen bonding is represented by aligning the peptide-group dipoles. These peptide groups are not contiguous in the amino acid sequence. The first dipoles to be aligned are those that are both sufficiently close in space to be arranged in approximately linear arrays termed dipole paths. The criteria used in the construction of dipole paths are: to assure good alignment of the greatest possible number of dipoles that are close in space; to optimize the electrostatic interactions between the dipoles that belong to different paths close in space; and to avoid locally unfavorable amino acid residue conformations. The equations for dipole alignment are solved separately for each path, and then the remaining single dipoles are aligned optimally with the electrostatic field from the dipoles that belong to the dipole-path network. A least-squares minimizer is used to keep the geometry of the alpha-carbon trace of the resulting backbone close to that of the input virtual-bond chain. This procedure is sufficient to convert the virtual-bond chain to a real chain; in applications to real systems, however, the final structure is obtained by minimizing the total ECEPP/2 (empirical conformational energy program for peptides) energy of the system, starting from the geometry resulting from the solution of the alignment equations. When applied to model alpha-helical and beta-sheet structures, the algorithm, followed by the ECEPP/2 energy minimization, resulted in an energy and backbone geometry characteristic of these alpha-helical and beta-sheet structures. Application to the alpha-carbon trace of the backbone of the crystallographic 5PTI structure of bovine pancreatic trypsin inhibitor, followed by ECEPP/2 energy minimization with C alpha-distance constraints, led to a structure with almost as low energy and root mean square deviation as the ECEPP/2 geometry analog of 5PTI, the best agreement between the crystal and reconstructed backbone being observed for the residues involved in the dipole-path network.     

A dynamic model for the structure of acyl carrier protein in solution

The determination of solution structures of proteins using two-dimensional NMR data is commonly based on the assumption that the structure can be represented by a single rigid conformer. We present here a procedure whereby this assumption can be relaxed and illustrate its application to acyl carrier protein from Escherichia coli, a small negatively charged protein with no internal disulfide bonds. The methodology rests on a model having two distinct conformers in dynamic equilibrium. Use of this two-state model results in a dramatic improvement in fit to cross-relaxation-derived distance constraints and a substantial lowering of molecular mechanics energies for individual conformers of acyl carrier protein. The two-state model retains the three-helix motif previously identified on the basis of a one-state structure, but substantial motion of loop regions and the C-terminal peptide, as well as partial disruption of the second helix, is suggested to occur. Support for the existence of these motions can be found in amide exchange rate and spin relaxation time data.     

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.     

Solution structure and backbone dynamics of the holo form of the frenolicin acyl carrier protein

During polyketide biosynthesis, acyl carrier proteins (ACPs) perform the central role of transferring polyketide intermediates between active sites of polyketide synthase. The 4'-phosphopantetheine prosthetic group of a holo-ACP is a long and flexible arm that can reach into different active sites and provide a terminal sulfhydryl group for the attachment of acyl groups through a thioester linkage. We have determined the solution structure and characterized backbone dynamics of the holo form of the frenolicin acyl carrier protein (fren holo-ACP) by nuclear magnetic resonance (NMR). Unambiguous assignments were made for 433 hydrogen atoms, 333 carbon atoms, and 84 nitrogen atoms, representing a total of 94.6% of the assignable atoms in this protein. From 879 meaningful NOEs and 45 angle constraints, a family of 24 structures has been calculated. The solution structure is composed of three major alpha-helices packed in a bundle with three additional short helices in intervening loops; one of the short helices slowly exchanges between two conformations. Superposition of the major helical regions on the mean structure yields average atomic rmsd values of 0.49 +/- 0.09 and 0.91 +/- 0.08 A for backbone and non-hydrogen atoms, respectively. Although the three-helix bundle fold is conserved among acyl carrier proteins involved in fatty acid synthases and polyketide synthases, a detailed comparison revealed that ACPs from polyketide biosynthetic pathways are more related to each other in tertiary fold than to their homologues from fatty acid biosynthetic pathways. Comparison of the free form of ACPs (NMR structures of fren ACP and the Bacillus subtilis ACP) with the substrate-bound form of ACP (crystal structure of butyryl-ACP from Escherichia coli) suggests that conformational exchange plays a role in substrate binding.     

The NMR solution structure of the pheromone Er-1 from the ciliated protozoan Euplotes raikovi.

The 3-dimensional structure of the pheromone Er-1 isolated from the ciliated protozoan Euplotes raikovi has been determined in aqueous solution by 1H NMR spectroscopy. The structure of this 40-residue protein was calculated with the distance geometry program DIANA on the basis of 503 upper distance constraints derived from nuclear Overhauser effects and 77 dihedral angle constraints derived from spin-spin coupling constants, and refined by restrained energy minimization with the program OPAL. The Er-1 solution structure is represented by a group of 20 conformers with an average RMS deviation relative to the mean structure of 0.55 A for the backbone atoms N, C alpha, and C', and 0.93 A for all heavy atoms of the complete polypeptide chain, residues 1-40. The molecular architecture is dominated by an up-down-up bundle of 3 alpha-helices formed by residues 2-9, 12-19, and 24-33. Although this core part coincides closely with the previously determined structure of the homologous pheromone Er-10, the C-terminal peptide segment adopts a novel conformation. This is of interest in view of previous suggestions, based on sequence comparisons, that this molecular region may be important for the different specificity of receptor recognition by different pheromones.     

Three-dimensional structure of the neurotoxin ATX Ia from Anemonia sulcata in aqueous solution determined by nuclear magnetic resonance spectroscopy

With the aid of 1H nuclear magnetic resonance (NMR) spectroscopy, the three-dimensional structure in aqueous solution was determined for ATX Ia, which is a 46 residue polypeptide neurotoxin of the sea anemone Anemonia sulcata. The input for the structure calculations consisted of 263 distance constraints from nuclear Overhauser effects (NOE) and 76 vicinal coupling constants. For the structure calculation several new or ammended programs were used in a revised strategy consisting of five successive computational steps. First, the program HABAS was used for a complete search of all backbone and chi 1 conformations that are compatible with the intraresidual and sequential NMR constraints. Second, using the program DISMAN, we extended this approach to pentapeptides by extensive sampling of all conformations that are consistent with the local and medium-range NMR constraints. Both steps resulted in the definition of additional dihedral angle constraints and in stereospecific assignments for a number of beta-methylene groups. In the next two steps DISMAN was used to obtain a group of eight conformers that contain no significant residual violations of the NMR constraints or van der Waals contacts. Finally, these structures were subjected to restrained energy refinement with a modified version of the molecular mechanics module of AMBER, which in addition to the energy force field includes potentials for the NOE distance constraints and the dihedral angle constraints. The average of the pairwise minimal RMS distances between the resulting refined conformers calculated for the well defined molecular core, which contains the backbone atoms of 35 residues and 20 interior side chains, is 1.5 +/- 0.3 A. This core is formed by a four-stranded beta-sheet connected by two well-defined loops, and there is an additional flexible loop consisting of the eleven residues 8-18. The core of the protein is stabilized by three disulfide bridges, which are surrounded by hydrophobic residues and shielded on one side by hydrophilic residues.     

Nuclear magnetic resonance solution structure of hirudin(1-51) and comparison with corresponding three-dimensional structures determined using the complete 65-residue hirudin polypeptide chain.

The three-dimensional structure of the N-terminal 51-residue domain of recombinant hirudin in aqueous solution was determined by 1H nuclear magnetic resonance (NMR) spectroscopy, and the resulting high-quality solution structure was compared with corresponding structures obtained from studies with the intact, 65-residue polypeptide chain of hirudin. On the basis of 580 distance constraints derived from nuclear Overhauser effects and 109 dihedral angle constraints, a group of 20 conformers representing the solution structure of hirudin(1-51) was computed with the program DIANA and energy-minimized with a modified version of the program AMBER. Residues 3 to 30 and 37 to 48 form a well-defined molecular core with two antiparallel beta-sheets composed of residues 14 to 16 and 20 to 22, and 27 to 31 and 36 to 40, and three reverse turns at residues 8 to 11 (type II), 17 to 20 (type II') and 23 to 26 (type II). The average root-mean-square deviation of the individual NMR conformers relative to their mean co-ordinates is 0.38 A for the backbone atoms and 0.77 A for all heavy atoms of these residues. Increased structural disorder was found for the N-terminal dipeptide segment, the loop at residues 31 to 36, and the C-terminal tripeptide segment. The solution structure of hirudin(1-51) has the same molecular architecture as the corresponding polypeptide segment in natural hirudin and recombinant desulfatohirudin. It is also closely similar to the crystal structure of the N-terminal 51-residue segment of hirudin in a hirudin-thrombin complex, with root-mean-square deviations of the crystal structure relative to the mean solution structure of 0.61 A for the backbone atoms and 0.91 A for all heavy atoms of residues 3 to 30 and 37 to 48. Further coincidence is found for the loop formed by residues 31 to 36, which shows increased structural disorder in all available solution structures of hirudin, and of which residues 32 to 35 are not observable in the electron density map of the thrombin complex. Significant local structural differences between hirudin(1-51) in solution and hirudin in the crystalline thrombin complex were identified mainly for the N-terminal tripeptide segment and residues 17 to 21. These are further analyzed in an accompanying paper.     

Three-dimensional solution structure of the HIV-1 protease complexed with DMP323, a novel cyclic urea-type inhibitor, determined by nuclear magnetic resonance spectroscopy.

The three-dimensional solution structure of the HIV-1 protease homodimer, MW 22.2 kDa, complexed to a potent, cyclic urea-based inhibitor, DMP323, is reported. This is the first solution structure of an HIV protease/inhibitor complex that has been elucidated. Multidimensional heteronuclear NMR spectra were used to assemble more than 4,200 distance and angle constraints. Using the constraints, together with a hybrid distance geometry/simulated annealing protocol, an ensemble of 28 NMR structures was calculated having no distance or angle violations greater than 0.3 A or 5 degrees, respectively. Neglecting residues in disordered loops, the RMS deviation (RMSD) for backbone atoms in the family of structures was 0.60 A relative to the average structure. The individual NMR structures had excellent covalent geometry and stereochemistry, as did the restrained minimized average structure. The latter structure is similar to the 1.8-A X-ray structure of the protease/DMP323 complex (Chang CH et al., 1995, Protein Science, submitted); the pairwise backbone RMSD calculated for the two structures is 1.22 A. As expected, the mismatch between the structures is greatest in the loops that are disordered and/or flexible. The flexibility of residues 37-42 and 50-51 may be important in facilitating substrate binding and product release, because these residues make up the respective hinges and tips of the protease flaps. Flexibility of residues 4-8 may play a role in protease regulation by facilitating autolysis.     

Solution conformation of thymosin beta 4: a nuclear magnetic resonance and simulated annealing study

The conformation of the polypeptide thymosin beta 4 in solutions of 60% (v/v) trifluoroethanol-d3 and 50% (v/v) hexafluoroisopropyl-d2 alcohol in water is investigated by nuclear magnetic resonance (NMR) spectroscopy. Under these conditions thymosin beta 4 adopts an ordered structure. By use of a combination of two-dimensional NMR techniques, the 1H NMR spectrum of thymosin beta 4 is assigned. A set of 180 approximate interproton distance constraints is derived from nuclear Overhauser enhancement (NOE) measurements. These, together with 33 phi constraints obtained for JNH alpha coupling data and the 23 psi dihedral angles identified on the basis of the pattern of short-range NOEs, form the basis of a three-dimensional structure determination by dynamical simulated annealing. The calculations are carried out starting from three initial structures, an alpha-helix, an extended beta-strand, and a mixed alpha/beta structure. Ten independent structures are computed from each starting structure by using different random number seeds for the assignments of the initial velocities. All 30 calculated structures satisfy the experimental constraints, display very small deviations from idealized covalent geometry, and possess good nonbonded contacts. Analysis of the 30 converged structures indicates that there are two helical regions extending from residues 4-16 and from residues 30-40, which are well defined both in terms of atomic root mean square differences and backbone torsion angles. For the two helical regions individually the average backbone rms difference between all pairs of structures is approximately 2 A. The two helices exhibit typical amino acid preferences for specific locations at the ends of helices.(ABSTRACT TRUNCATED AT 250 WORDS)     

The NMR solution structure of the pheromone Er-11 from the ciliated protozoan Euplotes raikovi.

The NMR solution structure of the pheromone Er-11, a 39-residue protein from the ciliated protozoan Euplotes raikovi, was calculated with the distance geometry program DIANA from 449 NOE upper distance constraints and 97 dihedral angle constraints, and the program OPAL was employed for structure refinement by molecular mechanics energy minimization in a water bath. For a group of 20 conformers used to characterize the solution structure, the average of the pairwise RMS deviations from the mean structure calculated for the backbone heavy atoms N, C alpha, and C' of residues 2-38 was 0.30 A. The molecular architecture is dominated by an up-down-up bundle of three short helices with residues 2-9, 12-19, and 22-32, which is closely similar to the previously determined structures of the homologous pheromones Er-1, Er-2, and Er-10. This finding provides structural evidence for the capability shown by these pheromones to compete with each other in binding reactions to their cell-surface receptors.     

A mammalian type I fatty acid synthase acyl carrier protein domain does not sequester acyl chains

The synthases that produce fatty acids in mammals (FASs) are arranged as large multidomain polypeptides. The growing fatty acid chain is bound covalently during chain elongation and reduction to the acyl carrier protein (ACP) domain that is then able to access each catalytic site. In this work we report the high-resolution nuclear magnetic resonance (NMR) solution structure of the isolated rat fatty acid synthase apoACP domain. The final ensemble of NMR structures and backbone (15)N relaxation studies show that apoACP adopts a single, well defined fold. On conversion to the holo form, several small chemical shift changes are observed on the ACP for residues surrounding the phosphopantetheine attachment site (as monitored by backbone (1)H-(15)N correlation experiments). However, there are negligible chemical shift changes when the holo form is modified to either the hexanoyl or palmitoyl forms. For further NMR analysis, a (13)C,(15)N-labeled hexanoyl-ACP sample was prepared and full chemical shift assignments completed. Analysis of two-dimensional F(2)-filtered and three-dimensional (13)C-edited nuclear Overhauser effect spectroscopy experiments revealed no detectable NOEs to the acyl chain. These experiments demonstrate that unlike other FAS ACPs studied, this Type I ACP does not sequester a covalently linked acyl moiety, although transient interactions cannot be ruled out. This is an important mechanistic difference between the ACPs from Type I and Type II FASs and may be significant for the modulation and regulation of these important mega-synthases.     

The NMR solution structure of the pheromone Er-2 from the ciliated protozoan Euplotes raikovi.

The NMR structure of the pheromone Er-2 from the ciliated protozoan Euplotes raikovi has been determined in aqueous solution. The structure of this 40-residue protein was calculated with the distance geometry program DIANA from 621 distance constraints and 89 dihedral angle constraints; the program OPAL was employed for the energy minimization. For a group of 20 conformers used to characterize the solution structure, the average pairwise RMS deviation from the mean structure calculated for the backbone heavy atoms N, C alpha, and C' of residues 3-37 was 0.31 A. The molecular architecture is dominated by an up-down-up bundle of 3 short helices of residues 5-11, 14-20, and 23-33, which is similar to the structures of the homologous pheromones Er-1 and Er-10. Novel structural features include a well-defined N-cap on the first helix, a 1-residue deletion in the second helix resulting in the formation of a 3(10)-helix rather than an alpha-helix as found in Er-1 and Er-10, and the simultaneous presence of 2 different conformations for the C-terminal tetrapeptide segment, i.e., a major conformation with the Leu 39-Pro 40 peptide bond in the trans form and a minor conformation with this peptide bond in the cis form.     

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.