Solution three-dimensional structure of surfactin: A cyclic lipopeptide studied by 1H-nmr,distance geometry,and molecular dynamics

The solution three-dimensional structure of the protonated [Leu7]-surfactin, an hepta-peptide extracted from Bacillus subtilis, has been determined from two-dimensional ¹Hnmr performed in ²H₆-dimethylsulfoxide and combined with molecular modeling. Experimental data included 9 coupling constants, 61 nuclear Overhauser effect derived distances, NH temperature coefficients, and ¹³C relaxation times. Two distance geometry (DISMAN) protocols converged toward models of the structure and the best of them were refined by restrained and unrestrained molecular dynamics (GROMOS). Two structures in accord with the set of experimental constraints are presented. Both are characterized by a “horse saddle” topology for ring atoms on which are attached the two polar Glu and Asp side chains showing an orientation clearly opposite to that of the C_11–13 aliphatic chain. Amphipathic and surface properties of surfactin are certainly related to the existence of such minor polar and a major hydrophobic domains. The particular “claw” configuration of acidic residues observed in surfactin gives important clues for the understanding of its cation binding and transporting ability. © 1994 John Wiley & Sons, Inc.     

Phosphoenolpyruvate carboxylase: three-dimensional structure and molecular mechanisms

Phosphoenolpyruvate carboxylase (PEPC; EC 4.1.1.31) catalyzes the irreversible carboxylation of phosphoenolpyruvate (PEP) to form oxaloacetate and Pi using Mg2+ or Mn2+ as a cofactor. PEPC plays a key role in photosynthesis by C4 and Crassulacean acid metabolism plants, in addition to its many anaplerotic functions. Recently, three-dimensional structures of PEPC from Escherichia coli and the C4 plant maize (Zea mays) were elucidated by X-ray crystallographic analysis. These structures reveal an overall square arrangement of the four identical subunits, making up a "dimer-of-dimers" and an eight-stranded beta barrel structure. At the C-terminal region of the beta barrel, the Mn2+ and a PEP analog interact with catalytically essential residues, confirmed by site-directed mutagenesis studies. At about 20A from the beta barrel, an allosteric inhibitor (aspartate) was found to be tightly bound to down-regulate the activity of the E. coli enzyme. In the case of maize C4-PEPC, the putative binding site for an allosteric activator (glucose 6-phosphate) was also revealed. Detailed comparison of the various structures of E. coli PEPC in its inactive state with maize PEPC in its active state shows that the relative orientations of the two subunits in the basal "dimer" are different, implicating an allosteric transition. Dynamic movements were observed for several loops due to the binding of either an allosteric inhibitor, a metal cofactor, a PEP analog, or a sulfate anion, indicating the functional significance of these mobile loops in catalysis and regulation. Information derived from these three-dimensional structures, combined with related biochemical studies, has established models for the reaction mechanism and allosteric regulation of this important C-fixing enzyme.     

Application of molecular dynamics with interproton distance restraints to three-dimensional protein structure determination. A model study of crambin

The applicability of restrained molecular dynamics for the determination of three-dimensional protein structures on the basis of short interproton distances (less than 4 A) that can be realistically determined from nuclear magnetic resonance measurements in solution is assessed. The model system used is the 1.2 A resolution crystal structure of the 46 residue protein crambin, from which a set of 240 approximate distance restraints, divided into three ranges (2.5 +/- 0.5, 3.0+0.5(-1.0) and 4 +/- 1 A), is derived. This interproton distance set comprises 159 short-range ([i-j] less than or equal to 5) and 56 ([i-j] greater than 5) long-range inter-residue distances and 25 intra-residue distances. Restrained molecular dynamics are carried out using a number of different protocols starting from two initial structures: a completely extended beta-strand; and an extended structure with two alpha-helices in the same positions as in the crystal structure (residues 7 to 19, and 23 to 30) and all other residues in the form of extended beta-strands. The root-mean-square (r.m.s.) atomic differences between these two initial structures and the crystal structure are 43 A and 23 A, respectively. It is shown that, provided protocols are used that permit the secondary structure elements to form at least partially prior to folding into a tertiary structure, convergence to the correct final structure, both globally and locally, is achieved. The r.m.s. atomic differences between the converged restrained dynamics structures and the crystal structure range from 1.5 to 2.2 A for the backbone atoms and from 2.0 to 2.8 A for all atoms. The r.m.s. atomic difference between the X-ray structure and the structure obtained by first averaging the co-ordinates of the converged restrained dynamics structures is even smaller: 1.0 A for the backbone atoms and 1.6 A for all atoms. These results provide a measure with which to judge future experimental results on proteins whose crystal structures are unknown. In addition, from an examination of the dynamics trajectories, it is shown that the convergence pathways followed by the various simulations are different.     

