Integrative structural studies of the SARS-CoV-2 spike protein during the fusion process (2022).

SARS-CoV-2 is the virus responsible for the COVID-19 pandemic and catastrophic, worldwide health and economic impacts. The spike protein on the viral surface is responsible for viral entry into the host cell. The binding of spike protein to the host cell receptor ACE2 is the first step leading to fusion of the host and viral membranes. Despite the vast amount of structure data that has been generated for the spike protein of SARS-CoV-2, many of the detailed structures of the spike protein in different stages of the fusion pathway are unknown, leaving a wealth of potential drug-target space unexplored. The atomic-scale structure of the complete S2 segment, as well as the complete fusion intermediate are also unknown and represent major gaps in our knowledge of the infectious pathway of SAR-CoV-2. The conformational changes of the spike protein during this process are similarly not well understood. Here we present structures of the spike protein at different stages of the fusion process. With the transitions being a necessary step before the receptor binding, we propose sites along the transition pathways as potential targets for drug development.

Modeling the Influenza A NP-vRNA-Polymerase Complex in Atomic Detail.

Seasonal flu is an acute respiratory disease that exacts a massive toll on human populations, healthcare systems and economies. The disease is caused by an enveloped Influenza virus containing eight ribonucleoprotein (RNP) complexes. Each RNP incorporates multiple copies of nucleoprotein (NP), a fragment of the viral genome (vRNA), and a viral RNA-dependent RNA polymerase (POL), and is responsible for packaging the viral genome and performing critical functions including replication and transcription. A complete model of an Influenza RNP in atomic detail can elucidate the structural basis for viral genome functions, and identify potential targets for viral therapeutics. In this work we construct a model of a complete Influenza A RNP complex in atomic detail using multiple sources of structural and sequence information and a series of homology-modeling techniques, including a motif-matching fragment assembly method. Our final model provides a rationale for experimentally-observed changes to viral polymerase activity in numerous mutational assays. Further, our model reveals specific interactions between the three primary structural components of the RNP, including potential targets for blocking POL-binding to the NP-vRNA complex. The methods developed in this work open the possibility of elucidating other functionally-relevant atomic-scale interactions in additional RNP structures and other biomolecular complexes.

Scaling molecular dynamics beyond 100,000 processor cores for large-scale biophysical simulations.

The growing interest in the complexity of biological interactions is continuously driving the need to increase system size in biophysical simulations, requiring not only powerful and advanced hardware but adaptable software that can accommodate a large number of atoms interacting through complex forcefields. To address this, we developed and implemented strategies in the GENESIS molecular dynamics package designed for large numbers of processors. Long-range electrostatic interactions were parallelized by minimizing the number of processes involved in communication. A novel algorithm was implemented for nonbonded interactions to increase single instruction multiple data (SIMD) performance, reducing memory usage for ultra large systems. Memory usage for neighbor searches in real-space nonbonded interactions was reduced by approximately 80%, leading to significant speedup. Using experimental data describing physical 3D chromatin interactions, we constructed the first atomistic model of an entire gene locus (GATA4). Taken together, these developments enabled the first billion-atom simulation of an intact biomolecular complex, achieving scaling to 65,000 processes (130,000 processor cores) with 1 ns/day performance. Published 2019. This article is a U.S. Government work and is in the public domain in the USA.

Histidine-Dependent Protein Methylation Is Required for Compartmentalization of CTP Synthase.

CTP synthase (CTPS) forms compartmentalized filaments in response to substrate availability and environmental nutrient status. However, the physiological role of filaments and mechanisms for filament assembly are not well understood. Here, we provide evidence that CTPS forms filaments in response to histidine influx during glutamine starvation. Tetramer conformation-based filament formation restricts CTPS enzymatic activity during nutrient deprivation. CTPS protein levels remain stable in the presence of histidine during nutrient deprivation, followed by rapid cell growth after stress relief. We demonstrate that filament formation is controlled by methylation and that histidine promotes re-methylation of homocysteine by donating one-carbon intermediates to the cytosolic folate cycle. Furthermore, we find that starvation stress and glutamine deficiency activate the GCN2/ATF4/MTHFD2 axis, which coordinates CTPS filament formation. CTPS filament formation induced by histidine-mediated methylation may be a strategy used by cancer cells to maintain homeostasis and ensure a growth advantage in adverse environments.

