Linghao Zhong

Sensitivity-enhanced solid-state NMR detection of expansin’s target in plant cell walls

Significance The protein expansin loosens the cell walls of plants for cell growth, but its carbohydrate target of binding has been elusive because of the difficulty of studying the noncrystalline plant cell wall by most structural biology techniques and the trace amount of expansin needed for wall loosening. We have now combined dynamic nuclear polarization sensitivity-enhanced solid-state NMR spectroscopy with 13C labeling of Arabidopsis thaliana plants and 13C, 15N labeling of expansin, to determine that expansin binds cellulose microfibrils to loosen the plant cell wall. The expansin binding site is enriched in the hemicellulose xyloglucan and has a different cellulose structure from bulk cellulose, shedding light on the mechanism of wall loosening. Structure determination of protein binding to noncrystalline macromolecular assemblies such as plant cell walls (CWs) poses a significant structural biology challenge. CWs are loosened during growth by expansin proteins, which weaken the noncovalent network formed by cellulose, hemicellulose, and pectins, but the CW target of expansins has remained elusive because of the minute amount of the protein required for activity and the complex nature of the CW. Using solid-state NMR spectroscopy, combined with sensitivity-enhancing dynamic nuclear polarization (DNP) and differential isotopic labeling of expansin and polysaccharides, we have now determined the functional binding target of expansin in the Arabidopsis thaliana CW. By transferring the electron polarization of a biradical dopant to the nuclei, DNP allowed selective detection of 13C spin diffusion from trace concentrations of 13C, 15N-labeled expansin in the CW to nearby polysaccharides. From the spin diffusion data of wild-type and mutant expansins, we conclude that to loosen the CW, expansin binds highly specific cellulose domains enriched in xyloglucan, whereas more abundant binding to pectins is unrelated to activity. Molecular dynamics simulations indicate short 13C-13C distances of 4–6 Å between a hydrophobic surface of the cellulose microfibril and an aromatic motif on the expansin surface, consistent with the observed NMR signals. DNP-enhanced 2D 13C correlation spectra further reveal that the expansin-bound cellulose has altered conformation and is enriched in xyloglucan, thus providing unique insight into the mechanism of CW loosening. DNP-enhanced NMR provides a powerful, generalizable approach for investigating protein binding to complex macromolecular targets.

Investigation of the effect of glycosylation on human prion protein by molecular dynamics.

Abstract Prion protein conformational isomerization, PrPc→PrPSc, has been attributed as the cause of TSE diseases such as mad-cow disease. The mechanism of such isomerization, however, is little known due the experimental difficulties in studying the scrapie form. Among factors that affect PrP isomerization, the role which glycosylation plays remains vague. The number of innumerous glycan species, together with their high flexibility, leads to ineffective structural characterization. In this research, we studied the effect of chitobiose glycosylation on human PrP, in both monomeric (huPrPmono) and dimeric (huPrPdimer) forms, by molecular dynamics (MD) simulations. Our results show that this glycosylation has minimal impact on the structure of huPrPmono. However, it affects the secondary structure of dimeric protein. An additional β-sheet strand is found while the glycosylation is absent in the huPrPdimer. Comparatively, when the protein is glycosylated with chitobiose, such β-sheet addition is not observed.

Exposure of hydrophobic core in human prion protein pathogenic mutant H187R.

Abstract Pathogenesis studies have revealed that H187R mutation of human prion protein (huPrP) is related to GSS type of TSE diseases. Its pathogenic mechanism is still unclear. We here studied the globular domain of this mutant protein by molecular dynamics simulations. Compared to the wide-type protein, the mutant has similar dynamics and stability profiles in our simulation. Conformational rearrangements are concentrated around the mutation site, due to the introduction the positively charged side chain of Arg187. The strong electrostatic repulsion between Arg156 and Arg187 drives both side chains away from their original positions, leaving its hydrophobic core to be solvent accessible. Such a unfavorable conformational change may destabilize the mutant protein and make it more susceptible to unfolding.

Hydrogen-bonding network and OH stretch vibration of cellulose: comparison of computational modeling with polarized IR and SFG spectra

