Protein Allostery Study in Cells Using NMR Spectroscopy.

Protein allostery is commonly observed in vitro. But how protein allostery behaves in cells is unknown. In this work, a protein monomer-dimer equilibrium system was built with the allosteric effect on the binding characterized using NMR spectroscopy through mutations away from the dimer interface. A chemical shift linear fitting method was developed that enabled us to accurately determine the dissociation constant. A total of 28 allosteric mutations were prepared and grouped to negative allosteric, nonallosteric, and positive allosteric modulators. ∼ 50% of mutations displayed the allosteric-state changes when moving from a buffered solution into cells. For example, there were no positive allosteric modulators in the buffered solution but eight in cells. The change in protein allostery is correlated with the interactions between the protein and the cellular environment. These interactions presumably drive the surrounding macromolecules in cells to transiently bind to the monomer and dimer mutational sites and change the free energies of the two species differently which generate new allosteric effects. These surrounding macromolecules create a new protein allostery pathway that is only present in cells.

Improving the Enzymatic Activity and Stability of a Lytic Polysaccharide Monooxygenase.

Lytic Polysaccharide Monooxygenases (LPMOs) are copper-dependent enzymes that play a pivotal role in the enzymatic conversion of the most recalcitrant polysaccharides, such as cellulose and chitin. Hence, protein engineering is highly required to enhance their catalytic efficiencies. To this effect, we optimized the protein sequence encoding for an LPMO from Bacillus amyloliquefaciens (BaLPMO10A) using the sequence consensus method. Enzyme activity was determined using the chromogenic substrate 2,6-Dimethoxyphenol (2,6-DMP). Compared with the wild type (WT), the variants exhibit up to a 93.7% increase in activity against 2,6-DMP. We also showed that BaLPMO10A can hydrolyze p-nitrophenyl-ß-D-cellobioside (PNPC), carboxymethylcellulose (CMC), and phosphoric acid-swollen cellulose (PASC). In addition to this, we investigated the degradation potential of BaLPMO10A against various substrates such as PASC, filter paper (FP), and Avicel, in synergy with the commercial cellulase, and it showed up to 2.7-, 2.0- and 1.9-fold increases in production with the substrates PASC, FP, and Avicel, respectively, compared to cellulase alone. Moreover, we examined the thermostability of BaLPMO10A. The mutants exhibited enhanced thermostability with an apparent melting temperature increase of up to 7.5 °C compared to the WT. The engineered BaLPMO10A with higher activity and thermal stability provides a better tool for cellulose depolymerization.

Engineering the Active Site Lid Dynamics to Improve the Catalytic Efficiency of Yeast Cytosine Deaminase.

Conformational dynamics is important for enzyme catalysis. However, engineering dynamics to achieve a higher catalytic efficiency is still challenging. In this work, we develop a new strategy to improve the activity of yeast cytosine deaminase (yCD) by engineering its conformational dynamics. Specifically, we increase the dynamics of the yCD C-terminal helix, an active site lid that controls the product release. The C-terminal is extended by a dynamical single α-helix (SAH), which improves the product release rate by up to ~8-fold, and the overall catalytic rate kcat by up to ~2-fold. It is also shown that the kcat increase is due to the favorable activation entropy change. The NMR H/D exchange data indicate that the conformational dynamics of the transition state analog complex increases as the helix is extended, elucidating the origin of the enhanced catalytic entropy. This study highlights a novel dynamics engineering strategy that can accelerate the overall catalysis through the entropy-driven mechanism.

A Tyrosine Phosphoproteome Analysis Approach Enabled by Selective Dephosphorylation with Protein Tyrosine Phosphatase.

Protein tyrosine phosphorylation (pTyr) plays a prominent role in signal transduction and regulation in all eukaryotic cells. In conventional immunoaffinity purification (IP) methods, phosphotyrosine peptides are isolated from the digest of cellular protein extracts with a phosphotyrosine-specific antibody and are identified by tandem mass spectrometry. However, low sensitivity, poor reproducibility, and high cost are universal concerns for IP approaches. In this study, we presented an antibody-free approach to identify phosphotyrosine peptides by using protein tyrosine phosphatase (PTP). It was found that most of the PTPs including PTP1B, TCPTP, and SHP1 can efficiently and selectively dephosphorylate phosphotyrosine peptides. We then designed a workflow by combining two Ti4+-IMAC-based phosphopeptide enrichment steps with PTP-catalyzed dephosphorylation for tyrosine phosphoproteomics analysis. This workflow was first validated by selective detection of phosphotyrosine peptides from semicomplex samples and then applied to analyze the tyrosine phosphoproteome of Jurkat T cells. Around 1000 putative former phosphotyrosine peptides were identified from less than 500 µg of cell lysate. The tyrosine phosphosites on the majority of these peptides could be unambiguously determined for over 70% of them possessing only one tyrosine residue. It was also found that the tyrosine sites identified by this method were highly complementary to those identified by the SH2 superbinder-based method. Therefore, the combination of Ti4+-IMAC enrichment with PTP dephosphorylation provides an alternative and cost-effective approach for tyrosine phosphoproteomics analysis.

