PDBx/mmCIF Ecosystem: Foundational Semantic Tools for Structural Biology

PDBx/mmCIF, Protein Data Bank Exchange (PDBx) macromolecular Crystallographic Information Framework (mmCIF), has become the data standard for structural biology. With its early roots in the domain of small-molecule crystallography, PDBx/mmCIF provides an extensible data representation that is used for deposition, archiving, remediation, and public dissemination of experimentally determined three-dimensional (3D) structures of biological macromolecules by the Worldwide Protein Data Bank (wwPDB, wwpdb.org). Extensions of PDBx/mmCIF are similarly used for computed structure models by ModelArchive (modelarchive.org), integrative/hybrid structures by PDB-Dev (pdb-dev.wwpdb.org), small angle scattering data by Small Angle Scattering Biological Data Bank SASBDB (sasbdb.org), and for models computed generated with the AlphaFold 2.0 deep learning software suite (alphafold.ebi.ac.uk). Community-driven development of PDBx/mmCIF spans three decades, involving contributions from researchers, software and methods developers in structural sciences, data repository providers, scientific publishers, and professional societies. Having a semantically rich and extensible data framework for representing a wide range of structural biology experimental and computational results, combined with expertly curated 3D biostructure data sets in public repositories, accelerates the pace of scientific discovery. Herein, we describe the architecture of the PDBx/mmCIF data standard, tools used to maintain representations of the data standard, governance, and processes by which data content standards are extended, plus community tools/software libraries available for processing and checking the integrity of PDBx/mmCIF data. Use cases exemplify how the members of the Worldwide Protein Data Bank have used PDBx/mmCIF as the foundation for its pipeline for delivering Findable, Accessible, Interoperable, and Reusable (FAIR) data to many millions of users worldwide.     

Dynamics and Conformational Changes in Human NEIL2 DNA Glycosylase Analyzed by Hydrogen/Deuterium Exchange Mass Spectrometry

Base excision DNA repair (BER) is necessary for removal of damaged nucleobases from the genome and their replacement with normal nucleobases. BER is initiated by DNA glycosylases, the enzymes that cleave the N-glycosidic bonds of damaged deoxynucleotides. Human endonuclease VIII-like protein 2 (hNEIL2), belonging to the helix–two-turn–helix structural superfamily of DNA glycosylases, is an enzyme uniquely specific for oxidized pyrimidines in non-canonical DNA substrates such as bubbles and loops. The structure of hNEIL2 has not been solved; its closest homologs with known structures are NEIL2 from opossum and from giant mimivirus. Here we analyze the conformational dynamics of free hNEIL2 using a combination of hydrogen/deuterium exchange mass spectrometry, homology modeling and molecular dynamics simulations. We show that a prominent feature of vertebrate NEIL2 – a large insert in its N-terminal domain absent from other DNA glycosylases – is unstructured in solution. It was suggested that helix–two-turn–helix DNA glycosylases undergo open–close transition upon DNA binding, with the large movement of their N- and C-terminal domains, but the open conformation has been elusive to capture. Our data point to the open conformation as favorable for free hNEIL2 in solution. Overall, our results are consistent with the view of hNEIL2 as a conformationally flexible protein, which may be due to its participation in the repair of non-canonical DNA structures and/or to the involvement in functional and regulatory protein–protein interactions.     

Methylglyoxal attenuates isoproterenol-induced increase in uncoupling protein 1 expression through activation of JNK signaling pathway in beige adipocytes

