Predicting the Structure of the Flavodoxin from Eschericia coli by Homology Modeling,Distance Geometry and Molecular Dynamics

Abstract

As part of our on-going development of a method, based upon distance geometry calculations, for predicting the structures of proteins from the known structures of their homologues, we have predicted the structure of the 176 residue Flavodoxin from Escherichia coli. This prediction was based upon the crystal structures of the homologous Flavodoxins from Anacystis nidulans, Chondrus crispus, Desulfovibrio vulgaris and Clostridium beijerinckii, whose sequence identities with Escherichia coli were 44%, 33%, 23% and 16%, respectively. A total of 13,043 distance constraints among the alpha-carbons of the Escherichia coli structure were derived from the sequence alignments with the known structures, together with 8,893 distance constraints among backbone and sidechain atoms of adjacent residues, 978 between the alpha-carbons and selected atoms of the flavin mononucleotide cofactor, 116 constraints to enforce conserved hydrogen bonds, and 452 constraints on the torsion angles in conserved residues. An ensemble of ten random Escherichia coli structures was computed from these constraints, with an average root mean square coordinate deviation (RMSD) among the alpha carbons of 0.85 Ångstroms (excluding the first 1 and last 6 residues, which have no corresponding residues in any of the homologues and hence were unconstrained); the corresponding average heavy-atom RMSD was 1.60 Å.

Since the distance geometry calculations were performed without hydrogen atoms, protons were added to the resulting structures and these structures embedded in a 50 × 50 × 40 Å solvent box with periodic boundary conditions. They were then subjected to a 20 picosecond dynamical simulated annealing procedure, starting at 300 K and gradually reduced to 10K, in which all the distance and torsion angle constraints were maintained by means of harmonic restraint functions. This was followed up by 1000 iterations of unrestrained conjugate gradients minimization. The goal of this energy refinement procedure was not to drastically modify the structures in an attempt at a priori prediction, but merely to improve upon the predictions obtained from the geometric constraints, particularly with regard to their local backbone and sidechain conformations and their hydrogen bonds. The resulting structures differed from the respective starting structures by an average of 1.52 Å in their heavy atom RMSD's, while the average RMSD among the heavy atoms of residues 2-170 increased slightly to 1.66 Å. We hope these structures will be good enough to enable the phase problem to be solved for the crystallographic data that is now being collected on this protein.     

Ab initio structure prediction of the antibody hypervariable H3 loop

Antibodies have the capability of binding a wide range of antigens due to the diversity of the six loops constituting the complementarity determining region (CDR). Among the six loops, the H3 loop is the most diverse in structure, length, and sequence identity. Prediction of the three‐dimensional structures of antibodies, especially the CDR loops, is an important step in the computational design and engineering of novel antibodies for improved affinity and specificity. Although it has been demonstrated that the conformation of the five non‐H3 loops can be accurately predicted by comparing their sequences against databases of canonical loop conformations, no such connection has been established for H3 loops. In this work, we present the results for ab initio structure prediction of the H3 loop using conformational sampling and energy calculations with the program Prime on a dataset of 53 loops ranging in length from 4 to 22 residues. When the prediction is performed in the crystal environment and including symmetry mates, the median backbone root mean square deviation (RMSD) is 0.5 Å to the crystal structure, with 91% of cases having an RMSD of less than 2.0 Å. When the prediction is performed in a noncrystallographic environment, where the scaffold is constructed by swapping the H3 loops between homologous antibodies, 70% of cases have an RMSD below 2.0 Å. These results show promise for ab initio loop predictions applied to modeling of antibodies. © 2012 Wiley Periodicals, Inc.     

The Genetic Algorithm and the Conformational Search of Polypeptides and Proteins

Abstract

The genetic algorithm is a technique of function optimization derived from the principles of evolutionary theory. We have adapted it to perform conformational search on polypeptides and proteins. The algorithm was first tested on several small polypeptides and the 46 amino acid protein crambin under the AMBER potential energy function. The probable global minimum conformations of the polypeptides were located 90% of the time and a non-native conformation of crambin was located that was 150kcal/mol lower in potential energy than the minimized crystal structure conformation. Next, we used a knowledge-based potential function to predict the structures of melittin, pancreatic polypeptide, and crambin. A 2.31 Å ΔRMS conformation of melittin and a 5.33 Å ΔRMS conformation of pancreatic polypeptide were located by genetic algorithm-based conformational search under the knowledge-based potential function. Although the ΔRMS of pancreatic polypeptide was somewhat high, most of the secondary structure was correct. The secondary structure of crambin was predicted correctly, but the potential failed to promote packing interactions. Finally, we tested the packing aspects of our potential function by attempting to predict the tertiary structure of cytochrome b ₅₆₂ given correct secondary structure as a constraint. The final predicted conformation of cytochrome b ₅₆₂ was an almost completely extended continuous helix which indicated that the knowledge-based potential was useless for tertiary structure prediction. This work serves as a warning against testing potential functions designed for tertiary structure prediction on small proteins.     

