NMRbox: A Resource for Biomolecular NMR Computation

Simple MATLAB and Python scripts for multi-exponential analysis

Multi-exponential decay is prevalent in magnetic resonance spectroscopy, relaxation, and imaging. This paper describes simple MATLAB and Python functions and scripts for regularized multi-exponential analysis methods for 1D and 2D data and example test problems and experiments. Regularized least-squares solutions provide production-quality outputs with robust stopping rules in ~5 and ~20 lines of code for 1D and 2D inversions, respectively. The software provides an open-architecture simple solution for transforming exponential decay data to the distribution of their decay lifetimes. Examples from magnetic resonance relaxation of a complex fluid, a Danish North Sea crude oil, and fluid mixtures in porous materials-brine/crude oil mixture in North Sea reservoir chalk-are presented. Developed codes may be incorporated in other software or directly used by other researchers, in magnetic resonance relaxation, diffusion, and imaging or other physical phenomena that require multi-exponential analysis.

Molecular Mechanism of pH-Induced Protrusion Configuration Switching in Piscine Betanodavirus Implies a Novel Antiviral Strategy

Many viruses contain surface spikes or protrusions that are essential for virus entry. These surface structures can thereby be targeted by antiviral drugs to treat viral infections. Nervous necrosis virus (NNV), a simple nonenveloped virus in the genus of betanodavirus, infects fish and damages aquaculture worldwide. NNV has 60 conspicuous surface protrusions, each comprising three protrusion domains (P-domain) of its capsid protein. NNV uses protrusions to bind to common receptors of sialic acids on the host cell surface to initiate its entry via the endocytic pathway. However, structural alterations of NNV in response to acidic conditions encountered during this pathway remain unknown, while detailed interactions of protrusions with receptors are unclear. Here, we used cryo-EM to discover that Grouper NNV protrusions undergo low-pH-induced compaction and resting. NMR and molecular dynamics (MD) simulations were employed to probe the atomic details. A solution structure of the P-domain at pH 7.0 revealed a long flexible loop (amino acids 311-330) and a pocket outlined by this loop. Molecular docking analysis showed that the N-terminal moiety of sialic acid inserted into this pocket to interact with conserved residues inside. MD simulations demonstrated that part of this loop converted to a β-strand under acidic conditions, allowing for P-domain trimerization and compaction. Additionally, a low-pH-favored conformation is attained for the linker connecting the P-domain to the NNV shell, conferring resting protrusions. Our findings uncover novel pH-dependent conformational switching mechanisms underlying NNV protrusion dynamics potentially utilized for facilitating NNV entry, providing new structural insights into complex NNV-host interactions with the identification of putative druggable hotspots on the protrusion.

The 1H, 15N, and 13C resonance assignments of a single-domain antibody against immunoglobulin G

Research on camelid-derived single-domain antibodies (sdAbs) has demonstrated their significant utility in diverse biotechnological applications, including therapy and diagnostic. This is largely due to their relative simplicity as monomeric proteins, ranging from 12 to 15 kDa, in contrast to immunoglobulin G (IgG) antibodies, which are glycosylated heterotetramers of 150-160 kDa. Single-domain antibodies exhibit high conformational stability and adopt the typical immunoglobulin domain fold, consisting of a two-layer sandwich of 7-9 antiparallel beta-strands. They contain three loops, known as complementary-determining regions (CDRs), which are assembled on the sdAb surface and are responsible for antigen recognition. The single-domain antibody examined in this study, sdAb-mrh-IgG, was engineered to recognize IgG from rats, mice, but it also weakly recognizes IgG from humans (Pleiner et al. 2018). A search of the Protein Data Bank revealed only one NMR structure of a single-domain antibody, which is unrelated to sdAb-mrh-IgG. The NMR chemical shift assignments of sdAb-mrh-IgG will be utilized to study its molecular dynamics and interactions with antigens in solution, which is fundamental for the rational design of novel single-domain antibodies.

The backbone NMR resonance assignments of the stabilized E. coli β clamp

The 81 kDa E. coli β clamp is a ring-shaped head-to-tail homodimer that encircles DNA and plays a central role in bacterial DNA replication by serving as a processivity factor for DNA polymerases and a binding platform for other DNA replication and repair proteins. Here we report the backbone 1H, 15N, and 13C NMR resonance assignments of the stabilized T45R/S107R β clamp variant obtained using standard TROSY-based triple-resonance experiments (BMRB 52548). The backbone assignments were aided by 13C and 15N edited NOESY experiments, allowing us to utilize our previously reported assignments of the β clamp ILV side-chain methyl groups (BMRB 51430, 51431). The backbone assignments of the T45R/S107R β clamp variant were transferred to the wild-type β clamp using a minimal set of TROSY-based 15N edited NOESY, NHCO and NHCA experiments (BMRB 52549). The reported backbone and previous ILV side-chain resonance assignments will enable NMR studies of the β clamp interactions and dynamics using amide and methyl groups as probes.

Mispacking of the F87 sidechain drives aggregation-promoting conformational fluctuations in the subunit interfaces of the transthyretin tetramer

Aberrant formation and deposition of human transthyretin (TTR) aggregates causes transthyretin amyloidosis. To initialize aggregation, transthyretin tetramers must first dissociate into monomers that partially unfold to promote entry into the aggregation pathway. The native TTR tetramer (T) is stabilized by docking of the F87 sidechain into an interfacial cavity enclosed by several hydrophobic residues including A120. We have previously shown that an alternative tetramer (T*) with mispacked F87 sidechains is more prone to dissociation and aggregation than the native T state. However, the molecular basis for the reduced stability in T* remains unclear. Here we report characterization of the A120L mutant, where steric hindrance is introduced into the F87 binding site. The x-ray structure of A120L shows that the F87 sidechain is displaced from its docking site across the subunit interface. In A120S, a naturally occurring pathogenic mutant that is less aggregation-prone than A120L, the F87 sidechain is correctly docked, as in the native TTR tetramer. Nevertheless, 19F-NMR aggregation assays show an elevated population of a monomeric aggregation intermediate in A120S relative to a control containing the native A120, due to accelerated tetramer dissociation and slowed monomer tetramerization. The mispacking of the F87 sidechain is associated with enhanced exchange dynamics for interfacial residues. At 298 K, the T* populations of various naturally occurring mutants fall between 4% and 7% (ΔG ~ 1.5-1.9 kcal/mol), consistent with the free energy change expected for undocking and solvent exposure of one of the four F87 sidechains in the tetramer (ΔG ~ 1.6 kcal/mol). Our data provide a molecular-level picture of the likely universal F87 sidechain mispacking in tetrameric TTR that promotes interfacial conformational dynamics and increases aggregation propensity.

Formerly degenerate seventh zinc finger domain from transcription factor ZNF711 rehabilitated by experimental NMR structure

Domain Z7 of nuclear transcription factor ZNF711 has the consensus last metal-ligand H23 found in odd-numbered zinc fingers of this protein replaced by a phenylalanine. Ever since the discovery of ZNF711, it has been thought that Z7 is probably non-functional because of the H23F substitution. The presence of H26 three positions downstream prompted us to examine if this histidine could substitute as the last metal-ligand. The Z7 domain adopts a stable tertiary structure upon metal-binding. The NMR structure of Zn2+-bound Z7 shows the classical ββα-fold of CCHH zinc fingers. Mutagenesis and pH titration experiments indicate that H26 is not involved in metal binding and that Z7 has a tridentate metal-binding site comprised of only residues C3, C6, and H19. By contrast, an F23H mutation that introduces a histidine in the consensus position forms a tetradentate ligand. The structure of the WT Z7 is stable causing restricted ring-flipping of phenylalanines 10 and 23. Dynamics are increased with either the H26A or F23H substitutions and aromatic ring rotation is no longer hindered in the two mutants. The mutations have only small effects on the Kd values for Zn2+ and Co2+ and retain the high thermal stability of the WT domain above 80°C. Like two previously reported designed zinc fingers with the last ligand replaced by water, the WT Z7 domain is catalytically active, hydrolyzing 4-nitrophenyl acetate. We discuss the implications of naturally occurring tridentate zinc fingers for cancer mutations and drug targeting of notoriously undruggable transcription factors.

