New strategies for the determination of macromolecular structure in solution

Non-crystallographic approaches to the determination of protein structure must solve the problem of insufficient and low information content experimental data. Most successful methods augment experimentation with theoretical constraints (for example, potential energy functions or optimization error metrics). We believe it is important to separate the contributions of experimentation and theory in the construction of protein structure. The PROTEAN system defines protein topology on the basis of experimental data alone. Its performance on three data sets, derived from the lac-repressor headpiece of E. coli, sperm whale myoglobin, and domain 1 of bacteriophage T4 lysozyme, indicates that there may be families of related conformations that are consistent with the experimental data. These conformations provide insight into the strengths and weaknesses in the data sets. They also provide a set of structures with which to begin theoretical refinements. We outline here a strategy which maintains a clear distinction between refinements based on theory and those based on experiment, and thus allows a careful analysis of the properties of such refinement methods.     

Crystal structure of human seminal diferric lactoferrin at 3.4 Angstrom resolution

Lactoferrin was purified from human seminal fluid obtained from the semen bank. The purified samples were saturated with Fe3+ and crystallized by microdialysis method. The crystals belong to orthorhombic space group P21212, with a = 55.9 Angstrom. b = 97.2 Angstrom, c = 156.1 Angstrom and Z = 4. The structure was determined with molecular replacement method and refined to an R factor of 18.7% for all the data to 3.4 Angstrom resolution. The overall structure of seminal lactoferrin is similar to human colostrum lactoferrin. The amino acid sequence of seminal lactoferrin shows that it has one amino acid less than human colostrum lactoferrin and the structure of its N-terminal region is far more ordered than other lactoferrins. The structure of the iron-binding site and its immediate surroundings indicate well defined features.     

Ensemble refinement of protein crystal structures: validation and application

X-ray crystallography typically uses a single set of coordinates and B factors to describe macromolecular conformations. Refinement of multiple copies of the entire structure has been previously used in specific cases as an alternative means of representing structural flexibility. Here, we systematically validate this method by using simulated diffraction data, and we find that ensemble refinement produces better representations of the distributions of atomic positions in the simulated structures than single-conformer refinements. Comparison of principal components calculated from the refined ensembles and simulations shows that concerted motions are captured locally, but that correlations dissipate over long distances. Ensemble refinement is also used on 50 experimental structures of varying resolution and leads to decreases in R(free) values, implying that improvements in the representation of flexibility observed for the simulated structures may apply to real structures. These gains are essentially independent of resolution or data-to-parameter ratio, suggesting that even structures at moderate resolution can benefit from ensemble refinement.     

Comparison of protein solution structures refined by molecular dynamics simulation in vacuum,with a generalized Born model,and with explicit water

The inclusion of explicit solvent water in molecular dynamics refinement of NMR structures ought to provide the most physically meaningful accounting for the effects of solvent on structure, but is computationally expensive. In order to evaluate the validity of commonly used vacuum refinements and of recently developed continuum solvent model methods, we have used three different methods to refine a set of NMR solution structures of a medium sized protein, Escherichia coliglutaredoxin 2, from starting structures calculated using the program DYANA. The three different refinement protocols used molecular dynamics simulated annealing with the program AMBER in vacuum (VAC), including a generalized Born (GB) solvent model, and a full calculation including explicit solvent water (WAT). The structures obtained using the three methods of refinements were very similar, a reflection of their generally well-determined nature. However, the structures refined with the generalized Born model were more similar to those from explicit water refinement than those refined in vacuum. Significant improvement was seen in the percentage of backbone dihedral angles in the most favored regions of , space and in hydrogen bond pattern for structures refined with the GB and WAT models, compared with the structures refined in vacuum. The explicit water calculation took an average of 200 h of CPU time per structure on an SGI cluster, compared to 15–90 h for the GB calculation (depending on the parameters used) and 2 h for the vacuum calculation. The generalized Born solvent model proved to be an excellent compromise between the vacuum and explicit water refinements, giving results comparable to those of the explicit water calculation. Some improvement for and angle distribution and hydrogen bond pattern can also be achieved by energy minimizing the vacuum structures with the GB model, which takes a much shorter time than MD simulations with the GB model.     