The three-dimensional structure of a Plasmodium falciparum cyclophilin in complex with the potent anti-malarial cyclosporin A

Cyclosporin A (CsA) is a potent anti-malarial compound in vitro and in vivo in mice though better known for its immunosuppressive properties in humans. Crystal structures of wild-type and a double mutant Plasmodium falciparum cyclophilin (PfCyP19 and mPfCyP19) complexed with CsA have been determined using diffraction terms to a resolution of 2.1 A (1 A=0.1 nm). The wild-type has a single PfCyP19/CsA complex per asymmetric unit in space group P1 and refined to an R-work of 0.15 and R-free of 0.19. An altered cyclophilin, with two accidental mutations, Phe120 to Leu in the CsA binding pocket and Leu171 to Trp at the C terminus, presents two complexes per asymmetric unit in the orthorhombic space group P2(1)2(1)2. This refined to an R-work of 0.18 and R-free 0.21. The mutations were identified from the crystallographic analysis and the C-terminal alteration helps to explain the different crystal forms obtained. PfCyP19 shares approximately 61 % sequence identity with human cyclophilin A (hCyPA) and the structures are similar, consisting of an eight-stranded antiparallel beta-barrel core capped by two alpha-helices. The fold creates a hydrophobic active-site, the floor of which is formed by side-chains of residues from four antiparallel beta-strands and the walls from loops and turns. We identified C-H.O hydrogen bonds between the drug and protein that may be an important feature of cyclophilins and suggest a general mode of interaction between hydrophobic molecules. Comparisons with cyclophilin-dipeptide complexes suggests that a specific C-H.O hydrogen bonding interaction may contribute to ligand binding. Residues Ser106, His99 and Asp130, located close to the active site and conserved in most cyclophilins, are arranged in a manner reminiscent of a serine protease catalytic triad. A Ser106Ala mutant was engineered to test the hypothesis that this triad contributes to CyP function. Mutant and wild-type enzymes were found to have similar catalytic properties.     

Induced peptide conformations in different antibody complexes: molecular modeling of the three-dimensional structure of peptide-antibody complexes using NMR-derived distance restraints.

Intramolecular interactions in bound cholera toxin peptide (CTP3) in three antibody complexes were studied by two-dimensional transferred NOE spectroscopy. These measurements together with previously recorded spectra that show intermolecular interactions in these complexes were used to obtain restraints on interproton distances in two of these complexes (TE32 and TE33). The NMR-derived distance restraints were used to dock the peptide into calculated models for the three-dimensional structure of the antibody combining site. It was found that TE32 and TE33 recognize a loop comprising the sequence VPGSQHID and a beta-turn formed by the sequence VPGS. The third antibody, TE34, recognizes a different epitope within the same peptide and a beta-turn formed by the sequence IDSQ. Neither of these two turns was observed in the free peptide. The formation of a beta-turn in the bound peptide gives a compact conformation that maximizes the contact with the antibody and that has greater conformational freedom than alpha-helix or beta-sheet secondary structure. A total of 15 antibody residues are involved in peptide contacts in the TE33 complex, and 73% of the contact area in the antibody combining site consists of the side chains of aromatic amino acids. A comparison of the NMR-derived models for CTP3 interacting with TE32 and TE33 with the previously derived model for TE34 reveals a relationship between amino acid sequence and combining site structure and function. (a) The three aromatic residues that interact with the peptide in TE32 and TE33 complexes, Tyr 32L, Tyr 32H, and Trp 50H, are invariant in all light chains sharing at least 65% identity with TE33 and TE32 and in all heavy chains sharing at least 75% identity with TE33. Although TE34 differs from TE32 and TE33 in its fine specificity, these aromatic residues are conserved in TE34 and interact with its antigen. Therefore, we conclude that the role of these three aromatic residues is to participate in nonspecific hydrophobic interactions with the antigen. (b) Residues 31, 31c, and 31e of CDR1 of the light chain interact with the antigen in all three antibodies that we have studied. The amino acids in these positions in TE34 differ from those in TE32 and TE33, and they are involved in specific polar interactions with the antigen. (c) CDR3 of the heavy chain varies considerably both in length and in sequence between TE34 and the two other anti-CTP3 antibodies. These changes modify the shape of the combining site and the hydrophobic and polar interactions of CDR3 with the peptide antigen.     