The reconstruction of side-chain conformations from protein Cα coordinates

A model of nine proteins including side-chain atoms have been built from the known C^α coordinates and amino acid sequences using a Monte Carlo Protein Building Annealing method. The Cartesian coordinates for the side-chain atoms were established with bond lengths and angles selected randomly from within previously determined ranges. A simulated annealing technique is used to generate some 300 structures with differing side-chain conformations. The atomic coordinates of the backbone atoms are fixed during the simulated annealing process. The coordinates of the side-chain atoms of 300 low energy conformations are averaged to obtain a mean structure that is minimized with the C^α atoms constrained to their position in the x-ray structure using the OPLS/AMBER force field with the GB/SA water model. The rms deviation of the main-chain atoms (without C^β) compared with the corresponding crystal structures is in the range 0.20–0.64 Å. The rms deviation of the side-chain atoms is between 1.72 and 2.71 Å and for all atoms is between 1.19 and 1.99 Å. The method is insensitive to random errors in the C^α positions and the computational requirement is modest. © 1997 John Wiley & Sons, Inc.     

Toward high‐resolution homology modeling of antibody Fv regions and application to antibody–antigen docking

High‐resolution homology models are useful in structure‐based protein engineering applications, especially when a crystallographic structure is unavailable. Here, we report the development and implementation of RosettaAntibody, a protocol for homology modeling of antibody variable regions. The protocol combines comparative modeling of canonical complementarity determining region (CDR) loop conformations and de novo loop modeling of CDR H3 conformation with simultaneous optimization of V_L‐V_H rigid‐body orientation and CDR backbone and side‐chain conformations. The protocol was tested on a benchmark of 54 antibody crystal structures. The median root mean square deviation (rmsd) of the antigen binding pocket comprised of all the CDR residues was 1.5 Å with 80% of the targets having an rmsd lower than 2.0 Å. The median backbone heavy atom global rmsd of the CDR H3 loop prediction was 1.6, 1.9, 2.4, 3.1, and 6.0 Å for very short (4–6 residues), short (7–9), medium (10–11), long (12–14) and very long (17–22) loops, respectively. When the set of ten top‐scoring antibody homology models are used in local ensemble docking to antigen, a moderate‐to‐high accuracy docking prediction was achieved in seven of fifteen targets. This success in computational docking with high‐resolution homology models is encouraging, but challenges still remain in modeling antibody structures for sequences with long H3 loops. This first large‐scale antibody–antigen docking study using homology models reveals the level of “functional accuracy” of these structural models toward protein engineering applications. Proteins 2009; 74:497–514. © 2008 Wiley‐Liss, Inc.     

Conformational variability of Gly-Gly segments in peptides: A comparison of the crystal structures of an acyclic pentapeptide and an octapeptide

The crystal structure of an acyclic pentapeptide, Boc-Gly-Gly-Leu-Aib-Val-OMe, reveals an extended conformation for the Gly-Gly segment, in contrast to the helical conformation determined earlier in the octapeptide Boc-Leu-Aib-Val-Gly-Gly-Leu-Aib-Val-OMe [I. L. Karle, A. Banerjee, S. Bhattacharjya, and P. Balaram [1996] Biopolymers, Vol. 38, pp. 515–526). The pentapeptide crystallizes in space group P2₁ with one molecule in the asymmetric unit. The cell parameters are: a = 10.979(2) Å, b = 9.625(2) Å, c = 14.141(2) Å, and β = 96.93(1)°, R = 6.7% for 2501 reflections (I > 3σ(I)). The Gly-Gly segment is extended (ϕ₁ = −92°, ψ₁ = −133°, ϕ₂ = 140°, ψ₂ = 170°), while the Leu-Aib segment adopts a type II β-turn conformation (ϕ₃ = −61°, ψ₃ = 130°, ϕ₄ = 71°, ψ₄ = 6°). The observed conformation for the pentapeptide permits rationalization of a structural transition observed for the octapeptide in solution. An analysis of Gly-Gly segments in peptide crystal structures shows a preference for either β-turn or extended conformations. © 1997 John Wiley & Sons, Inc.     

