Real-Time Observation of Conformational Changes and Translocation of Endogenous Cytochrome c within Intact Mitochondria.

Cytochrome c (cyt c) is a multifunctional protein with varying conformations. However, the conformation of cyt c in its native environment, mitochondria, is still unclear. Here, we applied NMR spectroscopy to investigate the conformation and location of endogenous cyt c within intact mitochondria at natural isotopic abundance, mainly using widespread methyl groups as probes. By monitoring time-dependent chemical shift perturbations, we observed that most cyt c is located in the inner mitochondrial membrane and partially unfolded, which is distinct from its native conformation in solution. When suffering oxidative stress, cyt c underwent oxidative modifications due to increasing reactive oxygen species (ROS), weakening electrostatic interactions with the membrane, and gradually translocating into the inner membrane spaces of mitochondria. Meanwhile, the lethality of oxidatively modified cyt c to cells was reduced compared with normal cyt c. Our findings significantly improve the understanding of the molecular mechanisms underlying the regulation of ROS by cyt c in mitochondria. Moreover, it highlights the potential of NMR to monitor high-concentration molecules at a natural isotopic abundance within intact cells or organelles.

Prediction of order parameters based on protein NMR structure ensemble and machine learning.

The fast motions of proteins at the picosecond to nanosecond timescale, known as fast dynamics, are closely related to protein conformational entropy and rearrangement, which in turn affect catalysis, ligand binding and protein allosteric effects. The most used NMR approach to study fast protein dynamics is the model free method, which uses order parameter S2 to describe the amplitude of the internal motion of local group. However, to obtain order parameter through NMR experiments is quite complex and lengthy. In this paper, we present a machine learning approach for predicting backbone 1H-15N order parameters based on protein NMR structure ensemble. A random forest model is used to learn the relationship between order parameters and structural features. Our method achieves high accuracy in predicting backbone 1H-15N order parameters for a test dataset of 10 proteins, with a Pearson correlation coefficient of 0.817 and a root-mean-square error of 0.131.

Critical role of minor eggcase silk component in promoting spidroin chain alignment and strong fiber formation.

Natural spider silk with extraordinary mechanical properties is typically spun from more than one type of spidroin. Although the main components of various spider silks have been widely studied, little is known about the molecular role of the minor silk components in spidroin self-assembly and fiber formation. Here, we show that the minor component of spider eggcase silk, TuSp2, not only accelerates self-assembly but remarkably promotes molecular chain alignment of spidroins upon physical shearing. NMR structure of the repetitive domain of TuSp2 reveals that its dimeric structure with unique charged surface serves as a platform to recruit different domains of the main eggcase component TuSp1. Artificial fiber spun from the complex between TuSp1 and TuSp2 minispidroins exhibits considerably higher strength and Young's modulus than its native counterpart. These results create a framework for rationally designing silk biomaterials based on distinct roles of silk components.

Effects of ligand binding on dynamics of fatty acid binding protein and interactions with membranes.

Intracellular transport of fatty acids involves binding of ligands to their carrier fatty acid binding proteins (FABPs) and interactions of ligand-free and -bound FABPs with membranes. Previous studies focused on ligand-free FABPs. Here, our amide hydrogen exchange data showed that oleic acid binding to human intestinal FABP (hIFABP) stabilizes the protein, most likely through enhancing the hydrogen-bonding network, and induces rearrangement of sidechains even far away from the ligand binding site. Using NMR relaxation techniques, we found that the ligand binding affects not only conformational exchanges between major and minor states but also the affinity of hIFABP to nanodiscs. Analyses of the relaxation and amide exchange data suggested that two minor native-like states existing in both ligand-free and -bound hIFABPs originate from global "breathing" motions, while one minor native-like state comes from local motions. The amide hydrogen exchange data also indicated that helix αII undergoes local unfolding through which ligands can exit from the binding cavity.

C-Terminal Domains of Spider Silk Proteins Having Divergent Structures but Conserved Functional Roles.

Spider silk is self-assembled from silk proteins or spidroins. C-terminal domains (CTDs) of various types of spidroins are relatively conserved in amino acid sequences and are suggested to adopt similar structures and perform similar functional roles in spidroin storage and silk formation. Here, we solved the structure of the CTD from a capture-spiral silk protein (CTDFl) and characterized its stability and fibril formation in the presence and absence of a reducing agent at different pH values. CTDFl adopts a dimeric structure with 8 helices, but the CTDs of other types of spidroins exist in a domain-swapped dimeric structure with 10 helices. Despite the structural differences, CTDFl is pH-responsive in stability and fibril formation, similar to the CTDs from minor and major ampullate spidroins. Thus, the functional role of CTDs in silk fiber formation seems conserved. Comparing wild-type CTDFl and its mutants, we found that the pH-responsive behavior results from the protonation of H76, which is conserved from different spider species. In addition, the fibril formation rate of CTDFl correlates with its instability, suggesting that structural changes are involved in fibril formation.