CD28 hinge used in chimeric antigen receptor (CAR) T-cells exhibits local structure and conformational exchange amidst global disorder

T-cell therapies based on chimeric antigen receptor (CAR) targeting of a tumor-specific antigen offer hope for patients with relapsed or refractory cancers. CAR hinge and transmembrane regions link antigen recognition domains to intracellular signal transduction domains. Here, we apply biophysical methods to characterize the structure and dynamic properties of the CD28 CAR hinge (CD28H) used in an FDA-approved CD19 CAR for the treatment of B-lineage leukemia/lymphoma. By using nuclear Overhauser effect spectroscopy (NOESY), which detects even transiently occupied structural motifs, we observed otherwise elusive local structural elements amidst overall disorder in CD28H, including a conformational switch from a native β-strand to a 310-helix and polyproline II helix-like structure. These local structural motifs contribute to an overall loosely formed extended geometry that could be captured by NOESY data. All FDA-approved CARs use prolines in the hinge region, which we find in CD28, and previously in CD8α, isomerize to promote structural plasticity and dynamics. These local structural elements may function in recognition and signaling events and constrain the spacing between the transmembrane and antigen recognition domains. Our study thus demonstrates a method for detecting local and transient structure within intrinsically disordered systems and moreover, our CD28H findings may inform future CAR design.

1H, 15N, and 13C resonance assignments of the N-terminal domain and ser-arg-rich intrinsically disordered region of the nucleocapsid protein of the SARS-CoV-2

The nucleocapsid (N) protein of SARS-CoV-2 is a multifunctional protein involved in nucleocapsid assembly and various regulatory functions. It is the most abundant protein during viral infection. Its functionality is closely related to its structure, which comprises two globular domains, the N-terminal domain (NTD) and the C-terminal domain (CTD), flanked by intrinsically disordered regions. The linker between the NTD and CTD includes a Serine-Arginine rich (SR) region, which is crucial for the regulation of the N protein's function. Here, we report the near-complete assignment of the construct containing the NTD followed by the SR region (NTD-SR). Additionally, we describe the dynamic nature of the SR region and compare it with all other available chemical shift assignments reported for the SR region.

Insights into the interaction between UGGT, the gatekeeper of folding in the ER, and its partner, the selenoprotein SEP15

The enzyme UDP-glucose: glycoprotein glucosyltransferase (UGGT) is the gatekeeper of protein folding within the endoplasmic reticulum (ER). One-third of the human proteome traverses the ER where folding and maturation are facilitated by a complex protein homeostasis network. Both glycan modifications and disulfide bonds are of key importance in the maturation of these ER proteins. The actions of UGGT are intimately linked to the glycan code for folding and maturation of secretory proteins in the ER. UGGT selectively glucosylates the N-linked glycan of misfolded proteins so that they can reenter the lectin-folding chaperone cycle and be retained within the ER for further attempts at folding. An intriguing aspect of UGGT function is its interaction with its poorly understood cochaperone, the 15 kDa selenoprotein known as SELENOF or SEP15. This small protein contains a rare selenocysteine residue proposed to act as an oxidoreductase toward UGGT substrates. AlphaFold2 predictions of the UGGT1/SEP15 complex provide insight into this complex at a structural level. The predicted UGGT1/SEP15 interaction interface was validated by mutagenesis and coimmunoprecipitation experiments. These results serve as a springboard for models of the integrated action of UGGT1 and SEP15.

Positive residues of the SARS-CoV-2 fusion domain are key contributors to the initiation of membrane fusion

SARS-CoV-2 is one of the most infectious viruses ever recorded. Despite a plethora of research over the last several years, the viral life cycle is still not well understood, particularly membrane fusion. This process is initiated by the fusion domain (FD), a highly conserved stretch of amino acids consisting of a fusion peptide (FP) and fusion loop (FL), which in synergy perturbs the target cells' lipid membrane to lower the energetic cost necessary for fusion. In this study, through a mutagenesis-based approach, we have investigated the basic residues within the FD (K825, K835, R847, K854) utilizing an in vitro fusion assay and 19F NMR, validated by traditional 13C 15N techniques. Alanine and charge-conserving mutants revealed every basic residue plays a highly specific role within the mechanism of initiating fusion. Intriguingly, K825A led to increased fusogenecity which was found to be correlated to the number of amino acids within helix one, further implicating the role of this specific helix within the FD's fusion mechanism. This work has found basic residues to be important within the FDs fusion mechanism and highlights K825A, a specific mutation made within the FD of the SARS-CoV-2 spike protein, as requiring further investigation due to its potential to contribute to a more virulent strain of SARS-CoV-2.

Oligomerization and a distinct tRNA-binding loop are important regulators of human arginyl-transferase function

The arginyl-transferase ATE1 is a tRNA-dependent enzyme that covalently attaches an arginine molecule to a protein substrate. Conserved from yeast to humans, ATE1 deficiency in mice correlates with defects in cardiovascular development and angiogenesis and results in embryonic lethality, while conditional knockouts exhibit reproductive, developmental, and neurological deficiencies. Despite the recent revelation of the tRNA binding mechanism and the catalytic cycle of yeast ATE1, the structure-function relationship of ATE1 in higher organisms is not well understood. In this study, we present the three-dimensional structure of human ATE1 in an apo-state and in complex with its tRNA cofactor and a peptide substrate. In contrast to its yeast counterpart, human ATE1 forms a symmetric homodimer, which dissociates upon binding of a substrate. Furthermore, human ATE1 includes a unique and extended loop that wraps around tRNAArg, creating extensive contacts with the T-arm of the tRNA cofactor. Substituting key residues identified in the substrate binding site of ATE1 abolishes enzymatic activity and results in the accumulation of ATE1 substrates in cells.

Insights into mechanisms of MALT1 allostery from NMR and AlphaFold dynamic analyses

Mucosa-associated lymphoid tissue lymphoma-translocation protein 1 (MALT1) is an attractive target for the development of modulatory compounds in the treatment of lymphoma and other cancers. While the three-dimensional structure of MALT1 has been previously determined through X-ray analysis, its dynamic behaviour in solution has remained unexplored. We present here dynamic analyses of the apo MALT1 form along with the E549A mutation. This investigation used NMR 15N relaxation and NOE measurements between side-chain methyl groups. Our findings confirm that MALT1 exists as a monomer in solution, and demonstrate that the domains display semi-independent movements in relation to each other. Our dynamic study, covering multiple time scales, along with the assessment of conformational populations by Molecular Dynamic simulations, Alpha Fold modelling and PCA analysis, put the side chain of residue W580 in an inward position, shedding light at potential mechanisms underlying the allosteric regulation of this enzyme.

PyINETA: Open-source platform for INADEQUATE-JRES integration in NMR metabolomics

Annotating compounds with high confidence is a critical element in metabolomics. 13C-detection NMR experiment INADEQUATE (incredible natural abundance double-quantum transfer experiment) stands out as a powerful tool for structural elucidation, whereas this valuable experiment is not often included in metabolomics studies. This is partly due to the lack of community platform that provides structural information based INADEQUATE. Also, it is often the case that a single study uses various NMR experiments synergistically to improve the quality of information or balance total NMR experiment time, but there is no public platform that can integrate the outputs of INADEQUATE and other NMR experiments either. Here, we introduce PyINETA, Python-based INADEQUATE network analysis. PyINETA is an open-source platform that provides structural information of molecules using INADEQUATE, conducts database search, and integrates information of INADEQUATE and a complementary NMR experiment 13C J-resolved experiment (13C-JRES). Those steps are carried out automatically, and PyINETA keeps track of all the pipeline parameters and outputs, ensuring the transparency of annotation in metabolomics. Our evaluation of PyINETA using a model mouse study showed that our pipeline successfully integrated INADEQUATE and 13C-JRES. The results showed that 13C-labeled amino acids that were fed to mice were transferred to different tissues, and, also, they were transformed to other metabolites. The distribution of those compounds was tissue-specific, showing enrichment of particular metabolites in liver, spleen, pancreas, muscle, or lung. The value of PyINETA was not limited to those known compounds; PyINETA also provided fragment information for unknown compounds. PyINETA is available on NMRbox.