The three-dimensional structure of alpha1-purothionin in solution: combined use of nuclear magnetic resonance, distance geometry and restrained molecular dynamics

The determination of the three-dimensional solution structure of α1-purothionin using a combination of metric matrix distance geometry and restrained molecular dynamics calculations based on n.m.r. data is presented. The experimental data comprise complete sequence-specific proton resonance assignments, a set of 310 approximate interproton distance restraints derived from nuclear Overhauser effects, 27 Ø backbone torsion angle restraints derived from vicinal coupling constants, 4 distance restraints from hydrogen bonds and 12 distance restraints from disulphide bridges. The average atomic rms difference between the final nine converged structures and the mean structure obtained by averaging their coordinates is 1.5 ± 0.1 å for the backbone atoms and 2.0 ± 0.1 å for all atoms. The overall shape of α1-purothionin is that of the capital letter L, similar to that of crambin, with the longer arm comprising two approximately parallel α-helices and the shorter arm a strand and a mini anti-parallel β sheet.     

Digestive lipases: from three-dimensional structure to physiology

Human gastric lipase (HGL) is a lipolytic enzyme that is secreted by the chief cells located in the fundic part of the stomach. HGL plays an important role in lipid digestion, since it promotes the subsequent hydrolytic action of pancreatic lipase in duodenal lumen. Physiological studies have shown that HGL is able of acting not only in the highly acid stomach environment but also in the duodenum in synergy with human pancreatic lipase (HPL). Recombinant HGL (r-HGL) was expressed in the baculovirus/insect cell system in the form of an active protein with a molecular mass of 45 kDa. The specific activities of r-HGL were found to be similar to that of the native enzyme when tested on various triacylglycerol (TG) substrates. The 3-D structure of r-HGL was the first solved within the mammalian acid lipase family. This globular enzyme (379 residues) shows a new feature, different from the other known lipases structures, which consists of a core domain having the alpha/beta hydrolase fold and a cap domain including a putative 'lid' of 30 residues covering the active site of the lipase (closed conformation). HPL is the major lipolytic enzyme involved in the digestion of dietary TG. HPL is a 50 kDa glycoprotein which is directly secreted as an active enzyme. HPL was the first mammalian lipase to be solved structurally, and it revealed the presence of two structural domains: a large N-terminal domain (residues 1-336) and a smaller C-terminal domain (residues 337-449). The large N-terminal domain belongs to the alpha/beta hydrolase fold and contains the active site. A surface loop called the lid domain (C237-C261) covers the active site in the closed conformation of the lipase. The 3-D structure of the lipase-procolipase complex illustrates how the procolipase might anchor the lipase at the interface in the presence of bile salts: procolipase binds to the C-terminal domain of HPL and exposes the hydrophobic tips of its fingers at the opposite site of its lipase-binding domain. These hydrophobic tips help to bring N-terminal domain into close conformation with the interface where the opening of the lid domain probably occurs. As a result of all these conformational changes, the open lid and the extremities of the procolipase form an impressive continuous hydrophobic plateau, extending over more than 50 A. This surface might able to interact strongly with a lipid-water interface. The biochemical, histochemical and clinical studies as well as the 3-D structures obtained will be a great help for a better understanding of the structure-function relationships of digestive lipases.     

Calculating three-dimensional changes in protein structure due to amino-acid substitutions: the variable region of immunoglobulins

A procedure (coupled perturbation procedure, CPP) is introduced as a specific method for calculating the detailed three-dimensional structure of a protein molecule which has a number of amino-acid substitutions relative to some previously determined "parent" protein structure. The accuracy of the procedure is tested by calculating the conformation of a region of the human immunoglobulin fragment Fab Kol based on the analogous region of the human immunoglobulin fragment Fab New. Both structures have previously been determined crystallographically. The calculated model is accurate to the extent that both of the sequence differences in the region are modeled correctly and that conformational changes in a number of nearby residues are correctly identified. CPP is shown to give better results than other commonly used modeling procedures when applied to the same problem.     