Structural Transition and Antibody Binding of EBOV GP and ZIKV E Proteins from Pre-Fusion to Fusion-Initiation State.

Membrane fusion proteins are responsible for viral entry into host cells—a crucial first step in viral infection. These proteins undergo large conformational changes from pre-fusion to fusion-initiation structures, and, despite differences in viral genomes and disease etiology, many fusion proteins are arranged as trimers. Structural information for both pre-fusion and fusion-initiation states is critical for understanding virus neutralization by the host immune system. In the case of Ebola virus glycoprotein (EBOV GP) and Zika virus envelope protein (ZIKV E), pre-fusion state structures have been identified experimentally, but only partial structures of fusion-initiation states have been described. While the fusion-initiation structure is in an energetically unfavorable state that is difficult to solve experimentally, the existing structural information combined with computational approaches enabled the modeling of fusion-initiation state structures of both proteins. These structural models provide an improved understanding of four different neutralizing antibodies in the prevention of viral host entry.

Identification and Characterization of a Class of MALAT1-like Genomic Loci.

The MALAT1 (Metastasis-Associated Lung Adenocarcinoma Transcript 1) gene encodes a noncoding RNA that is processed into a long nuclear retained transcript (MALAT1) and a small cytoplasmic tRNA-like transcript (mascRNA). Using an RNA sequence- and structure-based covariance model, we identified more than 130 genomic loci in vertebrate genomes containing the MALAT1 3' end triple-helix structure and its immediate downstream tRNA-like structure, including 44 in the green lizard Anolis carolinensis. Structural and computational analyses revealed a co-occurrence of components of the 3' end module. MALAT1-like genes in Anolis carolinensis are highly expressed in adult testis, thus we named them testis-abundant long noncoding RNAs (tancRNAs). MALAT1-like loci also produce multiple small RNA species, including PIWI-interacting RNAs (piRNAs), from the antisense strand. The 3' ends of tancRNAs serve as potential targets for the PIWI-piRNA complex. Thus, we have identified an evolutionarily conserved class of long noncoding RNAs (lncRNAs) with similar structural constraints, post-transcriptional processing, and subcellular localization and a distinct function in spermatocytes.

Functional and Structural Analysis of a Highly-Expressed Yersinia pestis Small RNA following Infection of Cultured Macrophages.

Non-coding small RNAs (sRNAs) are found in practically all bacterial genomes and play important roles in regulating gene expression to impact bacterial metabolism, growth, and virulence. We performed transcriptomics analysis to identify sRNAs that are differentially expressed in Yersinia pestis that invaded the human macrophage cell line THP-1, compared to pathogens that remained extracellular in the presence of host. Using ultra high-throughput sequencing, we identified 37 novel and 143 previously known sRNAs in Y. pestis. In particular, the sRNA Ysr170 was highly expressed in intracellular Yersinia and exhibited a log2 fold change ~3.6 higher levels compared to extracellular bacteria. We found that knock-down of Ysr170 expression attenuated infection efficiency in cell culture and growth rate in response to different stressors. In addition, we applied selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) analysis to determine the secondary structure of Ysr170 and observed structural changes resulting from interactions with the aminoglycoside antibiotic gentamycin and the RNA chaperone Hfq. Interestingly, gentamicin stabilized helix 4 of Ysr170, which structurally resembles the native gentamicin 16S ribosomal binding site. Finally, we modeled the tertiary structure of Ysr170 binding to gentamycin using RNA motif modeling. Integration of these experimental and structural methods can provide further insight into the design of small molecules that can inhibit function of sRNAs required for pathogen virulence.

