Estimation of the quality of refined protein crystal structures

Crystallographic R_work and R_free values, which are measures of the ability of the models of macromolecular structures to explain the crystallographic data on which they are based, are often used to assess structure quality. It is widely known, and confirmed here that both are sensitive to the methods used to compute them, and can be manipulated to improve the apparent quality of the model. As an alternative it is proposed here that the quality of crystallographic models should be assessed using a global goodness-of-fit metric R_O2A/R_work where R_O2A is the number of reflections used for refinement divided by the number of nonhydrogen atoms in the structure, and R_work is the working R-factor of the refined structure. Also, analysis of structures in the Protein Data Bank suggests that many data sets have been truncated at high resolution, thereby improving the R-factor statistics. To discourage this practice, it is proposed that the resolution of a dataset be defined as the resolution of the shell of data where <I/σ_I> falls to 1. The proposed goodness-of-fit metric encourages investigators to use all the data available rather than a truncated subset.     

Water-Refined Solution Structure of the Human Grb7-SH2 Domain in Complex with the erbB2 Receptor Peptide pY1139

We report a refinement in implicit water of the previously published solution structure of the Grb7-SH2 domain bound to the erbB2 receptor peptide pY1139. Structure quality measures indicate substantial improvement, with residues in the most favored regions of the Ramachandran plot increasing by 14 % and with WHAT IF statistics (Vriend, G. J. Mol. Graph., 1990, 8(1), 52-56) falling closer to expected values for well-refined structures.     

Molecular modelling of urease accessory interaction proteins of Helicobacter Pylori J 99 and predicting an interruption in interaction by Vigna radiata Defensins

Helicobacter pylori is the major causative agent of Gastric carcinoma. Significance of the urease accessory interaction proteins are emphasized in colonization of human gastric mucosa and efficient infection of H. pylori. Here an attempt is made to explore the structure and properties of urease accessory interaction proteins from Helicobacter pylori J 99. The proteins chosen for the study are ureH, ureI, nikR, groL and flgS based on the interaction map available from STRING database. The above mentioned proteins do not have a comprehensive three dimensional structure. Hence the models were generated using PSI-BLAST (Position Specific Iterative-Blast) and MODELLER 9V8. Physicochemical characterization encompasses pI, EC, AI, II and GRAVY. Secondary structure was predicted using PSI-PRED. Functional characterization was done by SOSUI and DISULFIND Servers and refinement of structure was done using Ramachandran plot analysis. RMS-Z values were calculated using Q-MEAN Server and CHIMERA was used for molecular simulation studies. Plant defensins from Vigna radiata are successfully docked to the modeled structures and thus interaction could be possibly prevented. These results will pave way for further selective inhibition of H. pylori colonization and in vivo survival by employing plant defensins from Vigna radiata (VrD1 & VrD2). The work will prove that plant defensins provides anticancer relief too.     

On the validation of crystallographic symmetry and the quality of structures

In 2008, Zwart and colleagues observed that the fraction of the structures deposited in the PDB alleged to have “pseudosymmetry” or “special noncrystallographic symmetry” (NCS) was about 6%, and that this percentage was rising annually. A few years later, Poon and colleagues found that 2% of all the crystal structures in the PDB belonged to higher symmetry space groups than those assigned to them. Here, I report an analysis of the X-ray diffraction data deposited for this class of structures, which shows that most of the “pseudosymmetry” and “special NCS” that has been reported is in fact true crystallographic symmetry (CS). This distinction is important because the credibility of crystal structures depends heavily on quality control statistics such as R_free that are unreliable when they are computed incorrectly, which they often are when CS is misidentified as “special NCS” or “pseudosymmetry”. When mistakes of this kind are made, artificially low values of R_free can give unjustified confidence in the accuracy of the reported structures.     

In silico model of DSF synthase RpfF protein from Xanthomonas oryzae pv. Oryzae: a novel target for bacterial blight of rice disease