Hydrogen bonds play critical roles in noncovalent directional interactions determining the crystal structure of cellulose. Although diffraction studies accurately determined the coordinates of carbon and oxygen atoms in crystalline cellulose, the structural information on hydrogen atoms involved in hydrogen-bonding is still elusive. This could be complemented by vibrational spectroscopy; but the assignment of the OH stretch peaks has been controversial. In this study, we performed calculations using density functional theory with dispersion corrections (DFT-D2) for the cellulose Iβ crystal lattices with the experimentally determined carbon and oxygen coordinates. DFT-D2 calculations revealed that the OH stretch vibrations of cellulose are highly coupled and delocalized through intra- and interchain hydrogen bonds involving all OH groups in the crystal. Additionally, molecular dynamics (MD) simulations of a single cellulose microfibril showed that the conformations of OH groups exposed at the microfibril surface are not well-defined. Comparison of the computation results with the experimentally determined IR dichroism of uniaxially aligned cellulose microfibrils and the peak positions of various cellulose crystals allowed unambiguous identification of OH stretch modes observed in the vibrational spectra of cellulose.

Cellulose Microfibril Twist, Mechanics, and Implication for Cellulose Biosynthesis

All-atom molecular dynamics simulations with explicit water solvent were used to investigate the microstructure and conformational dynamics of cellulose Iβ microfibrils as a function of microfibril length and cross-sectional size and shape. Cellulose microfibrils quickly develop a right-handed twist, which then remains stable over the entire 10 ns simulation time. The helical angle is independent of microfibril length and inversely proportional to its cross-sectional area, in accord with the expectations of continuum theory for an intrinsic chiral twist that is opposed by torsional shear. These calculations provide-to our knowledge-the first estimates of the shear modulus of a cellulose microfibril from MD simulations. The internal strains caused by this helical twist, propagated indefinitely along the microfibril axis, could be relaxed by periodic regions of amorphous structure along the axis of the cellulose microfibrils.

Peptide Mimotopes of Pneumococcal Capsular Polysaccharide of 6B Serotype: A Peptide Mimotope Can Bind to Two Unrelated Antibodies

Two groups of bacteriophage clones displaying the antigenic properties of serotype 6B pneumococcal capsular polysaccharide (PS) were obtained from different phage libraries expressing random heptameric peptides. One group, biopanned with a mouse mAb (Hyp6BM1), is comprised of 17 phage clones expressing 10 unique sequences of linear peptides. The other group, selected with another mAb (Hyp6BM8), contained six clones, all of which expressed the identical circular peptide. Phage clones expressing the linear peptides (e.g., PhaM1L3) bound only to Hyp6BM1, but not other 6B PS-specific mAb, and their binding could be inhibited with pneumococcal capsular type 6B PS only. In contrast, a phage clone expressing the circular peptide (PhaM8C1) cross-reacted with several other 6B PS-specific mAbs, and their binding could be inhibited with pneumococcal capsular PS of 6A and 6B serotypes. Two short peptides, PepM1L3 and PepM8C1, reflecting the peptide inserts of the corresponding phage clones, could inhibit the binding of the two clones to their respective mAb. Interestingly, the peptide insert in PhaM8C1 was identical to that in PhaB3C4, a previously reported mimotope of α(2→8) polysialic acid, Neisseria meningitidis group B PS. Indeed, PhaM8C1 bound to HmenB3 (a meningococcal Ab), and their association could be inhibited with α(2–8) polysialic acid, but not with 6B PS. Conversely, α(2–8) polysialic acid could not inhibit the binding of PhaM8C1 to Hyp6BM8. The two-dimensional nuclear magnetic resonance studies indicate that PepM8C1 peptide can assume several conformations in solution. The ability of this peptide to assume multiple conformations might account for its ability to mimic more than one Ag type.

Dynamic Multibody Protein Interactions Suggest Versatile Pathways for Copper Trafficking

As part of intracellular copper trafficking pathways, the human copper chaperone Hah1 delivers Cu(+) to the Wilsons Disease Protein (WDP) via weak and dynamic protein-protein interactions. WDP contains six homologous metal binding domains (MBDs) connected by flexible linkers, and these MBDs all can receive Cu(+) from Hah1. The functional roles of the MBD multiplicity in Cu(+) trafficking are not well understood. Building on our previous study of the dynamic interactions between Hah1 and the isolated fourth MBD of WDP, here we study how Hah1 interacts with MBD34, a double-domain WDP construct, using single-molecule fluorescence resonance energy transfer (smFRET) combined with vesicle trapping. By alternating the positions of the smFRET donor and acceptor, we systematically probed Hah1-MBD3, Hah1-MBD4, and MBD3-MBD4 interaction dynamics within the multidomain system. We found that the two interconverting interaction geometries were conserved in both intermolecular Hah1-MBD and intramolecular MBD-MBD interactions. The Hah1-MBD interactions within MBD34 are stabilized by an order of magnitude relative to the isolated single-MBDs, and thermodynamic and kinetic evidence suggest that Hah1 can interact with both MBDs simultaneously. The enhanced interaction stability of Hah1 with the multi-MBD system, the dynamic intramolecular MBD-MBD interactions, and the ability of Hah1 to interact with multiple MBDs simultaneously suggest an efficient and versatile mechanism for the Hah1-to-WDP pathway to transport Cu(+).