A Highly Conserved Yet Flexible Linker Is Part of a Polymorphic Protein-Binding Domain in Myosin-Binding Protein C.

The nuclear magnetic resonance (NMR) structure of the tri-helix bundle (THB) of the m-domain plus C2 (ΔmC2) of myosin-binding protein C (MyBP-C) has revealed a highly flexible seven-residue linker between the structured THB and C2. Bioinformatics shows significant patterns of conservation across the THB-linker sequence, with the linker containing a strictly conserved serine in all MyBP-C isoforms. Clinically linked mutations further support the functional significance of the THB-linker region. NMR, small-angle X-ray scattering, and binding studies show the THB-linker plus the first ten residues of C2 undergo dramatic changes when ΔmC2 binds Ca2+-calmodulin, with the linker and C2 N-terminal residues contributing significantly to the affinity. Modeling of all available experimental data indicates that the THB tertiary structure must be disrupted to form the complex. These results are discussed in the context of the THB-linker and the N-terminal residues of C2 forming a polymorphic binding domain that could accommodate multiple binding partners in the dynamic sarcomere.

Regulation of the mammalian elongation cycle by subunit rolling: a eukaryotic-specific ribosome rearrangement.

The extent to which bacterial ribosomes and the significantly larger eukaryotic ribosomes share the same mechanisms of ribosomal elongation is unknown. Here, we present subnanometer resolution cryoelectron microscopy maps of the mammalian 80S ribosome in the posttranslocational state and in complex with the eukaryotic eEF1A⋅Val-tRNA⋅GMPPNP ternary complex, revealing significant differences in the elongation mechanism between bacteria and mammals. Surprisingly, and in contrast to bacterial ribosomes, a rotation of the small subunit around its long axis and orthogonal to the well-known intersubunit rotation distinguishes the posttranslocational state from the classical pretranslocational state ribosome. We term this motion "subunit rolling." Correspondingly, a mammalian decoding complex visualized in substates before and after codon recognition reveals structural distinctions from the bacterial system. These findings suggest how codon recognition leads to GTPase activation in the mammalian system and demonstrate that in mammalia subunit rolling occurs during tRNA selection.

Cryo-EM structure of the small subunit of the mammalian mitochondrial ribosome.

The mammalian mitochondrial ribosomes (mitoribosomes) are responsible for synthesizing 13 membrane proteins that form essential components of the complexes involved in oxidative phosphorylation or ATP generation for the eukaryotic cell. The mammalian 55S mitoribosome contains significantly smaller rRNAs and a large mass of mitochondrial ribosomal proteins (MRPs), including large mito-specific amino acid extensions and insertions in MRPs that are homologous to bacterial ribosomal proteins and an additional 35 mito-specific MRPs. Here we present the cryo-EM structure analysis of the small (28S) subunit (SSU) of the 55S mitoribosome. We find that the mito-specific extensions in homologous MRPs generally are involved in inter-MRP contacts and in contacts with mito-specific MRPs, suggesting a stepwise evolution of the current architecture of the mitoribosome. Although most of the mito-specific MRPs and extensions of homologous MRPs are situated on the peripheral regions, they also contribute significantly to the formation of linings of the mRNA and tRNA paths, suggesting a tailor-made structural organization of the mito-SSU for the recruitment of mito-specific mRNAs, most of which do not possess a 5' leader sequence. In addition, docking of previously published coordinates of the large (39S) subunit (LSU) into the cryo-EM map of the 55S mitoribosome reveals that mito-specific MRPs of both the SSU and LSU are involved directly in the formation of six of the 15 intersubunit bridges.

Optimal aggregation of FcεRI with a structurally defined trivalent ligand overrides negative regulation driven by phosphatases.