Detecting the Hydrogen Bond Cooperativity in a Protein ß-Sheet by H/D Exchange.

The hydrogen bond (H-bond) cooperativity in the ß-sheet of GB3 is investigated by a NMR hydrogen/deuterium (H/D) exchange method. It is shown that the weakening of one backbone N-H…O=C H-bond between two ß-strands, ß1 and ß2, due to the exchange of NH to ND of the H-bond donor in ß1, perturbs the chemical shift of 13Cα, 13Cß, 1Hα, 1HN, and 15N of the H-bond acceptor and its following residue in ß2. Quantum mechanical calculations suggest that the -H-bond chemical shift isotope effect is caused by the structural reorganization in response to the H-bond weakening. This structural reorganization perturbs four neighboring H-bonds, with three being weaker and one being stronger, indicating that three H-bonds are cooperative and one is anticooperative with the perturbed H-bond. The sign of the cooperativity depends on the relative position of the H-bonds. This H-bond cooperativity, which contributes to ß-sheet stability overall, can be important for conformational coupling across the ß-sheet.

Quantifying Protein Electrostatic Interactions in Cells by Nuclear Magnetic Resonance Spectroscopy.

Most proteins perform their functions in cells. How the cellular environment modulates protein interactions is an important question. In this work, electrostatic interactions between protein charges were studied using in-cell nuclear magnetic resonance (NMR) spectroscopy. A total of eight charge pairs were introduced in protein GB3. Compared to the charge pair electrostatic interactions in a buffer, five charge pairs in cells displayed no apparent changes whereas three pairs had the interactions weakened by more than 70%. Further investigation suggests that the transfer free energy is responsible for the electrostatic interaction modulation. Both the transfer free energy of the folded state and that of the unfolded state can contribute to the cellular environmental effect on protein electrostatics, although the latter is generally larger (more negative) than the former. Our work highlights the importance of direct in-cell studies of protein interactions and thus protein function.

Structural Basis for Selective Oxidation of Phosphorylated Ethylphenols by Cytochrome P450 Monooxygenase CreJ.

Selective oxidation of C-H bonds in alkylphenols holds great significance for not only structural derivatization in pharma- and biomanufacturing but also biological degradation of these toxic chemicals in environmental protection. A unique chemomimetic biocatalytic system using enzymes from a p-cresol biodegradation pathway has recently been developed. As the central biocatalyst, the cytochrome P450 monooxygenase CreJ oxidizes diverse p- and m-alkylphenol phosphates with perfect stereoselectivity at different efficiencies. However, the mechanism of regio- and stereoselectivity of this chemomimetic biocatalytic system remained unclear. Here, using p- and m-ethylphenol substrates, we elucidate the CreJ-catalyzed key steps for selective oxidations. The crystal structure of CreJ in complex with m-ethylphenol phosphate was solved and compared with its complex structure with p-ethylphenol phosphate isomer. The results indicate that the conformational changes of substrate-binding residues are slight, while the substrate promiscuity is achieved mainly by the available space in the catalytic cavity. Moreover, the catalytic preferences of regio- and stereoselective hydroxylation for the two ethylphenol substrates is explored by molecular dynamics simulations. The ethyl groups in the complexes display different flexibilities, and the distances of the active oxygen to H pro-S and H pro-R of methylene agree with the experimental stereoselectivity. The regioselectivity can be explained by the distances and bond dissociation energy. These results provide not only the mechanistic insights into CreJ regio- and stereoselectivity but also the structural basis for further P450 enzyme design and engineering.IMPORTANCE The key cytochrome P450 monooxygenase CreJ showed excellent regio- and stereoselectivity in the oxidation of various alkylphenol substrates. C-H bond functionalization of these toxic alkylphenols holds great significance for both biological degradation of these environmental chemicals and production of value-added structural derivatives in pharmaceutical and biochemical industries. Our results, combined with in vitro enzymatic assays, crystal structure determination of enzyme-substrate complex, and molecular dynamics simulations, provide not only significant mechanism elucidation of the regio- and stereoselective catalyzation mediated by CreJ but also the promising directions for future engineering efforts of this enzyme toward more useful products. It also has great extendable potential to couple this multifunctional P450 enzyme with other biocatalysts (e.g., hydroxyl-based glycosylase) to access more alkylphenol-derived high-value chemicals through environment-friendly biocatalysis and biotransformation.

