Fine tuning of the catalytic activity of colicin E7 nuclease domain by systematic N‐terminal mutations

The nuclease domain of colicin E7 (NColE7) promotes the nonspecific cleavage of nucleic acids at its C‐terminal HNH motif. Interestingly, the deletion of four N‐terminal residues (446–449 NColE7 = KRNK) resulted in complete loss of the enzyme activity. R447A mutation was reported to decrease the nuclease activity, but a detailed analysis of the role of the highly positive and flexible N‐terminus is still missing. Here, we present the study of four mutants, with a decreased activity in the following order: NColE7  >> KGNK > KGNG ～ GGNK > GGNG. At the same time, the folding, the metal‐ion, and the DNA‐binding affinity were unaffected by the mutations as revealed by linear and circular dichroism spectroscopy, isothermal calorimetric titrations, and gel mobility shift experiments. Semiempirical quantum chemical calculations and molecular dynamics simulations revealed that K446, K449, and/or the N‐terminal amino group are able to approach the active centre in the absence of the other positively charged residues. The results suggested a complex role of the N‐terminus in the catalytic process that could be exploited in the design of a controlled nuclease.     

The role of the N-terminal loop in the function of the colicin E7 nuclease domain

Colicin E7 (ColE7) is a metallonuclease toxin of Escherichia coli belonging to the HNH superfamily of nucleases. It contains highly conserved amino acids in its HHX₁₄NX₈HX₃H ββα-type metal ion binding C-terminal active centre. However, the proximity of the arginine at the N-terminus of the nuclease domain of ColE7 (NColE7, 446–576) is necessary for the hydrolytic activity. This poses a possibility of allosteric activation control in this protein. To obtain more information on this phenomenon, two protein mutants were expressed, i.e. four and 25 N-terminal amino acids were removed from NColE7. The effect of the N-terminal truncation on the Zn²⁺ ion and DNA binding as well as on the activity was investigated in this study by mass spectrometry, synchrotron-radiation circular dichroism and fluorescence spectroscopy and agarose gel mobility shift assays. The dynamics of protein backbone movement was simulated by molecular dynamics. Semiempirical quantum chemical calculations were performed to obtain better insight into the structure of the active centre. The longer protein interacted with both Zn²⁺ ion and DNA more strongly than its shorter counterpart. The results were explained by the structural stabilization effect of the N-terminal amino acids on the catalytic centre. In agreement with this, the absence of the N-terminal sequences resulted in significantly increased movement of the backbone atoms compared with that in the native NColE7: in ΔN25-NColE7 the amino acid strings between residues 485–487, 511–515 and 570–571, and in ΔN4-NColE7 those between residues 467–468, 530–535 and 570–571.     

LALF32‐51‐E7, a HPV‐16 therapeutic vaccine candidate,forms protein body‐like structures when expressed in Nicotiana benthamiana leaves

High‐risk human papillomaviruses (HPVs) cause cervical cancer, and while there are good prophylactic vaccines on the market, these are ineffective against established infections, creating a clear need for therapeutic vaccines. The HPV E7 protein is one of the essential oncoproteins for the onset and maintenance of malignancy and is therefore an ideal therapeutic vaccine target. We fused the HPV‐16 E7 protein to the Limulus polyphemus antilipopolysaccharide factor (LALF_32‐51), a small hydrophobic peptide that can penetrate cell membranes and that has immunomodulatory properties. LALF_32‐51‐E7 was transiently expressed in Nicotiana benthamiana, and we previously determined that it accumulated better when targeted to chloroplasts compared to being localized in the cytoplasm. Subsequently, we aimed to prove whether LALF_32‐51‐E7 was indeed associated with the chloroplasts by determining its subcellular localization. The LALF_32‐51‐E7 gene was fused to one encoding enhanced GFP to generate a LG fusion protein, and localization was determined by confocal laser scanning microscopy and transmission electron microscopy (TEM). The fluorescence observed from chloroplast‐targeted LG was distinctively different from that of the cytoplasmic LG. Small spherical structures resembling protein bodies (PBs) were seen that clearly localized with the chloroplasts. Larger but less abundant PB‐like structures were also seen for the cytoplasmic LG. PB‐like structure formation was confirmed for both LG and LALF_32‐51‐E7 by TEM. LALF_32‐51‐E7 was indeed targeted to the chloroplasts by the chloroplast transit peptide used in this study, and it formed aggregated PB‐like structures. This study could open a new avenue for the use of LALF_32‐51 as a PB‐inducing peptide.     