To investigate why responses of mast cells to antigen-induced IgE receptor (FcεRI) aggregation depend nonlinearly on antigen dose, we characterized a new artificial ligand, DF3, through complementary modeling and experimentation. This ligand is a stable trimer of peptides derived from bacteriophage T4 fibritin, each conjugated to a hapten (DNP). We found low and high doses of DF3 at which degranulation of mast cells sensitized with DNP-specific IgE is minimal, but ligand-induced receptor aggregation is comparable to aggregation at an intermediate dose, optimal for degranulation. This finding makes DF3 an ideal reagent for studying the balance of negative and positive signaling in the FcεRI pathway. We find that the lipid phosphatase SHIP and the protein tyrosine phosphatase SHP-1 negatively regulate mast cell degranulation over all doses considered. In contrast, SHP-2 promotes degranulation. With high DF3 doses, relatively rapid recruitment of SHIP to the plasma membrane may explain the reduced degranulation response. Our results demonstrate that optimal secretory responses of mast cells depend on the formation of receptor aggregates that promote sufficient positive signaling by Syk to override phosphatase-mediated negative regulatory signals.

Fitness costs of rifampicin resistance in Mycobacterium tuberculosis are amplified under conditions of nutrient starvation and compensated by mutation in the ß' subunit of RNA polymerase.

Rifampicin resistance, a defining attribute of multidrug-resistant tuberculosis, is conferred by mutations in the ß subunit of RNA polymerase. Sequencing of rifampicin-resistant (RIF-R) clinical isolates of Mycobacterium tuberculosis revealed, in addition to RIF-R mutations, enrichment of potential compensatory mutations around the double-psi ß-barrel domain of the ß' subunit comprising the catalytic site and the exit tunnel for newly synthesized RNA. Sequential introduction of the resistance allele followed by the compensatory allele in isogenic Mycobacterium smegmatis showed that these mutations respectively caused and compensated a starvation enhanced growth defect by altering RNA polymerase activity. While specific combinations of resistance and compensatory alleles converged in divergent lineages, other combinations recurred among related isolates suggesting transmission of compensated RIF-R strains. These findings suggest nutrient poor growth conditions impose larger selective pressure on RIF-R organisms that results in the selection of compensatory mutations in a domain involved in catalysis and starvation control of RNA polymerase transcription.

Rule-based modeling: a computational approach for studying biomolecular site dynamics in cell signaling systems.

Rule-based modeling was developed to address the limitations of traditional approaches for modeling chemical kinetics in cell signaling systems. These systems consist of multiple interacting biomolecules (e.g., proteins), which themselves consist of multiple parts (e.g., domains, linear motifs, and sites of phosphorylation). Consequently, biomolecules that mediate information processing generally have the potential to interact in multiple ways, with the number of possible complexes and posttranslational modification states tending to grow exponentially with the number of binary interactions considered. As a result, only large reaction networks capture all possible consequences of the molecular interactions that occur in a cell signaling system, which is problematic because traditional modeling approaches for chemical kinetics (e.g., ordinary differential equations) require explicit network specification. This problem is circumvented through representation of interactions in terms of local rules. With this approach, network specification is implicit and model specification is concise. Concise representation results in a coarse graining of chemical kinetics, which is introduced because all reactions implied by a rule inherit the rate law associated with that rule. Coarse graining can be appropriate if interactions are modular, and the coarseness of a model can be adjusted as needed. Rules can be specified using specialized model-specification languages, and recently developed tools designed for specification of rule-based models allow one to leverage powerful software engineering capabilities. A rule-based model comprises a set of rules, which can be processed by general-purpose simulation and analysis tools to achieve different objectives (e.g., to perform either a deterministic or stochastic simulation).

Rise of the RNA machines: exploring the structure of long non-coding RNAs.