Structural basis of DNA binding to human YB-1 cold shock domain regulated by phosphorylation.

Human Y-box binding protein 1 (YB-1) is a multifunctional protein and overexpressed in many types of cancer. It specifically recognizes DNA/RNA through a cold shock domain (CSD) and regulates nucleic acid metabolism. The C-terminal extension of CSD and the phosphorylation of S102 are indispensable for YB-1 function. Until now, the roles of the C-terminal extension and phosphorylation in gene transcription and translation are still largely unknown. Here, we solved the structure of human YB-1 CSD with a C-terminal extension sequence (CSDex). The structure reveals that the extension interacts with several residues in the conventional CSD and adopts a rigid structure instead of being disordered. Either deletion of this extension or phosphorylation of S102 destabilizes the protein and results in partial unfolding. Structural characterization of CSDex in complex with a ssDNA heptamer shows that all the seven nucleotides are involved in DNA-protein interactions and the C-terminal extension provides a unique DNA binding site. Our DNA-binding study indicates that CSDex can recognize more DNA sequences than previously thought and the phosphorylation reduces its binding to ssDNA dramatically. Our results suggest that gene transcription and translation can be regulated by changing the affinity of CSDex binding to DNA and RNA through phosphorylation, respectively.

Easily Established and Multifunctional Synthetic Nanobody Libraries as Research Tools.

Nanobodies, or VHHs, refer to the antigen-binding domain of heavy-chain antibodies (HCAbs) from camelids. They have been widely used as research tools for protein purification and structure determination due to their small size, high specificity, and high stability, overcoming limitations with conventional antibody fragments. However, animal immunization and subsequent retrieval of antigen-specific nanobodies are expensive and complicated. Construction of synthetic nanobody libraries using DNA oligonucleotides is a cost-effective alternative for immunization libraries and shows great potential in identifying antigen-specific or even conformation-specific nanobodies. This review summarizes and analyses synthetic nanobody libraries in the current literature, including library design and biopanning methods, and further discusses applications of antigen-specific nanobodies obtained from synthetic libraries to research.

Conformational exchange of fatty acid binding protein induced by protein-nanodisc interactions.

Fatty acid binding proteins (FABPs) can facilitate the transfer of long-chain fatty acids between intracellular membranes across considerable distances. The transfer process involves fatty acids, their donor membrane and acceptor membrane, and FABPs, implying that potential protein-membrane interactions exist. Despite intensive studies on FABP-membrane interactions, the interaction mode remains elusive, and the protein-membrane association and dissociation rates are inconsistent. In this study, we used nanodiscs (NDs) as mimetic membranes to investigate FABP-membrane interactions. Our NMR experiments showed that human intestinal FABP interacts weakly with both negatively charged and neutral membranes, but it prefers the negatively charged one. Through simultaneous analysis of NMR relaxation in the rotating-frame (R1ρ), relaxation dispersion, chemical exchange saturation transfer, and dark-state exchange saturation transfer data, we estimated the affinity of the protein to negatively charged NDs, the dissociation rate, and apparent association rate. We further showed that the protein in the ND-bound state adopts a conformation different from the native structure and the second helix is very likely involved in interactions with NDs. We also found a membrane-induced FABP conformational state that exists only in the presence of NDs. This state is native-like, different from other conformational states in structure, unbound to NDs, and in dynamic equilibrium with the ND-bound state.

Dengue virus protease activity modulated by dynamics of protease cofactor.

The viral protease domain (NS3pro) of dengue virus is essential for virus replication, and its cofactor NS2B is indispensable for the proteolytic function. Although several NS3pro-NS2B complex structures have been obtained, the dynamic property of the complex remains poorly understood. Using NMR relaxation techniques, here we found that NS3pro-NS2B exists in both closed and open conformations that are in dynamic equilibrium on a submillisecond timescale in aqueous solution. Our structural information indicates that the C-terminal region of NS2B is disordered in the minor open conformation but folded in the major closed conformation. Using mutagenesis, we showed that the closed-open conformational equilibrium can be shifted by changing NS2B stability. Moreover, we revealed that the proteolytic activity of NS3pro-NS2B correlates well with the population of the closed conformation. Our results suggest that the closed-open conformational equilibrium can be used by both nature and humanity to control the replication of dengue virus.