Accuracy issues involved in modeling in vivo protein structures using PM7

Using the semiempirical method PM7, an attempt has been made to quantify the error in prediction of the in vivo structure of proteins relative to X‐ray structures. Three important contributory factors are the experimental limitations of X‐ray structures, the difference between the crystal and solution environments, and the errors due to PM7. The geometries of 19 proteins from the Protein Data Bank that had small R values, that is, high accuracy structures, were optimized and the resulting drop in heat of formation was calculated. Analysis of the changes showed that about 10% of this decrease in heat of formation was caused by faults in PM7, the balance being attributable to the X‐ray structure and the difference between the crystal and solution environments. A previously unknown fault in PM7 was revealed during tests to validate the geometries generated using PM7. Clashscores generated by the Molprobity molecular mechanics structure validation program showed that PM7 was predicting unrealistically close contacts between nonbonding atoms in regions where the local geometry is dominated by very weak noncovalent interactions. The origin of this fault was traced to an underestimation of the core‐core repulsion between atoms at distances smaller than the equilibrium distance. Proteins 2015; 83:1427–1435. © 2015 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Physical and functional interaction between hepatitis C virus NS5A protein and ovarian tumor protein deubiquitinase 7B

Hepatitis C virus (HCV) NS5A protein plays crucial roles in viral RNA replication, virus assembly, and viral pathogenesis. Although NS5A has no known enzymatic activity, it modulates various cellular pathways through interaction with cellular proteins. HCV NS5A (and other HCV proteins) are reportedly degraded through the ubiquitin–proteasome pathway; however, the physiological roles of ubiquitylation and deubiquitylation in HCV infection are largely unknown. To elucidate the role of deubiquitylation in HCV infection, an attempt was made to identify a deubiquitinase (DUB) that can interact with NS5A protein. An ovarian tumor protein (OTU), deubiquitinase 7B (OTUD7B), was identified as a novel NS5A‐binding protein. Co‐immunoprecipitation analyses showed that NS5A interacts with OTUD7B in both Huh‐7 and HCV RNA replicon cells. Immunofluorescence staining revealed that HCV NS5A protein colocalizes with OTUD7B in the cytoplasm. Moreover, HCV infection was found to enhance the nuclear localization of OTUD7B. The OTUD7B‐binding domain on NS5A was mapped using a series of NS5A deletion mutants. The present findings suggest that the domain I of NS5A is important and the region from amino acid 121 to 126 of NS5A essential for the interaction. Either V121A or V124A mutation in NS5A disrupts the NS5A‐OTUD7B interaction. The results of this in vivo ubiquitylation assay suggest that HCV NS5A enhances OTUD7B DUB activity. Taken together, these results suggest that HCV NS5A protein interacts with OTUD7B, thereby modulating its DUB activity.     

Conformational dynamics is more important than helical propensity for the folding of the all α‐helical protein Im7

Im7 folds via an on‐pathway intermediate that contains three of the four native α‐helices. The missing helix, helix III, is the shortest and its failure to be formed until late in the pathway is related to frustration in the structure. Im7H3M3, a 94‐residue variant of the 87‐residue Im7 in which helix III is the longest of the four native helices, also folds via an intermediate. To investigate the structural basis for this we calculated the frustration in the structure of Im7H3M3 and used NMR to investigate its dynamics. We found that the native state of Im7H3M3 is highly frustrated and in equilibrium with an intermediate state that lacks helix III, similar to Im7. Model‐free analysis identified residues with chemical exchange contributions to their relaxation that aligned with the residues predicted to have highly frustrated interactions, also like Im7. Finally, we determined properties of urea‐denatured Im7H3M3 and identified four clusters of interacting residues that corresponded to the α‐helices of the native protein. In Im7 the cluster sizes were related to the lengths of the α‐helices with cluster III being the smallest but in Im7H3M3 cluster III was also the smallest, despite this region forming the longest helix in the native state. These results suggest that the conformational properties of the urea‐denatured states promote formation of a three‐helix intermediate in which the residues that form helix III remain non‐helical. Thus it appears that features of the native structure are formed early in folding linked to collapse of the unfolded state.     