BACKGROUND: Rice plant diseases play a major role as biological constraints on production. One of such rice disease is bacterial leaf blight, caused by Xanthomonas oryzae pv. Oryzae (Xoo). The diffusible signal factor (DSF) synthesized by Xoo has a major role in virulence to rice plant. The DSF synthase RpfF protein, which is related to crotonase superfamily is responsible for the maintaining concentration of DSF. DSF-dependent quorum sensing (QS) system adopts protein- protein interaction mechanism to auto regulates the production of DSF. The antibacterial activity of pesticides against Xoo has not yet been completely understood. Three dimensional structure of RpfF protein was predicted using homology modeling method by MODELLER 9V9 software, SWISS MODEL and GENO3D online tools and structures were validated by Ramachandran plot, TM-Score and RMSD. 3D structure of RpfF (accession number AAL06345) was predicted using DSF synthase of Xanthomonas campestris pv. campestris (Xcc) (PDB ID: 3M6M) as a template. The stereo chemical check reveals the structure developed from the modeller was the best one and the potential ligand binding sites were identified by CASTp Server. The predicted RpfF model provides insight into its structure, active sites and aid in the development of novel inhibitors to control bacterial leaf blight in rice plant. DSF synthase RpfF protein could be used as a novel target to control infection.     

Atomic resolution crystal structures of HIV-1 protease and mutants V82A and I84V with saquinavir

Due to the limited distance data available from the experiments, the structures determined by NMR Spectroscopy may not always be as accurate as desired. Further refinement of the structures is often required and sometimes critical. With the increase of high quality protein structures determined and deposited in PDB Data Bank, commonly shared protein conformational properties can be extracted based on the statistical distributions of the properties in the structural database and used to improve the outcomes of the NMR-determined structures. Here we examine the distributions of protein interatomic distances in known protein structures. We show that based on these distributions, a set of mean-force potentials can be defined for proteins and employed to refine the NMR-determined structures. We report the test results on 70 NMR-determined structures and compare the potential energy, the Ramachandran plot, and the ensemble RMSD of the structures refined with and without using the derived mean-force potentials.     

Oligomerization of the SPP1 scaffolding protein

Viral scaffolding proteins direct polymerization of major capsid protein subunits into icosahedral procapsid structures. The scaffolding protein of bacteriophage SPP1 was engineered with a C-terminal hexahistidine tag (gp11-His₆) and purified. The protein is an α-helical-rich molecule with a very elongated shape as found for internal scaffolding proteins from other phages. It is a 3.3 S tetramer of 93.6 kDa at micromolar concentrations. Intersubunit cross-linking of these tetramers generated preferentially covalently bound dimers, revealing that gp11-His₆ is structurally a dimer of dimers. Incubation at temperatures above 37 °C correlated with a reduction of its α-helical content and a less effective intersubunit cross-linking. Complete loss of secondary structure was observed at temperatures above 60 °C. Refolding of gp11-His₆ thermally denatured at 65 °C led to reacquisition of the protein native ellipticity spectrum but the resulting population of molecules was heterogeneous. Its hydrodynamic behavior was compatible with a mix of 3.3 S elongated tetramers (∼ 90%) and a smaller fraction of 2.4 S dimers (∼ 10%). This population of gp11-His₆ was competent to direct polymerization of the SPP1 major capsid protein gp13 into procapsid-like structures in a newly developed assembly assay in vitro. Although native tetramers were active in assembly, refolded gp11-His₆ showed enhanced binding to gp13 revealing a more active species for interaction with the major capsid protein than native gp11-His₆.     

The high-resolution NMR structure of the early folding intermediate of the Thermus thermophilus ribonuclease H

Elucidation of the high-resolution structures of folding intermediates is a necessary but difficult step toward the ultimate understanding of the mechanism of protein folding. Here, using hydrogen-exchange-directed protein engineering, we populated the folding intermediate of the Thermus thermophilus ribonuclease H, which forms before the rate-limiting transition state, by removing the unfolded regions of the intermediate, including an α-helix and two β-strands (51 folded residues). Using multidimensional NMR, we solved the structure of this intermediate mimic to an atomic resolution (backbone rmsd, 0.51 Å). It has a native-like backbone topology and shows some local deviations from the native structure, revealing that the structure of the folded region of an early folding intermediate can be as well defined as the native structure. The topological parameters calculated from the structures of the intermediate mimic and the native state predict that the intermediate should fold on a millisecond time scale or less and form much faster than the native state. Other factors that may lead to the slow folding of the native state and the accumulation of the intermediate before the rate-limiting transition state are also discussed.     