Anti-Fungal Hevein-like Peptides Biosynthesized from Quinoa Cleavable Hololectins.

Chitin-binding hevein-like peptides (CB-HLPs) belong to a family of cysteine-rich peptides that play important roles in plant stress and defense mechanisms. CB-HLPs are ribosomally synthesized peptides that are known to be bioprocessed from the following two types of three-domain CB-HLP precursor architectures: cargo-carrying and non-cargo-carrying. Here, we report the identification and characterization of chenotides biosynthesized from the third type of precursors, which are cleavable hololectins of the quinoa (Chenopodium quinoa) family. Chenotides are 6-Cys-CB-HLPs of 29-31 amino acids, which have a third type of precursor architecture that encompasses a canonical chitin-binding domain that is involved in chitin binding and anti-fungal activities. Microbroth dilution assays and microscopic analyses showed that chenotides are effective against phyto-pathogenic fungi in the micromolar range. Structure determination revealed that chenotides are cystine knotted and highly compact, which could confer resistance against heat and proteolytic degradation. Importantly, chenotides are connected by a novel 18-residue Gly/Ala-rich linker that is a target for bioprocessing by cathepsin-like endopeptidases. Taken together, our findings reveal that chenotides are a new family of CB-HLPs from quinoa that are synthesized as a single multi-modular unit and bioprocessed to yield individual mature CB-HLPs. Importantly, such precursors constitute a new family of cleavable hololectins. This unusual feature could increase the biosynthetic efficiency of anti-fungal CB-HLPs, to provide an evolutionary advantage for plant survival and reproduction.

Ligand Entry into Fatty Acid Binding Protein via Local Unfolding Instead of Gap Widening.

Fatty acid binding proteins play an important role in the transportation of fatty acids. Despite intensive studies, how fatty acids enter the protein cavity for binding is still controversial. Here, a gap-closed variant of human intestinal fatty acid binding protein was generated by mutagenesis, in which the gap is locked by a disulfide bridge. According to its structure determined here by NMR, this variant has no obvious openings as the ligand entrance and the gap cannot be widened by internal dynamics. Nevertheless, it still takes up fatty acids and other ligands. NMR relaxation dispersion, chemical exchange saturation transfer, and hydrogen-deuterium exchange experiments show that the variant exists in a major native state, two minor native-like states, and two locally unfolded states in aqueous solution. Local unfolding of either ßB-ßD or helix 2 can generate an opening large enough for ligands to enter the protein cavity, but only the fast local unfolding of helix 2 is relevant to the ligand entry process.

Roseltide rT7 is a disulfide-rich, anionic, and cell-penetrating peptide that inhibits proteasomal degradation.

Disulfide-rich plant peptides with molecular masses of 2-6 kDa represent an expanding class of peptidyl-type natural products with diverse functions. They are structurally compact, hyperstable, and underexplored as cell-penetrating agents that inhibit intracellular functions. Here, we report the discovery of an anionic, 34-residue peptide, the disulfide-rich roseltide rT7 from Hibiscus sabdariffa (of the Malvaceae family) that penetrates cells and inhibits their proteasomal activities. Combined proteomics and NMR spectroscopy revealed that roseltide rT7 is a cystine-knotted, six-cysteine hevein-like cysteine-rich peptide. A pair-wise comparison indicated that roseltide rT7 is >100-fold more stable against protease degradation than its S-alkylated analog. Confocal microscopy studies and cell-based assays disclosed that after roseltide rT7 penetrates cells, it causes accumulation of ubiquitinated proteins, inhibits human 20S proteasomes, reduces tumor necrosis factor-induced IκBα degradation, and decreases expression levels of intercellular adhesion molecule-1. Structure-activity studies revealed that roseltide rT7 uses a canonical substrate-binding mechanism for proteasomal inhibition enabled by an IIML motif embedded in its proline-rich and exceptionally long intercysteine loop 4. Taken together, our results provide mechanistic insights into a novel disulfide-rich, anionic, and cell-penetrating peptide, representing a potential lead for further development as a proteasomal inhibitor in anti-cancer or anti-inflammatory therapies.