Structures of designed armadillo repeat proteins binding to peptides fused to globular domains

Designed armadillo repeat proteins (dArmRP) are α‐helical solenoid repeat proteins with an extended peptide binding groove that were engineered to develop a generic modular technology for peptide recognition. In this context, the term “peptide” not only denotes a short unstructured chain of amino acids, but also an unstructured region of a protein, as they occur in termini, loops, or linkers between folded domains. Here we report two crystal structures of dArmRPs, in complex with peptides fused either to the N‐terminus of Green Fluorescent Protein or to the C‐terminus of a phage lambda protein D. These structures demonstrate that dArmRPs bind unfolded peptides in the intended conformation also when they constitute unstructured parts of folded proteins, which greatly expands possible applications of the dArmRP technology. Nonetheless, the structures do not fully reflect the binding behavior in solution, that is, some binding sites remain unoccupied in the crystal and even unexpected peptide residues appear to be bound. We show how these differences can be explained by restrictions of the crystal lattice or the composition of the crystallization solution. This illustrates that crystal structures have to be interpreted with caution when protein–peptide interactions are characterized, and should always be correlated with measurements in solution.     

Direct Optical Microscopic Observation of the Microtubule Polymerization Intermediate Sheet Structure in the Presence of Gas7

The process of microtubule elongation is thought to consist of two stages—formation of a tubulin sheet structure and its closure into a tube. However, real-time observation of this process has been difficult. Here, by utilizing phospho-tau binding protein Gas7 (growth-arrest-specific protein 7), we visualized the polymer transformation process by dark-field microscopy. Upon elongation, thin and flexible structures, often similar to a curved hook, appeared at the end of microtubules. Electron microscopic observations supported the idea that these flexible structures are tubulin sheets. They maintained their length until they gradually became thick and rigid beginning in the central portion, resulting in straight microtubules. In the absence of Gas7, the sheet-like structure was rarely observed; moreover, when observed, it was fragile and engaged in typical dynamic instability. With Gas7, no catastrophe was observed. These results suggest that Gas7 enhances microtubule polymerization by stabilizing sheet intermediates and is a useful tool for analyzing microtubule transformation.     

Multiple solvent crystal structures of ribonuclease A: An assessment of the method

The multiple solvent crystal structures (MSCS) method uses organic solvents to map the surfaces of proteins. It identifies binding sites and allows for a more thorough examination of protein plasticity and hydration than could be achieved by a single structure. The crystal structures of bovine pancreatic ribonuclease A (RNAse A) soaked in the following organic solvents are presented: 50% dioxane, 50% dimethylformamide, 70% dimethylsulfoxide, 70% 1,6‐hexanediol, 70% isopropanol, 50% R,S,R‐bisfuran alcohol, 70% t‐butanol, 50% trifluoroethanol, or 1.0M trimethylamine‐N‐oxide. This set of structures is compared with four sets of crystal structures of RNAse A from the protein data bank (PDB) and with the solution NMR structure to assess the validity of previously untested assumptions associated with MSCS analysis. Plasticity from MSCS is the same as from PDB structures obtained in the same crystal form and deviates only at crystal contacts when compared to structures from a diverse set of crystal environments. Furthermore, there is a good correlation between plasticity as observed by MSCS and the dynamic regions seen by NMR. Conserved water binding sites are identified by MSCS to be those that are conserved in the sets of structures taken from the PDB. Comparison of the MSCS structures with inhibitor‐bound crystal structures of RNAse A reveals that the organic solvent molecules identify key interactions made by inhibitor molecules, highlighting ligand binding hot‐spots in the active site. The present work firmly establishes the relevance of information obtained by MSCS. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Design and structure of two new protein cages illustrate successes and ongoing challenges in protein engineering