Sugar-binding sites on the surface of the carbohydrate-binding module of CBH I from Trichoderma reesei.

Molecular dynamics simulations were carried out for a system consisting of the carbohydrate-binding module (CBM) of the cellulase CBH I from Trichoderma reesei (Hypocrea jecorina) in a concentrated solution of β-D-glucopyranose, to determine whether there is any tendency for the sugar molecules to bind to the CBM. In spite of the general tendency of glucose to behave as an osmolyte, a marked tendency for the sugar molecules to bind to the protein was observed. However, the glucose molecules tended to bind only to specific sites on the protein. As expected, the hydrophobic face of the sugar molecules, comprising the axial H1, H3, and H5 aliphatic protons, tended to adhere to the flat faces of the three tyrosine side chains on the planar binding surface of the CBM. However, a significant tendency to bind to a groove-like feature on the upper surface of the CBM was also observed. These results would not be inconsistent with a model of the mechanism for this globular domain in which the cellodextrin chain being removed from the surface of crystalline cellulose passes over the upper surface of the CBM, presumably then available for hydrolysis in the active site tunnel of this processive cellulase.

Quantum mechanical calculations on cellulose–water interactions: structures, energetics, vibrational frequencies and NMR chemical shifts for surfaces of Iα and Iβ cellulose

AbstractPeriodic and molecular cluster density functional theory calculations were performed on the Iα (001), Iα (021), Iβ (100), and Iβ (110) surfaces of cellulose with and without explicit H2O molecules of hydration. The energy-minimized H-bonding structures, water adsorption energies, vibrational spectra, and 13C NMR chemical shifts are discussed. The H-bonded structures and water adsorption energies (ΔEads) are used to distinguish hydrophobic and hydrophilic cellulose–water interactions. O–H stretching vibrational modes are assigned for hydrated and dry cellulose surfaces. Calculations of the 13C NMR chemical shifts for the C4 and C6 surface atoms demonstrate that these δ13C4 and δ13C6 values can be upfield shifted from the bulk values as observed without rotation of the hydroxymethyl groups from the bulk tg conformation to the gt conformation as previously assumed.

Insight into structural rearrangements and interdomain interactions related to electron transfer between flavin mononucleotide and heme in nitric oxide synthase: A molecular dynamics study.

Calmodulin (CaM) binding to nitric oxide synthase (NOS) enables a conformational change, in which the FMN domain shuttles between the FAD and heme domains to deliver electrons to the active site heme center. A clear understanding of this large conformational change is critical, since this step is the rate-limiting in NOS catalysis. Herein molecular dynamics simulations were conducted on a model of an oxygenase/FMN (oxyFMN) construct of human inducible NOS (iNOS). This is to investigate the structural rearrangements and the domain interactions related to the FMN-heme interdomain electron transfer (IET). We carried out simulations on the iNOS oxyFMN·CaM complex models in [Fe(III)][FMNH(-)] and [Fe(II)][FMNH] oxidation states, the pre- and post-IET states. The comparison of the dynamics and conformations of the iNOS construct at the two oxidation states has allowed us to identify key factors related to facilitating the FMN-heme IET process. The computational results demonstrated, for the first time, that the conformational change is redox-dependent. Predictions of the key interacting sites in optimal interdomain FMN/heme docking are well supported by experimental data in the literature. An intra-subunit pivot region is predicted to modulate the FMN domain motion and correlate with existence of a bottleneck in the conformational sampling that leads to the electron transfer-competent state. Interactions of the residues identified in this work are proposed to ensure that the FMN domain moves with appropriate degrees of freedom and docks to proper positions at the heme domain, resulting in efficient IET and nitric oxide production.