Exploring additivity effects of double mutations on the binding affinity of protein‐protein complexes

Additivity in binding affinity of protein‐protein complexes refers to the change in free energy of binding (ΔΔG_bind) for double (or multiple) mutations which is approximately equal to the sum of their corresponding single mutation ΔΔG_bind values. In this study, we have explored the additivity effect of double mutants, which shows a linear relationship between the binding affinity of double and sum of single mutants with a correlation of 0.90. However, the comparison of ΔΔG_bind values showed a mean absolute deviation of 0.86 kcal/mol, and 25.6% of the double mutants show a deviation of more than 1 kcal/mol, which are identified as non‐additive. The additivity effects have been analyzed based on the influence of structural features such as accessible surface area, long range order, binding propensity change, surrounding hydrophobicity, flexibility, atomic contacts between the mutations and distance between the 2 mutations. We found that non‐additive mutations tend to be closer to each other and have more contacts. We have also used machine learning methods to discriminate additive and non‐additive mutations using structure‐based features, which showed the accuracies in the range of 0.77–0.92 for protein‐protein complexes belonging to different functions. Further, we have compared the additivity effects of protein stability along with binding affinity and explored the similarities and differences between them. The results obtained in this study provide insights into the effects of various structural features on binding affinity of double mutants, and will aid the development of accurate methods to predict the binding affinity of double mutants.     

Euclidean sections of protein conformation space and their implications in dimensionality reduction

Dimensionality reduction is widely used in searching for the intrinsic reaction coordinates for protein conformational changes. We find the dimensionality?reduction methods using the pairwise root?mean?square deviation (RMSD) as the local distance metric face a challenge. We use Isomap as an example to illustrate the problem. We believe that there is an implied assumption for the dimensionality‐reduction approaches that aim to preserve the geometric relations between the objects: both the original space and the reduced space have the same kind of geometry, such as Euclidean geometry vs. Euclidean geometry or spherical geometry vs. spherical geometry. When the protein free energy landscape is mapped onto a 2D plane or 3D space, the reduced space is Euclidean, thus the original space should also be Euclidean. For a protein with N atoms, its conformation space is a subset of the 3N‐dimensional Euclidean space R^3N. We formally define the protein conformation space as the quotient space of R^3N by the equivalence relation of rigid motions. Whether the quotient space is Euclidean or not depends on how it is parameterized. When the pairwise RMSD is employed as the local distance metric, implicit representations are used for the protein conformation space, leading to no direct correspondence to a Euclidean set. We have demonstrated that an explicit Euclidean‐based representation of protein conformation space and the local distance metric associated to it improve the quality of dimensionality reduction in the tetra‐peptide and β‐hairpin systems. Proteins 2014; 82:2585–2596. © 2014 Wiley Periodicals, Inc.     

Effects of ligand binding on the stability of aldo–keto reductases: Implications for stabilizer or destabilizer chaperones

Ligands such as enzyme inhibitors stabilize the native conformation of a protein upon binding to the native state, but some compounds destabilize the native conformation upon binding to the non‐native state. The former ligands are termed “stabilizer chaperones” and the latter ones “destabilizer chaperones.” Because the stabilization effects are essential for the medical chaperone (MC) hypothesis, here we have formulated a thermodynamic system consisting of a ligand and a protein in its native‐ and non‐native state. Using the differential scanning fluorimetry and the circular dichroism varying the urea concentration and temperature, we found that when the coenzyme NADP⁺ was absent, inhibitors such as isolithocholic acid stabilized the aldo–keto reductase AKR1A1 upon binding, which showed actually the three‐state folding, but destabilized AKR1B10. In contrast, in the presence of NADP⁺, they destabilized AKR1A1 and stabilized AKR1B10. To explain these phenomena, we decomposed the free energy of stabilization (ΔΔG) into its enthalpy (ΔΔH) and entropy (ΔΔS) components. Then we found that in a relatively unstable protein showing the three‐state folding, native conformation was stabilized by the negative ΔΔH in association with the negative ΔΔS, suggesting that the stabilizer chaperon decreases the conformational fluctuation of the target protein or increase its hydration. However, in other cases, ΔΔG was essentially determined by the delicate balance between ΔΔH and ΔΔS. The proposed thermodynamic formalism is applicable to the system including multiple ligands with allosteric interactions. These findings would promote the development of screening strategies for MCs to regulate the target conformations.     