高分辨率高精度蛋白质晶体结构

自从六十年代初期应用X射线晶体衍射方法测定第一个蛋白质(肌红蛋白)的三维结构以来,我们关于蛋白质三维结构的知识已有了相当的积累.截止86年7月的统计,已有286个蛋白质分子的原子坐标存入国际蛋白质数据库.以此为基础,一些重要生命活动的结构机理已经在三维水平上得到精细阐明(参见),一个由分子生物学和X射线晶体学结合形成的新领域——蛋白质晶体学,应运而生.二十多年来,这一领域的发展已经历了二个阶段.第一阶段大体从六十年代初到七十年代中期,在这期间,分析方法和技术从突破发展到成熟和实用,并有近40个蛋白质结构测定出来,使我们获得了蛋白质结构知识的面面观.     

Molecular modeling studies on the ribosome

In the absence of a high resolution crystal structure for the ribosome, numerous research groups are carrying out low resolution structural studies using neutron diffraction, electron microscopy, fluorescence energy transfer, chemical crosslinking, chemical footprinting studies, and other methods. We have developed a computer-based refinement method for incorporating these data into low resolution three-dimensional models. The method is based on a molecular mechanics approach, with proteins represented by spherical particles of suitable diameter and the ribosomal RNA represented by a string of spherical pseudoatoms, one for each nucleotide. Experimental data are used to derive constraints that are introduced through a special force field (potential function). Models are refined by simulated annealing. Since every term in the force field is quadratic, any model that satisfies all of the input data has an energy of zero; higher energies indicate residual unsatisfied constraints. The residual energy provides a quantitative statement of model quality and can be used to identify conflicts in the experimental data. The method has been applied to the refinement of a low resolution model for the 30S subunit (the small subunit) of theE. coli ribosome. Since this is a very underdetermined system, the range of acceptable models has also been explored. This provides an estimate of the resolution of the structure, which is about 15 Å overall, with the uncertainty in position of individual nucleotides ranging from about 5 Å to 50 Å.     

A model-based parallel origin and orientation refinement algorithm for cryoTEM and its application to the study of virus structures

We present a model-based parallel algorithm for origin and orientation refinement for 3D reconstruction in cryoTEM. The algorithm is based upon the Projection Theorem of the Fourier Transform. Rather than projecting the current 3D model and searching for the best match between an experimental view and the calculated projections, the algorithm computes the Discrete Fourier Transform (DFT) of each projection and searches for the central section ("cut") of the 3D DFT that best matches the DFT of the projection. Factors that affect the efficiency of a parallel program are first reviewed and then the performance and limitations of the proposed algorithm are discussed. The parallel program that implements this algorithm, called PO(2)R, has been used for the refinement of several virus structures, including those of the 500 Angstroms diameter dengue virus (to 9.5 Angstroms resolution), the 850 Angstroms mammalian reovirus (to better than 7A), and the 1800 Angstroms paramecium bursaria chlorella virus (to 15 Angstroms).     

Relaxation matrix refinement of the solution structure of squash trypsin inhibitor