Thermodynamic coupling of the tandem RRM domains of hnRNP A1 underlie its pleiotropic RNA binding functions

The functional properties of RNA binding proteins (RBPs) require allosteric regulation through interdomain communication. Despite the importance of allostery to biological regulation, only a few studies have been conducted to describe the biophysical nature by which interdomain communication manifests in RBPs. Here, we show for hnRNP A1 that interdomain communication is vital for the unique stability of its amino-terminal domain, which consists of two RNA recognition motifs (RRMs). These RRMs exhibit drastically different stability under pressure. RRM2 unfolds as an individual domain but remains stable when appended to RRM1. Variants that disrupt interdomain communication between the tandem RRMs show a significant decrease in stability. Carrying these mutations over to the full-length protein for in vivo experiments revealed that the mutations affected the ability of the disordered carboxyl-terminal domain to engage in protein-protein interactions and influenced the protein's RNA binding capacity. Collectively, this work reveals that thermodynamic coupling between the tandem RRMs of hnRNP A1 accounts for its allosteric regulatory functions.

Structural transitions enable interleukin-18 maturation and signaling

Several interleukin-1 (IL-1) family members, including IL-1β and IL-18, require processing by inflammasome-associated caspases to unleash their activities. Here, we unveil, by cryoelectron microscopy (cryo-EM), two major conformations of the complex between caspase-1 and pro-IL-18. One conformation is similar to the complex of caspase-4 and pro-IL-18, with interactions at both the active site and an exosite (closed conformation), and the other only contains interactions at the active site (open conformation). Thus, pro-IL-18 recruitment and processing by caspase-1 is less dependent on the exosite than the active site, unlike caspase-4. Structure determination by nuclear magnetic resonance uncovers a compact fold of apo pro-IL-18, which is similar to caspase-1-bound pro-IL-18 but distinct from cleaved IL-18. Binding sites for IL-18 receptor and IL-18 binding protein are only formed upon conformational changes after pro-IL-18 cleavage. These studies show how pro-IL-18 is selected as a caspase-1 substrate, and why cleavage is necessary for its inflammatory activity.

Assignment of the Lassa virus transmembrane domain in the prefusion and postfusion states in detergent micelles

Lassa virus (LASV) is the most prevalent member of the arenavirus family and the causative agent of Lassa fever, a viral hemorrhagic fever. Although there are annual outbreaks in West Africa, and recently isolated cases worldwide, there are no current therapeutics or vaccines. As such, LASV poses a significant global public health threat. One of the key steps in LASV infection is delivering its genetic material by fusing its viral membrane with the host cell membrane. This process is facilitated by significant conformational changes within glycoprotein 2 (GP2), yielding distinct prefusion and postfusion structural states. However, structural information is missing to understand the changes that occur in the transmembrane domain (TM) during the fusion process. Previously, we showed that the TM undergoes pH-dependent structural changes that result in a helical extension. Here, we provide the 1H, 15N, and 13C assignment of the LASV TM backbone in the prefusion and postfusion states. We also provide the 1H, 15N, and 13C assignment of two mutants, G429P and D432P, which prevent this helical extension. These results will help understand the role the TM plays in membrane fusion and can lead to the design of therapeutics against LASV infection.

NMR 1H, 13C, 15N backbone resonance assignments of wild-type human K-Ras and its oncogenic mutants G12D and G12C bound to GTP

Human K-Ras protein, which is a member of the GTPase Ras family, hydrolyzes GTP to GDP and concomitantly converts from its active to its inactive state. It is a key oncoprotein, because several mutations, particularly those at residue position 12, occur with a high frequency in a wide range of human cancers. The K-Ras protein is therefore an important target for developing therapeutic anti-cancer agents. In this work we report the almost complete sequence-specific resonance assignments of wild-type and the oncogenic G12C and G12D mutants in the GTP-complexed active forms, including the functionally important Switch I and Switch II regions. These assignments serve as the basis for a comprehensive functional dynamics study of wild-type K-Ras and its G12 mutants.

Catalytic role of histidine-114 in the hydrolytic dehalogenation of chlorothalonil by Pseudomonas sp. CTN-3

Chlorothalonil (2,4,5,6-tetrachloroisophthalonitrile; TPN) is an environmentally persistent fungicide that sees heavy use in the USA and is highly toxic to aquatic species and birds, as well as a probable human carcinogen. The chlorothalonil dehalogenase from Pseudomonas sp. CTN-3 (Chd, UniProtKB C9EBR5) degrades TPN to its less toxic 4-OH-TPN analog making it an exciting candidate for the development of a bioremediation process for TPN; however, little is currently known about its catalytic mechanism. Therefore, an active site residue histidine-114 (His114) which forms a hydrogen bond with the Zn(II)-bound water/hydroxide and has been suggested to be the active site acid/base, was substituted by an Ala residue. Surprisingly, ChdH114A exhibited catalytic activity with a kcat value of 1.07 s-1, ~ 5% of wild-type (WT) Chd, and a KM of 32 µM. Thus, His114 is catalytically important but not essential. The electronic and structural aspects of the WT Chd and ChdH114A active sites were examined using UV-Vis and EPR spectroscopy on the catalytically competent Co(II)-substituted enzyme as well as all-atomistic molecular dynamics (MD) simulations. Combination of these data suggest His114 can quickly and reversibly move nearly 2 Å between one conformation that facilitates catalysis and another that enables product egress and active site recharge. In light of experimental and computational data on ChdH114A, Asn216 appears to play a role in substrate binding and preorganization of the transition-state while Asp116 likely facilitates the deprotonation of the Zn(II)-bound water in the absence of His114. Based on these data, an updated proposed catalytic mechanism for Chd is presented.

Effects of Xylanase A double mutation on substrate specificity and structural dynamics

While protein activity is traditionally studied with a major focus on the active site, the activity of enzymes has been hypothesized to be linked to the flexibility of adjacent regions, warranting more exploration into how the dynamics in these regions affects catalytic turnover. One such enzyme is Xylanase A (XylA), which cleaves hemicellulose xylan polymers by hydrolysis at internal β-1,4-xylosidic linkages. It contains a "thumb" region whose flexibility has been suggested to affect the activity. The double mutation D11F/R122D was previously found to affect activity and potentially bias the thumb region to a more open conformation. We find that the D11F/R122D double mutation shows substrate-dependent effects, increasing activity on the non-native substrate ONPX2 but decreasing activity on its native xylan substrate. To characterize how the double mutant causes these kinetics changes, nuclear magnetic resonance (NMR) and molecular dynamics (MD) simulations were used to probe structural and flexibility changes. NMR chemical shift perturbations revealed structural changes in the double mutant relative to the wild-type, specifically in the thumb and fingers regions. Increased slow-timescale dynamics in the fingers region was observed as intermediate-exchange line broadening. Lipari-Szabo order parameters show negligible changes in flexibility in the thumb region in the presence of the double mutation. To help understand if there is increased energetic accessibility to the open state upon mutation, alchemical free energy simulations were employed that indicated thumb opening is more favorable in the double mutant. These studies aid in further characterizing how flexibility in adjacent regions affects the function of XylA.

Structural basis of substrate recognition and allosteric activation of the proapoptotic mitochondrial HtrA2 protease