Comparative proteomic analyses reveal that the regulators of G‐protein signaling proteins regulate amino acid metabolism of the rice blast fungus Magnaporthe oryzae

The rice blast fungus Magnaporthe oryzae encodes eight regulators of G‐protein (GTP‐binding protein) signaling (RGS) proteins MoRgs1–MoRgs8 that orchestrate the growth, asexual/sexual production, appressorium differentiation, and pathogenicity. To address the mechanisms by which MoRgs proteins function, we conducted a 2DE proteome study and identified 82 differentially expressed proteins by comparing five ?Morgs mutants with wild‐type Guy11 strain. We found that the abundances of eight amino acid (AA) biosynthesis or degradation associated proteins were markedly altered in five ?Morgs mutants, indicating one of the main collective roles for the MoRgs proteins is to influence AA metabolism. We showed that MoRgs proteins have distinct roles in AA metabolism and nutrient responses from growth assays. In addition, we characterized MoLys20 (Lys is lysine), a homocitrate synthase, whose abundance was significantly decreased in the ?Morgs mutants. The ?Molys20 mutant is auxotrophic for lys and exogenous lys could partially rescue its auxotrophic defects. Deletion of MoLYS20 resulted in defects in conidiation and infection, as well as pathogenicity on rice. Overall, our results indicate that one of the critical roles for MoRgs proteins is to regulate AA metabolism, and that MoLys20 may be directly or indirectly regulated by MoRgs and participated in lys biosynthesis, thereby affecting fungal development and pathogenicity.     

Biochemical characteristics of guanine nucleotide binding protein alpha-subunit recombinant protein and three mutants: investigation of a domain motion involved in GDP-GTP exchange

Our previous studies using molecular dynamics have shown a hinge bending motion between the helical and the GTPase domains of G_T (Mello et al., 1998). The hypothesis that this motion is allowed by residues Gly56 and Gly179 and that this motion may affect the ligand exchange was tested in this work. Mutations of Gly 56 were carried out and the mutant proteins were expressed in Sf9 cells using the Baculovirus expression system. The recombinant proteins were purified using Ni-NTA affinity chromatography. The results for the (GDP/GTP) exchange assays showed that G56S and double mutants (D55G/G56S) proteins differ significantly from the wild type and D55G mutant forms. The K _d values for GTPS binding of those mutants have decreased by approximately 10-fold. No difference in the GTPase activity was detected for the mutants. Thus, the biochemical results obtained support the conclusions of the computational studies.     

Conformational analysis and design of cross-strand disulfides in antiparallel β-sheets

Cross‐strand disulfides bridge two cysteines in a registered pair of antiparallel β‐strands. A nonredundant data set comprising 5025 polypeptides containing 2311 disulfides was used to study cross‐strand disulfides. Seventy‐six cross‐strand disulfides were found of which 75 and 1 occurred at non‐hydrogen‐bonded (NHB) and hydrogen‐bonded (HB) registered pairs, respectively. Conformational analysis and modeling studies demonstrated that disulfide formation at HB pairs necessarily requires an extremely rare and positive χ¹ value for at least one of the cysteine residues. Disulfides at HB positions also have more unfavorable steric repulsion with the main chain. Thirteen pairs of disulfides were introduced in NHB and HB pairs in four model proteins: leucine binding protein (LBP), leucine, isoleucine, valine binding protein (LIVBP), maltose binding protein (MBP), and Top7. All mutants LIVBP T247C V331C showed disulfide formation either on purification, or on treatment with oxidants. Protein stability in both oxidized and reduced states of all mutants was measured. Relative to wild type, LBP and MBP mutants were destabilized with respect to chemical denaturation, although the sole exposed NHB LBP mutant showed an increase of 3.1°C in T _m . All Top7 mutants were characterized for stability through guanidinium thiocyanate chemical denaturation. Both exposed and two of the three buried NHB mutants were appreciably stabilized. All four HB Top7 mutants were destabilized (ΔΔG ⁰ = ?3.3 to ?6.7 kcal/mol). The data demonstrate that introduction of cross‐strand disulfides at exposed NHB pairs is a robust method of improving protein stability. All four exposed Top7 disulfide mutants showed mild redox activity. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Comparative modeling of the conformational stability of chymotrypsin inhibitor 2 protein mutants using amino acid sequence autocorrelation (AASA) and amino acid 3D autocorrelation (AA3DA) vectors and ensembles of Bayesian-regularized genetic neural networks

