Study of conformational rearrangement and refinement of structural homology models by the use of heteronuclear dipolar couplings

For an increasing fraction of proteins whose structures are being studied, sequence homology to known structures permits building of low resolution structural models. It is demonstrated that dipolar couplings, measured in a liquid crystalline medium, not only can validate such structural models, but also refine them. Here, experimental ¹H-¹⁵N, ¹H-¹³C, and ¹³C-¹³C dipolar couplings are shown to decrease the backbone rmsd between various homology models of calmodulin (CaM) and its crystal structure. Starting from a model of the Ca²⁺-saturated C-terminal domain of CaM, built from the structure of Ca²⁺-free recoverin on the basis of remote sequence homology, dipolar couplings are used to decrease the rmsd between the model and the crystal structure from 5.0 to 1.25 Å. A better starting model, built from the crystal structure of Ca²⁺-saturated parvalbumin, decreases in rmsd from 1.25 to 0.93 Å. Similarly, starting from the structure of the Ca²⁺-ligated CaM N-terminal domain, experimental dipolar couplings measured for the Ca²⁺-free form decrease the backbone rmsd relative to the refined solution structure of apo-CaM from 4.2 to 1.0 Å.     

Generation of native-like protein structures from limited NMR data,modern force fields and advanced conformational sampling

Determining an accurate initial native-like protein fold is one of the most important and time-consuming steps of de novo NMR structure determination. Here we demonstrate that high-quality native-like models can be rapidly generated from initial structures obtained using limited NOE assignments, through replica exchange molecular dynamics refinement with a generalized Born implicit solvent (REX/GB). Conventional structure calculations using an initial sparse NOE set were unable to identify a unique topology for the zinc-bound C-terminal domain of E. coli chaperone Hsp33, due to a lack of unambiguous long range NOEs. An accurate overall topology was eventually obtained through laborious hand identification of long range NOEs. However we were able to obtain high-quality models with backbone RMSD values of about 2 Å with respect to the final structures, using REX/GB refinement with the original limited set of initial NOE restraints. These models could then be used to make further assignments of ambiguous NOEs and thereby speed up the structure determination process. The ability to calculate accurate starting structures from the limited unambiguous NOE set available at the beginning of a structure calculation offers the potential of a much more rapid and automated process for NMR structure determination. Jianhan Chen: Authors contributed equally to this work.Hyung-Sik Won: Authors contributed equally to this work.     

RECOORD: a recalculated coordinate database of 500+ proteins from the PDB using restraints from the BioMagResBank

State-of-the-art methods based on CNS and CYANA were used to recalculate the nuclear magnetic resonance (NMR) solution structures of 500+ proteins for which coordinates and NMR restraints are available from the Protein Data Bank. Curated restraints were obtained from the BioMagResBank FRED database. Although the original NMR structures were determined by various methods, they all were recalculated by CNS and CYANA and refined subsequently by restrained molecular dynamics (CNS) in a hydrated environment. We present an extensive analysis of the results, in terms of various quality indicators generated by PROCHECK and WHAT_CHECK. On average, the quality indicators for packing and Ramachandran appearance moved one standard deviation closer to the mean of the reference database. The structural quality of the recalculated structures is discussed in relation to various parameters, including number of restraints per residue, NOE completeness and positional root mean square deviation (RMSD). Correlations between pairs of these quality indicators were generally low; for example, there is a weak correlation between the number of restraints per residue and the Ramachandran appearance according to WHAT_CHECK (r = 0.31). The set of recalculated coordinates constitutes a unified database of protein structures in which potential user- and software-dependent biases have been kept as small as possible. The database can be used by the structural biology community for further development of calculation protocols, validation tools, structure-based statistical approaches and modeling. The RECOORD database of recalculated structures is publicly available from http://www.ebi.ac.uk/msd/recoord.     

Towards atomic resolution with crystals grown in gel: the case of thaumatin seen at room temperature

One reason for introducing a gel in the crystallization medium of proteins is its ability to reduce convection in solution. This can lead to better nucleation and growth conditions, and to crystals having enhanced diffraction properties. We report here the X-ray characterization at room temperature of high-quality crystals of the intensely sweet thaumatin prepared in a sodium tartrate solution gelified with 0.15% (m/v) agarose. Using a synchrotron radiation, these crystals diffracted to a previously unachieved resolution. A diffraction dataset was collected from four crystals at a resolution of 1.2 A with a R(sym) of 3.6% and a completeness of 99%. Refinement was carried out to a final crystallographic R-factor of 12.0%. The quality of the electron density map allowed for the observation of fine structural details in the protein and its solvation shell. Crystallization in gel might be used more generally to improve the quality of macromolecular crystals. Advantages provided by the gelified medium in the frame of structural studies are emphasized.     

SESAME: A least-squares approach to the evaluation of protein structures computed from NMR data

Summary A method is proposed for defining a probability distribution on an ensemble of protein conformations from a 2D NOE spectrum, while at the same time back-calculating the experimental spectrum from the ensemble. This enables one to assess the relative quality and significance of the conformations, and to test the consistency of the ensemble as a whole with the experimental spectrum. The method eliminates the need to integrate the cross-peak intensities and is surprisingly insensitive to random noise in the spectrum. In this communication, these advantages are demonstrated by applying the method to simulated data, for which the correct result is already known.Abbreviations NOE nuclear Overhauser effect - BPTI basic pancreatic trypsin inhibitor - rmsd root-mean-square deviation     