Generating Five Independent Molecular Alignments for Simultaneous Protein Structure and Dynamics Determination Using Nuclear Magnetic Resonance Spectroscopy.

Residual dipolar couplings (RDCs) are commonly used in NMR for protein structure and dynamics studies, but it is challenging to generate five independent RDC data sets (required for simultaneous structure and dynamics determination) for most protein molecules in the magnetic field. In this work, a reporter protein with a lanthanide tag is introduced to create five independent alignments. This reporter protein is then attached to target proteins where five independent sets of RDCs are also obtained for the target proteins. The fitting of RDCs provides important information about the structure and dynamics of the target proteins. The method is simple and effective and, in principle, can be used to generate complete sets of RDCs for different protein molecules.

Selective oxidation of aliphatic C-H bonds in alkylphenols by a chemomimetic biocatalytic system.

Selective oxidation of aliphatic C-H bonds in alkylphenols serves significant roles not only in generation of functionalized intermediates that can be used to synthesize diverse downstream chemical products, but also in biological degradation of these environmentally hazardous compounds. Chemo-, regio-, and stereoselectivity; controllability; and environmental impact represent the major challenges for chemical oxidation of alkylphenols. Here, we report the development of a unique chemomimetic biocatalytic system originated from the Gram-positive bacterium Corynebacterium glutamicum The system consisting of CreHI (for installation of a phosphate directing/anchoring group), CreJEF/CreG/CreC (for oxidation of alkylphenols), and CreD (for directing/anchoring group offloading) is able to selectively oxidize the aliphatic C-H bonds of p- and m-alkylated phenols in a controllable manner. Moreover, the crystal structures of the central P450 biocatalyst CreJ in complex with two representative substrates provide significant structural insights into its substrate flexibility and reaction selectivity.

Combined genome-wide association study and transcriptome analysis reveal candidate genes for resistance to Fusarium ear rot in maize.

Fusarium ear rot, caused by Fusarium verticillioides, is a devastating fungal disease in maize that reduces yield and quality; moreover, F. verticillioides produces fumonisin mycotoxins, which pose serious threats to human and animal health. Here, we performed a genome-wide association study (GWAS) under three environmental conditions and identified 34 single-nucleotide polymorphisms (SNPs) that were significantly associated with Fusarium ear rot resistance. With reference to the maize B73 genome, 69 genes that overlapped with or were adjacent to the significant SNPs were identified as potential resistance genes to Fusarium ear rot. Comparing transcriptomes of the most resistant and most susceptible lines during the very early response to Fusarium ear rot, we detected many differentially expressed genes enriched for pathways related to plant immune responses, such as plant hormone signal transduction, phenylpropanoid biosynthesis, and cytochrome P450 metabolism. More than one-fourth of the potential resistance genes detected in the GWAS were differentially expressed in the transcriptome analysis, which allowed us to predict numbers of candidate genes for maize resistance to ear rot, including genes related to plant hormones, a MAP kinase, a PR5-like receptor kinase, and heat shock proteins. We propose that maize plants initiate early immune responses to Fusarium ear rot mainly by regulating the growth-defense balance and promoting biosynthesis of defense compounds.

Characterization of Residue Specific Protein Folding and Unfolding Dynamics in Cells.

In this work, we measured the millisecond residue specific protein folding and unfolding dynamics in E. coli cells for two protein GB3 mutants using NMR. The results show that the protein folding and unfolding dynamics in cells is different from that in buffer. Through a two-site exchange model, it is shown that both the population and the exchange rate are changed by the E. coli cellular environment. Further investigation suggests that the change is likely due to the quinary interaction with crowded molecules in the cell. Our work underlines the importance of cellular environment to protein folding kinetics and thermodynamics although this environmental effect may not be large enough to change the protein structure.

Arginine side chain stacking with peptide plane stabilizes the protein helix conformation in a cooperative way.