Predicting protein stability changes upon point mutation is important for understanding protein structure and designing new proteins. Autocorrelation vector formalism was extended to amino acid sequences and 3D conformations for encoding protein structural information with modeling purpose. Protein autocorrelation vectors were weighted by 48 amino acid/residue properties selected from the AAindex database. Ensembles of Bayesian-regularized genetic neural networks (BRGNNs) trained with amino acid sequence autocorrelation (AASA) vectors and amino acid 3D autocorrelation (AA3DA) vectors yielded predictive models of the change of unfolding Gibbs free energy change (ΔΔG) of chymotrypsin Inhibitor 2 protein mutants. The ensemble predictor described about 58 and 72% of the data variances in test sets for AASA and AA3DA models, respectively. Optimum sequence and 3D-based ensembles exhibit high effects on relevant structural (volume, solvent-accessible surface area), physico-chemical (hydrophilicity/hydrophobicity-related) and thermodynamic (hydration parameters) properties.     

Alterations by a defect in a rice G protein α subunit in probenazole and pathogen-induced responses

The rice dwarf1 (d1) mutant, which is deficient in an α subunit (Gα) of heterotrimeric G protein, was used to obtain specific evidence on the functions of Gα protein in defence signalling in rice. Using proteome analysis, a probenazole‐inducible protein (PBZ1) was detected in the cytosolic fraction of leaf blade of the wild type, but not the d1 mutant. After treatment with probenazol, PBZ1 reached maximal levels at 72 h in the wild type but 96 h in the d1 mutant. The induction of PBZ1 by probenazole treatment was inhibited by protein kinase inhibitors. A 48‐kDa putative mitogen‐activated protein kinase (MAPK) and a 55‐kDa putative Ca²⁺‐dependent protein kinase (CDPK) showed lower activities in the cytosolic fraction of the d1 mutant than that of the wild type. The activities of these protein kinases were enhanced at 24 h in the wild type and 48 h in the d1 mutant after probenazole treatment. Although the d1 mutant responded to the rice blast fungus similarly to the wild type, the d1 mutant developed rice blight symptoms earlier than the wild type when infected with Xoo. In addition, the blight symptoms were more severe on the mutant than on the wild type, and wilting was frequently observed in the d1 mutant. Furthermore, induction by the bacterial infection of the 48‐kDa putative MAPK and PBZ1 was delayed by 2 and 4 d, respectively, in the d1 mutant compared with the wild type. These results indicate that the Gα protein plays a role in the induction of PBZ1 and protein kinases by probenazole and Xoo, and suggest that the 48‐kDa putative MAPK may be involved in a signalling pathway for resistance to bacterial infection.     

Low expression of lipid-linked oligosaccharide due to a functionally altered Dol-P-Man synthase reduces protein glycosylation in cAMP-dependent protein kinase deficient Chinese hamster ovary cells