The structure of the small squash trypsin inhibitor CMTI-I is refined by directly minimizing the difference between the observed two-dimensional nuclear Overhauser enhancement (NOE) intensities and those calculated by the full relaxation matrix approach. To achieve this, a term proportional to this difference was added to the potential energy function of the molecular dynamics program X-PLOR. Derivatives with respect to atomic co-ordinates are calculated analytically. Spin diffusion effects are thus accounted for fully during the refinement. Initial structures for the refinement were those determined recently by solution nuclear magnetic resonance using the isolated two-spin approximation to derive distance range estimates. The fits to the nuclear magnetic resonance data improve significantly with only small shifts in the refined structures during a few cycles of conjugate gradient minimization. However, larger changes (approximately 1 A) in the conformation occur during simulated annealing, which is accompanied by a further reduction of the difference between experimental and calculated two-dimensional NOE intensities. The refined structures are closer to the X-ray structure of the inhibitor complexed with trypsin than the initial structures. The root-mean-square difference for backbone atoms between the initial structures and the X-ray structure is 0.96 A, and that between the refined structures and the X-ray structure 0.61 A.     

Refinement of 3D models of horseradish peroxidase isoenzyme C: Predictions of 2D NMR assignments and substrate binding sites

In this study, two alternative three-dimensional (3D) models of horseradish peroxidase (HRP-C)—differing mainly in the structure of a long untemplated insertion—were refined, systematically assessed, and used to make predictions that can both guide and be tested by future experimental studies. A key first step in the model-building process was a procedure for multiple sequence alignment based on structurally conserved regions and key conserved residues, including those side chains providing ligands to the two Ca²⁺ binding sites. The model refinements reported here include (1) optimization of side-chain conformations; (3) addition of structural waters using a template-independent procedure; (2) structural refinement of the untemplated 34 amino acid insertion located between the F and G helices, using both energy criteria and NMR data; (4) unconstrained energy optimization of the refined models. Using these procedures, two refined structures of HRP-C were obtained, differing mainly in the conformation of this long insertion. The presence of residues in this insertion that could potentially interact with bound substrates suggests a functional role that may be related to the general ability of class III peroxidases to form stable 1:1 complexes with a variety of substrates. The structural validity of the models was systematically assessed by a variety of criteria. Most notably, the ProsaII z scores and Profiles 3D scores of the two HRP-C models indicated that they are significantly better than would be obtained by simple amino acid replacement, using any of the known structures as a template. These two 3D HRP-C models, were then used to predict candidate residues for the assignment of NOESY cross-peaks previously noted in 2D-NMR studies. Specifically, the residues known as Ile X, Phe A, Phe B, aliphatic residue Q, and Ile T. Candidate substrate binding sites were also identified and compared with experimentally based predictions. This work is timely because new X-ray structures are anticipated that will facilitate the validation of these procedures. © 1996 Wiley-Liss, Inc.     

Computational insights into the O<Subscript>2</Subscript>-evolving complex of photosystem II

Mechanistic investigations of the water-splitting reaction of the oxygen-evolving complex (OEC) of photosystem II (PSII) are fundamentally informed by structural studies. Many physical techniques have provided important insights into the OEC structure and function, including X-ray diffraction (XRD) and extended X-ray absorption fine structure (EXAFS) spectroscopy as well as mass spectrometry (MS), electron paramagnetic resonance (EPR) spectroscopy, and Fourier transform infrared spectroscopy applied in conjunction with mutagenesis studies. However, experimental studies have yet to yield consensus as to the exact configuration of the catalytic metal cluster and its ligation scheme. Computational modeling studies, including density functional (DFT) theory combined with quantum mechanics/molecular mechanics (QM/MM) hybrid methods for explicitly including the influence of the surrounding protein, have proposed chemically satisfactory models of the fully ligated OEC within PSII that are maximally consistent with experimental results. The inorganic core of these models is similar to the crystallographic model upon which they were based, but comprises important modifications due to structural refinement, hydration, and proteinaceous ligation which improve agreement with a wide range of experimental data. The computational models are useful for rationalizing spectroscopic and crystallographic results and for building a complete structure-based mechanism of water-splitting in PSII as described by the intermediate oxidation states of the OEC. This review summarizes these recent advances in QM/MM modeling of PSII within the context of recent experimental studies.     

Outside‐binding site mutations modify the active site's shapes in neuraminidase from influenza A H1N1