Novel, profound and unexpected roles of long non-coding RNAs (lncRNAs) are emerging in critical aspects of gene regulation. Thousands of lncRNAs have been recently discovered in a wide range of mammalian systems, related to development, epigenetics, cancer, brain function and hereditary disease. The structural biology of these lncRNAs presents a brave new RNA world, which may contain a diverse zoo of new architectures and mechanisms. While structural studies of lncRNAs are in their infancy, we describe existing structural data for lncRNAs, as well as crystallographic studies of other RNA machines and their implications for lncRNAs. We also discuss the importance of dynamics in RNA machine mechanism. Determining commonalities between lncRNA systems will help elucidate the evolution and mechanistic role of lncRNAs in disease, creating a structural framework necessary to pursue lncRNA-based therapeutics.

A structural model of the E. coli PhoB dimer in the transcription initiation complex.

BACKGROUND: There exist > 78,000 proteins and/or nucleic acids structures that were determined experimentally. Only a small portion of these structures corresponds to those of protein complexes. While homology modeling is able to exploit knowledge-based potentials of side-chain rotomers and backbone motifs to infer structures for new proteins, no such general method exists to extend our understanding of protein interaction motifs to novel protein complexes. RESULTS: We use a Motif Binding Geometries (MBG) approach, to infer the structure of a protein complex from the database of complexes of homologous proteins taken from other contexts (such as the helix-turn-helix motif binding double stranded DNA), and demonstrate its utility on one of the more important regulatory complexes in biology, that of the RNA polymerase initiating transcription under conditions of phosphate starvation. The modeled PhoB/RNAP/σ-factor/DNA complex is stereo-chemically reasonable, has sufficient interfacial Solvent Excluded Surface Areas (SESAs) to provide adequate binding strength, is physically meaningful for transcription regulation, and is consistent with a variety of known experimental constraints. CONCLUSIONS: Based on a straightforward and easy to comprehend concept, "proteins and protein domains that fold similarly could interact similarly", a structural model of the PhoB dimer in the transcription initiation complex has been developed. This approach could be extended to enable structural modeling and prediction of other bio-molecular complexes. Just as models of individual proteins provide insight into molecular recognition, catalytic mechanism, and substrate specificity, models of protein complexes will provide understanding into the combinatorial rules of cellular regulation and signaling.

Global structure of a three-way junction in a phi29 packaging RNA dimer determined using site-directed spin labeling.

The condensation of bacteriophage phi29 genomic DNA into its preformed procapsid requires the DNA packaging motor, which is the strongest known biological motor. The packaging motor is an intricate ring-shaped protein/RNA complex, and its function requires an RNA component called packaging RNA (pRNA). Current structural information on pRNA is limited, which hinders studies of motor function. Here, we used site-directed spin labeling to map the conformation of a pRNA three-way junction that bridges binding sites for the motor ATPase and the procapsid. The studies were carried out on a pRNA dimer, which is the simplest ring-shaped pRNA complex and serves as a functional intermediate during motor assembly. Using a nucleotide-independent labeling scheme, stable nitroxide radicals were attached to eight specific pRNA sites without perturbing RNA folding and dimer formation, and a total of 17 internitroxide distances spanning the three-way junction were measured using Double Electron-Electron Resonance spectroscopy. The measured distances, together with steric chemical constraints, were used to select 3662 viable three-way junction models from a pool of 65 billion. The results reveal a similar conformation among the viable models, with two of the helices (H(T) and H(L)) adopting an acute bend. This is in contrast to a recently reported pRNA tetramer crystal structure, in which H(T) and H(L) stack onto each other linearly. The studies establish a new method for mapping global structures of complex RNA molecules, and provide information on pRNA conformation that aids investigations of phi29 packaging motor and developments of pRNA-based nanomedicine and nanomaterial.

The Brucella TIR-like protein TcpB interacts with the death domain of MyD88.