Recovering normal networks from shortest inter-taxa distance information

Phylogenetic networks are a type of leaf-labelled, acyclic, directed graph used by biologists to represent the evolutionary history of species whose past includes reticulation events. A phylogenetic network is tree–child if each non-leaf vertex is the parent of a tree vertex or a leaf. Up to a certain equivalence, it has been recently shown that, under two different types of weightings, edge-weighted tree–child networks are determined by their collection of distances between each pair of taxa. However, the size of these collections can be exponential in the size of the taxa set. In this paper, we show that, if we have no “shortcuts”, that is, the networks are normal, the same results are obtained with only a quadratic number of inter-taxa distances by using the shortest distance between each pair of taxa. The proofs are constructive and give cubic-time algorithms in the size of the taxa sets for building such weighted networks.

     

Ancient weapons: the three-dimensional structure of amoebapore A

The recently solved three-dimensional structure of amoebapore A, the major pore-forming protein of Entamoeba histolytica, represents the first tertiary structure determined from a parasitic toxin. The implications derived from this solved structure, together with biochemical data, paint a picture of a unique activation mechanism and reveal that a histidine-mediated dimerization of the protein acts as the molecular switch for the formation of oligomeric pores in target cell membranes.     

The determination of the three-dimensional structure of barley serine proteinase inhibitor 2 by nuclear magnetic resonance, distance geometry and restrained molecular dynamics

The solution structure of the 64 residue structured domain (residues 20-83) of barley serine proteinase inhibitor 2 (BSPI-2) is determined on the basis of 403 interproton distance, 34 phi backbone torsion angle and 26 hydrogen bonding restraints derived from n.m.r. measurements. A total of 11 converged structures were computed using a metric matrix distance geometry algorithm and refined by restrained molecular dynamics. The average rms difference between the final 11 structures and the mean structure obtained by averaging their coordinates is 1.4 +/- 0.2 A for the backbone atoms and 2.1 +/- 0.1 A for all atoms. The overall structure, which is almost identical to that found by X-ray crystallography, is disc shaped and consists of a central four component mixed parallel and antiparallel beta-sheet flanked by a 13 residue alpha-helix on one side and the reactive site loop on the other.     

A molecular ruler for measuring quantitative distance distributions

We report a novel molecular ruler for measurement of distances and distance distributions with accurate external calibration. Using solution X-ray scattering we determine the scattering interference between two gold nanocrystal probes attached site-specifically to a macromolecule of interest. Fourier transformation of the interference pattern provides a model-independent probability distribution for the distances between the probe centers-of-mass. To test the approach, we measure end-to-end distances for a variety of DNA structures. We demonstrate that measurements with independently prepared samples and using different X-ray sources are highly reproducible, we demonstrate the quantitative accuracy of the first and second moments of the distance distributions, and we demonstrate that the technique recovers complex distribution shapes. Distances measured with the solution scattering-interference ruler match the corresponding crystallographic values, but differ from distances measured previously with alternate ruler techniques. The X-ray scattering interference ruler should be a powerful tool for relating crystal structures to solution structures and for studying molecular fluctuations.     

Use of restrained molecular dynamics in water to determine three-dimensional protein structure: prediction of the three-dimensional structure of Ecballium elaterium trypsin inhibitor II

Refinement of distance geometry (DG) structures of EETI-II (Heitz et al.: Biochemistry 28:2392-2398, 1989), a member of the squash family trypsin inhibitor, have been carried out by restrained molecular dynamics (RMD) in water. The resulting models show better side chain apolar/polar surface ratio and estimated solvation free energy than structures refined "in vacuo." The consistent lower values of residual NMR constraint violations, apolar/polar surface ratio, and solvation free energy for one of these refined structures allowed prediction of the 3D folding and disulfide connectivity of EETI-II. Except for the few first residues for which no NMR constraints were available, this computer model fully agreed with X-ray structures of CMTI-I (Bode et al.: FEBS Lett. 242:285-292, 1989) and EETI-II complexed with trypsin that appeared after the RMD simulation was completed. Restrained molecular dynamics in water is thus proved to be highly valuable for refinement of DG structures. Also, the successful use of apolar/polar surface ratio and of solvation free energy reinforce the analysis of Novotny et al. (Proteins 4:19-30, 1988) and shows that these criteria are useful indicators of correct versus misfolded models.     