Prediction of the folding of short polypeptide segments by uniform conformational sampling

A procedure, CONGEN, for uniformly sampling the conformational spaceof short polypeptide segments in proteins has been implemented. Because thetime required for this sampling grows exponentially with the number of residues, parameters are introduced to limit the conformational space that has to be explored. This is done by the use of the empirical energy function ofCHARMM [B. R. Brooks, R. E. Bruccoleri, B. D. Olafson, D. J. States, S. Swaminathan and M. Karplus (1983) J. Comput. Chem. 4 , 187-217] and truncating the search when conformations of grossly unfavorable energy are sampled. Tests are made to determine control parameters that optimize the search without excluding important configurations. When applied to known protein structures, the resulting procedure is generally capable of generating conformations where the lowest energy conformation matches the known structure within a rms deviation of 1 Å.     

Learning about protein folding via potential functions

Over the last few years we have developed an empirical potential function that solves the protein structure recognition problem: given the sequence for an n-residue globular protein and a collection of plausible protein conformations, including the native conformation for that sequence, identify the correct, native conformation. Having determined this potential on the basis of only some 6500 native/nonnative pairs of structures for 58 proteins, we find it recognizes the native conformation for essentially all compact, soluble, globular proteins having known native conformations in comparisons with 10⁴ to 10⁶ reasonable alternative conformations apiece. In this sense, the potential encodes nearly all the essential features of globular protein conformational preference. In addition it “knows” about many additional factors in protein folding, such as the stabilization of multimeric proteins, quaternary structure, the role of disulfide bridges and ligands, proproteins vs. processed proteins, and minimal strand lengths in globular proteins. Comparisons are made with other sorts of protein folding problems, and applications in protein conformational determination and prediction are discussed. © 1994 Wiley-Liss, Inc.     

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

Abstract