Chinese hamster ovary cells express a wide variety of glycoproteins with M_r ranging from 15,000 to 200,000 dalton and higher. Glycosylation of these proteins was much less in cAMP-dependent protein kinase (PKA)-deficient mutants which expressed either (i) a defective C-subunit with altered substrate specificity and having no detectable type II kinase (mutant 10215); or (ii) an altered RI subunit and having no detectable type II kinase (mutant 10248); or (iii) exhibited the lowest level of total kinase with no detectable type I kinase but having a small amount of type II kinase (mutant 10260). Addition of 8Br-cAMP enhanced protein glycosylation index in wild type cells 10001 by 120% but only 7 to 23% in the mutant cells. The rate of lipid-linked oligosaccharide (LLO) biosynthesis was linear for 1 h in all cell types, but the total amount of LLO expressed was much less in PKA-deficient mutants. Pulse-chase experiments indicated that the t_1/2 for LLO turnover was also twice as high in PKA-deficient cells as in the wild type. Size exclusion chromatography of the mild-acid released oligosaccharide confirmed that both wild type and the mutant cells synthesized Glc₃Man₉GlcNAc₂-PP-Dol as the most predominating species with no accumulation of Man₅GlcNAc₂-PP-Dol in the mutants. Kinetic studies exhibited a reduced mannosylphosphodolichol synthase (DPMS) activity in mutant cells with a K_m for GDP-mannose 160 to 400% higher than that of the wild type. In addition, the k_cat for DPMS was also reduced 2 to 4-fold in these mutant cells. Exogenously added Dol-P failed to rescue the k_cat for DPMS in CHO cell mutants; however, in vitro protein phosphorylation with a cAMP-dependent protein kinase restored their kinetic activity to the level of the wild type. Published in 2004.     

Analysis of the functions of TaGW2 homoeologs in wheat grain weight and protein content traits

GW2 is emerging as a key genetic determinant of grain weight in cereal crops; it has three homoeologs (TaGW2‐A1, ‐B1 and ‐D1) in hexaploid common wheat (Triticum aestivum L.). Here, by analyzing the gene editing mutants that lack one (B1 or D1), two (B1 and D1) or all three (A1, B1 and D1) homoeologs of TaGW2, several insights are gained into the functions of TaGW2‐B1 and ‐D1 in common wheat grain traits. First, both TaGW2‐B1 and ‐D1 affect thousand‐grain weight (TGW) by influencing grain width and length, but the effect conferred by TaGW2‐B1 is stronger than that of TaGW2‐D1. Second, there exists functional interaction between TaGW2 homoeologs because the TGW increase shown by a double mutant (lacking B1 and D1) was substantially larger than that of their single mutants. Third, both TaGW2‐B1 and ‐D1 modulate cell number and length in the outer pericarp of developing grains, with TaGW2‐B1 being more potent. Finally, TaGW2 homoeologs also affect grain protein content as this parameter was generally increased in the mutants, especially in the lines lacking two or three homoeologs. Consistent with this finding, two wheat end‐use quality‐related parameters, flour protein content and gluten strength, were considerably elevated in the mutants. Collectively, our data shed light on functional difference between and additive interaction of TaGW2 homoeologs in the genetic control of grain weight and protein content traits in common wheat, which may accelerate further research on this important gene and its application in wheat improvement.     

A gibberellin-regulated protein phosphorylated by a putative Ca2+-dependent protein kinase is G-protein mediated in rice root

The GA-signal transduction pathways downstream to the Gα protein in rice seedling root were investigated using in-gel kinase assay and in vitro protein phosphorylation techniques with a Gα protein defective mutant, d1. A 50-kDa protein kinase was detected downstream to Gα protein in the membrane fraction of rice seedling roots using an in-gel kinase assay with histone III-S as a substrate. The activity of a 50-kDa protein kinase increased in the wild-type rice by gibberellin (GA₃) treatment, but did not change in the d1 mutant. This protein kinase activity was inhibited by the Ca²⁺ chelator ethyleneglycol-bis-(beta-aminoethylether)-N,N,N ¹,N ¹-tetraacetic acid (EGTA), protein kinase inhibitors, staurosporine and H7, and calmodulin antagonist, trifluoperazine, suggesting that the 50-kDa protein kinase is a putative plant Ca²⁺-dependent protein kinase (CDPK). The activity of the 50-kDa putative CDPK reached its highest level at 3 h after GA₃ treatment and then gradually declined with time. In order to identify the endogenous substrate for 50-kDa putative CDPK, two-dimensional polyacrylamide gel electrophoresis followed by in vitro protein phosphorylation was carried out. The phosphorylation activity of an endogenous protein PP30, identified as an unknown protein having molecular weight 30 kDa and isoelectric point 5.8 was increased in the wild-type rice by GA₃ treatment, compared with the d1 mutant. The addition of GA₃ treated membrane fraction, which predominantly represent a 50-kDa putative CDPK further increased the phosphorylation of PP30. Almost similar to GA₃ treatment, phosphorylation activity of PP30 was also increased by the treatment with cholera toxin in the wild-type rice but not in d1 mutant. These results suggest that the 50-kDa putative CDPK and an unknown protein, PP30 promoted by GA₃ treatment are G-protein mediated in rice seedling roots.     