A combined experimental and computational study is performed for arginine side chain stacking with the protein α-helix. Theremostability measurements of Aristaless homeodomain, a helical protein, suggest that mutating the arginine residue R106, R137 or R141, which has the guanidino side chain stacking with the peptide plane, to alanine, destabilizes the protein. The R-PP stacking has an energy of ∼0.2-0.4 kcal/mol. This stacking interaction mainly comes from dispersion and electrostatics, based on MP2 calculations with the energy decomposition analysis. The calculations also suggest that the stacking stabilizes 2 backbone-backbone h-bonds (i→i-4 and i-3→i-7) in a cooperative way. Desolvation and electrostatic polarization are responsible for cooperativity with the i→i-4 and i-3→i-7 h-bonds, respectively. This cooperativity is supported by a protein α-helices h-bond survey in the pdb databank where stacking shortens the corresponding h-bond distances.

Direct Observation of CH/CH van der Waals Interactions in Proteins by NMR.

van der Waals interactions are important to protein stability and function. These interactions are usually identified empirically based on protein 3D structures. In this work, we performed a solution nuclear magnetic resonance (NMR) spectroscopy study of van der Waals interactions by detecting the through-space vdw JCC-coupling between protein aliphatic side chain groups. Specifically, vdw JCC-coupling values up to ∼0.5 Hz were obtained between the methyl and nearby aliphatic groups in protein GB3, providing direct experimental evidence for the van der Waals interactions. Quantum mechanical calculations suggest that the J-coupling is correlated with the exchange-repulsion term of van der Waals interaction. NMR detection of vdw JCC-coupling offers a new tool to characterize such interactions in proteins.

Processive Degradation of Crystalline Cellulose by a Multimodular Endoglucanase via a Wirewalking Mode.

Processive hydrolysis of crystalline cellulose by cellulases is a critical step for lignocellulose deconstruction. The classic Trichoderma reesei exoglucanase TrCel7A, which has a closed active-site tunnel, starts each processive run by threading the tunnel with a cellulose chain. Loop regions are necessary for tunnel conformation, resulting in weak thermostability of fungal exoglucanases. However, endoglucanase CcCel9A, from the thermophilic bacterium Clostridium cellulosi, comprises a glycoside hydrolase (GH) family 9 module with an open cleft and five carbohydrate-binding modules (CBMs) and hydrolyzes crystalline cellulose processively. How CcCel9A and other similar GH9 enzymes bind to the smooth surface of crystalline cellulose to achieve processivity is still unknown. Our results demonstrate that the C-terminal CBM3b and three CBMX2s enhance productive adsorption to cellulose, while the CBM3c adjacent to the GH9 is tightly bound to 11 glucosyl units, thereby extending the catalytic cleft to 17 subsites, which facilitates decrystallization by forming a supramodular binding surface. In the open cleft, the strong interaction forces between substrate-binding subsites and glucosyl rings enable cleavage of the hydrogen bonds and extraction of a single cellulose chain. In addition, subsite -4 is capable of drawing the chain to its favored location. Cellotetraose is released from the open cleft as the initial product to achieve high processivity, which is further hydrolyzed to cellotriose, cellobiose and glucose by the catalytic cleft of the endoglucanase. On this basis, we propose a wirewalking mode for processive degradation of crystalline cellulose by an endoglucanase, which provides insights for rational design of industrial cellulases.

Structural insights into the substrate specificity of a glycoside hydrolase family 5 lichenase from Caldicellulosiruptor sp. F32.

Glycoside hydrolase (GH) family 5 is one of the largest GH families with various GH activities including lichenase, but the structural basis of the GH5 lichenase activity is still unknown. A novel thermostable lichenase F32EG5 belonging to GH5 was identified from an extremely thermophilic bacterium Caldicellulosiruptor sp. F32. F32EG5 is a bi-functional cellulose and a lichenan-degrading enzyme, and exhibited a high activity on ß-1,3-1,4-glucan but side activity on cellulose. Thin-layer chromatography and NMR analyses indicated that F32EG5 cleaved the ß-1,4 linkage or the ß-1,3 linkage while a 4-O-substitued glucose residue linked to a glucose residue through a ß-1,3 linkage, which is completely different from extensively studied GH16 lichenase that catalyses strict endo-hydrolysis of the ß-1,4-glycosidic linkage adjacent to a 3-O-substitued glucose residue in the mixed-linked ß-glucans. The crystal structure of F32EG5 was determined to 2.8 Šresolution, and the crystal structure of the complex of F32EG5 E193Q mutant and cellotetraose was determined to 1.7 Šresolution, which revealed that the exit subsites of substrate-binding sites contribute to both thermostability and substrate specificity of F32EG5. The sugar chain showed a sharp bend in the complex structure, suggesting that a substrate cleft fitting to the bent sugar chains in lichenan is a common feature of GH5 lichenases. The mechanism of thermostability and substrate selectivity of F32EG5 was further demonstrated by molecular dynamics simulation and site-directed mutagenesis. These results provide biochemical and structural insights into thermostability and substrate selectivity of GH5 lichenases, which have potential in industrial processes.