In recent years, new protein engineering methods have produced more than a dozen symmetric, self‐assembling protein cages whose structures have been validated to match their design models with near‐atomic accuracy. However, many protein cage designs that are tested in the lab do not form the desired assembly, and improving the success rate of design has been a point of recent emphasis. Here we present two protein structures solved by X‐ray crystallography of designed protein oligomers that form two‐component cages with tetrahedral symmetry. To improve on the past tendency toward poorly soluble protein, we used a computational protocol that favors the formation of hydrogen‐bonding networks over exclusively hydrophobic interactions to stabilize the designed protein–protein interfaces. Preliminary characterization showed highly soluble expression, and solution studies indicated successful cage formation by both designed proteins. For one of the designs, a crystal structure confirmed at high resolution that the intended tetrahedral cage was formed, though several flipped amino acid side chain rotamers resulted in an interface that deviates from the precise hydrogen‐bonding pattern that was intended. A structure of the other designed cage showed that, under the conditions where crystals were obtained, a noncage structure was formed wherein a porous 3D protein network in space group I2₁3 is generated by an off‐target twofold homomeric interface. These results illustrate some of the ongoing challenges of developing computational methods for polar interface design, and add two potentially valuable new entries to the growing list of engineered protein materials for downstream applications.     

Solution structures of polcalcin Phl p 7 in three ligation states: Apo‐, hemi‐Mg2+‐bound,and fully Ca2+‐bound

Polcalcins are small EF‐hand proteins believed to assist in regulating pollen‐tube growth. Phl p 7, from timothy grass (Phleum pratense), crystallizes as a domain‐swapped dimer at low pH. This study describes the solution structures of the recombinant protein in buffered saline at pH 6.0, containing either 5.0 mM EDTA, 5.0 mM Mg²⁺, or 100 μM Ca²⁺. Phl p 7 is monomeric in all three ligation states. In the apo‐form, both EF‐hand motifs reside in the closed conformation, with roughly antiparallel N‐ and C‐terminal helical segments. In 5.0 mM Mg²⁺, the divalent ion is bound by EF‐hand 2, perturbing interhelical angles and imposing more regular helical structure. The structure of Ca²⁺‐bound Phl p 7 resembles that previously reported for Bet v 4—likewise exposing apolar surface to the solvent. Occluded in the apo‐ and Mg²⁺‐bound forms, this surface presumably provides the docking site for Phl p 7 targets. Unlike Bet v 4, EF‐hand 2 in Phl p 7 includes five potential anionic ligands, due to replacement of the consensus serine residue at –x (residue 55 in Phl p 7) with aspartate. In the Phl p 7 crystal structure, D55 functions as a helix cap for helix D. In solution, however, D55 apparently serves as a ligand to the bound Ca²⁺. When Mg²⁺ resides in site 2, the D55 carboxylate withdraws to a distance consistent with a role as an outer‐sphere ligand. ¹⁵N relaxation data, collected at 600 MHz, indicate that backbone mobility is limited in all three ligation states. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Representing and comparing protein folds and fold families using three‐dimensional shape‐density representations

The question of how best to compare and classify the (three‐dimensional) structures of proteins is one of the most important unsolved problems in computational biology. To help tackle this problem, we have developed a novel shape‐density superposition algorithm called 3D‐Blast which represents and superposes the shapes of protein backbone folds using the spherical polar Fourier correlation technique originally developed by us for protein docking. The utility of this approach is compared with several well‐known protein structure alignment algorithms using receiver‐operator‐characteristic plots of queries against the “gold standard” CATH database. Despite being completely independent of protein sequences and using no information about the internal geometry of proteins, our results from searching the CATH database show that 3D‐Blast is highly competitive compared to current state‐of‐the‐art protein structure alignment algorithms. A novel and potentially very useful feature of our approach is that it allows an average or “consensus” fold to be calculated easily for a given group of protein structures. We find that using consensus shapes to represent entire fold families also gives very good database query performance. We propose that using the notion of consensus fold shapes could provide a powerful new way to index existing protein structure databases, and that it offers an objective way to cluster and classify all of the currently known folds in the protein universe. Proteins 2012. © 2011 Wiley Periodicals, Inc.     

