Structure of porcine heart cytoplasmic malate dehydrogenase: combining X-ray diffraction and chemical sequence data in structural studies

The amino acid sequence of cytoplasmic malate dehydrogenase (sMDH) has been determined by a combination of X-ray crystallographic and chemical sequencing methods. The initial molecular model incorporated an "X-ray amino acid sequence" that was derived primarily from an evaluation of a multiple isomorphous replacement phased electron density map calculated at 2.5-A resolution. Following restrained least-squares crystallographic refinement, difference electron density maps were calculated from model phases, and attempts were made to upgrade the X-ray amino acid sequence. The method used to find the positions of peptides in the X-ray structure was similar to those used for studying protein homology and was shown to be successful for large fragments. For sMDH, X-ray methods by themselves were insufficient to derive a complete amino acid sequence, even with partial chemical sequence data. However, for this relatively large molecule at medium resolution, the electron density maps were of considerable help in determining the linear position of peptide fragments. The N-acetylated polypeptide chain of sMDH has 331 amino acids and has been crystallographically refined to an R factor of 19% for 2.5-A resolution diffraction data.     

Crystallographic structure refinement of Chromatium high potential iron protein at two Angstroms resolution.

The structure of Chromatium high potential iron protein (HiPIP) has been refined by semiautomatic Fo-Fc (observed minus calculated structure amplitude Fourier methods to a convential R index, R=sum of the absolute value of Fo-Fc divided by the sum of Fo, of 24.7% for a model in which bond distances and angles are constrained to standard values. Bond length and angle constraints were applied only intermittenly during the computations. At a late stage of the refinement, atomic parameters for only the Fe4S4 cluster plus the 4 associated cystein S-gamma atoms were adjusted by least squares methods and kept fixed during the rest of the refinement. The refined model consists of 625 of the 632 nonhydrogen atoms in the protein plus 75 water molecules. Seven side chain atoms could not be located in the final electron density map. A computer program rather than visual inspection was used wherever possible in the refinement: for locating water molecules, for removing water molecules that too closely approach other atoms, for deleting atoms that lay in regions of low electron density, and for evaluating the progress of refinement. Fo-Fc Fourier refinement is sufficiently economical to be applied routinely in protein crystal structure determinations. The complete HiPIP refinement required approximately 12 hours of CDC 3600 computer time and cost less than $3000 starting from a "trial structure," based upon multipe isomorphoous replacement phases, which gave an R of 43%...     

X-ray crystal structure of the bovine alpha-chymotrypsin/eglin c complex at 2.6 A resolution

The crystal structure of the molecular complex formed by bovine alpha-chymotrypsin and the recombinant serine proteinase inhibitor eglin c from Hirudo medicinalis has been solved using monoclinic crystals of the complex, reported previously. Four circle diffractometer data at 3.0 A resolution were employed to determine the structure by molecular replacement techniques. Bovine alpha-chymotrypsin alone was used as the search model; it allowed us to correctly orient and translate the enzyme in the unit cell and to obtain sufficient electron density for positioning the eglin c molecule. After independent rigid body refinement of the two complex components, the molecular model yielded a crystallographic R factor of 0.39. Five iterative cycles of restrained crystallographic refinement and model building were conducted, gradually increasing resolution. The current R factor at 2.6 A resolution (diffractometer data) is 0.18. The model includes 56 solvent molecules. Eglin c binds to bovine alpha-chymotrypsin in a manner consistent with other known serine proteinase/inhibitor complex structures. The reactive site loop shows the expected conformation for productive binding and is in tight contact with bovine alpha-chymotrypsin between subsites P3 and P'2; Leu 451 acts as the P1 residue, located in the primary specificity S1 site of the enzyme. Hydrogen bonds equivalent to those observed in complexes of trypsin(ogen) with the pancreatic basic- and secretory-inhibitors are found around the scissile peptide bond.     

Crystal structure of the trigonal form of bovine beta-lactoglobulin and of its complex with retinol at 2.5 A resolution

The structure of the trigonal crystal form of bovine beta-lactoglobulin has been determined by X-ray diffraction methods. An electron density map, calculated with phases obtained by the multiple isomorphous replacement method, served as a starting point for alternate cycles of model building and restrained least-squares refinement. The model of the molecule fitted to the initial Fourier map was the one built for the orthorhombic crystal form of beta-lactoglobulin, solved at 2.8 A resolution (1 A = 0.1 nm). The final R factor for 1456 atoms (1276 non-hydrogen protein atoms and 180 solvent atoms) is 0.22, including 5245 reflections from 6.0 to 2.5 A. The molecule shows significant differences in the two crystal forms mentioned, mainly due to different packing. In the trigonal form, the species crystallized does not appear to be dimeric, but a linear polymer with tight intermolecular contacts. A difference electron density map between the complex of beta-lactoglobulin with retinol and the native protein shows no significant peaks in the cavity which, in the similar retinol-binding protein, binds the chromophore. Instead, differences are found at a surface pocket, which is limited almost completely by hydrophobic residues.     

Molecular dynamics study of a complex between the human histocompatibility antigen HLA-A2 and the IMP58-66 nonapeptide from influenza virus matrix protein.

The structure of the influenza-virus-matrix-protein (IMP) 58-66 nonapeptide, bound to the major-histocompatibility-complex-encoded human leukocyte antigen (HLA) A2 protein was studied by molecular dynamics simulation. Starting from the extra electron density map of peptides co-crystallized with HLA-A2, the nonapeptide IMP58-66 was docked residue by residue in the protein binding cleft. The complex was simulated for 100 ps in a shell of 1372 water molecules. The averaged simulated HLA-A2 conformation was found to be similar to the crystal structure (0.182 nm RMS deviation, for the backbone atoms of the alpha 1-alpha 2 domain). Nine out of the 14 hydrogen bonds observed in the antigen-binding site were reproduced in the simulation. The IMP58-66 peptide exhibits an extended conformation with kinks at positions 3 and 5. The side chains of residues 2, 3 and 9 develop van der Waals' interactions with hydrophobic pockets of HLA-A2, corresponding to polymorphic residues of the major-histocompatibility-complex-encoded proteins. Both the N-terminus and C-terminus of the nonapeptide were anchored in the antigen-binding groove by hydrogen bonds with conserved amino acids. The N-terminus was more flexible and contacts four HLA-A2 conserved tyrosines (Tyr7, Tyr59, Tyr159 and Tyr171) and Glu63 by direct or water-relayed hydrogen bonds. Water intercalation occurred only around the N-terminus of the peptide, the C-terminal carboxylate forming strong hydrogen bonds with polar residues (Tyr84 and Thr143) and a salt bridge with Lys146 all over the molecular dynamics simulation. This model is fully compatible with the recently published crystal structure of the HLA-B27 protein, complexed by a mixture of self nonapeptides.     

Structure-activity relationships of methionyl-tRNA synthetase: graphics modelling and genetic engineering

The 3D structure of the methionyl-tRNA synthetase from E. coli has been investigated using X-ray analysis^1,2 at a resolution 1.8Å. 90% of the molecule is now well defined and the zinc atom has been identified in a buried region of the molecule, close to the active site. At the same time, the refinement of the complex ATP-MetRS at 2.5 Å has been carried out. The crystallographic R factor has been assigned a value of 25% at 2.5Å with an overall temperature factor of 9Ų and 22% when the individual temperature factors are refined. A Fourier difference map clearly reveals the electron density of the bound ATP, showing the phosphate groups deeply plunging into the active site. In parallel, the synthetase gene has been used to probe some of the enzyme structure-activity relationships. A series of 60 modified enzymes truncated at the C-terminus have been constructed in vitro and assayed for activity. In agreement with the graphics model, the results show that a minimum of 534 residues is necessary to sustain the aminoacylation reaction. A programme of site-directed mutagenesis is in progress: residues thought to be important for the catalytic activity, the metal coordination and tRNA interaction are being modified. Preliminary results are discussed in the light of the crystallographic model.     

Modeling conformational changes in cyclosporin A.

NMR and X-ray structures for the immunosuppressant cyclosporin A (CsA) reveal a remarkable difference between the unbound (free) conformation in organic solvents and the conformation bound to cyclophilin. We have performed computer simulations of the molecular dynamics of CsA under a variety of conditions and confirmed the stability of these two conformations at room temperature in water and in vacuum. However, when the free conformation was modeled in vacuum at 600 K, a transition pathway leading to the bound conformation was observed. This involved a change in the cis MeLeu-9 peptide bond to a trans conformation and the movement of the side chains forming the dominant hydrophobic cluster (residues MeBmt-1, MeLeu-4, MeLeu-6, and MeLeu-10) to the opposite side of the plane formed by the backbone atoms in the molecular ring. The final conformation had a backbone RMS deviation from the bound conformation of 0.53 A and was as stable in dynamics simulations as the bound conformation. Our calculations allowed us to make a detailed analysis of a transition pathway between the free and the bound conformations of CsA and to identify two distinct regions of coordinated movement in CsA, both of which underwent transitions independently.     

Membrane partitioning of peptide aggregates: coarse-grained molecular dynamics simulations

Abstract

Coarse-grained molecular dynamics (CGMD) simulation technique (MARTINI force field) is applied to monitor the aggregation of helical peptides representing the transmembrane sequence and its extension of bone marrow stromal cell antigen 2 (BST-2). One of the peptides is coupled with a protein transducing domain (PTD) of nine arginine residues (R9) at its N-terminal side as well as a peptide, pep11**, which has been shown to bind to human papilloma virus 16 (HPV16) E6 oncoprotein. A short hydrophobic stretch of the transmembrane domain (TMD) of BST-2 aggregates the fastest and inserts into a lipid membrane. An aggregate of R9-pep11** attaches to the membrane via simultaneous contact of many arginine residues. Monomers from the aggregates of the shortest of the hydrophobic TMDs dissolve into the opposing leaflet when the aggregate spans the bilayer. A ‘flipping’ of the individual monomeric peptides is not observed.

Communicated by Ramaswamy H. Sarma     

Normal mode refinement: crystallographic refinement of protein dynamic structure. II. Application to human lysozyme.

The dynamic structure of a protein, human lysozyme, is determined by the normal mode refinement of X-ray crystal structure. This method uses the normal modes of both internal and external motions to distinguish the real internal dynamics from the external terms such as lattice disorder, and gives an anisotropic and concerted picture of atomic fluctuations. The refinement is carried out with diffraction data of 5.0 to 1.8 A resolution, which are collected on an imaging plate. The results of the refinement show: (1) Debye-Waller factor consists of two parts, highly anisotropic internal fluctuations and almost isotropic external terms. The former is smaller than the latter by a factor of 0.72 in the scale of B-factor. Therefore, the internal dynamics cannot be recognized directly from the apparent electron density distribution. (2) The internal fluctuations show basically similar features as those predicted by the normal mode analysis, with almost the same amplitude and a similar level of anisotropy. (3) Correlations of fluctuations are detected between two lobes forming the active site cleft, which move simultaneously in opposite directions. This corresponds to the hinge-bending motion of lysozyme.     

Structure of wheat serine carboxypeptidase II at 3.5-A resolution. A new class of serine proteinase

The structure of serine carboxypeptidase II from wheat bran has been determined to 3.5-A resolution by multiple isomorphous replacement, solvent flattening, and crystallographic refinement. The amino acid sequence has been fit to the electron density map and the model refined to a conventional crystallographic R factor of 20.9%. The molecule is an alpha + beta protein and contains a "catalytic triad" (Asp338, His397, and Ser146) similar in arrangement to those in chymotrypsin and subtilisin. The -fold of the polypeptide backbone is, however, completely different from those enzymes. This suggests that this is a third example of convergent evolution to a common enzymatic mechanism. The -fold is, on the other hand, surprisingly similar to that of the zinc proteinase carboxypeptidase A.     

Structure prediction and R115866 binding study of human CYP26A1: homology modelling,fold recognition,molecular docking and MD simulations

Two three-dimensional (3D) models of human cytochrome P450 26A1 (CYP26A1) were constructed using the programs Modeller and Sybyl-GeneFold, respectively. After refinement by molecular mechanics and molecular dynamics (MD) simulations, the two models were validated by structure analysis-validation online server. Subsequently, a flexible docking study was performed on the model constructed by GeneFold with the potent and specific inhibitor R115866 to examine the enzyme–inhibitor interactions. From the docking results, we can see R115866 interacts with amino acid residues at the active site by multiple hydrophobic interactions including the side chains of His111, Trp112, Ser115, Val116, Leu125, Ser126, Leu221, Phe222, Glu296, Phe299, Gly300, Glu303, Thr304, Pro371 and the cofactor heme. Trp112 and Thr304 form hydrogen bonds with R115866 and play important roles in stabilising the complex. This constructed CYP26A1 model may provide an opportunity to understand the action mode of the enzyme and could be useful in designing novel retinoic acid metabolism blocking agents (RAMBAs).     

Modeling and simulation of the human delta opioid receptor

A model for the human delta opioid receptor has been generated via sequence alignment, structure building using the crystal structure of bovine rhodopsin as a template, and refinement by molecular dynamics simulation. The model building suggested that, in addition to the previously postulated interaction between D128 and Y308, an internal salt bridge also exists between residues D128 and R192, both of which are conserved in all the opioid receptors. The model and salt bridge were then shown to be stable during a 20-nsec simulation in a lipid bilayer. It is therefore proposed that both of these interactions play a role in stabilizing the inactive state of the receptor. The model is also used in an effort to rationalize many of the mutational studies performed on delta opioid receptors, and to suggest a plausible explanation for the differences between known delta opioid agonists and antagonists.     

Crystallographic analysis of AcrB

A His-tagged derivative of the multidrug efflux pump AcrB could be crystallized in three different space groups (R3, R32 and P321). Experimental MAD-phasing maps from R32 AcrB(His) crystals were obtained to a resolution of 3.5 A. Datasets of native and substrate soaked AcrB(His) crystals were collected at the Swiss Light Source X06SA beamline up to a resolution of 2.7 A and refinement of these data provided good quality electron density maps, which allowed us to complement the published AcrB structure (PDB code 1iwg). Introduction of amino acids 860-865 and 868 lacking in the 1iwg structure and deletion of a highly disordered region (amino acids 669-678) improved R(free) and average B factors in the 2.7 A model. We could not identify significant densities indicating specific antibiotic binding sites in the AcrB R32 space group datasets under the soaking conditions tested.     

Three-dimensional structure of murine anti-p-azophenylarsonate Fab 36-71. 1. X-ray crystallography, site-directed mutagenesis, and modeling of the complex with hapten

The structure of the antigen-binding fragment (Fab) of an anti-p-azophenylarsonate monoclonal antibody, 36-71, bearing a major cross-reactive idiotype of A/J mice has been refined to an R factor of 24.8% at a resolution of 1.85 A. The previously solved partial structure of this Fab at a resolution of 2.9 A (Rose et al., 1990) was used as an initial model for refinement against the high-resolution data. The complex with hapten has been modeled by docking the small-molecule crystal structure of phenylarsonic acid into the structure of the native Fab on the basis of a low-resolution electron density map of the complex. In this model, residue Arg-96 in the light chain and residues Asn-35, Trp-47, and Ser-99 in the heavy chain contact the arsonate moiety of the hapten; an additional bond is found between the arsonate group and a tightly bound water molecule. The phenyl moiety of the hapten packs against two tyrosine side chains at positions 50 and 106 in the heavy chain. Residue Arg-96 in the light chain had been implicated as involved in hapten binding on the basis of previous experiments, and indeed, this residue appears to play a crucial role in this model. Experiments employing site-directed mutagenesis directly support this conclusion. The heavy-chain complementarity-determining regions have novel conformations not previously observed in immunoglobulins except for the recently solved anti-p-azophenylarsonate Fab R 19.9 (Lascombe et al., 1989).     

Simulated annealing coupled replica exchange molecular dynamics—An efficient conformational sampling method

Molecular dynamics simulated annealing (SA-MD) simulations are frequently used for refinement and optimization of peptide and protein structures. Depending on the simulation conditions and simulation length SA-MD simulations can be trapped in locally stable conformations far from the global optimum. As an alternative replica exchange molecular dynamics (RexMD) simulations can be used which allow exchanges between high and low simulation temperatures at all stages of the simulation. A significant drawback of RexMD simulations is, however, the rapid increase of the replica number with increasing system size to cover a desired temperature range. A combined SA-MD and RexMD approach termed SA-RexMD is suggested that employs a small number of replicas (4) and starts initially with a set of high simulation temperatures followed by gradual cooling of the set of temperatures until a target temperature has been reached. The protocol has been applied for the folding of several peptide systems and for the refinement of protein model structures. In all the cases, the SA-RexMD method turned out to be significantly more efficient in reaching low energy structures and also structures close to experiment compared to continuous MD simulations at the target temperature and to SA-MD simulations at the same computational demand. The approach is well suited for applications in structure refinement and for systematic force field improvement.     

On using time-averaging restraints in molecular dynamics simulation

Introducing experimental values as restraints into molecular dynamics (MD) simulations to bias the values of particular molecular properties, such as nuclear Overhauser effect intensities or distances, 3J coupling constants, chemical shifts or crystallographic structure factors, towards experimental values is a widely used structure refinement method. To account for the averaging of experimentally derived quantities inherent in the experimental techniques, time-averaging restraining methods may be used. In the case of structure refinement using 3J coupling constants from NMR experiments, time-averaging methods previously proposed can suffer from large artificially induced structural fluctuations. A modified time-averaged restraining potential energy function is proposed which overcomes this problem. The different possible approaches are compared using stochastic dynamics simulations of antamanide, a cyclic peptide of ten residues.     

The high-resolution crystal structure of porcine pepsinogen.

The structure of porcine pepsinogen at pH 6.1 has been refined to an R-factor of 0.173 for data extending to 1.65 A. The final model contains 180 solvent molecules and lacks density for residues 157-161. The structure of this aspartic proteinase zymogen possesses many of the characteristics of pepsin, the mature enzyme. The secondary structure of the zymogen consists predominantly of beta-sheet, with an approximate 2-fold axis of symmetry. The activation peptide packs into the active site cleft, and the N-terminus (1P-9P) occupies the position of the mature N-terminus (1-9). Thus changes upon activation include excision of the activation peptide and proper relocation of the mature N-terminus. The activation peptide or residues of the displaced mature N-terminus make specific interactions with the substrate binding subsites. The active site of pepsinogen is intact; thus the lack of activity of pepsinogen is not due to a deformation of the active site. Nine ion pairs in pepsinogen may be important in the advent of activation and involve the activation peptide or regions of the mature N-terminus which are relocated in the mature enzyme. The activation peptide-pepsin junction, 44P-1, is characterized by high thermal parameters and weak density, indicating a flexible structure which would be accessible to cleavage. Pepsinogen is an appropriate model for the structures of other zymogens in the aspartic proteinase family.     

Cross-beta order and diversity in nanocrystals of an amyloid-forming peptide

The seven-residue peptide GNNQQNY from the N-terminal region of the yeast prion protein Sup35, which forms amyloid fibers, colloidal aggregates and highly ordered nanocrystals, provides a model system for characterizing the elusively protean cross-beta conformation. Depending on preparative conditions, orthorhombic and monoclinic crystals with similar lath-shaped morphology have been obtained. Ultra high-resolution (<0.5A spacing) electron diffraction patterns from single nanocrystals show that the peptide chains pack in parallel cross-beta columns with approximately 4.86A axial spacing. Mosaic striations 20-50 nm wide observed by electron microscopy indicate lateral size-limiting crystal growth related to amyloid fiber formation. Frequently obtained orthorhombic forms, with apparent space group symmetry P2(1)2(1)2(1), have cell dimensions ranging from /a/=22.7-21.2A, /b/=39.9-39.3A, /c/=4.89-4.86A for wet to dried states. Electron diffraction data from single nanocrystals, recorded in tilt series of still frames, have been mapped in reciprocal space. However, reliable integrated intensities cannot be obtained from these series, and dynamical electron diffraction effects present problems in data analysis. The diversity of ordered structures formed under similar conditions has made it difficult to obtain reproducible X-ray diffraction data from powder specimens; and overlapping Bragg reflections in the powder patterns preclude separated structure factor measurements for these data. Model protofilaments, consisting of tightly paired, half-staggered beta strands related by a screw axis, can be fit in the crystal lattices, but model refinement will require accurate structure factor measurements. Nearly anhydrous packing of this hydrophilic peptide can account for the insolubility of the crystals, since the activation energy for rehydration may be extremely high. Water-excluding packing of paired cross-beta peptide segments in thin protofilaments may be characteristic of the wide variety of anomalously stable amyloid aggregates.     

Low-resolution structure refinement in electron microscopy

A real-space structure refinement method, originally developed for macromolecular X-ray crystallography, has been applied to protein structure analysis by electron microscopy (EM). This method simultaneously optimizes the fit of an atomic model to a density map and the stereo-chemical properties of the model by minimizing an energy function. The performance of this method is characterized at different resolution and signal-to-noise ratio conditions typical for EM electron density maps. A multi-resolution scheme is devised to improve the convergence of the refinement on the global energy minimum. Applications of the method to various model systems are demonstrated here. The first case is the arrangement of FlgE molecules in the helical filament of flagellar hook, in which refinement with segmented rigid bodies improves the density correlation and reduces severe van der Waals contacts among the symmetry-related subunits. The second case is a conformational analysis of the NSF AAA ATPase in which a multi-conformer model is used in the refinement to investigate the arrangement of the two ATPase domains in the molecule. The third case is a docking simulation in which the crystal structure of actin and the NOE data from NMR experiments on the dematin headpiece are combined with a low-resolution EM density map to generate an atomic model of the F-actin-dematin headpiece structure.     

Modeling the structure of Pyrococcus furiosus rubredoxin by homology to other X-ray structures.

The three-dimensional structure of rubredoxin from the hyperthermophilic archaebacterium, Pyrococcus furiosus, has been modeled from the X-ray crystal structures of three homologous proteins from Clostridium pasteurianum, Desulfovibrio gigas, and Desulfovibrio vulgaris. All three homology models are similar. When comparing the positions of all heavy atoms and essential hydrogen atoms to the recently solved crystal structure (Day, M. W., et al., 1992, Protein Sci. 1, 1494-1507) of the same protein, the homology model differ from the X-ray structure by 2.09 A root mean square (RMS). The X-ray and the zinc-substituted NMR structures (Blake, P. R., et al., 1992b, Protein Sci. 1, 1508-1521) show a similar level of difference (2.05 A RMS). On average, the homology models are closer to the X-ray structure than to the NMR structures (2.09 vs. 2.42 A RMS).