Methylglyoxal (MG) is a metabolite derived from glycolysis whose levels in the blood and tissues of patients with diabetes are higher than those of healthy individuals, suggesting that MG is associated with the development of diabetic complications. However, it remains unknown whether high levels of MG are a cause or consequence of diabetes. Here, we show that MG negatively affects the expression of uncoupling protein 1 (UCP1), which is involved in thermogenesis and the regulation of systemic metabolism. Decreased Ucp1 expression is associated with obesity and type 2 diabetes. We found that MG attenuated the increase in Ucp1 expression following treatment with isoproterenol in beige adipocytes. However, MG did not affect protein kinase A signaling, the core coordinator of isoproterenol-induced Ucp1 expression. Instead, MG activated c-Jun N-terminal kinase (JNK) and p38 mitogen-activated protein kinases. We found that JNK inhibition, but not p38, recovered isoproterenol-stimulated Ucp1 expression under MG treatment. Altogether, these results suggest an inhibitory role of MG on the thermogenic function of beige adipocytes through the JNK signaling pathway.     

Genome modeling: From chromatin fibers to genes

The intricacies of the 3D hierarchical organization of the genome have been approached by many creative modeling studies. The specific model/simulation technique combination defines and restricts the system and phenomena that can be investigated. We present the latest modeling developments and studies of the genome, involving models ranging from nucleosome systems and small polynucleosome arrays to chromatin fibers in the kb-range, chromosomes, and whole genomes, while emphasizing gene folding from first principles. Clever combinations allow the exploration of many interesting phenomena involved in gene regulation, such as nucleosome structure and dynamics, nucleosome-nucleosome stacking, polynucleosome array folding, protein regulation of chromatin architecture, mechanisms of gene folding, loop formation, compartmentalization, and structural transitions at the chromosome and genome levels. Gene-level modeling with full details on nucleosome positions, epigenetic factors, and protein binding, in particular, can in principle be scaled up to model chromosomes and cells to study fundamental biological regulation.     

Cryo-EM structures of cardiac thin filaments reveal the 3D architecture of troponin

Troponin is an essential component of striated muscle and it regulates the sliding of actomyosin system in a calcium-dependent manner. Despite its importance, the structure of troponin has been elusive due to its high structural heterogeneity. In this study, we analyzed the 3D structures of murine cardiac thin filaments using a cryo-electron microscope equipped with a Volta phase plate (VPP). Contrast enhancement by a VPP enabled us to reconstruct the entire repeat of the thin filament. We determined the orientation of troponin relative to F-actin and tropomyosin, and characterized the interactions between troponin and tropomyosin. This study provides a structural basis for understanding the molecular mechanism of actomyosin system.     

Structure-based analyses of gut microbiome-related proteins by neural networks and molecular dynamics simulations

An imbalance in the gut microbiome is linked to immune disorders, such as autoimmune, allergic, and chronic inflammatory disorders. Elucidation of disease mechanisms is a matter of urgency. It requires precise elucidation of the structure-based mechanisms of protein interactions involved in disease onset. In addition, an understanding of the protein dynamics is vital because these fluctuations affect the function and interaction of disease-associated proteins. Experimental evaluation of not only protein interactions, functions, and structures but also the dynamics are time-consuming; therefore, computational predictions are necessary to elucidate disease mechanisms. Here, we introduce recent studies on structure-based analyses of proteins using computational approaches, particularly artificial intelligence (AI) and molecular dynamics (MD) simulations.     

Smart de novo Macromolecular Structure Modeling from Cryo-EM Maps

The study of macromolecular structures has expanded our understanding of the amazing cell machinery and such knowledge has changed how the pharmaceutical industry develops new vaccines in recent years. Traditionally, X-ray crystallography has been the main method for structure determination, however, cryogenic electron microscopy (cryo-EM) has increasingly become more popular due to recent advancements in hardware and software. The number of cryo-EM maps deposited in the EMDataResource (formerly EMDatabase) since 2002 has been dramatically increasing and it continues to do so. De novo macromolecular complex modeling is a labor-intensive process, therefore, it is highly desirable to develop software that can automate this process. Here we discuss our automated, data-driven, and artificial intelligence approaches including map processing, feature extraction, modeling building, and target identification. Recently, we have enabled DNA/RNA modeling in our deep learning-based prediction tool, DeepTracer. We have also developed DeepTracer-ID, a tool that can identify proteins solely based on the cryo-EM map. In this paper, we will present our accumulated experiences in developing deep learning-based methods surrounding macromolecule modeling applications.     