Folding of Top7 in unbiased all‐atom Monte Carlo simulations

For computational studies of protein folding, proteins with both helical and β‐sheet secondary structure elements are very challenging, as they expose subtle biases of the physical models. Here, we present reproducible folding of a 92 residue α/β protein (residues 3–94 of Top7, PDB ID: 1QYS) in computer simulations starting from random initial conformations using a transferable physical model which has been previously shown to describe the folding and thermodynamic properties of about 20 other smaller proteins of different folds. Top7 is a de novo designed protein with two α‐helices and a five stranded β‐sheet. Experimentally, it is known to be unusually stable for its size, and its folding transition distinctly deviates from the two‐state behavior commonly seen in natural single domain proteins. In our all‐atom implicit solvent parallel tempering Monte Carlo simulations, Top7 shows a rapid transition to a group of states with high native‐like secondary structure, and a much slower subsequent transition to the native state with a root mean square deviation of about 3.5 Å from the experimentally determined structure. Consistent with experiments, we find Top7 to be thermally extremely stable, although the simulations also find a large number of very stable non‐native states with high native‐like secondary structure. Proteins 2013; 81:1446–1456. © 2013 Wiley Periodicals, Inc.     

The role of BMP‐7 in chondrogenic and osteogenic differentiation of human bone marrow multipotent mesenchymal stromal cells in vitro

This study addresses the role of bone morphogenetic protein‐7 (BMP‐7) in chondrogenic and osteogenic differentiation of human bone marrow multipotent mesenchymal stromal cells (BM MSCs) in vitro. BM MSCs were expanded and differentiated in the presence or absence of BMP‐7 in monolayer and three‐dimensional cultures. After 3 days of stimulation, BMP‐7 significantly inhibited MSC growth in expansion cultures. When supplemented in commonly used induction media for 7–21 days, BMP‐7 facilitated both chondrogenic and osteogenic differentiation of MSCs. This was evident by specific gene and protein expression analyses using real‐time PCR, Western blot, histological, and immunohistochemical staining. BMP‐7 supplementation appeared to enhance upregulation of lineage‐specific markers, such as type II and type IX collagens (COL2A1, COL9A1) in chondrogenic and secreted phosphoprotein 1 (SPP1), osteocalcin (BGLAP), and osterix (SP7) in osteogenic differentiation. BMP‐7 in the presence of TGF‐β3 induced superior chondrocytic proteoglycan accumulation, type II collagen, and SOX9 protein expression in alginate and pellet cultures compared to either factor alone. BMP‐7 increased alkaline phosphatase activity and dose‐dependently accelerated calcium mineralization of osteogenic differentiated MSCs. The potential of BMP‐7 to promote adipogenesis of MSCs was restricted under osteogenic conditions, despite upregulation of adipocyte gene expression. These data suggest that BMP‐7 is not a singular lineage determinant, rather it promotes both chondrogenic and osteogenic differentiation of MSCs by co‐ordinating with initial lineage‐specific signals to accelerate cell fate determination. BMP‐7 may be a useful enhancer of in vitro differentiation of BM MSCs for cell‐based tissue repair. J. Cell. Biochem. 109: 406–416, 2010. © 2009 Wiley‐Liss, Inc.     

First structural evidence of sequestration of mRNA cap structures by type 1 ribosome inactivating protein from Momordica balsamina