The mitochondrial serine protease HtrA2 is a human homolog of the Escherichia coli Deg-proteins exhibiting chaperone and proteolytic roles. HtrA2 is involved in both apoptotic regulation via its ability to degrade inhibitor-of-apoptosis proteins (IAPs), as well as in cellular maintenance as part of the cellular protein quality control machinery, by preventing the possible toxic accumulation of aggregated proteins. In this study, we use advanced solution NMR spectroscopy methods combined with biophysical characterization and biochemical assays to elucidate the crucial role of the substrate recognizing PDZ domain. This domain regulates the protease activity of HtrA2 by triggering an intricate allosteric network involving the regulatory loops of the protease domain. We further show that divalent metal ions can both positively and negatively modulate the activity of HtrA2, leading to a refined model of HtrA2 regulation within the apoptotic pathway.

Structural basis of nearest-neighbor cooperativity in the ring-shaped gene regulatory protein TRAP from protein engineering and cryo-EM

Homotropic cooperativity is widespread in biological regulation. The homo-oligomeric ring-shaped trp RNA binding attenuation protein (TRAP) from bacillus binds multiple tryptophan ligands (Trp) and becomes activated to bind a specific sequence in the 5' leader region of the trp operon mRNA. Ligand-activated binding to this specific RNA sequence regulates downstream biosynthesis of Trp in a feedback loop. Characterized TRAP variants form 11- or 12-mer rings and bind Trp at the interface between adjacent subunits. Various studies have shown that a pair of loops that gate each Trp binding site is flexible in the absence of the ligand and become ordered upon ligand binding. Thermodynamic measurements of Trp binding have revealed a range of cooperative behavior for different TRAP variants, even if the averaged apparent affinities for Trp have been found to be similar. Proximity between the ligand binding sites, and the ligand-coupled disorder-to-order transition has implicated nearest-neighbor interactions in cooperativity. To establish a solid basis for describing nearest-neighbor cooperativity we engineered dodecameric (12-mer) TRAP variants constructed with two subunits connected by a flexible linker (dTRAP). We mutated one of the protomers such that only every other site was competent for Trp binding. Thermodynamic and structural studies using native mass spectrometry, NMR spectroscopy, and cryo-EM provided unprecedented detail into the thermodynamic and structural basis for the observed ligand binding cooperativity. Such insights can be useful for understanding allosteric control networks and for the development of new ones with defined ligand sensitivity and regulatory control.

NMR characterization of RNA binding property of the DEAD-box RNA helicase DDX3X and its implications for helicase activity

The DEAD-box RNA helicase (DDX) plays a central role in many aspects of RNA metabolism by remodeling the defined structure of RNA molecules. While a number of structural studies have revealed the atomistic details of the interaction between DDX and RNA ligands, the molecular mechanism of how this molecule unwinds a structured RNA into an unstructured single-stranded RNA (ssRNA) has largely remained elusive. This is due to challenges in structurally characterizing the unwinding intermediate state and the lack of thermodynamic details underlying this process. In this study, we use solution nuclear magnetic resonance (NMR) spectroscopy to characterize the interaction of human DDX3X, a member of the DDX family, with various RNA ligands. Our results show that the inherent binding affinity of DDX3X for ssRNA is significantly higher than that for structured RNA elements. This preferential binding, accompanied by the formation of a domain-closed conformation in complex with ssRNA, effectively stabilizes the denatured ssRNA state and thus underlies the unwinding activity of DDX3X. Our results provide a thermodynamic and structural basis for the DDX function, whereby DDX can recognize and remodel a distinct set of structured RNAs to participate in a wide range of physiological processes.

Toward More Selective Antibiotic Inhibitors: A Structural View of the Complexed Binding Pocket of E. coli Peptide Deformylase

Peptide deformylase (PDF) is involved in bacterial protein maturation processes. Originating from the interest in a new antibiotic, tremendous effort was put into the refinement of PDF inhibitors (PDFIs) and their selectivity. We obtained a full NMR backbone assignment the emergent additional protein backbone resonances of ecPDF 1-147 in complex with 2-(5-bromo-1H-indol-3-yl)-N-hydroxyacetamide (2), a potential new structural scaffold for more selective PDFIs. We also determined the complex crystal structures of E. coli PDF (ecPDF fl) and 2. Our structure suggests an alternative ligand conformation within the protein, a possible starting point for further selectivity optimization. The orientation of the second ligand conformation in the crystal structure points toward a small region of the S1' pocket, which differs between bacterial PDFs and human PDF. Moreover, we analyzed the binding mode of 2 via NMR TITAN line shape analysis, revealing an induced fit mechanism.

Glutamylation of Npm2 and Nap1 acidic disordered regions increases DNA mimicry and histone chaperone efficiency

Histone chaperones-structurally diverse, non-catalytic proteins enriched with acidic intrinsically disordered regions (IDRs)-protect histones from spurious nucleic acid interactions and guide their deposition into and out of nucleosomes. Despite their conservation and ubiquity, the function of the chaperone acidic IDRs remains unclear. Here, we show that the Xenopus laevis Npm2 and Nap1 acidic IDRs are substrates for TTLL4 (Tubulin Tyrosine Ligase Like 4)-catalyzed post-translational glutamate-glutamylation. We demonstrate that to bind, stabilize, and deposit histones into nucleosomes, chaperone acidic IDRs function as DNA mimetics. Our biochemical, computational, and biophysical studies reveal that glutamylation of these chaperone polyelectrolyte acidic stretches functions to enhance DNA electrostatic mimicry, promoting the binding and stabilization of H2A/H2B heterodimers and facilitating nucleosome assembly. This discovery provides insights into both the previously unclear function of the acidic IDRs and the regulatory role of post-translational modifications in chromatin dynamics.

Formerly degenerate seventh zinc finger domain from transcription factor ZNF711 rehabilitated by experimental NMR structure

Domain Z7 of nuclear transcription factor ZNF711 has the consensus last metal-ligand H23 found in odd-numbered zinc-fingers of this protein replaced by a phenylalanine. Ever since the discovery of ZNF711 it has been thought that Z7 is probably non-functional because of the H23F substitution. The presence of H26 three positions downstream prompted us to examine if this histidine could substitute as the last metal ligand. The Z7 domain adopts a stable tertiary structure upon metal binding. The NMR structure of Zn2+-bound Z7 shows the classical ββα-fold of CCHH zinc fingers. Mutagenesis and pH titration experiments indicate that H26 is not involved in metal binding and that Z7 has a tridentate metal-binding site comprised of only residues C3, C6, and H19. By contrast, an F23H mutation that introduces a histidine in the consensus position forms a tetradentate ligand. The structure of the WT Z7 is stable causing restricted ring-flipping of phenyalanines 10 and 23. Dynamics are increased with either the H26A or F23H substitutions and aromatic ring rotation is no longer hindered in the two mutants. The mutations have only small effects on the Kd values for Zn2+ and Co2+ and retain the high thermal stability of the WT domain above 80 °C. Like two previously reported designed zinc fingers with the last ligand replaced by water, the WT Z7 domain is catalytically active, hydrolyzing 4-nitophenyl acetate. We discuss the implications of naturally occurring tridentate zinc fingers for cancer mutations and drug targeting of notoriously undruggable transcription factors. Our findings that Z7 can fold with only a subset of three metal ligands suggests the recent view that most everything about protein structure can be predicted through homology modeling might be premature for at least the resilient and versatile zinc-finger motif.

Cancer-associated polybromo-1 bromodomain 4 missense variants variably impact bromodomain ligand binding and cell growth suppression