SP3-Enabled Rapid and High Coverage Chemoproteomic Identification of Cell-State–Dependent Redox-Sensitive Cysteines

Proteinaceous cysteine residues act as privileged sensors of oxidative stress. As reactive oxygen and nitrogen species have been implicated in numerous pathophysiological processes, deciphering which cysteines are sensitive to oxidative modification and the specific nature of these modifications is essential to understanding protein and cellular function in health and disease. While established mass spectrometry-based proteomic platforms have improved our understanding of the redox proteome, the widespread adoption of these methods is often hindered by complex sample preparation workflows, prohibitive cost of isotopic labeling reagents, and requirements for custom data analysis workflows. Here, we present the SP3-Rox redox proteomics method that combines tailored low cost isotopically labeled capture reagents with SP3 sample cleanup to achieve high throughput and high coverage proteome-wide identification of redox-sensitive cysteines. By implementing a customized workflow in the free FragPipe computational pipeline, we achieve accurate MS1-based quantitation, including for peptides containing multiple cysteine residues. Application of the SP3-Rox method to cellular proteomes identified cysteines sensitive to the oxidative stressor GSNO and cysteine oxidation state changes that occur during T cell activation.     

Deep-lipidotyping by mass spectrometry: recent technical advances and applications

In-depth structural characterization of lipids is an essential component of lipidomics. There has been a rapid expansion of mass spectrometry methods that are capable of resolving lipid isomers at various structural levels over the past decade. These developments finally make deep-lipidotyping possible, which provides new means to study lipid metabolism and discover new lipid biomarkers. In this review, we discuss recent advancements in tandem mass spectrometry (MS/MS) methods for identification of complex lipids beyond the species (known headgroup information) and molecular species (known chain composition) levels. These include identification at the levels of carbon-carbon double bond (C=C) location and sn-position, as well as characterization of acyl chain modifications. We also discuss the integration of isomer-resolving MS/MS methods with different lipid analysis workflows and their applications in lipidomics. The results showcase the distinct capabilities of deep-lipidotyping in untangling the metabolism of individual isomers and sensitive phenotyping by using relative fractional quantitation of the isomers.     

A structural metagenomics pipeline for examining the gut microbiome

Metagenomic sequencing data provide a rich resource from which to expand our understanding of differential protein functions involved in human health. Here, we outline a pipeline that combines microbial whole genome sequencing with protein structure data to yield a structural metagenomics-informed atlas of microbial enzyme families of interest. Visualizing metagenomics data through a structural lens facilitates downstream studies including targeted inhibition and probe-based proteomics to define at the molecular level how different enzyme orthologs impact in vivo function. Application of this pipeline to gut microbial enzymes like glucuronidases, TMA lyases, and bile salt hydrolases is expected to pinpoint their involvement in health and disease and may aid in the development of therapeutics that target specific enzymes within the microbiome.     

Principles of Alternating Access in Multidrug and Toxin Extrusion (MATE) Transporters

The multidrug and toxin extrusion (MATE) transporters catalyze active efflux of a broad range of chemically- and structurally-diverse compounds including antimicrobials and chemotherapeutics, thus contributing to multidrug resistance in pathogenic bacteria and cancers. Multiple methodological approaches have been taken to investigate the structural basis of energy transduction and substrate translocation in MATE transporters. Crystal structures representing members from all three MATE subfamilies have been interpreted within the context of an alternating access mechanism that postulates occupation of distinct structural intermediates in a conformational cycle powered by electrochemical ion gradients. Here we review the structural biology of MATE transporters, integrating the crystallographic models with biophysical and computational studies to define the molecular determinants that shape the transport energy landscape. This holistic analysis highlights both shared and disparate structural and functional features within the MATE family, which underpin an emerging theme of mechanistic diversity within the framework of a conserved structural scaffold.     