Special Issue "Selected Papers from the 8th Asia-Pacific NMR (APNMR) Symposium: Recent Advances in NMR Spectroscopy".

Asia-Pacific NMR (APNMR) has been an important scientific event in the region, engaging a large number of NMR scientists from academia and industries [...].

Structural Basis of Oligomerization of N-Terminal Domain of Spider Aciniform Silk Protein.

Spider silk is self-assembled from water-soluble silk proteins through changes in the environment, including pH, salt concentrations, and shear force. The N-terminal domains of major and minor ampullate silk proteins have been found to play an important role in the assembly process through salt- and pH-dependent dimerization. Here, we identified the sequences of the N-terminal domains of aciniform silk protein (AcSpN) and major ampullate silk protein (MaSpN) from Nephila antipodiana (NA). Different from MaSpN, our biophysical characterization indicated that AcSpN assembles to form large oligomers, instead of a dimer, upon condition changes from neutral to acidic pH and/or from a high to low salt concentration. Our structural studies, by nuclear magnetic resonance spectroscopy and homology modelling, revealed that AcSpN and MaSpN monomers adopt similar overall structures, but have very different charge distributions contributing to the differential self-association features. The intermolecular interaction interfaces for AcSp oligomers were identified using hydrogen-deuterium exchange mass spectrometry and mutagenesis. On the basis of the monomeric structure and identified interfaces, the oligomeric structures of AcSpN were modelled. The structural information obtained will facilitate an understanding of silk fiber formation mechanisms for aciniform silk protein.

Atomistic Insights into the Functional Instability of the Second Helix of Fatty Acid Binding Protein.

Structural dynamics of fatty acid binding proteins (FABPs), which accommodate poorly soluble ligands in the internalized binding cavities, are intimately related to their function. Recently, local unfolding of the α-helical cap in a variant of human intestinal FABP (IFABP) has been shown to correlate with the kinetics of ligand association, shedding light on the nature of the critical conformational reorganization. Yet, the physical origin and mechanism of the functionally relevant transient unfolding remain elusive. Here, we investigate the intrinsic structural instability of the second helix (αII) of IFABP in comparison with other segments of the protein using hydrogen-exchange NMR spectroscopy, microsecond molecular dynamics simulations, and enhanced sampling techniques. Although tertiary interactions positively contribute to the stability of helices in IFABP, the intrinsic unfolding tendency of αII is encoded in its primary sequence and can be described by the Lifson-Roig theory in the absence of tertiary interactions. The unfolding pathway of αII in intact proteins involves an on-pathway intermediate state that is characterized with the fraying of the last helical turn, captured by independent enhanced sampling methods. The simulations in this work, combined with hydrogen-exchange NMR data, provide new, to our knowledge, atomistic insights into the functional local unfolding of FABPs.

Investigating energy-based pool structure selection in the structure ensemble modeling with experimental distance constraints: The example from a multidomain protein Pub1.

The structural variations of multidomain proteins with flexible parts mediate many biological processes, and a structure ensemble can be determined by selecting a weighted combination of representative structures from a simulated structure pool, producing the best fit to experimental constraints such as interatomic distance. In this study, a hybrid structure-based and physics-based atomistic force field with an efficient sampling strategy is adopted to simulate a model di-domain protein against experimental paramagnetic relaxation enhancement (PRE) data that correspond to distance constraints. The molecular dynamics simulations produce a wide range of conformations depicted on a protein energy landscape. Subsequently, a conformational ensemble recovered with low-energy structures and the minimum-size restraint is identified in good agreement with experimental PRE rates, and the result is also supported by chemical shift perturbations and small-angle X-ray scattering data. It is illustrated that the regularizations of energy and ensemble-size prevent an arbitrary interpretation of protein conformations. Moreover, energy is found to serve as a critical control to refine the structure pool and prevent data overfitting, because the absence of energy regularization exposes ensemble construction to the noise from high-energy structures and causes a more ambiguous representation of protein conformations. Finally, we perform structure-ensemble optimizations with a topology-based structure pool, to enhance the understanding on the ensemble results from different sources of pool candidates.

A Method for Determining Structure Ensemble of Large Disordered Protein: Application to a Mechanosensing Protein.