The polybromo, brahma-related gene 1-associated factors (PBAF) chromatin remodeling complex subunit polybromo-1 (PBRM1) contains six bromodomains that recognize and bind acetylated lysine residues on histone tails and other nuclear proteins. PBRM1 bromodomains thus provide a link between epigenetic posttranslational modifications and PBAF modulation of chromatin accessibility and transcription. As a putative tumor suppressor in several cancers, PBRM1 protein expression is often abrogated by truncations and deletions. However, ∼33% of PBRM1 mutations in cancer are missense and cluster within its bromodomains. Such mutations may generate full-length PBRM1 variant proteins with undetermined structural and functional characteristics. Here, we employed computational, biophysical, and cellular assays to interrogate the effects of PBRM1 bromodomain missense variants on bromodomain stability and function. Since mutations in the fourth bromodomain of PBRM1 (PBRM1-BD4) comprise nearly 20% of all cancer-associated PBRM1 missense mutations, we focused our analysis on PBRM1-BD4 missense protein variants. Selecting 16 potentially deleterious PBRM1-BD4 missense protein variants for further study based on high residue mutational frequency and/or conservation, we show that cancer-associated PBRM1-BD4 missense variants exhibit varied bromodomain stability and ability to bind acetylated histones. Our results demonstrate the effectiveness of identifying the unique impacts of individual PBRM1-BD4 missense variants on protein structure and function, based on affected residue location within the bromodomain. This knowledge provides a foundation for drawing correlations between specific cancer-associated PBRM1 missense variants and distinct alterations in PBRM1 function, informing future cancer personalized medicine approaches.

Efficient and economic protein labeling for NMR in mammalian expression systems: Application to a preT-cell and T-cell receptor protein

Protein nuclear magnetic resonance (NMR) spectroscopy relies on the ability to isotopically label polypeptides, which is achieved through heterologous expression in various host organisms. Most commonly, Escherichia coli is employed by leveraging isotopically substituted ammonium and glucose to uniformly label proteins with 15N and 13C, respectively. Moreover, E. coli can grow and express proteins in uniformly deuterium-substituted water (D2O), a strategy useful for experiments targeting high molecular weight proteins. Unfortunately, many proteins, particularly those requiring specific posttranslational modifications like disulfide bonding or glycosylation for proper folding and/or function, cannot be readily expressed in their functional forms using E. coli-based expression systems. One such class of proteins includes T-cell receptors and their related preT-cell receptors. In this study, we present an expression system for isotopic labeling of proteins using a nonadherent human embryonic kidney cell line, Expi293F, and a specially designed media. We demonstrate the application of this platform to the β subunit common to both receptors. In addition, we show that this expression system and media can be used to specifically label amino acids Phe, Ile, Val, and Leu in this system, utilizing an amino acid-specific labeling protocol that allows targeted incorporation at high efficiency without significant isotopic scrambling. We demonstrate that this system can also be used to express proteins with fluorinated amino acids. We were routinely able to obtain an NMR sample with a concentration of 200 μM from 30 mL of culture media, utilizing less than 20 mg of the labeled amino acids.

Structure of alpha-synuclein fibrils derived from human Lewy body dementia tissue

The defining feature of Parkinson disease (PD) and Lewy body dementia (LBD) is the accumulation of alpha-synuclein (Asyn) fibrils in Lewy bodies and Lewy neurites. Here we develop and validate a method to amplify Asyn fibrils extracted from LBD postmortem tissue samples and use solid state nuclear magnetic resonance (SSNMR) studies to determine atomic resolution structure. Amplified LBD Asyn fibrils comprise a mixture of single protofilament and two protofilament fibrils with very low twist. The protofilament fold is highly similar to the fold determined by a recent cryo-electron microscopy study for a minority population of twisted single protofilament fibrils extracted from LBD tissue. These results expand the structural characterization of LBD Asyn fibrils and approaches for studying disease mechanisms, imaging agents and therapeutics targeting Asyn.

Perspective: use and reuse of NMR-based metabolomics data: what works and what remains challenging

BACKGROUND

The National Cancer Institute issued a Request for Information (RFI; NOT-CA-23-007) in October 2022, soliciting input on using and reusing metabolomics data. This RFI aimed to gather input on best practices for metabolomics data storage, management, and use/reuse.

AIM OF REVIEW

The nuclear magnetic resonance (NMR) Interest Group within the Metabolomics Association of North America (MANA) prepared a set of recommendations regarding the deposition, archiving, use, and reuse of NMR-based and, to a lesser extent, mass spectrometry (MS)-based metabolomics datasets. These recommendations were built on the collective experiences of metabolomics researchers within MANA who are generating, handling, and analyzing diverse metabolomics datasets spanning experimental (sample handling and preparation, NMR/MS metabolomics data acquisition, processing, and spectral analyses) to computational (automation of spectral processing, univariate and multivariate statistical analysis, metabolite prediction and identification, multi-omics data integration, etc.) studies.

KEY SCIENTIFIC CONCEPTS OF REVIEW

We provide a synopsis of our collective view regarding the use and reuse of metabolomics data and articulate several recommendations regarding best practices, which are aimed at encouraging researchers to strengthen efforts toward maximizing the utility of metabolomics data, multi-omics data integration, and enhancing the overall scientific impact of metabolomics studies.

The structural properties of full-length annexin A11

Annexin A11 (ANXA11) is a calcium-dependent phospholipid-binding protein belonging to the annexin protein family and implicated in the neurodegenerative amyotrophic lateral sclerosis. Structurally, ANXA11 contains a conserved calcium-binding C-terminal domain common to all annexins and a putative intrinsically unfolded N-terminus specific for ANXA11. Little is known about the structure and functions of this region of the protein. By analogy with annexin A1, it was suggested that residues 38 to 59 within the ANXA11 N-terminus could form a helical region that would be involved in interactions. Interestingly, this region contains residues that, when mutated, may lead to clinical manifestations. In the present study, we have studied the structural features of the full-length protein with special attention to the N-terminal region using a combination of biophysical techniques which include nuclear magnetic resonance and small angle X-ray scattering. We show that the N-terminus is intrinsically disordered and that the overall features of the protein are not markedly affected by the presence of calcium. We also analyzed the 38-59 helix hypothesis using synthetic peptides spanning both the wild-type sequence and clinically relevant mutations. We show that the peptides have a remarkable character typical of a native helix and that mutations do not alter the behaviour suggesting that they are required for interactions rather than being structurally important. Our work paves the way to a more thorough understanding of the ANXA11 functions.

Mispacking of the F87 sidechain drives aggregation-promoting conformational fluctuations in the subunit interfaces of the transthyretin tetramer

Aberrant formation and deposition of human transthyretin (TTR) aggregates causes transthyretin amyloidosis. To initialize aggregation, transthyretin tetramers must first dissociate into monomers that partially unfold to promote entry into the aggregation pathway. The native TTR tetramer (T) is stabilized by docking of the F87 sidechain into an interfacial cavity enclosed by several hydrophobic residues including A120. We have previously shown that an alternative tetramer (T*) with mispacked F87 sidechains is more prone to dissociation and aggregation than the native T state. However, the molecular basis for the reduced stability in T* remains unclear. Here we report characterization of the A120L mutant, where steric hindrance is introduced into the F87 binding site. The X-ray structure of A120L shows that the F87 sidechain is displaced from its docking site across the subunit interface. In A120S, a naturally occurring pathogenic mutant that is less aggregation-prone than A120L, the F87 sidechain is correctly docked, as in the native TTR tetramer. Nevertheless, 19F-NMR aggregation assays show an elevated population of a monomeric aggregation intermediate in A120S relative to a control containing the native A120, due to accelerated tetramer dissociation and slowed monomer tetramerization. The mispacking of the F87 sidechain is associated with enhanced exchange dynamics for interfacial residues. At 298 K, the T* populations of various naturally occurring mutants fall between 4-7% (ΔG ~ 1.5-1.9 kcal/mol), consistent with the free energy change expected for undocking and solvent exposure of one of the four F87 sidechains in the tetramer (ΔG ~ 1.6 kcal/mol). Our data provide a molecular-level picture of the likely universal F87 sidechain mispacking in tetrameric TTR that promotes interfacial conformational dynamics and increases aggregation propensity.

SAND: Automated Time-Domain Modeling of NMR Spectra Applied to Metabolite Quantification

