pH‐selective mutagenesis of protein–protein interfaces: In silico design of therapeutic antibodies with prolonged half‐life

Understanding the effects of mutation on pH‐dependent protein binding affinity is important in protein design, especially in the area of protein therapeutics. We propose a novel method for fast in silico mutagenesis of protein–protein complexes to calculate the effect of mutation as a function of pH. The free energy differences between the wild type and mutants are evaluated from a molecular mechanics model, combined with calculations of the equilibria of proton binding. The predicted pH‐dependent energy profiles demonstrate excellent agreement with experimentally measured pH‐dependency of the effect of mutations on the dissociation constants for the complex of turkey ovomucoid third domain (OMTKY3) and proteinase B. The virtual scanning mutagenesis identifies all hotspots responsible for pH‐dependent binding of immunoglobulin G (IgG) to neonatal Fc receptor (FcRn) and the results support the current understanding of the salvage mechanism of the antibody by FcRn based on pH‐selective binding. The method can be used to select mutations that change the pH‐dependent binding profiles of proteins and guide the time consuming and expensive protein engineering experiments. As an application of this method, we propose a computational strategy to search for mutations that can alter the pH‐dependent binding behavior of IgG to FcRn with the aim of improving the half‐life of therapeutic antibodies in the target organism. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Investigation on bioactive protection of the amino acids derived from LEA protein on insulin by molecular simulation

Nowadays heat-sensitive protein medicines are increasingly showing their importance in the treatment of various diseases. Their popularisation and application are meeting a great challenge because of their heat lability. In this study, human insulin as a heat-sensitive protein medicine and 66 amino acids derived from a Group 3 late embryogenesis abundant protein fragment as a complex bioactive protectant, were chosen to be investigated to determine whether these amino acids can be used to protect the insulin from denaturation due to drying. The experiments were carried out by using a replica exchange molecular dynamics (REMD) simulation and GROMACS software with Gromos96 (53a6) force field. The REMD results indicate that those amino acids can effectively prevent the reversal between hydrophilic and hydrophobic surface. Both the configurations and secondary structures of the protected insulin were preserved very well. The H-bonding and electrostatic interactions between the insulin and the protectant play key roles in the bioactive protection of insulin. These results agree well with the water replacement hypothesis. All the results prove that these amino acids are a perfect bioactive protectant for heat-sensitive protein medicines.     

Molecular dynamics simulation on the decomposition of type SII hydrogen hydrate and the performance of tetrahydrofuran as a stabiliser

Molecular dynamics simulation is used to study the decomposition and stability of SII hydrogen and hydrogen/tetrahydrofuran (THF) hydrates at 150 K, 220 K and 100 bar. The modelling of the microscopic decomposition process of hydrogen hydrate indicates that the decomposition of hydrogen hydrate is led by the diffusive behaviour of H₂ molecules. The hydrogen/THF hydrate presents higher stability, by comparing the distributions of the tetrahedral angle of H₂O molecules, radial distribution functions of H₂O molecules and mean square displacements or diffusion coefficients of H₂O and H₂ molecules in hydrogen hydrate with those in hydrogen/THF hydrate. It is also found that the resistance of the diffusion behaviour of H₂O and H₂ molecules can be enhanced by encaging THF molecules in the (5¹²6⁴) cavities. Additionally, the motion of THF molecules is restricted due to its high interaction energy barrier. Accordingly, THF, as a stabiliser, is helpful in increasing the stability of hydrogen hydrate.     

ATM pathway is essential for ionizing radiation-induced autophagy

BackgroundATM plays an important role in response to DNA damage, while the roles of ATM in radiation-induced autophagy are still unclear in cervical cancer cells.MethodsHuman cervical cancer cells, Hela, were used, and cell models with ATM^?/? and MAPK14^?/? were established by gene engineering. Western blot was implemented to detect protein expression. MDC staining and GFP-LC3 relocalization were used to detect autophagy. CCK-8 was used to detect cell viability. Radiosensitivity was analyzed by colony formation assays. Co-immunoprecipitation was used to detect the interaction between different proteins, and apoptosis was detected by flow cytometry.ResultsAfter radiation autophagy was induced, illustrated by the increase of MAPLC3-II/MAPLC3-I ratio and decrease of p62, and phosphorylation of ATM simultaneously increased. ATM^?/? cells displayed hypersensitivity but had no influence on IR-induced apoptosis. Then inhibitor of ATM, KU55933, ATM and MAPK14 silencing were used, and autophagy was induced by IR more than 200% in control, and only by 35.72%, 53.18% and 24.76% in KU55933-treated cells, ATM^?/? and MAPK14^?/? cells, respectively. KU55933 inhibited IR-induced autophagy by activating mTOR pathways. ATM silencing decreased the expression of MAPK14 and mTOR signals significantly. Beclin's bond to PI3KIII and their interaction increased after IR, while in ATM^?/? and MAPK14^?/? cells this interaction decreased after IR. Both ATM and MAPK14 interacted with Beclin, while ATM^?/? and MAPK14^?/? cells showed no interaction.ConclusionsATM could promote IR-induced autophagy via the MAPK14 pathway, the mTOR pathway, and Beclin/PI3KIII complexes, which contributed to the effect of ATM on radiosensitivity.     

Protein dynamics and motions in relation to their functions: several case studies and the underlying mechanisms