A Novel Spectral Annotation Strategy Streamlines Reporting of Mono-ADP-ribosylated Peptides Derived from Mouse Liver and Spleen in Response to IFN-γ

Mass-spectrometry-enabled ADP-ribosylation workflows are developing rapidly, providing researchers a variety of ADP-ribosylome enrichment strategies and mass spectrometric acquisition options. Despite the growth spurt in upstream technologies, systematic ADP-ribosyl (ADPr) peptide mass spectral annotation methods are lacking. HCD-dependent ADP-ribosylome studies are common, but the resulting MS2 spectra are complex, owing to a mixture of b/y-ions and the m/p-ion peaks representing one or more dissociation events of the ADPr moiety (m-ion) and peptide (p-ion). In particular, p-ions that dissociate further into one or more fragment ions can dominate HCD spectra but are not recognized by standard spectral annotation workflows. As a result, annotation strategies that are solely reliant upon the b/y-ions result in lower spectral scores that in turn reduce the number of reportable ADPr peptides. To improve the confidence of spectral assignments, we implemented an ADPr peptide annotation and scoring strategy. All MS2 spectra are scored for the ADPr m-ions, but once spectra are assigned as an ADPr peptide, they are further annotated and scored for the p-ions. We implemented this novel workflow to ADPr peptides enriched from the liver and spleen isolated from mice post 4 h exposure to systemic IFN-γ. HCD collision energy experiments were first performed on the Orbitrap Fusion Lumos and the Q Exactive, with notable ADPr peptide dissociation properties verified with CID (Lumos). The m-ion and p-ion series score distributions revealed that ADPr peptide dissociation properties vary markedly between instruments and within instrument collision energy settings, with consequences on ADPr peptide reporting and amino acid localization. Consequentially, we increased the number of reportable ADPr peptides by 25% (liver) and 17% (spleen) by validation and the inclusion of lower confidence ADPr peptide spectra. This systematic annotation strategy will streamline future reporting of ADPr peptides that have been sequenced using any HCD/CID-based method.     

Principles and Methods in Computational Membrane Protein Design

After decades of progress in computational protein design, the design of proteins folding and functioning in lipid membranes appears today as the next frontier. Some notable successes in the de novo design of simplified model membrane protein systems have helped articulate fundamental principles of protein folding, architecture and interaction in the hydrophobic lipid environment. These principles are reviewed here, together with the computational methods and approaches that were used to identify them. We provide an overview of the methodological innovations in the generation of new protein structures and functions and in the development of membrane-specific energy functions. We highlight the opportunities offered by new machine learning approaches applied to protein design, and by new experimental characterization techniques applied to membrane proteins. Although membrane protein design is in its infancy, it appears more reachable than previously thought.     

First 3D-Structural Data of Full-Length Guanylyl Cyclase 1 in Rod-Outer-Segment Preparations of Bovine Retina by Cross-Linking/Mass Spectrometry

The rod-outer-segment guanylyl cyclase 1 (ROS-GC1) is a key transmembrane protein for retinal phototransduction. Mutations of ROS-GC1 correlate with different retinal diseases that often lead to blindness. No structural data are available for ROS-GC1 so far. We performed a 3D-structural analysis of native ROS-GC1 from bovine retina by cross-linking/mass spectrometry (XL-MS) and computational modeling. Absolute quantification and activity measurements of native ROS-GC1 were performed by MS-based assays directly in bovine retina samples. Our data present the first 3D-structural analysis of active, full-length ROS-GC1 derived from bovine retina. We propose a novel domain organization for the intracellular domain ROS-GC1. Our XL-MS data of native ROS-GC1 from rod-outer-segment preparations of bovine retina agree with a dimeric architecture. Our integrated approach can serve as a blueprint for conducting 3D-structural studies of membrane proteins in their native environment.     

Probing the sORF-Encoded Peptides of Deinococcus radiodurans in Response to Extreme Stress