Dual-Functional Small Molecules for Generating an Efficient Cytochrome P450BM3 Peroxygenase.

We report a unique strategy for the development of a H2 O2 -dependent cytochrome P450BM3 system, which catalyzes the monooxygenation of non-native substrates with the assistance of dual-functional small molecules (DFSMs), such as N-(ω-imidazolyl fatty acyl)-l-amino acids. The acyl amino acid group of DFSM is responsible for bounding to enzyme as an anchoring group, while the imidazolyl group plays the role of general acid-base catalyst in the activation of H2 O2 . This system affords the best peroxygenase activity for the epoxidation of styrene, sulfoxidation of thioanisole, and hydroxylation of ethylbenzene among those P450-H2 O2 system previously reported. This work provides the first example of the activation of the normally H2 O2 -inert P450s through the introduction of an exogenous small molecule. This approach improves the potential use of P450s in organic synthesis as it avoids the expensive consumption of the reduced nicotinamide cofactor NAD(P)H and its dependent electron transport system. This introduces a promising approach for exploiting enzyme activity and function based on direct chemical intervention in the catalytic process.

Quinary Interactions Weaken the Electric Field Generated by Protein Side-Chain Charges in the Cell-like Environment.

The intramolecular electric field (e-field) generated by protein GB3 side-chain charges K/E10, K/E19, and D/K40 was measured in the absence or presence of macromolecular crowding. The e-field responds differently to different crowding agents-dextran, Ficoll, BSA, and E. coli cell lysate. Dextran and Ficoll have no effect on the e-field. The lysate generally weakens the e-field but the amplitude of weakening varies greatly. For example, the e-field by K19 is reduced by 67% in the presence of 90 g/L lysate, corresponding to a charge change from 0.9 to 0.3 e for K19, whereas the e-fields by D/K40 are weakened only by ∼7% under the same lysate concentration. The extent of the e-field weakening by BSA is in between that by Ficoll (dextran) and lysate. Further investigations suggest that the e-field weakening mechanism by lysate is similar to that by NaCl. That is, the e-field generated by a protein surface charge affects the distribution of lysate which creates a reaction field and weakens the protein e-field. Our study indicates that the protein electrostatic property can be changed significantly due to quinary interaction with the cell environment.

qRfg3, a novel quantitative resistance locus against Gibberella stalk rot in maize.

KEY MESSAGE: A quantitative trait locus  qRfg3 imparts recessive resistance to maize Gibberella stalk rot. qRfg3 has been mapped into a 350-kb interval and could reduce the disease severity index by ~26.6%. Gibberella stalk rot, caused by the fungal pathogen Fusarium graminearum, severely affects maize yield and grain quality worldwide. To identify more resistance quantitative trait loci (QTLs) against this disease, we analyzed a recombinant inbred line (RIL) population derived from a cross between resistant H127R and susceptible C7-2 inbred lines. Within this population, maize resistance to Gibberella stalk rot had high broad-sense heritability. A major QTL, qRfg3, on chromosome 3 was consistently detected across three field trials, accounting for 10.7-19.4% of the total phenotypic variation. Using a progeny-based sequential fine-mapping strategy, we narrowed qRfg3 down to an interval of ~350 kb. We further demonstrated that qRfg3 is a recessive resistance locus to Gibberella stalk rot that reduced the disease severity index by ~26.6%. Both the gene location and recessive genetic mode distinguish qRfg3 from other stalk rot resistance loci. Hence, qRfg3 is valuable as a complement to existing resistance QTLs to improve maize resistance to Gibberella stalk rot.

Entropy Drives the Formation of Salt Bridges in the Protein GB3.

Salt bridges are very common in proteins. But what drives the formation of protein salt bridges is not clear. In this work, we determined the strength of four salt bridges in the protein GB3 by measuring the ΔpKa values of the basic residues that constitute the salt bridges with a highly accurate NMR titration method at different temperatures. The results show that the ΔpKa values increase with temperature, thus indicating that the salt bridges are stronger at higher temperatures. Fitting of ΔpKa values to the van't Hoff equation yields positive ΔH and ΔS values, thus indicating that entropy drives salt-bridge formation. Molecular dynamics simulations show that the protein and solvent make opposite contributions to ΔH and ΔS. Specifically, the enthalpic gain contributed from the protein is more than offset by the enthalpic loss contributed from the solvent, whereas the entropic gain originates from the desolvation effect.