Developments in untargeted nuclear magnetic resonance (NMR) metabolomics enable the profiling of thousands of biological samples. The exploitation of this rich source of information requires a detailed quantification of spectral features. However, the development of a consistent and automatic workflow has been challenging because of extensive signal overlap. To address this challenge, we introduce the software Spectral Automated NMR Decomposition (SAND). SAND follows on from the previous success of time-domain modeling and automatically quantifies entire spectra without manual interaction. The SAND approach uses hybrid optimization with Markov chain Monte Carlo methods, employing subsampling in both time and frequency domains. In particular, SAND randomly divides the time-domain data into training and validation sets to help avoid overfitting. We demonstrate the accuracy of SAND, which provides a correlation of ∼0.9 with ground truth on cases including highly overlapped simulated data sets, a two-compound mixture, and a urine sample spiked with different amounts of a four-compound mixture. We further demonstrate an automated annotation using correlation networks derived from SAND decomposed peaks, and on average, 74% of peaks for each compound can be recovered in single clusters. SAND is available in NMRbox, the cloud computing environment for NMR software hosted by the Network for Advanced NMR (NAN). Since the SAND method uses time-domain subsampling (i.e., random subset of time-domain points), it has the potential to be extended to a higher dimensionality and nonuniformly sampled data.

Structure-function relationships in ShKT domain peptides: ShKT-Ts1 from the sea anemone Telmatactis stephensoni

Diverse structural scaffolds have been described in peptides from sea anemones, with the ShKT domain being a common scaffold first identified in ShK toxin from Stichodactyla helianthus. ShK is a potent blocker of voltage-gated potassium channels (KV 1.x), and an analog, ShK-186 (dalazatide), has completed Phase 1 clinical trials in plaque psoriasis. The ShKT domain has been found in numerous other species, but only a tiny fraction of ShKT domains has been characterized functionally. Despite adopting the canonical ShK fold, some ShKT peptides from sea anemones inhibit KV 1.x, while others do not. Mutagenesis studies have shown that a Lys-Tyr (KY) dyad plays a key role in KV 1.x blockade, although a cationic residue followed by a hydrophobic residue may also suffice. Nevertheless, ShKT peptides displaying an ShK-like fold and containing a KY dyad do not necessarily block potassium channels, so additional criteria are needed to determine whether new ShKT peptides might show activity against potassium channels. In this study, we used a combination of NMR and molecular dynamics (MD) simulations to assess the potential activity of a new ShKT peptide. We determined the structure of ShKT-Ts1, from the sea anemone Telmatactis stephensoni, examined its tissue localization, and investigated its activity against a range of ion channels. As ShKT-Ts1 showed no activity against KV 1.x channels, we used MD simulations to investigate whether solvent exposure of the dyad residues may be informative in rationalizing and potentially predicting the ability of ShKT peptides to block KV 1.x channels. We show that either a buried dyad that does not become exposed during MD simulations, or a partially exposed dyad that becomes buried during MD simulations, correlates with weak or absent activity against KV 1.x channels. Therefore, structure determination coupled with MD simulations, may be used to predict whether new sequences belonging to the ShKT family may act as potassium channel blockers.

Interpreting the effect of mutations to protein binding sites from large-scale genomic screens

Predicting the functionality of missense mutations is extremely difficult. Large-scale genomic screens are commonly performed to identify mutational correlates or drivers of disease and treatment resistance, but interpretation of how these mutations impact protein function is limited. One such consequence of mutations to a protein is to impact its ability to bind and interact with partners or small molecules such as ATP, thereby modulating its function. Multiple methods exist for predicting the impact of a single mutation on protein-protein binding energy, but it is difficult in the context of a genomic screen to understand if these mutations with large impacts on binding are more common than statistically expected. We present a methodology for taking mutational data from large-scale genomic screens and generating functional and statistical insights into their role in the binding of proteins both with each other and their small molecule ligands. This allows a quantitative and statistical analysis to determine whether mutations impacting protein binding or ligand interactions are occurring more or less frequently than expected by chance. We achieve this by calculating the potential impact of any possible mutation and comparing an expected distribution to the observed mutations. This method is applied to examples demonstrating its ability to interpret mutations involved in protein-protein binding, protein-DNA interactions, and the evolution of therapeutic resistance.

Visualizing a two-state conformational ensemble in stem-loop 3 of the transcriptional regulator 7SK RNA

Structural plasticity is integral to RNA function; however, there are currently few methods to quantitatively resolve RNAs that have multiple structural states. NMR spectroscopy is a powerful approach for resolving conformational ensembles but is size-limited. Chemical probing is well-suited for large RNAs but provides limited structural and kinetics information. Here, we integrate the two approaches to visualize a two-state conformational ensemble for the central stem-loop 3 (SL3) of 7SK RNA, a critical element for 7SK RNA function in transcription regulation. We find that the SL3 distal end exchanges between two equally populated yet structurally distinct states in both isolated SL3 constructs and full-length 7SK RNA. We rationally designed constructs that lock SL3 into a single state and demonstrate that both chemical probing and NMR data fit to a linear combination of the two states. Comparison of vertebrate 7SK RNA sequences shows either or both states are highly conserved. These results provide new insights into 7SK RNA structural dynamics and demonstrate the utility of integrating chemical probing with NMR spectroscopy to gain quantitative insights into RNA conformational ensembles.

Toward the mechanism of jarastatin (rJast) inhibition of the integrin αVβ3

Disintegrins are a family of cysteine-rich small proteins that were first identified in snake venom. The high divergence of disintegrins gave rise to a plethora of functions, all related to the interaction with integrins. Disintegrins evolved to interact selectively with different integrins, eliciting many physiological outcomes and being promising candidates for the therapy of many pathologies. We used NMR to determine the structure and dynamics of the recombinant disintegrin jarastatin (rJast) and its interaction with the cancer-related integrin αVβ3. rJast displayed the canonical fold of a medium-sized disintegrin and showed complex dynamic in multiple timescales. We used NMR experiments to map the interaction of rJast with αVβ3, and molecular docking followed by molecular dynamics (MD) simulation to describe the first structural model of a disintegrin/integrin complex. We showed that not only the RGD loop participates in the interaction, but also the N-terminal domain. rJast plasticity was essential for the interaction with αVβ3 and correlated with the main modes of motion depicted in the MD trajectories. In summary, our study provides novel structural insights that enhance our comprehension of the mechanisms underlying disintegrin functionality.

pH-dependent conformational change within the Lassa virus transmembrane domain elicits efficient membrane fusion

Lassa virus (LASV) is the most prevalent member of the arenavirus family and the causative agent of Lassa fever, a viral hemorrhagic fever. Although there are annual outbreaks in West Africa, and recently isolated cases worldwide, there are no current therapeutics or vaccines, which poses LASV as a significant global public health threat. One of the key steps in LASV infection is the delivery of its genetic material by fusing its viral membrane with the host cell membrane. This process is facilitated by significant conformational changes within glycoprotein 2 (GP2), yielding distinct prefusion and postfusion structural states. However, structural information is missing to understand the changes that occur in the transmembrane domain during the fusion process. Here, we used CD and NMR spectroscopy to show that the transmembrane domain has pH-dependent conformational changes that result in an extension of the alpha helix at the N-terminal end. Proline mutants of key residues in that region prevent the helical extension, as seen in CD and NMR. We developed a modified lipid mixing assay to study the importance of this extension on the function of GP2. Our assay shows that membrane fusion efficiency is optimal at low pH values but introducing the proline mutants results in lower fusion efficiency. These results indicate that these pH-dependent conformational changes are important to the fusion mechanism. This information can be used to design therapeutics to combat Lassa virus infections and prevent its potential spread.

Structure, dynamics, and stability of the smallest and most complex 71 protein knot