Organisms have developed different mechanisms to respond to stresses. However, the roles of small ORF–encoded peptides (SEPs) in these regulatory systems remain elusive, which is partially because of the lack of comprehensive knowledge regarding these biomolecules. We chose the extremophile Deinococcus radiodurans R1 as a model species and conducted large-scale profiling of the SEPs related to the stress response. The integrated workflow consisting of multiple omics approaches for SEP identification was streamlined, and an SEPome of D. radiodurans containing 109 novel and high-confidence SEPs was drafted. Forty-four percent of these SEPs were predicted to function as antimicrobial peptides. Quantitative peptidomics analysis indicated that the expression of SEP068184 was upregulated upon oxidative treatment and gamma irradiation of the bacteria. SEP068184 was conserved in Deinococcus and exhibited negative regulation of oxidative stress resistance in a comparative phenotypic assay of its mutants. Further quantitative and interactive proteomics analyses suggested that SEP068184 might function through metabolic pathways and interact with cytoplasmic proteins. Collectively, our findings demonstrate that SEPs are involved in the regulation of oxidative resistance, and the SEPome dataset provides a rich resource for research on the molecular mechanisms of the response to extreme stress in organisms.     

Comparative Analysis of T-Cell Spatial Proteomics and the Influence of HIV Expression

As systems biology approaches to virology have become more tractable, highly studied viruses such as HIV can now be analyzed in new unbiased ways, including spatial proteomics. We employed here a differential centrifugation protocol to fractionate Jurkat T cells for proteomic analysis by mass spectrometry; these cells contain inducible HIV-1 genomes, enabling us to look for changes in the spatial proteome induced by viral gene expression. Using these proteomics data, we evaluated the merits of several reported machine learning pipelines for classification of the spatial proteome and identification of protein translocations. From these analyses, we found that classifier performance in this system was organelle dependent, with Bayesian t-augmented Gaussian mixture modeling outperforming support vector machine learning for mitochondrial and endoplasmic reticulum proteins but underperforming on cytosolic, nuclear, and plasma membrane proteins by QSep analysis. We also observed a generally higher performance for protein translocation identification using a Bayesian model, Bayesian analysis of differential localization experiments, on row-normalized data. Comparative Bayesian analysis of differential localization experiment analysis of cells induced to express the WT viral genome versus cells induced to express a genome unable to express the accessory protein Nef identified known Nef-dependent interactors such as T-cell receptor signaling components and coatomer complex. Finally, we found that support vector machine classification showed higher consistency and was less sensitive to HIV-dependent noise. These findings illustrate important considerations for studies of the spatial proteome following viral infection or viral gene expression and provide a reference for future studies of HIV-gene-dropout viruses.     

Electron microscopy of cellular ultrastructure in three dimensions

Electron microscopy in three dimensions (3D) of cells and tissues can be essential for understanding the ultrastructural aspects of biological processes. The quest for 3D information reveals challenges at many stages of the workflow, from sample preparation, to imaging, data analysis and segmentation. Here, we outline several available methods, including volume SEM imaging, cryo-TEM and cryo-STEM tomography, each one occupying a different domain in the basic tradeoff between field-of-view and resolution. We discuss the considerations for choosing a suitable method depending on research needs and highlight recent developments that are essential for making 3D volume imaging of cells and tissues a standard tool for cellular and structural biologists.     

The ugly,bad, and good stories of large-scale biomolecular simulations

Molecular modeling of large biomolecular assemblies exemplifies a disruptive area holding both promises and contentions. Propelled by peta and exascale computing, several simulation methodologies have now matured into user-friendly tools that are successfully employed for modeling viruses, membranous nano-constructs, and key pieces of the genetic machinery. We present three unifying biophysical themes that emanate from some of the most recent multi-million atom simulation endeavors. Despite connecting molecular changes with phenotypic outcomes, the quality measures of these simulations remain questionable. We discuss the existing and upcoming strategies for constructing representative ensembles of large systems, how new computing technologies will boost this area, and make a point that integrative modeling guided by experimental data is the future of biomolecular computations.