The recent occurrence of 2009 influenza A (H1N1) pandemic as well as others has raised concern of a far more dangerous outcome should this virus becomes resistant to current drug therapies. The number of clinical cases that are resistant to oseltamivir (Tamiflu®) is larger than the limited number of neuraminidase (NA) mutations (H275Y, N295S, and I223R) that have been identified at the active site and that are associated to oseltamivir resistance. In this study, we have performed a comparative analysis between a set of NAs that have the most representative mutations located outside the active site. The recently crystallized NA‐oseltamivir complex (PDB ID: 3NSS) was used as a wild‐type structure. After selecting the target NA sequences, their three‐dimensional (3D) structure was built using 3NSS as a template by homology modeling. The 3D NA models were refined by molecular dynamics (MD) simulations. The refined models were used to perform a docking study, using oseltamivir as a ligand. Furthermore, the docking results were refined by free‐energy analysis using the MM‐PBSA method. The analysis of the MD simulation results showed that the NA models reached convergence during the first 10 ns. Visual inspection and structural measures showed that the mutated NA active sites show structural variations. The docking and MM‐PBSA results from the complexes showed different binding modes and free energy values. These results suggest that distant mutations located outside the active site of NA affect its structure and could be considered to be a new source of resistance to oseltamivir, which agrees with reports in the clinical literature. © 2012 Wiley Periodicals, Inc.     

Refinement of the solution structure of the B DNA hexamer 5'd(C-G-T-A-C-G)2 on the basis of inter-proton distance data

A restrained least-squares refinement of the solution structure of the self-complementary B DNA hexamer 5'd(C-G-T-A-C-G)2 is presented. The structure is refined on the basis of 190 inter-proton distances determined by pre-steady-state nuclear Overhauser enhancement measurements. Two refinements were carried out starting from two initial B DNA structures differing by an overall root-mean-square (r.m.s.) difference of 0.32 A. In both cases, the final r.m.s. difference between the experimental and calculated inter-proton distances was 0.12 A compared to 0.61 A and 0.58 A for the two initial structures. The difference between the two refined structures is small, with an overall r.m.s. difference of 0.16 A, and represents the error in the refined co-ordinates. The refined structures have a B-type conformation with local structural variations in backbone and glycosidic bond torsion angles, and base-pair propellor twist, base roll, base tilt and local helical twist angles.     

The crystal structure of poplar apoplastocyanin at 1.8-A resolution. The geometry of the copper-binding site is created by the polypeptide

The three-dimensional structure of apoplastocyanin from poplar leaves (Populus nigra var. italica) has been determined by x-ray diffraction at 1.8-A resolution. The structure closely resembles that of the holoprotein. In particular, the positions of the copper-binding residues in the apo- and holoproteins differ by only 0.1-0.3 A. This indicates that the irregular geometry of the "type 1" copper site is imposed upon the metal atom by the polypeptide moiety. A 180 degrees rotation of one solvent-exposed histidine imidazole ring about C beta-C gamma appears to facilitate access to the copper site. The close structural similarity between apo-, Cu-(II)-, and Cu(I)-plastocyanin was initially demonstrated by means of electron density difference maps. Two series of restrained least squares refinement calculations for apoplastocyanin, originating from different sets of atomic positional parameters, were carried out in parallel. Both refinements converged to the same model which, when fully refined, had a residual R = 0.16. Forty-two water molecules were located during the refinement.     

Protein structure determination using a database of interatomic distance probabilities