Structure prediction and evolution of a halo-acid dehalogenase of Burkholderia mallei

Environmental pollutants containing halogenated organic compounds e.g. haloacid, can cause a plethora of health problems. The structural and functional analyses of the gene responsible of their degradation are an important aspect for environmental studies and are important to human well-being. It has been shown that some haloacids are toxic and mutagenic. Microorganisms capable of degrading these haloacids can be found in the natural environment. One of these, a soil-borne Burkholderia mallei posses the ability to grow on monobromoacetate (MBA). This bacterium produces a haloacid dehalogenase that allows the cell to grow on MBA, a highly toxic and mutagenic environmental pollutant. For the structural and functional analysis, a 346 amino acid encoding protein sequence of haloacid dehalogenase is retrieve from NCBI data base. Primary and secondary structure analysis suggested that the high percentage of helices in the structure makes the protein more flexible for folding, which might increase protein interactions. The consensus protein sub-cellular localization predictions suggest that dehalogenase protein is a periplasmic protein 3D2GO server, suggesting that it is mainly employed in metabolic process followed by hydrolase activity and catalytic activity. The tertiary structure of protein was predicted by homology modeling. The result suggests that the protein is an unstable protein which is also an important characteristic of active enzyme enabling them to bind various cofactors and substrate for proper functioning. Validation of 3D structure was done using Ramachandran plot ProsA-web and RMSD score. This predicted information will help in better understanding of mechanism underlying haloacid dehalogenase encoding protein and its evolutionary relationship.     

Structure of a premicellar complex of alkyl sulfates with the interfacial binding surfaces of four subunits of phospholipase A2