Conformational energy range of ligands in protein crystal structures: The difficult quest for accurate understanding

In this review, we address a fundamental question: What is the range of conformational energies seen in ligands in protein‐ligand crystal structures? This value is important biophysically, for better understanding the protein‐ligand binding process; and practically, for providing a parameter to be used in many computational drug design methods such as docking and pharmacophore searches. We synthesize a selection of previously reported conflicting results from computational studies of this issue and conclude that high ligand conformational energies really are present in some crystal structures. The main source of disagreement between different analyses appears to be due to divergent treatments of electrostatics and solvation. At the same time, however, for many ligands, a high conformational energy is in error, due to either crystal structure inaccuracies or incorrect determination of the reference state. Aside from simple chemistry mistakes, we argue that crystal structure error may mainly be because of the heuristic weighting of ligand stereochemical restraints relative to the fit of the structure to the electron density. This problem cannot be fixed with improvements to electron density fitting or with simple ligand geometry checks, though better metrics are needed for evaluating ligand and binding site chemistry in addition to geometry during structure refinement. The ultimate solution for accurately determining ligand conformational energies lies in ultrahigh‐resolution crystal structures that can be refined without restraints.     

Neutron and X-ray crystallographic analysis of Achromobacter protease I at pD 8.0: Protonation states and hydration structure in the free-form

The structure of the free-form of Achromobacter protease I (API) at pD 8.0 was refined by simultaneous use of single crystal X-ray and neutron diffraction data sets to investigate the protonation states of key catalytic residues of the serine protease. Occupancy refinement of the catalytic triad in the active site of API free-form showed that ca. 30% of the imidazole ring of H57 and ca. 70% of the hydroxyl group of S194 were deuterated. This observation indicates that a major fraction of S194 is protonated in the absence of a substrate. The protonation state of the catalytic triad in API was compared with the bovine β-trypsin–BPTI complex. The comparison led to the hypothesis that close contact of a substrate with S194 could lower the acidity of its hydroxyl group, thereby allowing H57 to extract the hydrogen from the hydroxyl group of S194. H210, which is a residue specific to API, does not form a hydrogen bond with the catalytic triad residue D113. Instead, H210 forms a hydrogen bond network with S176, H177 and a water molecule. The close proximity of the bulky, hydrophobic residue W169 may protect this hydrogen bond network, and this protection may stabilize the function of API over a wide pH range.     

Evidence of colorectal cancer risk associated variant Lys25Ser in the proximity of human bone morphogenetic protein 2

Colorectal cancer (CRC) is the third most prevalent cancer and fourth leading cause of cancer-related deaths globally. It has been shown that the nsSNP variants play an important role in diseases, however it remained unclear how these variants are associated with the disease. Recently, several CRC risk associated SNPs have been discovered, however rs961253 (Lys25Arg at 20p12.3) located in the proximity of bone morphogenetic protein 2 (Bmp2) and fermitin family homolog 1 Fermt1 genes have been reported to be highly associated with the CRC risk. Here we provide evidence for the first time in silico biological functional and structural implications of non-synonymous (nsSNPs) CRC disease-associated variant Lys25Arg via molecular dynamic (MD) simulation. Protein structural analysis was performed with a particular variant allele (A/C, Lys25Arg) and compared with the predicted native protein structure. Our results showed that this nsSNP will cause changes in the protein structure and as a result is associated with the disease. In addition to the native and mutant 3D structures of CRC associated risk allele protein domain (CRAPD), they were also analyzed using solvent accessibility models for further protein stability confirmation. Taken together, this study confirmed that this variant has functional effect and structural impact on the CRAPD and may play an important role in CRC disease progression; hence it could be a reasonable approach for studying the effect of other deleterious variants in future studies.     

Sampling and scoring: A marriage made in heaven

Most structure prediction algorithms consist of initial sampling of the conformational space, followed by rescoring and possibly refinement of a number of selected structures. Here we focus on protein docking, and show that while decoupling sampling and scoring facilitates method development, integration of the two steps can lead to substantial improvements in docking results. Since decoupling is usually achieved by generating a decoy set containing both non‐native and near‐native docked structures, which can be then used for scoring function construction, we first review the roles and potential pitfalls of decoys in protein–protein docking, and show that some type of decoys are better than others for method development. We then describe three case studies showing that complete decoupling of scoring from sampling is not the best choice for solving realistic docking problems. Although some of the examples are based on our own experience, the results of the CAPRI docking and scoring experiments also show that performing both sampling and scoring generally yields better results than scoring the structures generated by all predictors. Next we investigate how the selection of training and decoy sets affects the performance of the scoring functions obtained. Finally, we discuss pathways to better alignment of the two steps, and show some algorithms that achieve a certain level of integration. Although we focus on protein–protein docking, our observations most likely also apply to other conformational search problems, including protein structure prediction and the docking of small molecules to proteins.Proteins 2013; 81:1874–1884. © 2013 Wiley Periodicals, Inc.     

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.