The solution structure of motilin from NMR distance constraints, distance geometry, molecular dynamics, and an iterative full relaxation matrix refinement

A model of the structure of the 22 amino acid residue gastrointestinal peptide hormone motilin in 30% hexafluoro-2-propanol has been obtained by using distance constraints obtained from two-dimensional nuclear Overhauser enhancements. A set of initial structures have been generated by using the distance geometry program DIANA, and 10 of these structures have been refined by using restrained molecular dynamics (AMBER). The resulting structures are virtually indistinguishable in terms of constraint violations and energies and display less than 0.5-A root mean square deviations (RMSD) of the backbone atom positions from Tyr7 to Lys20. A comparison of back-calculated and experimental NOE intensities indicates that RMSD's are not the best indicators of the goodness of fit or of the precision with which the structure is defined. The structure was further refined by fitting the experimental NOE data using an iterative full relaxation matrix analysis. The mean error between the observed and calculated backbone NOE intensities for the final refined structure was 0.23 for the full length of the molecule, 0.18 for the region from Glu9 to Lys20, and 0.29 for the region from Phe1 to Gly8. R factors for the same regions were 0.27, 0.19, and 0.43, respectively. All of the NOE-determined structures consistently display an alpha-helix which extends from Glu9 to Lys20. Considerable lack of definition of structure exists at the amino and carboxyl ends of the molecule and also in the vicinity of Thr6-Tyr7-Gly8. A tendency to form a wide turn appears to exist over the sequence Pro3-Ile4-Phe5-Thr6, but the structure in this region is not well defined by the NOE data.     

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.     

Simultaneous determination of disulphide bridge topology and three-dimensional structure using ambiguous intersulphur distance restraints: Possibilities and limitations

Knowledge of the native disulphide bridge topology allows the introduction of conformational restraints between remote parts of the peptide chain. This information is therefore of great importance for the successful determination of the three-dimensional structure of cysteine-rich proteins by NMR spectroscopy. In this paper we investigate the limitations of using ambiguous intersulphur restraints [Nilges, M. (1995) J. Mol. Biol., 245, 645–660] associated with NMR experimental information to determine the native disulphide bridge pattern. Using these restraints in a simulated annealing protocol we have determined the correct topology of numerous examples, including a protein with seven disulphide bridges (phospholipase A₂) and a protein in which 25% of the total number of residues are cysteines (-conotoxin GIIIB). We have also characterised the behaviour of the method when only limited experimental data is available, and find that the proposed protocol permits disulphide bridge determination even with a small number of restraints (around 5 NOEs – including a long-range restraint – per residue). In addition, we have shown that under these conditions the use of a reduced penalty function allows the identification of misassigned NOE restraints. These results indicate that the use of ambiguous intersulphur distances with the proposed simulated annealing protocol is a general method for the determination of disulphide bridge topology, particularly interesting in the first steps of NMR study of cysteine-rich proteins. Comparison with previously proposed protocols indicates that the presented method is more reliable and the interpretation of results is straightforward.     

Oligonucleotide arrays: information from replication and spatial structure

MOTIVATION: The introduction of oligonucleotide DNA arrays has resulted in much debate concerning appropriate models for the measurement of gene expression. By contrast, little account has been taken of the possibility of identifying the physical imperfections in the raw data. RESULTS: This paper demonstrates that, with the use of replicates and an awareness of the spatial structure, deficiencies in the data can be identified, the possibility of their correction can be ascertained and correction can be effected (by use of local scaling) where possible. The procedures were motivated by data from replicates of Arabidopsis thaliana using the GeneChip ATH1-121501 microarray. Similar problems are illustrated for GeneChip Human Genome U133 arrays and for the newer and larger GeneChip Wheat Genome microarray. AVAILABILITY: R code is freely available on request.     

The three-dimensional structure of porin from Rhodobacter capsulatus at 3 A resolution

The crystal structure of porin from Rhodobacter capsulatus strain 37b4 has been solved at 3.0 A (1 A = 0.1 nm) resolution by multiple isomorphous replacement and solvent-flattening. The three pores of the trimer are well defined in the electron density map. Each pore consists of a 16-stranded beta-barrel which traverses the membrane as a tube. Near its center the tube is narrowed by chain segments protruding from the inner wall of the barrel that form an eye-let with an irregular cross-section of about 6 A by 10 A. The eye-let has an axial length of about 10 A; it defines the exclusion limit for diffusing particles.