E75、R78和D82是影响大肠埃希氏菌FtsZ自身组装及其与MreB相互作用的重要氨基酸

【目的】探索大肠埃希氏菌Escherichia coli FtsZ突变体FtsZ~(E75A)、FtsZ~(R78G)和FtsZ~(D82A)对FtsZ自身组装和FtsZ-MreB相互作用的影响。【方法】利用常规分子克隆和定点突变技术,构建FtsZ及其突变体表达载体,亲和纯化得到相应的目标蛋白;通过同源重组构建QN6(ftsZ::yfp-cat)、QN7(ftsZ~(E75A)::yfp-cat)、QN8(ftsZ~(R78G)::yfp-cat)和QN9(ftsZ~(D82A)::yfp-cat)菌株;利用活细胞成像技术观察FtsZ及其突变体的胞内定位模式;免疫沉淀和细菌双杂交实验检测FtsZ/FtsZ*-FtsZ*或FtsZ/FtsZ*-MreB间的相互作用;光扫描检测定点突变对FtsZ组装特性的影响。【结果】FtsZ~(E75A)、FtsZ~(R78G)和FtsZ~(D82A)突变体的功能活性降低、各突变体在E.coli内不能正确的定位和形成功能性Z环;FtsZ/FtsZ*-FtsZ*单体间的相互作用减弱或消失,FtsZ*-MreB相互作用破坏;FtsZ突变体体外聚合效率降低。【结论】FtsZ E75、R78和D82是影响FtsZ正确组装和功能及FtsZ-MreB相互作用的重要氨基酸。     

Reversibility of DnaA protein activity in the'irreversible'dnaA204 mutant of Escherichia coli

The dnaA204 mutant, one of the so-called irreversible dnaA mutants which cannot reinitiate chromosome replication upon a shift from non-permissive to permissive growth temperature in the absence of protein synthesis, was reinvestigated using flow cytometry and marker frequency analysis. In a temperature downshift experiment and in the presence of protein synthesis the dnaA204 mutant reinitiates chromosome replication very fast. Using a lac promoter-controlled wild type or a dnaA204 mutant gene carried on a plasmid, we have observed instantaneous initiation of replication when synthesis of DnaA protein is induced in the dnaA204 mutant at 42δC. The data indicate that the dnaA204 mutant after a shift to 42δC still contains functional DnaA protein, but that the activity level is below the initiation threshold. Thus, after synthesis of very small amounts of additional DnaA protein, initiation occurs very fast both after a shift to 30δC, and after induction of DnaA protein synthesis at 42 C. A model describing the processing of DnaA protein in mutants and in the wild type Is presented.     

Functional analysis of the PsbX protein by deletion of the corresponding gene in Synechocystis sp. PCC 6803

The psbX gene (sml0002) coding for a 4.1 kDa protein in Photosystem II of plants and cyanobacteria was deleted in both wild type and in a Photosystem I-less mutant of the cyanobacterium Synechocystis sp. PCC 6803. Polymerase chain reaction and sequencing analysis showed that the mutants had completely segregated. Deletion of the PsbX protein does not seem to influence growth rate, electron transport or water oxidation ability. Whereas a high light induction of the psbX mRNA could be observed in wild type, deletion of the gene did not lead to high light sensibility. Light saturation measurements and 77K fluorescence measurements indicated a minor disconnection of the antenna in the deletion mutant. Furthermore, fluorescence induction measurements as well as immuno-staining of the D1 protein showed that the amount of Photosystem II complexes in the mutants was reduced by 30%. Therefore, PsbX does not seem to be necessary for the Photosystem II electron transport, but directly or indirectly involved in the regulation of the amount of functionally active Photosystem II centres in Synechocystis sp. PCC 6803.