Proteins are dynamic entities in cellular solution with functions governed essentially by their dynamic personalities. We review several dynamics studies on serine protease proteinase K and HIV-1 gp120 envelope glycoprotein to demonstrate the importance of investigating the dynamic behaviors and molecular motions for a complete understanding of their structure–function relationships. Using computer simulations and essential dynamic (ED) analysis approaches, the dynamics data obtained revealed that: (i) proteinase K has highly flexible substrate-binding site, thus supporting the induced-fit or conformational selection mechanism of substrate binding; (ii) Ca²⁺ removal from proteinase K increases the global conformational flexibility, decreases the local flexibility of substrate-binding region, and does not influence the thermal motion of catalytic triad, thus explaining the experimentally determined decreased thermal stability, reduced substrate affinity, and almost unchanged catalytic activity upon Ca²⁺ removal; (iii) substrate binding affects the large concerted motions of proteinase K, and the resulting dynamic pocket can be connected to substrate binding, orientation, and product release; (iv) amino acid mutations 375 S/W and 423 I/P of HIV-1 gp120 have distinct effects on molecular motions of gp120, facilitating 375 S/W mutant to assume the CD4-bound conformation, while 423 I/P mutant to prefer for CD4-unliganded state. The mechanisms underlying protein dynamics and protein–ligand binding, including the concept of the free energy landscape (FEL) of the protein–solvent system, how the ruggedness and variability of FEL determine protein’s dynamics, and how the three ligand-binding models, the lock-and-key, induced-fit, and conformational selection are rationalized based on the FEL theory are discussed in depth.     

Protein Kinase LTRPK1 Influences Cold Adaptation and Microtubule Stability in Rice

Rice LTRPK1, which encodes a member of the casein kinase I family, has been reported to be involved in root development, hormone response, and metabolic processes. Here we further show that LTRPK1 participates in stress resistance by regulating cytoskeleton rearrangement and formation of cold tolerance and adaptation. Semiquantitative RT-PCR analysis revealed enhanced expression of LTRPK1 in plants subject to low-temperature stress at 4 °C, suggesting a role in low-temperature-related cell responses and signal transduction pathways. Further analysis of LTRPK1-deficient transgenic plants showed that under low-temperature treatment, the growth rate of transgenic plant primary roots, which is commonly used as an indicator for cold stress response abilities, was less inhibited than that of control plants. Moreover, damage to the plasma membrane of root cells in LTRPK1-deficient plants was greater than that of controls as measured by relative electrical conductivity (REC). The malondialdehyde (MDA) content of LTRPK1-deficient plants also increased over that of the control, indicating increased plasma membrane permeability. Further immunofluorescence localization observations indicated that microtubules of transgenic plants subject to low temperature disassembled more rapidly, whereas the control plant microtubules in most cells of the root elongation zone kept their normal habitus, which suggested that LTRPK1-deficient plants had reduced capacity to resist low-temperature stress through regulation of microtubule assembly. These results demonstrate involvement of LTRPK1 in low-temperature stress and provide new insight for rice breeding and germplasm innovation to improve crop cold tolerance.     

Proposal of Ensifer psoraleae sp. nov., Ensifer sesbaniae sp. nov., Ensifer morelense comb. nov. and Ensifer americanum comb. nov

In a survey of rhizobia associated with the native legumes in Yunnan Province, China, seven and nine strains isolated from the root nodules of Psoralea corylifolia, Sesbania cannabina and Medicago lupulina were respectively classified into the novel genomic species groups I and II in the genus Ensifer (former Sinorhizobium) based on the sequence analyses of the 16S rRNA gene. Analyses of concatenated housekeeping genes (atpD, recA and glnII) further revealed that they were distinct lineages in the genus, and group I was most similar to Ensifer terangae and Ensifer garamanticus (both with 94.2% similarity), while group II was most similar to Ensifer adhaerens (94.0%). These groups could be distinguished from closely related species by DNA–DNA relatedness, MALID-TOF MS, cellular fatty acid profiles and a series of phenotypic characters. Therefore, two novel species were proposed: Ensifer psoraleae sp. nov. (seven strains, type strain CCBAU 65732^T = LMG 26835^T = HAMBI 3286^T) and Ensifer sesbaniae sp. nov. (nine strains, type strain CCBAU 65729^T = LMG 26833^T = HAMBI 3287^T). They had a DNA G + C mol% (T_m) of 58.9 and 60.4, respectively. Both of the type strains formed effective nodules on common bean (Phaseolus vulgaris) and their hosts of origin. In addition, the previously described species Sinorhizobium morelense and Sinorhizobium americanum were renamed as Ensifer morelense comb. nov. and Ensifer americanum comb. nov. according to the accumulated data from different studies.     

Serum proteomic profiling for the early diagnosis of colorectal cancer

No ideal serum biomarker currently exists for the early diagnosis of colorectal cancer (CRC). Magnetic bead‐based fractionation coupled with MALDI‐TOF MS was used to screen serum samples from CRC patients, healthy controls, and other cancer patients. A diagnostic model with five proteomic features (m/z 1778.97, 1866.16, 1934.65, 2022.46, and 4588.53) was generated using Fisher algorithm with best performance. The Fisher‐based model could discriminate CRC patients from the controls with 100% (46/46) sensitivity and 100% (35/35) specificity in the training set, 95.6% (43/45) sensitivity and 83.3% (35/42) specificity in the test set. We further validated the model with 94.4% (254/269) sensitivity and 75.5% (83/110) specificity in the external independent group. In other cancers group, the Fisher‐based model classified 25 of 46 samples (54.3%) as positive and the other 21 as negative. With FT‐ICR‐MS, the proteomic features of m/z 1778.97, 1866.16, 1934.65, and 2022.46, of which intensities decreased significantly in CRC, were identified as fragments of complement C3f. Therefore, the Fisher‐based model containing five proteomic features was able to effectively differentiate CRC patients from healthy controls and other cancers with a high sensitivity and specificity, and may be CRC‐specific. Serum complement C3f, which was significantly decreased in CRC group, may be relevant to the incidence of CRC. J. Cell. Biochem. 114: 448–455, 2013. © 2012 Wiley Periodicals, Inc.