Structure characterization of intrinsically disordered proteins (IDPs) remains a key obstacle in understanding their functional mechanisms. Due to the highly dynamic feature of IDPs, structure ensembles instead of static unique structures are often derived from experimental data. Several state-of-the-art computational methods have been developed to select an optimal ensemble from a pregenerated structure pool, but they suffer from low efficiency for large IDPs. Here we present a matching pursuit genetic algorithm (MPGA) for structure ensemble determination, which takes advantages from both matching pursuit (MP) to reduce the search space and genetic algorithm (GA) to reduce the restriction on constraint types. The MPGA method is validated using a reference ensemble with predefined structures. In comparison with the conventional GA, MPGA takes much less computational time for large IDPs. The utility of the method is demonstrated by application to structure ensemble determination of a mechanosensing protein domain with 306 amino acids. The structure ensemble determined reveals that the N-terminal region 1-240 is more compact than the C-terminal region 240-306. The unique structural feature explains why only a small portion of YXXP tyrosine residues can be phosphorylated easily by kinases in the absence of extension force and why the phosphorylation is force-dependent.

Pin1At encoding a peptidyl-prolyl cis/trans isomerase regulates flowering time in Arabidopsis.

Floral transition in plants is regulated by an integrated network of flowering genetic pathways. We show that an Arabidopsis PIN1-type parvulin 1, Pin1At, controls floral transition by accelerating cis/trans isomerization of the phosphorylated Ser/Thr-Pro motifs in two MADS-domain transcription factors, SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1) and AGAMOUS-LIKE 24 (AGL24). Pin1At regulates flowering, which is genetically mediated by AGL24 and SOC1. Pin1At interacts with the phosphorylated AGL24 and SOC1 in vitro and with AGL24 and SOC1 in vivo and accelerates the cis/trans conformational change of phosphorylated Ser/Thr-Pro motifs of AGL24 and SOC1. We further demonstrate that these Ser/Thr-Pro motifs are important for Pin1At function in promoting flowering through AGL24 and SOC1 and that the interaction between Pin1At and AGL24 mediates the AGL24 stability in the nucleus. Taken together, we propose that phosphorylation-dependent prolyl cis/trans isomerization of key transcription factors is an important flowering regulatory mechanism.

A New Method for Determining Structure Ensemble: Application to a RNA Binding Di-Domain Protein.

Structure ensemble determination is the basis of understanding the structure-function relationship of a multidomain protein with weak domain-domain interactions. Paramagnetic relaxation enhancement has been proven a powerful tool in the study of structure ensembles, but there exist a number of challenges such as spin-label flexibility, domain dynamics, and overfitting. Here we propose a new (to our knowledge) method to describe structure ensembles using a minimal number of conformers. In this method, individual domains are considered rigid; the position of each spin-label conformer and the structure of each protein conformer are defined by three and six orthogonal parameters, respectively. First, the spin-label ensemble is determined by optimizing the positions and populations of spin-label conformers against intradomain paramagnetic relaxation enhancements with a genetic algorithm. Subsequently, the protein structure ensemble is optimized using a more efficient genetic algorithm-based approach and an overfitting indicator, both of which were established in this work. The method was validated using a reference ensemble with a set of conformers whose populations and structures are known. This method was also applied to study the structure ensemble of the tandem di-domain of a poly (U) binding protein. The determined ensemble was supported by small-angle x-ray scattering and nuclear magnetic resonance relaxation data. The ensemble obtained suggests an induced fit mechanism for recognition of target RNA by the protein.

Molecular effects of cancer-associated somatic mutations on the structural and target recognition properties of Keap1.

Kelch-like ECH-associated protein 1 (Keap1) plays an important regulatory role in the nuclear factor erythroid 2-related factor 2 (Nrf2)-dependent oxidative stress response pathway. It functions as a repressor of Nrf2, a key transcription factor that initiates the expression of cytoprotective enzymes during oxidative stress to protect cells from damage caused by reactive oxygen species. Recent studies show that mutations of Keap1 can lead to aberrant activation of the antioxidant pathway, which is associated with different types of cancers. To gain a mechanistic understanding of the links between Keap1 mutations and cancer pathogenesis, we have investigated the molecular effects of a series of mutations (G333C, G350S, G364C, G379D, R413L, R415G, A427V, G430C and G476R) on the structural and target recognition properties of Keap1 by using nuclear magnetic resonance (NMR) spectroscopy, circular dichroism (CD) and isothermal titration calorimetry (ITC). Depending on their locations in the protein, these mutations are found to exert differential effects on the protein stability and target binding. Together with the proposed hinge-and-latch mechanism of Nrf2-Keap1 binding in the literature, our results provide important insight into the molecular affect of different somatic mutations on Keap1's function as an Nrf2 repressor.