This is the first structural evidence of recognition of mRNA cap structures by a ribosome inactivating protein. It is well known that a unique cap structure is formed at the 5′ end of mRNA for carrying out various processes including mRNA maturation, translation initiation, and RNA turnover. The binding studies and crystal structure determinations of type 1 ribosome inactivating protein (RIP‐1) from Momordica balsamina (MbRIP‐1) were carried out with mRNA cap structures including (i) N7‐methyl guanine (m7G), (ii) N7‐methyl guanosine diphosphate (m7GDP), and (iii) N7‐methyl guanosine triphosphate (m7GTP). These compounds showed affinities to MbRIP‐1 at nanomolar concentrations. The structure determinations of the complexes of MbRIP‐1 with m7G, m7GDP, and m7GTP at 2.65, 1.77, and 1.75 Å resolutions revealed that all the three compounds bound to MbRIP‐1 in the substrate binding site at the positions which are slightly shifted towards Glu85 as compared to those of rRNA substrates. In this position, Glu85 forms several hydrogen bonds with guanine moiety while N‐7 methyl group forms van der Waals contacts. However, the guanine rings are poorly stacked in these complexes. Thus, the mode of binding by MbRIP‐1 to mRNA cap structures is different which results in the inhibition of depurination. Since some viruses are known to exploit the capping property of the host, this action of MbRIP‐1 may have implications for the antiviral activity of this protein in vivo. The understanding of the mode of binding of MbRIP‐1 to cap structures may also assist in the design of anti‐viral agents. Proteins 2012. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Comparison of NMR and crystal structures of membrane proteins and computational refinement to improve model quality

Membrane proteins are challenging to study and restraints for structure determination are typically sparse or of low resolution because the membrane environment that surrounds them leads to a variety of experimental challenges. When membrane protein structures are determined by different techniques in different environments, a natural question is “which structure is most biologically relevant?” Towards answering this question, we compiled a dataset of membrane proteins with known structures determined by both solution NMR and X‐ray crystallography. By investigating differences between the structures, we found that RMSDs between crystal and NMR structures are below 5 Å in the membrane region, NMR ensembles have a higher convergence in the membrane region, crystal structures typically have a straighter transmembrane region, have higher stereo‐chemical correctness, and are more tightly packed. After quantifying these differences, we used high‐resolution refinement of the NMR structures to mitigate them, which paves the way for identifying and improving the structural quality of membrane proteins.     

The role of BRD7 in embryo development and glucose metabolism

Bromodomain‐containing protein 7 (BRD7) is a member of bromodomain‐containing protein family and its function has been implicated in several diseases. We have previously shown that BRD7 plays a role in metabolic processes. However, the effect of BRD7 deficiency in glucose metabolism and its role in in vivo have not been fully revealed. Here, we report the essential role of BRD7 during embryo development. Mice homozygous for BRD7 led to embryonic lethality at mid‐gestation. Homozygous BRD7 knockout (KO) mice showed retardation in development, and eventually all BRD7 KO embryos died in utero prior to E16.5. Partial knockdown of Brd7 gene displayed mild changes in glucose metabolism.     

Protein models docking benchmark 2

Structural characterization of protein–protein interactions is essential for our ability to understand life processes. However, only a fraction of known proteins have experimentally determined structures. Such structures provide templates for modeling of a large part of the proteome, where individual proteins can be docked by template‐free or template‐based techniques. Still, the sensitivity of the docking methods to the inherent inaccuracies of protein models, as opposed to the experimentally determined high‐resolution structures, remains largely untested, primarily due to the absence of appropriate benchmark set(s). Structures in such a set should have predefined inaccuracy levels and, at the same time, resemble actual protein models in terms of structural motifs/packing. The set should also be large enough to ensure statistical reliability of the benchmarking results. We present a major update of the previously developed benchmark set of protein models. For each interactor, six models were generated with the model‐to‐native C^α RMSD in the 1 to 6 Å range. The models in the set were generated by a new approach, which corresponds to the actual modeling of new protein structures in the “real case scenario,” as opposed to the previous set, where a significant number of structures were model‐like only. In addition, the larger number of complexes (165 vs. 63 in the previous set) increases the statistical reliability of the benchmarking. We estimated the highest accuracy of the predicted complexes (according to CAPRI criteria), which can be attained using the benchmark structures. The set is available at http://dockground.bioinformatics.ku.edu . Proteins 2015; 83:891–897. © 2015 Wiley Periodicals, Inc.