Proteins can spontaneously tie a variety of intricate topological knots through twisting and threading of the polypeptide chains. Recently developed artificial intelligence algorithms have predicted several new classes of topological knotted proteins, but the predictions remain to be authenticated experimentally. Here, we showed by X-ray crystallography and solution-state NMR spectroscopy that Q9PR55, an 89-residue protein from Ureaplasma urealyticum, possesses a novel 71 knotted topology that is accurately predicted by AlphaFold 2, except for the flexible N terminus. Q9PR55 is monomeric in solution, making it the smallest and most complex knotted protein known to date. In addition to its exceptional chemical stability against urea-induced unfolding, Q9PR55 is remarkably robust to resist the mechanical unfolding-coupled proteolysis by a bacterial proteasome, ClpXP. Our results suggest that the mechanical resistance against pulling-induced unfolding is determined by the complexity of the knotted topology rather than the size of the molecule.

Structural and functional characterisation of Tst2, a novel TRPV1 inhibitory peptide from the Australian sea anemone Telmatactis stephensoni

Sea anemone venoms are complex mixtures of biologically active compounds, including disulfide-rich peptides, some of which have found applications as research tools, and others as therapeutic leads. Our recent transcriptomic and proteomic studies of the Australian sea anemone Telmatactis stephensoni identified a transcript for a peptide designated Tst2. Tst2 is a 38-residue peptide showing sequence similarity to peptide toxins known to interact with a range of ion channels (NaV, TRPV1, KV and CaV). Recombinant Tst2 (rTst2, which contains an additional Gly at the N-terminus) was produced by periplasmic expression in Escherichia coli, enabling the production of both unlabelled and uniformly 13C,15N-labelled peptide for functional assays and structural studies. The LC-MS profile of the recombinant Tst2 showed a pure peak with molecular mass 6 Da less than that of the reduced form of the peptide, indicating the successful formation of three disulfide bonds from its six cysteine residues. The solution structure of rTst2 was determined using multidimensional NMR spectroscopy and revealed that rTst2 adopts an inhibitor cystine knot (ICK) structure. rTst2 was screened using various functional assays, including patch-clamp electrophysiological and cytotoxicity assays. rTst2 was inactive against voltage-gated sodium channels (NaV) and the human voltage-gated proton (hHv1) channel. rTst2 also did not possess cytotoxic activity when assessed against Drosophila melanogaster flies. However, the recombinant peptide at 100 nM showed >50% inhibition of the transient receptor potential subfamily V member 1 (TRPV1) and slight (∼10%) inhibition of transient receptor potential subfamily A member 1 (TRPA1). Tst2 is thus a novel ICK inhibitor of the TRPV1 channel.

REDEN: Interactive multi-fitting decomposition-based NMR peak picking assistant

We present a new program REDEN (Residual Decomposition of NMR peaks) designed to perform identification of peaks in NMR spectra. This integrated, cross-platform, open-source software visually assists with explicit peak picking through decomposition of NMR peaks on the frequency domain data. It provides a distinctive interactive workflow with iPick due to its integration with the POKY suite, providing users with a seamless and efficient experience. The decomposition of peaks operates in a chosen region of an NMR spectrum by multi-fitting simulated peaks with four lineshape fitting options as support, Gaussian, Lorentzian, a fast/optimized Lorentzian, and Pseudo-Voigt. Furthermore, REDEN provides a way to fine-tune for the users in two operating modes (Basic and Advanced). REDEN is pre-built in the POKY suite, which is available from https://poky.clas.ucdenver.edu.

New insights into the structure and function of the complex between the Escherichia coli Hsp70, DnaK, and its nucleotide-exchange factor, GrpE

The 70 kDa heat shock proteins (Hsp70s) play a pivotal role in many cellular functions using allosteric communication between their nucleotide-binding domain (NBD) and substrate-binding domain, mediated by an interdomain linker, to modulate their affinity for protein clients. Critical to modulation of the Hsp70 allosteric cycle, nucleotide-exchange factors (NEFs) act by a conserved mechanism involving binding to the ADP-bound NBD and opening of the nucleotide-binding cleft to accelerate the release of ADP and binding of ATP. The crystal structure of the complex between the NBD of the Escherichia coli Hsp70, DnaK, and its NEF, GrpE, was reported previously, but the GrpE in the complex carried a point mutation (G122D). Both the functional impact of this mutation and its location on the NEF led us to revisit the DnaK NBD/GrpE complex structurally using AlphaFold modeling and validation by solution methods that report on protein conformation and mutagenesis. This work resulted in a new model for the DnaK NBD in complex with GrpE in which subdomain IIB of the NBD rotates more than in the crystal structure, resulting in an open conformation of the nucleotide-binding cleft, which now resembles more closely what is seen in other Hsp/NEF complexes. Moreover, the new model is consistent with the increased ADP off-rate accompanying GrpE binding. Excitingly, our findings point to an interdomain allosteric signal in DnaK triggered by GrpE binding.

Predicting multiple conformations via sequence clustering and AlphaFold2

AlphaFold2 (ref. 1) has revolutionized structural biology by accurately predicting single structures of proteins. However, a protein's biological function often depends on multiple conformational substates2, and disease-causing point mutations often cause population changes within these substates3,4. We demonstrate that clustering a multiple-sequence alignment by sequence similarity enables AlphaFold2 to sample alternative states of known metamorphic proteins with high confidence. Using this method, named AF-Cluster, we investigated the evolutionary distribution of predicted structures for the metamorphic protein KaiB5 and found that predictions of both conformations were distributed in clusters across the KaiB family. We used nuclear magnetic resonance spectroscopy to confirm an AF-Cluster prediction: a cyanobacteria KaiB variant is stabilized in the opposite state compared with the more widely studied variant. To test AF-Cluster's sensitivity to point mutations, we designed and experimentally verified a set of three mutations predicted to flip KaiB from Rhodobacter sphaeroides from the ground to the fold-switched state. Finally, screening for alternative states in protein families without known fold switching identified a putative alternative state for the oxidoreductase Mpt53 in Mycobacterium tuberculosis. Further development of such bioinformatic methods in tandem with experiments will probably have a considerable impact on predicting protein energy landscapes, essential for illuminating biological function.

Solution Structure of Poly(UG) RNA

Poly(UG) or "pUG" RNAs are UG or GU dinucleotide repeat sequences which are highly abundant in eukaryotes. Post-transcriptional addition of pUGs to RNA 3' ends marks mRNAs as vectors for gene silencing in C. elegans. We previously determined the crystal structure of pUG RNA bound to the ligand N-methyl mesoporphyrin IX (NMM), but the structure of free pUG RNA is unknown. Here we report the solution structure of the free pUG RNA (GU)12, as determined by nuclear magnetic resonance spectroscopy and small and wide-angle x-ray scattering (NMR-SAXS-WAXS). The low complexity sequence and 4-fold symmetry of the structure result in overlapped NMR signals that complicate chemical shift assignment. We therefore utilized single site-specific deoxyribose modifications which did not perturb the structure and introduced well-resolved methylene signals that are easily identified in NMR spectra. The solution structure ensemble has a root mean squared deviation (RMSD) of 0.62 Å and is a compact, left-handed quadruplex with a Z-form backbone, or "pUG fold." Overall, the structure agrees with the crystal structure of (GU)12 bound to NMM, indicating the pUG fold is unaltered by docking of the NMM ligand. The solution structure reveals conformational details that could not be resolved by x-ray crystallography, which explain how the pUG fold can form within longer RNAs.

Legionella pneumophila macrophage infectivity potentiator protein appendage domains modulate protein dynamics and inhibitor binding