An implementation of the variable-target-function procedure, first introduced by Braun and Gō [W. Braun and N. Gō, 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 Å for the α-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 Å for the α-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.     

Crystal structure,conformational analysis,and molecular dynamics of tetra-O-methyl-(+) -catechin

The structure of tetra-O-methyl- (+) -catechin has been determined in the crystalline state. Two independent molecules, denoted structure A and structure B, exist in the unit cell. Crystals are triclinic, space group P1, a = 4.8125(2) Å, b = 12.9148(8) Å, c = 13.8862(11) Å, α = 86.962(6) °, β = 89.120(5)°, γ = 88.044(5)°, Z = 2, D_c = 1.336 g cm^?3, R = 0.033 for 6830 observations. The heterocyclic rings of the crystal structures are compared to previous results for 8-bromotetra-O-methyl-(+)-catechin, penta-O-acetyl-(+)-catechin, and (?) -epicatechin. One of the two molecules has a heterocyclic ring conformation similar to that observed previously for (?)-epicatechin, and the other has a heterocyclic ring conformation similar to one predicted earlier in a theoretical analysis of dimers of (+)-catechin and (?) -epicatechin. Both structure A and structure B in the crystal have heterocyclic ring conformations that place the dimethoxyphenyl substituent at C(2) in the equatorial position. However, this heterocyclic ring conformation does not explain the proton nmr coupling constant measured in solution. Molecular dynamics simulations show an equatorial ? axial interconversion of the heterocyclic ring, which can explain the nmr results. © 1993 John Wiley & Sons, Inc.     

Protein conformational landscapes: Energy minimization and clustering of a long molecular dynamics trajectory

Using energy minimization and cluster analysis, we have analyzed a 1020 ps molecular dynamics trajectory of solvated bovine pancreatic trypsin inhibitor. Elucidation of conformational sub states in this way both illustrates the degree of conformational convergence in the simulation and reduces the structural data to a tractable subset. The relative movement of structures upon energy minimization was used to estimate the sizes of features on the protein potential energy surface. The structures were analyzed using their pairwise root-mean-square C_α deviations, which gave a global measure of conformational changes that would not be apparent by monitoring single degrees of freedom. At time scales of 0.1 ps, energy minimization detected sharp transitions between energy minima separated by 0.1 Å rms deviation. Larger conformational clusters containing these smaller minima and separated by 0.25 Å were seen at 1 ps time scales. Both of these small features of the conformational landscape were characterized by movements in loop regions associated with small, correlated backbone dihedral angle shifts. On a nanosecond time scale, the main features of the protein energy landscape were clusters separated by over 0.7 Å rms deviation, with only seven of these sub states visited over the 1 ns trajectory. These substates, discernible both before and after energy minimization, differ mainly in a monotonic pivot of the loop residues 11–18 over the course of the simulation. This loop contains lysine 17, which specifically binds to trypsin in the active site. The trajectory did not return to previously visited clusters, indicating that this trajectory has not been shown to have completely sampled the conformational substates available to it. Because the apparent convergence to a single region of conformation space depends on both the time scale of observation and the size of the conformational features examined, convergence must be operationally defined within the context of the simulation. © 1995 Wiley-Liss, Inc.     

Ab initio structure prediction for small polypeptides and protein fragments using genetic algorithms

Ab initio folding simulations have been performed on three peptides, using a genetic algorithm-based search method which operates on a full atom representation. Conformations are evaluated with an empirical force field parameterized by a potential of mean force analysis of experimental structures. The dominant terms in the force field are local and nonlocal main chain electrostatics and the hydrophobic effect. Two of the simulated structures were for fragments of complete proteins (eosinophil-derived neurotoxin (EDN) and the subtilisin propeptide) that were identified as being likely initiation sites for folding. The experimental structure of one of these (EDN) was subsequently found to be consistent with that prediction (using local hydrophobic burial as the determinant for independent folding). The simulations of the structures of these two peptides were only partly successful. The most successful folding simulation was that of a 22-residue peptide corresponding to the membrane binding domain of blood coagulation factor VIII (Membind). Three simulations were performed on this peptide and the lowest energy conformation was found to be the most similar to the experimental structure. The conformation of this peptide was determined with a C_α rms deviation of 4.4 Å. Although these simulations were partly successful there are still many unresolved problems, which we expect to be able to address in the next structure prediction experiment. © 1995 Wiley-Liss, Inc.     

Monte Carlo loop refinement and virtual screening of the thyroid-stimulating hormone receptor transmembrane domain

Metropolis Monte Carlo (MMC) loop refinement has been performed on the three extracellular loops (ECLs) of rhodopsin and opsin-based homology models of the thyroid-stimulating hormone receptor transmembrane domain, a class A type G protein-coupled receptor. The Monte Carlo sampling technique, employing torsion angles of amino acid side chains and local moves for the six consecutive backbone torsion angles, has previously reproduced the conformation of several loops with known crystal structures with accuracy consistently less than 2?Å. A grid-based potential map, which includes van der Waals, electrostatics, hydrophobic as well as hydrogen-bond potentials for bulk protein environment and the solvation effect, has been used to significantly reduce the computational cost of energy evaluation. A modified sigmoidal distance-dependent dielectric function has been implemented in conjunction with the desolvation and hydrogen-bonding terms. A long high-temperature simulation with 2?kcal/mol repulsion potential resulted in extensive sampling of the conformational space. The slow annealing leading to the low-energy structures predicted secondary structure by the MMC technique. Molecular docking with the reported agonist reproduced the binding site within 1.5?Å. Virtual screening performed on the three lowest structures showed that the ligand-binding mode in the inter-helical region is dependent on the ECL conformations.     

Analysis of six protein structures predicted by comparative modeling techniques

The protein structures of six comparative modeling targets were predicted in a procedure that relied on improved energy minimization, without empirical rules, to position all new atoms. The structures of human nucleoside diphosphate kinase NM23-H2, HPr from Mycoplasma capricolum, 2Fe-2S ferredoxin from Haloarcula marismortui, eosinophil-derived neurotoxin (EDN), mouse cellular retinoic acid protein I (CRABP1), and P450eryf were predicted with root mean square deviations on Cα atoms of 0.69, 0.73, 1.11, 1.48, 1.69, and 1.73 Å, respectively, compared to the target crystal structures. These differences increased as the sequence similarity between the target and parent proteins decreased from about 60 to 20% identity. More residues were predicted than form the common region shared by the two crystal structures. In most cases insertions or deletions between the target and the related protein of known structure were not correctly positioned. One two residue insertion in CRABP1 was predicted in the correct conformation, while a nine residue insertion in EDN was predicted in the correct spatial region, although not in the correct conformation. The positions of common cofactors and their binding sites were predicted correctly, even when overall sequence similarity was low. © 1995 Wiley-Liss, Inc.     

Recognition of native structure from complete enumeration of low-resolution models with constraints

Complete sets of low-resolution conformations are generated for eight small proteins by rotating the C_α-C_α virtual bonds at selected flexible regions, while the remaining structural elements are assumed to move in rigid blocks. Several filtering criteria are used to reduce the ensemble size and to ensure the sampling of well-constructed conformations. These filters, based on structure and energy constraints deduced from knowledge-based studies, include the excluded volume requirement, the radius of gyration constraint, and the occurrence of sufficiently strong attractive inter-residue potentials to stabilize compact forms. About 8,000 well-constructed decoys or “probable folds” (PFs) are constructed for each protein. A correlation between root-mean-square (rms) deviations from X-ray structure and total energies is observed, revealing a decrease in energy as the rms deviation decreases. The conformation with the lowest energy exhibits an rms deviation smaller than 3.0 Å, in most of the proteins considered. The results are highly sensitive to the choice of flexible regions. A strong tendency to assume native state rotational angles is revealed for some flexible bonds from the analysis of the distributions of dihedral angles in the PFs, suggesting the formation of foldons near these locally stable regions at early folding pathway. Proteins 32:211–222, 1998. © 1998 Wiley-Liss, Inc.     

Conformational Preferences of Heterochiral Peptides. Crystal Structures of Heterochiral Peptides Boc-(D) Val-(D) Ala-Leu-Ala-OMe and Boc-Val-Ala-Leu-(D) Ala-OMe-Enhanced Stability of β-sheet Through C-H…O Hydrogen Bonds

Abstract

The crystal structures of Boc-(D) Val-(D) Ala-Leu-Ala-OMe (vaLA) and Boc-Val-Ala-Leu-(D) Ala-OMe (VALa) have been determined. vaLA crystallises in space group P2₁2₁2₁ with a = 9.401 (4), b = 17.253 (5), c = 36.276 (9)Å, V = 5884 (3) ų, Z = 8, R = 0.086. VALa crystallises in space group P2₁ with a = 9.683 (9), b = 17.355 (7), c = 18.187 (9) Å, β = 95.84 (8)°, V = 3040(4) ų, Z = 4, R = 0.125. There are two molecules in the asymmetric unit in antiparallel β-sheet arrangement in both the structures. Several of the C^α hydrogens are in hydrogen bonding contact with the carbonyl oxygen in the adjacent strand.

An analysis of the observed conformational feature of D-chiral amino acid residues in oligopeptides, using coordinates of 123 crystal structures selected from the 1998 release of CSD has been carried out. This shows that all the residues except D-isoleucine prefer both extended and α_L conformation though the frequence of occurence may not be equal. In addition to this, D-leucine, valine, proline and phenylalanine have assumed α_R conformations in solid state. D-leucine has a strong preference for helical conformation in linear peptides whereas they prefer an extended conformation in cyclic peptides.     

Antibodies as a model system for comparative model refinement

Predicting the conformations of loops is a critical aspect of protein comparative (homology) modeling. Despite considerable advances in developing loop prediction algorithms, refining loops in homology models remains challenging. In this work, we use antibodies as a model system to investigate strategies for more robustly predicting loop conformations when the protein model contains errors in the conformations of side chains and protein backbone surrounding the loop in question. Specifically, our test system consists of partial models of antibodies in which the “scaffold” (i.e., the portion other than the complementarity determining region, CDR, loops) retains native backbone conformation, whereas the CDR loops are predicted using a combination of knowledge‐based modeling (H1, H2, L1, L2, and L3) and ab initio loop prediction (H3). H3 is the most variable of the CDRs. Using a previously published method, a test set of 10 shorter H3 loops (5–7 residues) are predicted to an average backbone (N? Cα? C? O) RMSD of 2.7 Å while 11 longer loops (8–9 residues) are predicted to 5.1 Å, thus recapitulating the difficulties in refining loops in models. By contrast, in control calculations predicting the same loops in crystal structures, the same method reconstructs the loops to an average of 0.5 and 1.4 Å for the shorter and longer loops, respectively. We modify the loop prediction method to improve the ability to sample near‐native loop conformations in the models, primarily by reducing the sensitivity of the sampling to the loop surroundings, and allowing the other CDR loops to optimize with the H3 loop. The new method improves the average accuracy significantly to 1.3 Å RMSD and 3.1 Å RMSD for the shorter and longer loops, respectively. Finally, we present results predicting 8–10 residue loops within complete comparative models of five nonantibody proteins. While anecdotal, these mixed, full‐model results suggest our approach is a promising step toward more accurately predicting loops in homology models. Furthermore, while significant challenges remain, our method is a potentially useful tool for predicting antibody structures based on a known Fv scaffold. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Molecular Conformation of 2′-Deoxy-2′-methylidene-cytidine: A Potent Antineoplastic Nucleoside

Abstract

2′-Deoxy-2′-methylidenecytidine (DMDC), a potent inhibitor of the growth of tumor cells, was crystallized with two different forms. One is dihydrated (DMDC·2H₂O) and the other is its hydrochloride salt (DMDC·HCLl). Both crystal and molecular structures have been determined by the X-ray diffraction method. In both forms the glycosidic and sugar conformations are anti and C(4′)-exo, respectively, whereas the conformation about the exocyclic bond is trans for DMDC·2H₂O and gauche ⁺ for DMDC·HCl. Proton nuclear magnetic resonance data of DMDC indicate a preference for the anti C(4′)-exo conformation found in the solid state. These molecular conformations were compared with the related pyrimidine nucleosides. When the cytosine bases are brought into coincidence, DMDC displays the exocyclic C(4′)-C(5′) bond located on the very close position to those of pyrimidine nucleosides with typical overall conformations. On the other hand, the hydroxyl O(3′)-H groups are separated by ca. 3 Å in the cases of DMDC and other pyrimidine nucleosides which have the C(2′)-endo sugar conformation. This result may be useful for the implication about the mechanism of the biological activity of DMDC.     

Variants of human thymidylate synthase with loop 181–197 stabilized in the inactive conformation

Loop 181–197 of human thymidylate synthase (hTS) populates two major conformations, essentially corresponding to the loop flipped by 180°. In one of the conformations, the catalytic Cys195 residue lies distant from the active site making the enzyme inactive. Ligands stabilizing this inactive conformation may function as allosteric inhibitors. To facilitate the search for such inhibitors, we have expressed and characterized several mutants designed to shift the equilibrium toward the inactive conformer. In most cases, the catalytic efficiency of the mutants was only somewhat impaired with values of k_cat/K_m reduced by factors in a 2–12 range. One of the mutants, M190K, is however unique in having the value of k_cat/K_m smaller by a factor of ～7500 than the wild type. The crystal structure of this mutant is similar to that of the wt hTS with loop 181–197 in the inactive conformation. However, the direct vicinity of the mutation, residues 188–194 of this loop, assumes a different conformation with the positions of C_α shifted up to 7.2 Å. This affects region 116–128, which became ordered in M190K while it is disordered in wt. The conformation of 116–128 is however different than that observed in hTS in the active conformation. The side chain of Lys190 does not form contacts and is in solvent region. The very low activity of M190K as compared to another mutant with a charged residue in this position, M190E, suggests that the protein is trapped in an inactive state that does not equilibrate easily with the active conformer.     

Crystal and molecular structure of benzyloxycarbonyl-α-aminoisobutyryl-L-prolyl methylamide: The observation of an X2-Pro3 Type III β-Turn

The crystal and molecular structure of N-benzyloxycarbonyl-α-aminoisobutyryl-L -prolyl methylamide, the amino terminal dipeptide fragment of alamethicin, has been determined using direct methods. The compound crystallizes in the orthorhombic system with the space group P2₁2₁2₁. Cell dimensions are a = 7.705 Å, b = 11.365 Å, and c = 21.904 Å. The structure has been refined using conventional procedures to a final R factor of 0.054. The molecular structure possesses a 4 → 1 intramolecular N-H—O hydrogen bond formed between the CO group of the urethane moiety and the NH group of the methylamide function. The peptide backbone adopts the type III β-turn conformation, with ?₂ = ?51.0°, ψ₂ = ?39.7°, ?₃ = ?65.0°, ψ₃ = ?25.4°. An unusual feature is the occurrence of the proline residue at position 3 of the β-turn. The observed structure supports the view that Aib residues initiate the formation of type III β-turn conformations. The pyrrolidine ring is puckered in C^γ-exo fashion.