The properties of three discrete premicellar complexes (E₁^#, E₂^#, E₃^#) of pig pancreatic group-IB secreted phospholipase A2 (sPLA₂) with monodisperse alkyl sulfates have been characterized [Berg, O. G. et al., Biochemistry 43, 7999–8013, 2004]. Here we have solved the 2.7 Å crystal structure of group-IB sPLA₂ complexed with 12 molecules of octyl sulfate (C₈S) in a form consistent with a tetrameric oligomeric that exists during the E₁^# phase of premicellar complexes. The alkyl tails of the C₈S molecules are centered in the middle of the tetrameric cluster of sPLA₂ subunits. Three of the four sPLA₂ subunits also contain a C₈S molecule in the active site pocket. The sulfate oxygen of a C₈S ligand is complexed to the active site calcium in three of the four protein active sites. The interactions of the alkyl sulfate head group with Arg-6 and Lys-10, as well as the backbone amide of Met-20, are analogous to those observed in the previously solved sPLA₂ crystal structures with bound phosphate and sulfate anions. The cluster of three anions found in the present structure is postulated to be the site for nucleating the binding of anionic amphiphiles to the interfacial surface of the protein, and therefore this binding interaction has implications for interfacial activation of the enzyme.     

Modeling and phylogenetic analysis of cytosolic ascorbate peroxidase (OsAPX1) from rice reveal signature motifs that may play a role in stress tolerance

Ascorbate peroxidase (APX) is a crucial, haeme-containing enzyme of the ascorbate glutathione cycle that detoxifies reactive oxygen species in plants by catalyzing the conversion of hydrogen peroxide to water using ascorbate as a specific electron donor. Different APX isoforms are present in discrete subcellular compartments in rice and their expression is stress regulated. We revealed the homology model of OsAPX1 protein using the crystal structure of soybean GmAPX1 (PDB ID: 2XIF) as template by Modeller 9.12. The resultant OsAPX1 model structure was refined by PROCHECK, ProSA, Verify3D and RMSD that indicated the model structure is reliable with 83 % amino acid sequence identity with template, RMSD (1.4 Å), Verify3D (86.06 %), Zscores (-8.44) and Ramachandran plot analysis showed that conformations for 94.6% of amino acid residues are within the most favoured regions. Investigation revealed two conserved signatures for haeme ligand binding and peroxidase activity in the alpha helical region that may play a significant role during stress.     

Solvent-induced differentiation of protein backbone hydrogen bonds in calmodulin

In apo and holoCaM, almost half of the hydrogen bonds (H-bonds) at the protein backbone expected from the corresponding NMR or X-ray structures were not detected by h3JNC' couplings. The paucity of the h3JNC' couplings was considered in terms of dynamic features of these structures. We examined a set of seven proteins and found that protein-backbone H-bonds form two groups according to the h3JNC' couplings measured in solution. H-bonds that have h3JNC' couplings above the threshold of 0.2 Hz show distance/angle correlation among the H-bond geometrical parameters, and appear to be supported by the backbone dynamics in solution. The other H-bonds have no such correlation; they populate the water-exposed and flexible regions of proteins, including many of the CaM helices. The observed differentiation in a dynamical behavior of backbone H-bonds in apo and holoCaM appears to be related to protein functions.     

Insight from the structural molecular model of cytidylate kinase from Mycobacterium tuberculosis

Mycobacterium tuberculosis is a gram-positive bacterium causes tuberculosis in human. H37Rv strain is a pathogenic strain utilized for tuberculosis research. The cytidylate mono-phosphate (CMP) kinase of Mycobacterium tuberculosis belongs to the family nucleoside mono-phosphate kinase (NMK), this enzyme is required for the bacterial growth. Therefore, it is important to study the structural and functional features of this enzyme in the control of the disease. Hence, we developed the structural molecular model of the CMP kinase protein from Mycobacterium tuberculosis by homology modeling using the software MODELLER (9v10). Based on sequence similarity with protein of known structure (template) of Mycobacterium smegmatis (PDB ID: 3R20) was chosen from protein databank (PDB) by using BLASTp. The energy of constructed models was minimized and the qualities of the models were evaluated by PROCHECK and VERRIFY-3D. Resulted Ramachandran plot analysis showed that conformations for 100.00% of amino acids residues are within the most favored regions. A possible homologous deep cleft active site was identified in the Model using CASTp program. Amino acid composition and polarity of that protein was observed by CLC-Protein Workbench tool. Expasy''s Prot-param server and CYC_REC tool were used for physiochemical and functional characterization of the protein. Studied of secondary structure of that protein was carried out by computational program, ProFunc. The structure is finally submitted in Protein Model Database. The predicted model permits initial inferences about the unexplored 3D structure of the CMP kinase and may be promote in relational designing of molecules for structure-function studies.     

Crystal Structures of the CERT START Domain with Inhibitors Provide Insights into the Mechanism of Ceramide Transfer

The cytosolic protein CERT transfers ceramide from the endoplasmic reticulum to the Golgi apparatus where ceramide is converted to SM. The C-terminal START (steroidogenic acute regulatory protein-related lipid transfer) domain of CERT binds one ceramide molecule in its central amphiphilic cavity. (1R,3R)-N-(3-Hydroxy-1-hydroxymethyl-3-phenylpropyl)alkanamide (HPA), a synthesized analogue of ceramide, inhibits ceramide transfer by CERT. Here we report crystal structures of the CERT START domain in complex with HPAs of varying acyl chain lengths. In these structures, one HPA molecule is buried in the amphiphilic cavity where the amide and hydroxyl groups of HPA form a hydrogen-bond network with specific amino acid residues. The Ω1 loop, which has been suggested to function as a gate of the cavity, adopts a different conformation when bound to HPA than when bound to ceramide. In the Ω1 loop region, Trp473 shows the largest difference between these two structures. This residue exists inside of the cavity in HPA-bound structures, while it is exposed to the outside of the protein in the apo-form and ceramide-bound complex structures. Surface plasmon resonance experiments confirmed that Trp473 is important for interaction with membranes. These results provide insights into not only the molecular mechanism of inhibition by HPAs but also possible mechanisms by which CERT interacts with ceramide.     

Stabilization of yeast hexokinase A by polyol osmolytes: correlation with the physicochemical properties of aqueous solutions

Osmolytes of the polyol series are known to accumulate in biological systems under stress and stabilize the structures of a wide variety of proteins. While increased surface tension of aqueous solutions has been considered an important factor in protein stabilization effect, glycerol is an exception, lowering the surface tension of water. To clarify this anomalous effect, the effect of a series of polyols on the thermal stability of a highly thermolabile two domain protein yeast hexokinase A has been investigated by differential scanning calorimetry and by monitoring loss in the biological activity of the enzyme as a function of time. A larger increase in the T(m) of domain 1 compared with that of domain 2, varying linearly with the number of hydroxyl groups in polyols, has been observed, sorbitol being the best stabilizer against both thermal as well as urea denaturation. Polyols help retain the activity of the enzyme considerably and a good correlation of the increase in T(m) (DeltaT(m)) and the retention of activity with the increase in the surface tension of polyol solutions, except glycerol, which breaks this trend, has been observed. However, the DeltaT(m) values show a linear correlation with apparent molal heat capacity and volume of aqueous polyol solutions including glycerol. These results suggest that while bulk solution properties contribute significantly to protein stabilization, interfacial properties are not always a good indicator of the stabilizing effect. A subtle balance of various weak binding and exclusion effects of the osmolytes mediated by water further regulates the stabilizing effect. Understanding these aspects is critical in the rational design of stable protein formulations.     

Molecular cloning and polymorphism analysis of the prion protein gene in Tan sheep of Ningxia, China

The resistance or susceptibility of sheep to scrapie is associated with polymorphisms of the prion protein gene (PRNP), particularly, single nucleotide polymorphisms (SNPs) in amino acid positions 136, 154 and 171. The prion protein (PrP) gene sequence and the deduced amino acid alignment of prion protein in Tan sheep, a local Chinese sheep breed traditionally raised in Ningxia, northwestern China, were determined and variability of the PrP amino acids sequence was analyzed in this study. The PrP nucleic acids and amino acids sequences of 112 Tan sheep were highly homogenous, although polymorphism of the PrP gene was detected at several sites, particularly codons 106, 154, and 171. The analysis of both sequences revealed that the most predominant allele at codons 136, 154 and 171 in Tan sheep was ARQ, which was known to be associated with high susceptibility to scrapie in sheep. The result suggests that Tan sheep is potentially susceptible to scrapie. Our findings provide valuable information for future breeding projects to scrapie resistance in Tan sheep.     

Structural basis for bile salt inhibition of pancreatic phospholipase A2

Bile salt interactions with phospholipid monolayers of fat emulsions are known to regulate the actions of gastrointestinal lipolytic enzymes in order to control the uptake of dietary fat. Specifically, on the lipid/aqueous interface of fat emulsions, the anionic portions of amphipathic bile salts have been thought to interact with and activate the enzyme group-IB phospholipase A2 (PLA2) derived from the pancreas. To explore this regulatory process, we have determined the crystal structures of the complexes of pancreatic PLA2 with the naturally occurring bile salts: cholate, glycocholate, taurocholate, glycochenodeoxycholate, and taurochenodeoxycholate. The five PLA2-bile salt complexes each result in a partly occluded active site, and the resulting ligand binding displays specific hydrogen bonding interactions and extensive hydrophobic packing. The amphipathic bile salts are bound to PLA2 with their polar hydroxyl and sulfate/carboxy groups oriented away from the enzyme's hydrophobic core. The impaired catalytic and interface binding functions implied by these structures provide a basis for the previous numerous observations of a biphasic dependence of the rate of PLA2 catalyzed hydrolysis of zwitterionic glycerophospholipids in the presence of bile salts. The rising or activation phase is consistent with enhanced binding and activation of the bound PLA2 by the bile salt induced anionic charge in a zwitterionic interface. The falling or inhibitory phase can be explained by the formation of a catalytically inert stoichiometric complex between PLA2 and any bile salts in which it forms a stable complex. The model provides new insight into the regulatory role that specific PLA2-bile salt interactions are likely to play in fat metabolism.     

Structural and functional characterization of the Geobacillus copper nitrite reductase: Involvement of the unique N-terminal region in the interprotein electron transfer with its redox partner

The crystal structures of copper-containing nitrite reductase (CuNiR) from the thermophilic Gram-positive bacterium Geobacillus kaustophilus HTA426 and the amino (N)-terminal 68 residue-deleted mutant were determined at resolutions of 1.3 Å and 1.8 Å, respectively. Both structures show a striking resemblance with the overall structure of the well-known CuNiRs composed of two Greek key β-barrel domains; however, a remarkable structural difference was found in the N-terminal region. The unique region has one β-strand and one α-helix extended to the northern surface of the type-1 copper site. The superposition of the Geobacillus CuNiR model on the electron-transfer complex structure of CuNiR with the redox partner cytochrome c₅₅₁ in other denitrifier system led us to infer that this region contributes to the transient binding with the partner protein during the interprotein electron transfer reaction in the Geobacillus system. Furthermore, electron-transfer kinetics experiments using N-terminal residue-deleted mutant and the redox partner, Geobacillus cytochrome c₅₅₁, were carried out. These structural and kinetics studies demonstrate that the region is directly involved in the specific partner recognition.     

Applying an empirical hydropathic forcefield in refinement may improve low-resolution protein X-ray crystal structures

Background

The quality of X-ray crystallographic models for biomacromolecules refined from data obtained at high-resolution is assured by the data itself. However, at low-resolution, >3.0 Å, additional information is supplied by a forcefield coupled with an associated refinement protocol. These resulting structures are often of lower quality and thus unsuitable for downstream activities like structure-based drug discovery.

Methodology

An X-ray crystallography refinement protocol that enhances standard methodology by incorporating energy terms from the HINT (Hydropathic INTeractions) empirical forcefield is described. This protocol was tested by refining synthetic low-resolution structural data derived from 25 diverse high-resolution structures, and referencing the resulting models to these structures. The models were also evaluated with global structural quality metrics, e.g., Ramachandran score and MolProbity clashscore. Three additional structures, for which only low-resolution data are available, were also re-refined with this methodology.

Results

The enhanced refinement protocol is most beneficial for reflection data at resolutions of 3.0 Å or worse. At the low-resolution limit, ≥4.0 Å, the new protocol generated models with Cα positions that have RMSDs that are 0.18 Å more similar to the reference high-resolution structure, Ramachandran scores improved by 13%, and clashscores improved by 51%, all in comparison to models generated with the standard refinement protocol. The hydropathic forcefield terms are at least as effective as Coulombic electrostatic terms in maintaining polar interaction networks, and significantly more effective in maintaining hydrophobic networks, as synthetic resolution is decremented. Even at resolutions ≥4.0 Å, these latter networks are generally native-like, as measured with a hydropathic interactions scoring tool.     

Three-dimensional (3D) structure prediction of the American and African oil-palms β-ketoacyl-[ACP] synthase-II protein by comparative modelling

Background: The fatty-acid profile of the vegetable oils determines its properties and nutritional value. Palm-oil obtained from the African oil-palm [Elaeis guineensis Jacq. (Tenera)] contains 44% palmitic acid (C_16:0), but, palm-oil obtained from the American oilpalm [Elaeis oleifera] contains only 25% C16:0. In part, the b-ketoacyl-[ACP] synthase II (KASII) [EC: 2.3.1.179] protein is responsible for the high level of C_16:0 in palm-oil derived from the African oil-palm. To understand more about E. guineensis KASII (EgKASII) and E. oleifera KASII (EoKASII) proteins, it is essential to know its structures. Hence, this study was undertaken. Objective: The objective of this study was to predict three-dimensional (3D) structure of EgKASII and EoKASII proteins using molecular modelling tools. Materials and Methods: The amino-acid sequences for KASII proteins were retrieved from the protein database of National Center for Biotechnology Information (NCBI), USA. The 3D structures were predicted for both proteins using homology modelling and ab-initio technique approach of protein structure prediction. The molecular dynamics (MD) simulation was performed to refine the predicted structures. The predicted structure models were evaluated and root mean square deviation (RMSD) and root mean square fluctuation (RMSF) values were calculated. Results: The homology modelling showed that EgKASII and EoKASII proteins are 78% and 74% similar with Streptococcus pneumonia KASII and Brucella melitensis KASII, respectively. The EgKASII and EoKASII structures predicted by using ab-initio technique approach shows 6% and 9% deviation to its structures predicted by homology modelling, respectively. The structure refinement and validation confirmed that the predicted structures are accurate. Conclusion: The 3D structures for EgKASII and EoKASII proteins were predicted. However, further research is essential to understand the interaction of EgKASII and EoKASII proteins with its substrates.