Macrophage infectivity potentiator (MIP) proteins are widespread in human pathogens including Legionella pneumophila, the causative agent of Legionnaires' disease and protozoans such as Trypanosoma cruzi. All MIP proteins contain a FKBP (FK506 binding protein)-like prolyl-cis/trans-isomerase domain that hence presents an attractive drug target. Some MIPs such as the Legionella pneumophila protein (LpMIP) have additional appendage domains of mostly unknown function. In full-length, homodimeric LpMIP, the N-terminal dimerization domain is linked to the FKBP-like domain via a long, free-standing stalk helix. Combining X-ray crystallography, NMR and EPR spectroscopy and SAXS, we elucidated the importance of the stalk helix for protein dynamics and inhibitor binding to the FKBP-like domain and bidirectional crosstalk between the different protein regions. The first comparison of a microbial MIP and a human FKBP in complex with the same synthetic inhibitor was made possible by high-resolution structures of LpMIP with a [4.3.1]-aza-bicyclic sulfonamide and provides a basis for designing pathogen-selective inhibitors. Through stereospecific methylation, the affinity of inhibitors to L. pneumophila and T. cruzi MIP was greatly improved. The resulting X-ray inhibitor-complex structures of LpMIP and TcMIP at 1.49 and 1.34 Å, respectively, provide a starting point for developing potent inhibitors against MIPs from multiple pathogenic microorganisms.

Backbone and sidechain NMR assignments of residues 1-81 from yeast Sis1 in complex with an Hsp70 C-terminal EEVD peptide

Molecular chaperones aid proteins to fold and assemble without modifying their final structure, requiring, in several folding processes, the interplay between members of the Hsp70 and Hsp40 families. Here, we report the NMR chemical shift assignments for 1 H, 15 N, and 13 C nuclei of the backbone and side chains of the J-domain of the class B Hsp40 from Saccharomyces cerevisiae, Sis1, complexed with the C-terminal EEVD motif of Hsp70. The data revealed information on the structure and backbone dynamics that add significantly to the understanding of the J-domain-Hsp70-EEVD mechanism of interaction.

Human Fis1 directly interacts with Drp1 in an evolutionarily conserved manner to promote mitochondrial fission

Mitochondrial fission protein 1 (Fis1) and dynamin-related protein 1 (Drp1) are the only two proteins evolutionarily conserved for mitochondrial fission, and directly interact in Saccharomyces cerevisiae to facilitate membrane scission. However, it remains unclear if a direct interaction is conserved in higher eukaryotes as other Drp1 recruiters, not present in yeast, are known. Using NMR, differential scanning fluorimetry, and microscale thermophoresis, we determined that human Fis1 directly interacts with human Drp1 (KD = 12-68 μM), and appears to prevent Drp1 assembly, but not GTP hydrolysis. Similar to yeast, the Fis1-Drp1 interaction appears governed by two structural features of Fis1: its N-terminal arm and a conserved surface. Alanine scanning mutagenesis of the arm identified both loss-of-function and gain-of-function alleles with mitochondrial morphologies ranging from highly elongated (N6A) to highly fragmented (E7A), demonstrating a profound ability of Fis1 to govern morphology in human cells. An integrated analysis identified a conserved Fis1 residue, Y76, that upon substitution to alanine, but not phenylalanine, also caused highly fragmented mitochondria. The similar phenotypic effects of the E7A and Y76A substitutions, along with NMR data, support that intramolecular interactions occur between the arm and a conserved surface on Fis1 to promote Drp1-mediated fission as in S. cerevisiae. These findings indicate that some aspects of Drp1-mediated fission in humans derive from direct Fis1-Drp1 interactions that are conserved across eukaryotes.

Homomeric interactions of the MPZ Ig domain and their relation to Charcot-Marie-Tooth disease

Mutations in MPZ (myelin protein zero) can cause demyelinating early-onset Charcot-Marie-Tooth type 1B disease or later onset type 2I/J disease characterized by axonal degeneration, reflecting the diverse roles of MPZ in Schwann cells. MPZ holds apposing membranes of the myelin sheath together, with the adhesion role fulfilled by its extracellular immunoglobulin-like domain (IgMPZ), which oligomerizes. Models for how the IgMPZ might form oligomeric assemblies has been extrapolated from a protein crystal structure in which individual rat IgMPZ subunits are packed together under artificial conditions, forming three weak interfaces. One interface organizes the IgMPZ into tetramers, a second 'dimer' interface links tetramers together across the intraperiod line, and a third hydrophobic interface that mediates binding to lipid bilayers or the same hydrophobic surface on another IgMPZ domain. Presently, there are no data confirming whether the proposed IgMPZ interfaces actually mediate oligomerization in solution, whether they are required for the adhesion activity of MPZ, whether they are important for myelination, or whether their loss results in disease. We performed nuclear magnetic resonance spectroscopy and small angle X-ray scattering analysis of wild-type IgMPZ as well as mutant forms with amino acid substitutions designed to interrupt its presumptive oligomerization interfaces. Here, we confirm the interface that mediates IgMPZ tetramerization, but find that dimerization is mediated by a distinct interface that has yet to be identified. We next correlated different types of Charcot-Marie-Tooth disease symptoms to subregions within IgMPZ tetramers. Variants causing axonal late-onset disease (CMT2I/J) map to surface residues of IgMPZ proximal to the transmembrane domain. Variants causing early-onset demyelinating disease (CMT1B) segregate into two groups: one is described by variants that disrupt the stability of the Ig-fold itself and are largely located within the core of the IgMPZ domain; whereas another describes a region on the surface of IgMPZ tetramers, accessible to protein interactions. Computational docking studies predict that this latter disease-relevant subregion may potentially mediate dimerization of IgMPZ tetramers.

In-Silico Selection of Wound-Healing Plant Secondary Molecules and Their Pro-Healing Activities on Experimental Models

This study explored the potential of plant-derived molecules (PDMs) as a medicinal treatment for skin wounds. To assess their healing properties, 34 potential drug molecules (PDMs) and ten therapeutic targets were subjected to molecular docking and dynamics analysis, with allantoin used as a standard compound. Although aristolochic acid had the most potent inhibitory effect, its toxicity made it unsuitable for testing on cells and mice. Therefore, β-caryophyllene (BC) and caryophyllene oxide (BCoxide) were chosen for further testing. The results showed that BC-treated HaCat cells had significantly improved scratch area closure, and both BC and BCoxide treatment produced positive effects such as reduced dermal cellularity and mast cells, decreased levels of inflammation markers IL-6 and TNF-α, and an increase in collagen deposition in mice tissues. However, these treatments did not accelerate wound healing. This study suggests that the PDMs selected based on in-silico results have significant potential for pro-healing abilities. It is essential to conduct further research to confirm our findings.

Potential inhibitors of FemC to combat Staphylococcus aureus: virtual screening, molecular docking, dynamics simulation, and MM-PBSA analysis

FemC is a methicillin resistance factor involved in the alterations of peptidoglycan and glutamine synthesis in Staphylococcus aureus. To identify the potent antibacterial agents, antibacterial molecules were screened against the predicted and validated FemC model. Based on docking scores, presence of essential interactions with active site residues of FemC, pharmacokinetic, and ADMET properties, six candidates were shortlisted and subjected to molecular dynamics to evaluate the stability of FemC-ligand complexes. Further, per residue decomposition analysis and Molecular Mechanics/Poisson-Boltzmann Surface Area (MMPBSA) analysis confirmed that S15, M16, S17, R31, R43, Q47, K48 and R49 of FemC played a vital role in the formation of lower energy stable FemC-inhibitor(s) complexes. Therefore, in the present study, the reported six molecules (Z317461228, Z92241701, Z30923155, Z30202349, Z2609517102 and Z92470167) may pave the path to design the scaffold of novel potent antimicrobials against S. aureus.Communicated by Ramaswamy H. Sarma.

13C and 15N resonance assignments of alpha synuclein fibrils amplified from Lewy Body Dementia tissue

Fibrils of the protein α-synuclein (Asyn) are implicated in the pathogenesis of Parkinson Disease, Lewy Body Dementia, and Multiple System Atrophy. Numerous forms of Asyn fibrils have been studied by solid-state NMR and resonance assignments have been reported. Here, we report a new set of 13C, 15N assignments that are unique to fibrils obtained by amplification from postmortem brain tissue of a patient diagnosed with Lewy Body Dementia.