The accelerated pace of genomic sequencing has increased the demand for structural models of gene products. Improved quantitative methods are needed to study the many systems (e.g., macromolecular assemblies) for which data are scarce. Here, we describe a new molecular dynamics method for protein structure determination and molecular modeling. An energy function, or database potential, is derived from distributions of interatomic distances obtained from a database of known structures. X-ray crystal structures are refined by molecular dynamics with the new energy function replacing the Van der Waals potential. Compared to standard methods, this method improved the atomic positions, interatomic distances, and side-chain dihedral angles of structures randomized to mimic the early stages of refinement. The greatest enhancement in side-chain placement was observed for groups that are characteristically buried. More accurate calculated model phases will follow from improved interatomic distances. Details usually seen only in high-resolution refinements were improved, as is shown by an R-factor analysis. The improvements were greatest when refinements were carried out using X-ray data truncated at 3.5 A. The database potential should therefore be a valuable tool for determining X-ray structures, especially when only low-resolution data are available.     

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%...     

Paramagnetism-based refinement strategy for the solution structure of human alpha-parvalbumin

In the frame of a research aimed at the detailed structural characterization of human calcium-binding proteins of the EF-hand family, the solution structure of human alpha-parvalbumin has been solved by NMR and refined with the help of substitution of the Ca(2+) ion in the EF site with the paramagnetic Dy(3+) ion. A simple (1)H-(15)N HSQC spectrum allowed the NH assignments based on the properties of Dy(3+). This allowed us to exploit pseudocontact shifts and residual dipolar couplings for solution structure refinement. The backbone and heavy atom RMSD are 0.55 +/- 0.08 and 1.02 +/- 0.08 A, respectively, and decrease to 0.39 +/- 0.05 and 0.90 +/- 0.06 A upon refinement with paramagnetism-based restraints. The RMSD for the metal itself in the EF site in the refined structure is 0.26 +/- 0.12 A. Backbone NH R(1), R(2), and NOE measured at two temperatures show the protein to be relatively rigid. The NH orientations are well determined by the paramagnetism-based restraints. This allows us to detect small but significant local structural differences with the orthologue protein from rat, whose X-ray structure is available at 2.0 A resolution. All differences are related to local changes in the amino acidic composition.     

Structural refinement of protein segments containing secondary structure elements: Local sampling, knowledge-based potentials, and clustering

In this article, we present an iterative, modular optimization (IMO) protocol for the local structure refinement of protein segments containing secondary structure elements (SSEs). The protocol is based on three modules: a torsion-space local sampling algorithm, a knowledge-based potential, and a conformational clustering algorithm. Alternative methods are tested for each module in the protocol. For each segment, random initial conformations were constructed by perturbing the native dihedral angles of loops (and SSEs) of the segment to be refined while keeping the protein body fixed. Two refinement procedures based on molecular mechanics force fields - using either energy minimization or molecular dynamics - were also tested but were found to be less successful than the IMO protocol. We found that DFIRE is a particularly effective knowledge-based potential and that clustering algorithms that are biased by the DFIRE energies improve the overall results. Results were further improved by adding an energy minimization step to the conformations generated with the IMO procedure, suggesting that hybrid strategies that combine both knowledge-based and physical effective energy functions may prove to be particularly effective in future applications.     

Molecular and crystal structure of galactinol dihydrate [1-O-(alpha-D-galactopyranosyl)-myo-inositol dihydrate

The crystal structure of galactinol dihydrate has been determined by X-ray diffraction. The crystal belongs to the orthorhombic system, space group P2(1)2(1)2, a = 15.898(6), b = 19.357(5), c = 5.104(4) A, and Z = 4. The structure was refined to R = 0.044 for 1818 observed structure amplitudes. The primary hydroxyl group exhibits twofold orientational disorder. The linkage conformation is close to those of alpha-(1 --> 4) linkages in methyl alpha-maltotrioside tetrahydrate and erlose trihydrate. Although there is no interring hydrogen bond in galactinol, an indirect interring hydrogen bond including a water molecule is present. The observed conformation is additionally stabilized by the indirect interring hydrogen bond. The global minimum in the relaxed-residue energy map based on the MM3(92) force-field is close to the observed conformation in the crystal structure. All hydroxyl, ring and water oxygen atoms are involved in a complex three-dimensional hydrogen-bonding network.