The pathogen Brucella melitensis secretes a Toll/interleukin-1 receptor (TIR) domain containing protein that abrogates host innate immune responses. In this study, we have characterized the biochemical interactions of Brucella TIR-like protein TcpB with host innate immune adaptor proteins. Using protein-fragment complementation assays based on Gaussia luciferase and green fluorescent protein, we find that TcpB interacts directly with MyD88 and that this interaction is significantly stronger than the interaction of TcpB with TIRAP, the only other adaptor protein that detectably interacts with TcpB. Surprisingly, the TcpB-MyD88 interaction depends on the death domain (DD) of MyD88, and TcpB does not interact with the isolated TIR domain of MyD88. TcpB disrupts MyD88(DD)-MyD88(DD), MyD88(DD)-MyD88(TIR) and MyD88(DD)-MyD88 interactions but not MyD88-MyD88 or MyD88(TIR)-MyD88(TIR) interactions. Structural models consistent with these results suggest how TcpB might inhibit TLR signaling by targeting MyD88 via a DD-TIR domain interface.

A fast dynamic mode of the EF-G-bound ribosome.

A key intermediate in translocation is an 'unlocked state' of the pre-translocation ribosome in which the P-site tRNA adopts the P/E hybrid state, the L1 stalk domain closes and ribosomal subunits adopt a ratcheted configuration. Here, through two- and three-colour smFRET imaging from multiple structural perspectives, EF-G is shown to accelerate structural and kinetic pathways in the ribosome, leading to this transition. The EF-G-bound ribosome remains highly dynamic in nature, wherein, the unlocked state is transiently and reversibly formed. The P/E hybrid state is energetically favoured, but exchange with the classical P/P configuration persists; the L1 stalk adopts a fast dynamic mode characterized by rapid cycles of closure and opening. These data support a model in which P/E hybrid state formation, L1 stalk closure and subunit ratcheting are loosely coupled, independent processes that must converge to achieve the unlocked state. The highly dynamic nature of these motions, and their sensitivity to conformational and compositional changes in the ribosome, suggests that regulating the formation of this intermediate may present an effective avenue for translational control.

Spontaneous formation of the unlocked state of the ribosome is a multistep process.

The mechanism of substrate translocation through the ribosome is central to the rapid and faithful translation of mRNA into proteins. The rate-limiting step in translocation is an unlocking process that includes the formation of an "unlocked" intermediate state, which requires the convergence of large-scale conformational events within the ribosome including tRNA hybrid states formation, closure of the ribosomal L1 stalk domain, and subunit ratcheting. Here, by imaging of the pretranslocation ribosome complex from multiple structural perspectives using two- and three-color single-molecule fluorescence resonance energy transfer, we observe that tRNA hybrid states formation and L1 stalk closure, events central to the unlocking mechanism, are not tightly coupled. These findings reveal that the unlocked state is achieved through a stochastic-multistep process, where the extent of conformational coupling depends on the nature of tRNA substrates. These data suggest that cellular mechanisms affecting the coupling of conformational processes on the ribosome may regulate the process of translation elongation.

Mismatched dNTP incorporation by DNA polymerase beta does not proceed via globally different conformational pathways.

Understanding how DNA polymerases control fidelity requires elucidation of the mechanisms of matched and mismatched dNTP incorporations. Little is known about the latter because mismatched complexes do not crystallize readily. In this report, we employed small-angle X-ray scattering (SAXS) and structural modeling to probe the conformations of different intermediate states of mammalian DNA polymerase beta (Pol beta) in its wild-type and an error-prone variant, I260Q. Our structural results indicate that the mismatched ternary complex lies in-between the open and the closed forms, but more closely resembles the open form for WT and the closed form for I260Q. On the basis of molecular modeling, this over-stabilization of mismatched ternary complex of I260Q is likely caused by formation of a hydrogen bonding network between the side chains of Gln(260), Tyr(296), Glu(295) and Arg(258), freeing up Asp(192) to coordinate MgdNTP. These results argue against recent reports suggesting that mismatched dNTP incorporations follow a conformational path distinctly different from that of matched dNTP incorporation, or that its conformational closing is a major contributor to fidelity.