Advances, interactions, and future developments in the CNS, Phenix, and Rosetta structural biology software systems

Structural characterization of macro domain-containing Thoeris antiphage defense systems

Thoeris defense systems protect bacteria from infection by phages via abortive infection. In these systems, ThsB proteins serve as sensors of infection and generate signaling nucleotides that activate ThsA effectors. Silent information regulator and SMF/DprA-LOG (SIR2-SLOG) containing ThsA effectors are activated by cyclic ADP-ribose (ADPR) isomers 2'cADPR and 3'cADPR, triggering abortive infection via nicotinamide adenine dinucleotide (NAD+) depletion. Here, we characterize Thoeris systems with transmembrane and macro domain (TM-macro)-containing ThsA effectors. We demonstrate that ThsA macro domains bind ADPR and imidazole adenine dinucleotide (IAD), but not 2'cADPR or 3'cADPR. Combining crystallography, in silico predictions, and site-directed mutagenesis, we show that ThsA macro domains form nucleotide-induced higher-order oligomers, enabling TM domain clustering. We demonstrate that ThsB can produce both ADPR and IAD, and we identify a ThsA TM-macro-specific ThsB subfamily with an active site resembling deoxy-nucleotide and deoxy-nucleoside processing enzymes. Collectively, our study demonstrates that Thoeris systems with SIR2-SLOG and TM-macro ThsA effectors trigger abortive infection via distinct mechanisms.

Farnesyltransferase inhibitor lonafarnib suppresses respiratory syncytial virus infection by blocking conformational change of fusion glycoprotein

Respiratory syncytial virus (RSV) is the major cause of bronchiolitis and pneumonia in young children and the elderly. There are currently no approved RSV-specific therapeutic small molecules available. Using high-throughput antiviral screening, we identified an oral drug, the prenylation inhibitor lonafarnib, which showed potent inhibition of the RSV fusion process. Lonafarnib exhibited antiviral activity against both the RSV A and B genotypes and showed low cytotoxicity in HEp-2 and human primary bronchial epithelial cells (HBEC). Time-of-addition and pseudovirus assays demonstrated that lonafarnib inhibits RSV entry, but has farnesyltransferase-independent antiviral efficacy. Cryo-electron microscopy revealed that lonafarnib binds to a triple-symmetric pocket within the central cavity of the RSV F metastable pre-fusion conformation. Mutants at the RSV F sites interacting with lonafarnib showed resistance to lonafarnib but remained fully sensitive to the neutralizing monoclonal antibody palivizumab. Furthermore, lonafarnib dose-dependently reduced the replication of RSV in BALB/c mice. Collectively, lonafarnib could be a potential fusion inhibitor for RSV infection.

An adeno-associated virus variant enabling efficient ocular-directed gene delivery across species

Recombinant adeno-associated viruses (rAAVs) have emerged as promising gene therapy vectors due to their proven efficacy and safety in clinical applications. In non-human primates (NHPs), rAAVs are administered via suprachoroidal injection at a higher dose. However, high doses of rAAVs tend to increase additional safety risks. Here, we present a novel AAV capsid (AAVv128), which exhibits significantly enhanced transduction efficiency for photoreceptors and retinal pigment epithelial (RPE) cells, along with a broader distribution across the layers of retinal tissues in different animal models (mice, rabbits, and NHPs) following intraocular injection. Notably, the suprachoroidal delivery of AAVv128-anti-VEGF vector completely suppresses the Grade IV lesions in a laser-induced choroidal neovascularization (CNV) NHP model for neovascular age-related macular degeneration (nAMD). Furthermore, cryo-EM analysis at 2.1 Å resolution reveals that the critical residues of AAVv128 exhibit a more robust advantage in AAV binding, the nuclear uptake and endosome escaping. Collectively, our findings highlight the potential of AAVv128 as a next generation ocular gene therapy vector, particularly using the suprachoroidal delivery route.

Error estimates in atom coordinates and B factors in macromolecular crystallography

The overall diffraction precision index (DPI) of a biological macromolecule crystal structure was first described by Cruickshank in 1999. This topical review proceeds from this point and describes the subsequent elaboration of the index to individual atom coordinates. Additional developments were introduced by the availability of a webserver, which provides a transformed PDB entry with individual atom coordinate errors derived from applying the DPI method using the parameters provided by the authors and then subsequently added to the PDB file. This webserver has been extensively used and harnessed in describing non-covalent distance error estimates as well as assessing the significance, or otherwise, of atom movements in a variety of studies. The standard uncertainties on a biological macromolecule's atomic displacement parameters (the 'B factors') has been an entirely different challenge but is obviously important since the crystallographic community has developed the habit of quoting B factors to a false precision in papers. This can convey a false certainty in the dynamics of a structure. A method involving parallelisation of workflows for diffraction image data processing does however offer estimates of the precision of B factors.

Structure of the human ATP synthase

Biological energy currency ATP is produced by F₁F_o-ATP synthase. However, the molecular mechanism for human ATP synthase action remains unknown. Here, we present snapshot images for three main rotational states and one substate of human ATP synthase using cryoelectron microscopy. These structures reveal that the release of ADP occurs when the β subunit of F₁F_o-ATP synthase is in the open conformation, showing how ADP binding is coordinated during synthesis. The accommodation of the symmetry mismatch between F₁ and F_o motors is resolved by the torsional flexing of the entire complex, especially the γ subunit, and the rotational substep of the c subunit. Water molecules are identified in the inlet and outlet half-channels, suggesting that the proton transfer in these two half-channels proceed via a Grotthus mechanism. Clinically relevant mutations are mapped to the structure, showing that they are mainly located at the subunit-subunit interfaces, thus causing instability of the complex.

The dynamic nature of netrin-1 and the structural basis for glycosaminoglycan fragment-induced filament formation

Netrin-1 is a bifunctional chemotropic guidance cue that plays key roles in diverse cellular processes including axon pathfinding, cell migration, adhesion, differentiation, and survival. Here, we present a molecular understanding of netrin-1 mediated interactions with glycosaminoglycan chains of diverse heparan sulfate proteoglycans (HSPGs) and short heparin oligosaccharides. Whereas interactions with HSPGs act as platform to co-localise netrin-1 close to the cell surface, heparin oligosaccharides have a significant impact on the highly dynamic behaviour of netrin-1. Remarkably, the monomer-dimer equilibrium of netrin-1 in solution is abolished in the presence of heparin oligosaccharides and replaced with highly hierarchical and distinct super assemblies leading to unique, yet unknown netrin-1 filament formation. In our integrated approach we provide a molecular mechanism for the filament assembly which opens fresh paths towards a molecular understanding of netrin-1 functions.

Domain-Based Protein Docking with Extremely Large Conformational Changes

Proteins are key components in many processes in living cells, and physical interactions with other proteins and nucleic acids often form key parts of their functions. In many cases, large flexibility of proteins as they interact is key to their function. To understand the mechanisms of these processes, it is necessary to consider the 3D structures of such protein complexes. When such structures are not yet experimentally determined, protein docking has long been present to computationally generate useful structure models. However, protein docking has long had the limitation that the consideration of flexibility is usually limited to very small movements or very small structures. Methods have been developed which handle minor flexibility via normal mode or other structure sampling, but new methods are required to model ordered proteins which undergo large-scale conformational changes to elucidate their function at the molecular level. Here, we present Flex-LZerD, a framework for docking such complexes. Via partial assembly multidomain docking and an iterative normal mode analysis admitting curvilinear motions, we demonstrate the ability to model the assembly of a variety of protein-protein and protein-nucleic acid complexes.

Cyclic ADP ribose isomers: Production, chemical structures, and immune signaling

Cyclic adenosine diphosphate (ADP)-ribose (cADPR) isomers are signaling molecules produced by bacterial and plant Toll/interleukin-1 receptor (TIR) domains via nicotinamide adenine dinucleotide (oxidized form) (NAD+) hydrolysis. We show that v-cADPR (2'cADPR) and v2-cADPR (3'cADPR) isomers are cyclized by O-glycosidic bond formation between the ribose moieties in ADPR. Structures of 2'cADPR-producing TIR domains reveal conformational changes that lead to an active assembly that resembles those of Toll-like receptor adaptor TIR domains. Mutagenesis reveals a conserved tryptophan that is essential for cyclization. We show that 3'cADPR is an activator of ThsA effector proteins from the bacterial antiphage defense system termed Thoeris and a suppressor of plant immunity when produced by the effector HopAM1. Collectively, our results reveal the molecular basis of cADPR isomer production and establish 3'cADPR in bacteria as an antiviral and plant immunity-suppressing signaling molecule.

Complementarity of neutron, XFEL and synchrotron crystallography for defining the structures of metalloenzymes at room temperature

Room-temperature macromolecular crystallography allows protein structures to be determined under close-to-physiological conditions, permits dynamic freedom in protein motions and enables time-resolved studies. In the case of metalloenzymes that are highly sensitive to radiation damage, such room-temperature experiments can present challenges, including increased rates of X-ray reduction of metal centres and site-specific radiation-damage artefacts, as well as in devising appropriate sample-delivery and data-collection methods. It can also be problematic to compare structures measured using different crystal sizes and light sources. In this study, structures of a multifunctional globin, dehaloperoxidase B (DHP-B), obtained using several methods of room-temperature crystallographic structure determination are described and compared. Here, data were measured from large single crystals and multiple microcrystals using neutrons, X-ray free-electron laser pulses, monochromatic synchrotron radiation and polychromatic (Laue) radiation light sources. These approaches span a range of 18 orders of magnitude in measurement time per diffraction pattern and four orders of magnitude in crystal volume. The first room-temperature neutron structures of DHP-B are also presented, allowing the explicit identification of the hydrogen positions. The neutron data proved to be complementary to the serial femtosecond crystallography data, with both methods providing structures free of the effects of X-ray radiation damage when compared with standard cryo-crystallography. Comparison of these room-temperature methods demonstrated the large differences in sample requirements, data-collection time and the potential for radiation damage between them. With regard to the structure and function of DHP-B, despite the results being partly limited by differences in the underlying structures, new information was gained on the protonation states of active-site residues which may guide future studies of DHP-B.

Limits and potential of combined folding and docking

MOTIVATION

In the last decade, de novo protein structure prediction accuracy for individual proteins has improved significantly by utilising deep learning (DL) methods for harvesting the co-evolution information from large multiple sequence alignments (MSAs). The same approach can, in principle, also be used to extract information about evolutionary-based contacts across protein-protein interfaces. However, most earlier studies have not used the latest DL methods for inter-chain contact distance prediction. This article introduces a fold-and-dock method based on predicted residue-residue distances with trRosetta.

RESULTS

The method can simultaneously predict the tertiary and quaternary structure of a protein pair, even when the structures of the monomers are not known. The straightforward application of this method to a standard dataset for protein-protein docking yielded limited success. However, using alternative methods for generating MSAs allowed us to dock accurately significantly more proteins. We also introduced a novel scoring function, PconsDock, that accurately separates 98% of correctly and incorrectly folded and docked proteins. The average performance of the method is comparable to the use of traditional, template-based or ab initio shape-complementarity-only docking methods. Moreover, the results of conventional and fold-and-dock approaches are complementary, and thus a combined docking pipeline could increase overall docking success significantly. This methodology contributed to the best model for one of the CASP14 oligomeric targets, H1065.

AVAILABILITY AND IMPLEMENTATION

All scripts for predictions and analysis are available from https://github.com/ElofssonLab/bioinfo-toolbox/ and https://gitlab.com/ElofssonLab/benchmark5/. All models joined alignments, and evaluation results are available from the following figshare repository https://doi.org/10.6084/m9.figshare.14654886.v2.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Remdesivir overcomes the S861 roadblock in SARS-CoV-2 polymerase elongation complex

Remdesivir (RDV), a nucleotide analog with broad-spectrum features, has exhibited effectiveness in COVID-19 treatment. However, the precise working mechanism of RDV when targeting the viral RNA-dependent RNA polymerase (RdRP) has not been fully elucidated. Here, we solve a 3.0-Å structure of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RdRP elongation complex (EC) and assess RDV intervention in polymerase elongation phase. Although RDV could induce an “i+3” delayed termination in meta-stable complexes, only pausing and subsequent elongation are observed in the EC. A comparative investigation using an enterovirus RdRP further confirms similar delayed intervention and demonstrates that steric hindrance of the RDV-characteristic 1′-cyano at the −4 position is responsible for the “i+3” intervention, although two representative Flaviviridae RdRPs do not exhibit similar behavior. A comparison of representative viral RdRP catalytic complex structures indicates that the product RNA backbone encounters highly conserved structural elements, highlighting the broad-spectrum intervention potential of 1′-modified nucleotide analogs in anti-RNA virus drug development.

Getting to the bottom of lncRNA mechanism: structure-function relationships

While long non-coding RNAs are known to play key roles in disease and development, relatively few structural studies have been performed for this important class of RNAs. Here, we review functional studies of long non-coding RNAs and expose the need for high-resolution 3-D structural studies, discussing the roles of long non-coding RNAs in the cell and how structure–function relationships might be used to elucidate further understanding. We then describe structural studies of other classes of RNAs using chemical probing, nuclear magnetic resonance, small-angle X-ray scattering, X-ray crystallography, and cryogenic electron microscopy (cryo-EM). Next, we review early structural studies of long non-coding RNAs to date and describe the way forward for the structural biology of long non-coding RNAs in terms of cryo-EM.

Cryo-EM structure of mycobacterial cytochrome bd reveals two oxygen access channels

Cytochromes bd are ubiquitous amongst prokaryotes including many human-pathogenic bacteria. Such complexes are targets for the development of antimicrobial drugs. However, an understanding of the relationship between the structure and functional mechanisms of these oxidases is incomplete. Here, we have determined the 2.8 Å structure of Mycobacterium smegmatis cytochrome bd by single-particle cryo-electron microscopy. This bd oxidase consists of two subunits CydA and CydB, that adopt a pseudo two-fold symmetrical arrangement. The structural topology of its Q-loop domain, whose function is to bind the substrate, quinol, is significantly different compared to the C-terminal region reported for cytochromes bd from Geobacillus thermodenitrificans (G. th) and Escherichia coli (E. coli). In addition, we have identified two potential oxygen access channels in the structure and shown that similar tunnels also exist in G. th and E. coli cytochromes bd. This study provides insights to develop a framework for the rational design of antituberculosis compounds that block the oxygen access channels of this oxidase.

Cytochromes bd oxidase (Cyt-bd) catalyzes the reduction of oxygen to water and is the terminal oxidase in the respiratory chain of prokaryotes. Here, the authors present the 2.8 Å cryo-EM structure of Mycobacterium smegmatis Cyt-bd and identify two potential oxygen access channels in the structure, which is of interest for the development of novel antituberculosis drugs.

Plant "helper" immune receptors are Ca2+-permeable nonselective cation channels

Plant nucleotide-binding leucine-rich repeat receptors (NLRs) regulate immunity and cell death. In Arabidopsis, a subfamily of "helper" NLRs is required by many "sensor" NLRs. Active NRG1.1 oligomerized, was enriched in plasma membrane puncta, and conferred cytoplasmic calcium ion (Ca2+) influx in plant and human cells. NRG1.1-dependent Ca2+ influx and cell death were sensitive to Ca2+ channel blockers and were suppressed by mutations affecting oligomerization or plasma membrane enrichment. Ca2+ influx and cell death mediated by NRG1.1 and ACTIVATED DISEASE RESISTANCE 1 (ADR1), another helper NLR, required conserved negatively charged N-terminal residues. Whole-cell voltage-clamp recordings demonstrated that Arabidopsis helper NLRs form Ca2+-permeable cation channels to directly regulate cytoplasmic Ca2+ levels and consequent cell death. Thus, helper NLRs transduce cell death signals directly.

Catalytically impaired TrpA subunit of tryptophan synthase from Chlamydia trachomatis is an allosteric regulator of TrpB

Intracellular growth and pathogenesis of Chlamydia species is controlled by the availability of tryptophan, yet the complete biosynthetic pathway for l‐Trp is absent among members of the genus. Some representatives, however, preserve genes encoding tryptophan synthase, TrpAB – a bifunctional enzyme catalyzing the last two steps in l‐Trp synthesis. TrpA (subunit α) converts indole‐3‐glycerol phosphate into indole and glyceraldehyde‐3‐phosphate (α reaction). The former compound is subsequently used by TrpB (subunit β) to produce l‐Trp in the presence of l‐Ser and a pyridoxal 5′‐phosphate cofactor (β reaction). Previous studies have indicated that in Chlamydia, TrpA has lost its catalytic activity yet remains associated with TrpB to support the β reaction. Here, we provide detailed analysis of the TrpAB from C. trachomatis D/UW‐3/CX, confirming that accumulation of mutations in the active site of TrpA renders it enzymatically inactive, despite the conservation of the catalytic residues. We also show that TrpA remains a functional component of the TrpAB complex, increasing the activity of TrpB by four‐fold. The side chain of non‐conserved βArg267 functions as cation effector, potentially rendering the enzyme less susceptible to the solvent ion composition. The observed structural and functional changes detected herein were placed in a broader evolutionary and genomic context, allowing identification of these mutations in relation to their trp gene contexts in which they occur. Moreover, in agreement with the in vitro data, partial relaxation of purifying selection for TrpA, but not for TrpB, was detected, reinforcing a partial loss of TrpA functions during the course of evolution.

PDB Code(s): 6V82;

Architecture of the mycobacterial succinate dehydrogenase with a membrane-embedded Rieske FeS cluster

Targeting energy metabolism in Mycobacterium tuberculosis has emerged as a new paradigm in antituberculosis drug discovery. Succinate dehydrogenase is considered the regulator of respiration in M. tuberculosis. Mycobacteria contains two different succinate dehydrogenase enzymes designated Sdh1 and Sdh2. Sdh1 has recently been identified as a new class of succinate dehydrogenase. In this study, we have determined M. smegmatis Sdh1 structures alone and in the presence of ubiquinone-1, revealing that Sdh1 has a novel electron transfer pathway and a unique substrate-binding site. These data show that the structure of M. tuberculosis Sdh1 is significantly different by comparison with the human counterpart making a good antituberculosis drug target.

Complex II, also known as succinate dehydrogenase (SQR) or fumarate reductase (QFR), is an enzyme involved in both the Krebs cycle and oxidative phosphorylation. Mycobacterial Sdh1 has recently been identified as a new class of respiratory complex II (type F) but with an unknown electron transfer mechanism. Here, using cryoelectron microscopy, we have determined the structure of Mycobacterium smegmatis Sdh1 in the presence and absence of the substrate, ubiquinone-1, at 2.53-Å and 2.88-Å resolution, respectively. Sdh1 comprises three subunits, two that are water soluble, SdhA and SdhB, and one that is membrane spanning, SdhC. Within these subunits we identified a quinone-binding site and a rarely observed Rieske-type [2Fe-2S] cluster, the latter being embedded in the transmembrane region. A mutant, where two His ligands of the Rieske-type [2Fe-2S] were changed to alanine, abolished the quinone reduction activity of the Sdh1. Our structures allow the proposal of an electron transfer pathway that connects the substrate-binding and quinone-binding sites. Given the unique features of Sdh1 and its essential role in Mycobacteria, these structures will facilitate antituberculosis drug discovery efforts that specifically target this complex.

Cryo-EM structure of trimeric Mycobacterium smegmatis succinate dehydrogenase with a membrane-anchor SdhF

Diheme-containing succinate:menaquinone oxidoreductases (Sdh) are widespread in Gram-positive bacteria but little is known about the catalytic mechanisms they employ for succinate oxidation by menaquinone. Here, we present the 2.8 Å cryo-electron microscopy structure of a Mycobacterium smegmatis Sdh, which forms a trimer. We identified the membrane-anchored SdhF as a subunit of the complex. The 3 kDa SdhF forms a single transmembrane helix and this helix plays a role in blocking the canonically proximal quinone-binding site. We also identified two distal quinone-binding sites with bound quinones. One distal binding site is formed by neighboring subunits of the complex. Our structure further reveals the electron/proton transfer pathway for succinate oxidation by menaquinone. Moreover, this study provides further structural insights into the physiological significance of a trimeric respiratory complex II. The structure of the menaquinone binding site could provide a framework for the development of Sdh-selective anti-mycobacterial drugs.

Diheme-containing succinate:menaquinone oxidoreductases (Sdh) are members of the complex II superfamily. Here, the authors present the 2.8 Å cryo-EM structure of Mycobacterium smegmatis Sdh2, which reveals membrane-anchored SdhF as a component of the complex and they discuss the electron/proton transfer pathway in the Sdh2 trimer.

Benzylaminoethyureido-Tailed Benzenesulfonamides: Design, Synthesis, Kinetic and X-ray Investigations on Human Carbonic Anhydrases

A drug design strategy of carbonic anhydrase inhibitors (CAIs) belonging to sulfonamides incorporating ureidoethylaminobenzyl tails is presented. A variety of substitution patterns on the ring and the tails, located on para- or meta- positions with respect to the sulfonamide warheads were incorporated in the new compounds. Inhibition of human carbonic anhydrases (hCA) isoforms I, II, IX and XII, involving various pathologies, was assessed with the new compounds. Selective inhibitory profile towards hCA II was observed, the most active compounds being low nM inhibitors (KIs of 2.8-9.2 nM, respectively). Extensive X-ray crystallographic analysis of several sulfonamides in an adduct with hCA I allowed an in-depth understanding of their binding mode and to lay a detailed structure-activity relationship.

Structures of SALSA/DMBT1 SRCR domains reveal the conserved ligand-binding mechanism of the ancient SRCR fold

The structures of SALSA SRCR domains 1 and 8 reveal a cation-dependent mechanism for ligand recognition, contributing to important roles in the immune system and cellular signalling. The cation-binding sites are conserved across all SRCR domains, suggesting conserved functional mechanisms.

The scavenger receptor cysteine-rich (SRCR) family of proteins comprises more than 20 membrane-associated and secreted molecules. Characterised by the presence of one or more copies of the ∼110 amino-acid SRCR domain, this class of proteins have widespread functions as antimicrobial molecules, scavenger receptors, and signalling receptors. Despite the high level of structural conservation of SRCR domains, no unifying mechanism for ligand interaction has been described. The SRCR protein SALSA, also known as DMBT1/gp340, is a key player in mucosal immunology. Based on detailed structural data of SALSA SRCR domains 1 and 8, we here reveal a novel universal ligand-binding mechanism for SALSA ligands. The binding interface incorporates a dual cation-binding site, which is highly conserved across the SRCR superfamily. Along with the well-described cation dependency on most SRCR domain–ligand interactions, our data suggest that the binding mechanism described for the SALSA SRCR domains is applicable to all SRCR domains. We thus propose to have identified in SALSA a conserved functional mechanism for the SRCR class of proteins.

Structural basis of keto acid utilization in nonribosomal depsipeptide synthesis

Nonribosomal depsipeptides are natural products composed of amino and hydroxy acid residues. The hydroxy acid residues often derive from α-keto acids, reduced by ketoreductase domains in the depsipeptide synthetases. Biochemistry and structures reveal the mechanism of discrimination for α-keto acids and a remarkable architecture: flanking intact adenylation and ketoreductase domains are sequences separated by >1,100 residues that form a split 'pseudoA_sub' domain, structurally important for the depsipeptide module's synthetic cycle.

The Crystal Structure of the Manganese Superoxide Dismutase from Geobacillus stearothermophilus: Parker and Blake (1988) Revisited

Superoxide dismutase (SOD) is an almost ubiquitous metalloenzyme in aerobic organisms that catalyses the disproportionation of superoxide. Geobacillus stearothermophilus MnSOD is the only published MnSOD structure that does not have its coordinates publicly available, yet it is one of the more cited structures in the SOD literature. The structure has now been refined with modern programs, yielding a significantly improved structure which has been deposited in the Protein Data Bank. Importantly, the further refined structure reveals the presence of a catalytically important fifth ligand, water, to the metal centre, as observed in other SOD structures.

Passenger sequences can promote interlaced dimers in a common variant of the maltose-binding protein

The maltose-binding protein (MBP) is one of the most frequently used protein tags due to its capacity to stabilize, solubilize and even crystallize recombinant proteins that are fused to it. Given that MBP is thought to be a highly stable monomeric protein with known characteristics, fused passenger proteins are often studied without being cleaved from MBP. Here we report that a commonly used engineered MBP version (mutated to lower its surface entropy) can form interlaced dimers when fused to short protein sequences derived from the focal adhesion kinase (FAK) or the homologous protein tyrosine kinase 2 (PYK2). These MBP dimers still bind maltose and can interconvert with monomeric forms in vitro under standard conditions despite a contact surface of more than 11,000 Ų. We demonstrate that both the mutations in MBP and the fused protein sequences were required for dimer formation. The FAK and PYK2 sequences are less than 40% identical, monomeric, and did not show specific interactions with MBP, suggesting that a variety of sequences can promote this MBP dimerization. MBP dimerization was abrogated by reverting two of the eight mutations introduced in the engineered MBP. Our results provide an extreme example for induced reversible domain-swapping, with implications for protein folding dynamics. Our observations caution that passenger-promoted MBP dimerization might mislead experimental characterization of the fused protein sequences, but also suggest a simple mutation to stop this phenomenon.

An inhibitor of complement C5 provides structural insights into activation

The complement system is a crucial antimicrobial system in the human body. However, controlling its regulation is essential, and failure to do so is implicated in a number of clinical inflammatory pathologies leading to great interest in therapeutic complement inhibition. We have identified and characterized a class of complement inhibitors from biting ticks. Utilizing both cryoelectron microscopy and X-ray crystallography we provide a comprehensive understanding of their mechanism of inhibition at the level of the terminal pathway of complement. We present a high-resolution cryoelectron microscopy structure of complement C5, the molecule targeted by the major therapeutic Eculizumab. In addition, we reveal the fold of the CirpT family of tick inhibitors and their unique mode of inhibition.

The complement system is a crucial part of innate immune defenses against invading pathogens. The blood-meal of the tick Rhipicephalus pulchellus lasts for days, and the tick must therefore rely on inhibitors to counter complement activation. We have identified a class of inhibitors from tick saliva, the CirpT family, and generated detailed structural data revealing their mechanism of action. We show direct binding of a CirpT to complement C5 and have determined the structure of the C5–CirpT complex by cryoelectron microscopy. This reveals an interaction with the peripheral macro globulin domain 4 (C5_MG4) of C5. To achieve higher resolution detail, the structure of the C5_MG4–CirpT complex was solved by X-ray crystallography (at 2.7 Å). We thus present the fold of the CirpT protein family, and provide detailed mechanistic insights into its inhibitory function. Analysis of the binding interface reveals a mechanism of C5 inhibition, and provides information to expand our biological understanding of the activation of C5, and thus the terminal complement pathway.

Molstack: A platform for interactive presentations of electron density and cryo-EM maps and their interpretations

In the Special Issue on Tools for Protein Science in 2018, we presented Molstack: a concept of a cloud-based platform for sharing electron density maps and their interpretations. Molstack is a web platform that allows the interactive visualization of density maps through the simultaneous presentation of multiple datasets and models in a way that allows for easy pairwise comparison. We anticipated that the users of this conceptually simple platform would find many different uses for their projects, and we were not mistaken. We have observed researchers use Molstack to present experimental evidence for their models in the form of electron density maps, omit maps, and anomalous difference density maps. Users also employed Molstack to present alternative interpretations of densities, including rerefinements and speculative interpretations. While we anticipated these types of projects to be the main use cases, we were pleased to see Molstack used to display superpositions of different models, as a tool for story-driven presentations, and for collaboration as well. Here, we present developments in the platform that were driven by user feedback, highlight several cases that used Molstack to enhance the publication, and discuss possible directions for the platform.

The Molecular Mechanism of Cellular Attachment for an Archaeal Virus

Sulfolobus turreted icosahedral virus (STIV) is a model archaeal virus and member of the PRD1-adenovirus lineage. Although STIV employs pyramidal lysis structures to exit the host, knowledge of the viral entry process is lacking. We therefore initiated studies on STIV attachment and entry. Negative stain and cryoelectron micrographs showed virion attachment to pili-like structures emanating from the Sulfolobus host. Tomographic reconstruction and sub-tomogram averaging revealed pili recognition by the STIV C381 turret protein. Specifically, the triple jelly roll structure of C381 determined by X-ray crystallography shows that pilus recognition is mediated by conserved surface residues in the second and third domains. In addition, the STIV petal protein (C557), when present, occludes the pili binding site, suggesting that it functions as a maturation protein. Combined, these results demonstrate a role for the namesake STIV turrets in initial cellular attachment and provide the first molecular model for viral attachment in the archaeal domain of life.

NAD+ cleavage activity by animal and plant TIR domains in cell death pathways

SARM1 (sterile alpha and TIR motif containing 1) is responsible for depletion of nicotinamide adenine dinucleotide in its oxidized form (NAD⁺) during Wallerian degeneration associated with neuropathies. Plant nucleotide-binding leucine-rich repeat (NLR) immune receptors recognize pathogen effector proteins and trigger localized cell death to restrict pathogen infection. Both processes depend on closely related Toll/interleukin-1 receptor (TIR) domains in these proteins, which, as we show, feature self-association-dependent NAD⁺ cleavage activity associated with cell death signaling. We further show that SARM1 SAM (sterile alpha motif) domains form an octamer essential for axon degeneration that contributes to TIR domain enzymatic activity. The crystal structures of ribose and NADP⁺ (the oxidized form of nicotinamide adenine dinucleotide phosphate) complexes of SARM1 and plant NLR RUN1 TIR domains, respectively, reveal a conserved substrate binding site. NAD⁺ cleavage by TIR domains is therefore a conserved feature of animal and plant cell death signaling pathways.

Conservation of the structure and function of bacterial tryptophan synthases

Tryptophan biosynthesis is one of the most characterized processes in bacteria, in which the enzymes from Salmonella typhimurium and Escherichia coli serve as model systems. Tryptophan synthase (TrpAB) catalyzes the final two steps of tryptophan biosynthesis in plants, fungi and bacteria. This pyridoxal 5'-phosphate (PLP)-dependent enzyme consists of two protein chains, α (TrpA) and β (TrpB), functioning as a linear αββα heterotetrameric complex containing two TrpAB units. The reaction has a complicated, multistep mechanism resulting in the β-replacement of the hydroxyl group of l-serine with an indole moiety. Recent studies have shown that functional TrpAB is required for the survival of pathogenic bacteria in macrophages and for evading host defense. Therefore, TrpAB is a promising target for drug discovery, as its orthologs include enzymes from the important human pathogens Streptococcus pneumoniae, Legionella pneumophila and Francisella tularensis, the causative agents of pneumonia, legionnaires' disease and tularemia, respectively. However, specific biochemical and structural properties of the TrpABs from these organisms have not been investigated. To fill the important phylogenetic gaps in the understanding of TrpABs and to uncover unique features of TrpAB orthologs to spearhead future drug-discovery efforts, the TrpABs from L. pneumophila, F. tularensis and S. pneumoniae have been characterized. In addition to kinetic properties and inhibitor-sensitivity data, structural information gathered using X-ray crystallo-graphy is presented. The enzymes show remarkable structural conservation, but at the same time display local differences in both their catalytic and allosteric sites that may be responsible for the observed differences in catalysis and inhibitor binding. This functional dissimilarity may be exploited in the design of species-specific enzyme inhibitors.

Structural insight into TRPV5 channel function and modulation

TRPV5 (transient receptor potential vanilloid 5) is a unique calcium-selective TRP channel essential for calcium homeostasis. Unlike other TRPV channels, TRPV5 and its close homolog, TRPV6, do not exhibit thermosensitivity or ligand-dependent activation but are constitutively open at physiological membrane potentials and modulated by calmodulin (CaM) in a calcium-dependent manner. Here we report high-resolution electron cryomicroscopy structures of truncated and full-length TRPV5 in lipid nanodiscs, as well as of a TRPV5 W583A mutant and TRPV5 in complex with CaM. These structures highlight the mechanism of calcium regulation and reveal a flexible stoichiometry of CaM binding to TRPV5.

A model of molecular vector machine of proteins

Based on the analysis of N_iH⋯O_i-4C bond region of the proteins pentafragments using the quantum chemistry representations, a group of 20 vectors of action has been isolated. These vectors of action are reconstructed by a group of 20 side chains of amino acids being the irreducible representations of vectors. As a result, a model of molecular vector machine of proteins (MVM) has been proposed. This machine includes proteins pentafragments, a system of 20 vectors inside dodecahedron coming out of the O_i-4 atom as from the center and tetrahedral α_i-carbon atom with the side chains of amino acids that change as the polypeptide chain grows. A model with the structured canonical set of amino acids based on the principles of antisymmetry was transferred onto the dodecahedron as an MVM element (Karasev et al., 2005; Karasev et al., 2007). The perspective of MVM application for predicting and designing the secondary structure of proteins has been discussed.

Structure of Calcarisporiella thermophila Hsp104 Disaggregase that Antagonizes Diverse Proteotoxic Misfolding Events

Hsp104 is an AAA+ protein disaggregase with powerful amyloid-remodeling activity. All nonmetazoan eukaryotes express Hsp104 while eubacteria express an Hsp104 ortholog, ClpB. However, most studies have focused on Hsp104 from Saccharomyces cerevisiae and ClpB orthologs from two eubacterial species. Thus, the natural spectrum of Hsp104/ClpB molecular architectures and protein-remodeling activities remains largely unexplored. Here, we report two structures of Hsp104 from the thermophilic fungus Calcarisporiella thermophila (CtHsp104), a 2.70Å crystal structure and 4.0Å cryo-electron microscopy structure. Both structures reveal left-handed, helical assemblies with all domains clearly resolved. We thus provide the highest resolution and most complete view of Hsp104 hexamers to date. We also establish that CtHsp104 antagonizes several toxic protein-misfolding events in vivo where S. cerevisiae Hsp104 is ineffective, including rescue of TDP-43, polyglutamine, and α-synuclein toxicity. We suggest that natural Hsp104 variation is an invaluable, untapped resource for illuminating therapeutic disaggregases for fatal neurodegenerative diseases.

Combining Evolutionary Covariance and NMR Data for Protein Structure Determination

Accurate protein structure determination by solution-state NMR is challenging for proteins greater than about 20kDa, for which extensive perdeuteration is generally required, providing experimental data that are incomplete (sparse) and ambiguous. However, the massive increase in evolutionary sequence information coupled with advances in methods for sequence covariance analysis can provide reliable residue-residue contact information for a protein from sequence data alone. These "evolutionary couplings (ECs)" can be combined with sparse NMR data to determine accurate 3D protein structures. This hybrid "EC-NMR" method has been developed using NMR data for several soluble proteins and validated by comparison with corresponding reference structures determined by X-ray crystallography and/or conventional NMR methods. For small proteins, only backbone resonance assignments are utilized, while for larger proteins both backbone and some sidechain methyl resonance assignments are generally required. ECs can be combined with sparse NMR data obtained on deuterated, selectively protonated protein samples to provide structures that are more accurate and complete than those obtained using such sparse NMR data alone. EC-NMR also has significant potential for analysis of protein structures from solid-state NMR data and for studies of integral membrane proteins. The requirement that ECs are consistent with NMR data recorded on a specific member of a protein family, under specific conditions, also allows identification of ECs that reflect alternative allosteric or excited states of the protein structure.

The Protozoan Trichomonas vaginalis Targets Bacteria with Laterally Acquired NlpC/P60 Peptidoglycan Hydrolases

The human eukaryotic pathogen Trichomonas vaginalis causes trichomoniasis, a prevalent sexually transmitted infection. This extracellular protozoan is intimately associated with the human vaginal mucosa and microbiota, but key aspects of the complex interactions between the parasite and the vaginal bacteria remain elusive. We report that T. vaginalis has acquired, by lateral gene transfer from bacteria, genes encoding peptidoglycan hydrolases of the NlpC/P60 family. Two of the T. vaginalis enzymes were active against bacterial peptidoglycan, retaining the active-site fold and specificity as dl-endopeptidases. The endogenous NlpC/P60 genes are transcriptionally upregulated in T. vaginalis in the presence of bacteria. The overexpression of an exogenous copy enables the parasite to outcompete bacteria from mixed cultures, consistent with the biochemical activity of the enzyme. Our study results highlight the relevance of the interactions of this eukaryotic pathogen with bacteria, a poorly understood aspect of the biology of this important human parasite.IMPORTANCETrichomonas vaginalis is a parasitic protozoan of the human urogenital tract that causes trichomoniasis, a very common sexually transmitted disease. Despite residing extracellularly and in close association with the vaginal bacteria (i.e., the microbiota), very little is known about the nature of the parasite-bacterium interactions. Our study showed that this parasite had acquired genes from bacteria which retained their original function. They produce active enzymes capable of degrading peptidoglycan, a unique polymer of the bacterial cell envelope, helping the parasite to outcompete bacteria in mixed cultures. This study was the first to show that a laterally acquired group of genes enables a eukaryotic mucosal pathogen to control bacterial population. We highlight the importance of understanding the interactions between pathogens and microbiota, as the outcomes of these interactions are increasingly understood to have important implications on health and disease.

Conformations of an RNA Helix-Junction-Helix Construct Revealed by SAXS Refinement of MD Simulations

RNA is involved in a broad range of biological processes that extend far beyond translation. Many of RNA's recently discovered functions rely on folding to a specific conformation or transitioning between conformations. The RNA structure contains rigid, short basepaired regions connected by more flexible linkers. Studies of model constructs such as small helix-junction-helix (HJH) motifs are useful in understanding how these elements work together to determine RNA conformation. Here, we reveal the full ensemble of solution structures assumed by a model RNA HJH. We apply small-angle x-ray scattering and an ensemble optimization method to selectively refine models generated by all-atom molecular dynamics simulations. The expectation of a broad distribution of helix orientations, at and above physiological ionic strength, is not met. Instead, this analysis shows that the HJH structures are dominated by two distinct conformations at moderate to high ionic strength. Atomic structures, selected from the molecular dynamics simulations, reveal strong base-base interactions in the junction that critically constrain the conformational space available to the HJH molecule and lead to a surprising re-extension at high salt. These results are corroborated by comparison with previous single-molecule fluorescence resonance energy transfer experiments on the same constructs.

A helical LC3-interacting region mediates the interaction between the retroviral restriction factor Trim5α and mammalian autophagy-related ATG8 proteins

The retroviral restriction factor tripartite motif–containing 5α (Trim5α) acts during the early postentry stages of the retroviral life cycle to block infection by a broad range of retroviruses, disrupting reverse transcription and integration. The mechanism of this restriction is poorly understood, but it has recently been suggested to involve recruitment of components of the autophagy machinery, including members of the mammalian autophagy-related 8 (ATG8) family involved in targeting proteins to the autophagosome. To better understand the molecular details of this interaction, here we utilized analytical ultracentrifugation to characterize the binding of six ATG8 isoforms and determined the crystal structure of the Trim5α Bbox coiled-coil region in complex with one member of the mammalian ATG8 proteins, autophagy-related protein LC3 B (LC3B). We found that Trim5α binds all mammalian ATG8s and that, unlike the typical LC3-interacting region (LIR) that binds to mammalian ATG8s through a β-strand motif comprising approximately six residues, LC3B binds to Trim5α via the α-helical coiled-coil region. The orientation of the structure demonstrated that LC3B could be accommodated within a Trim5α assembly that can bind the retroviral capsid. However, mutation of the binding interface does not affect retroviral restriction. Comparison of the typical linear β-strand LIR with our atypical helical LIR reveals a conservation of the presentation of residues that are required for the interaction with LC3B. This observation expands the range of LC3B-binding proteins to include helical binding motifs and demonstrates a link between Trim5α and components of the autophagosome.

Model validation: local diagnosis, correction and when to quit

Traditionally, validation was considered to be a final gatekeeping function, but refinement is smoother and results are better if model validation actively guides corrections throughout structure solution. This shifts emphasis from global to local measures: primarily geometry, conformations and sterics. A fit into the wrong local minimum conformation usually produces outliers in multiple measures. Moving to the right local minimum should be prioritized, rather than small shifts across arbitrary borderlines. Steric criteria work best with all explicit H atoms. `Backrub' motions should be used for side chains and `P-perp' diagnostics to correct ribose puckers. A `water' may actually be an ion, a relic of misfitting or an unmodeled alternate. Beware of wishful thinking in modeling ligands. At high resolution, internally consistent alternate conformations should be modeled and geometry in poor density should not be downweighted. At low resolution, CaBLAM should be used to diagnose protein secondary structure and ERRASER to correct RNA backbone. All atoms should not be forced inside density, beware of sequence misalignment, and very rare conformations such as cis-non-Pro peptides should be avoided. Automation continues to improve, but the crystallographer still must look at each outlier, in the context of density, and correct most of them. For the valid few with unambiguous density and something that is holding them in place, a functional reason should be sought. The expectation is a few outliers, not zero.

The ancestral retinoic acid receptor was a low-affinity sensor triggering neuronal differentiation

Retinoic acid (RA) is an important intercellular signaling molecule in vertebrate development, with a well-established role in the regulation of hox genes during hindbrain patterning and in neurogenesis. However, the evolutionary origin of the RA signaling pathway remains elusive. To elucidate the evolution of the RA signaling system, we characterized RA metabolism and signaling in the marine annelid Platynereis dumerilii, a powerful model for evolution, development, and neurobiology. Binding assays and crystal structure analyses show that the annelid retinoic acid receptor (RAR) binds RA and activates transcription just as vertebrate RARs, yet with a different ligand-binding pocket and lower binding affinity, suggesting a permissive rather than instructive role of RA signaling. RAR knockdown and RA treatment of swimming annelid larvae further reveal that the RA signal is locally received in the medial neuroectoderm, where it controls neurogenesis and axon outgrowth, whereas the spatial colinear hox gene expression in the neuroectoderm remains unaffected. These findings suggest that one early role of the new RAR in bilaterian evolution was to control the spatially restricted onset of motor and interneuron differentiation in the developing ventral nerve cord and to indicate that the regulation of hox-controlled anterior-posterior patterning arose only at the base of the chordates, concomitant with a high-affinity RAR needed for the interpretation of a complex RA gradient.

Theoretical restrictions on longest implicit time scales in Markov state models of biomolecular dynamics

Markov state models (MSMs) have been widely used to analyze computer simulations of various biomolecular systems. They can capture conformational transitions much slower than an average or maximal length of a single molecular dynamics (MD) trajectory from the set of trajectories used to build the MSM. A rule of thumb claiming that the slowest implicit time scale captured by an MSM should be comparable by the order of magnitude to the aggregate duration of all MD trajectories used to build this MSM has been known in the field. However, this rule has never been formally proved. In this work, we present analytical results for the slowest time scale in several types of MSMs, supporting the above rule. We conclude that the slowest implicit time scale equals the product of the aggregate sampling and four factors that quantify: (1) how much statistics on the conformational transitions corresponding to the longest implicit time scale is available, (2) how good the sampling of the destination Markov state is, (3) the gain in statistics from using a sliding window for counting transitions between Markov states, and (4) a bias in the estimate of the implicit time scale arising from finite sampling of the conformational transitions. We demonstrate that in many practically important cases all these four factors are on the order of unity, and we analyze possible scenarios that could lead to their significant deviation from unity. Overall, we provide for the first time analytical results on the slowest time scales captured by MSMs. These results can guide further practical applications of MSMs to biomolecular dynamics and allow for higher computational efficiency of simulations.

MolProbity: More and better reference data for improved all-atom structure validation

This paper describes the current update on macromolecular model validation services that are provided at the MolProbity website, emphasizing changes and additions since the previous review in 2010. There have been many infrastructure improvements, including rewrite of previous Java utilities to now use existing or newly written Python utilities in the open-source CCTBX portion of the Phenix software system. This improves long-term maintainability and enhances the thorough integration of MolProbity-style validation within Phenix. There is now a complete MolProbity mirror site at http://molprobity.manchester.ac.uk. GitHub serves our open-source code, reference datasets, and the resulting multi-dimensional distributions that define most validation criteria. Coordinate output after Asn/Gln/His "flip" correction is now more idealized, since the post-refinement step has apparently often been skipped in the past. Two distinct sets of heavy-atom-to-hydrogen distances and accompanying van der Waals radii have been researched and improved in accuracy, one for the electron-cloud-center positions suitable for X-ray crystallography and one for nuclear positions. New validations include messages at input about problem-causing format irregularities, updates of Ramachandran and rotamer criteria from the million quality-filtered residues in a new reference dataset, the CaBLAM Cα-CO virtual-angle analysis of backbone and secondary structure for cryoEM or low-resolution X-ray, and flagging of the very rare cis-nonProline and twisted peptides which have recently been greatly overused. Due to wide application of MolProbity validation and corrections by the research community, in Phenix, and at the worldwide Protein Data Bank, newly deposited structures have continued to improve greatly as measured by MolProbity's unique all-atom clashscore.

Structure of a novel antibacterial toxin that exploits elongation factor Tu to cleave specific transfer RNAs

Contact-dependent growth inhibition (CDI) is a mechanism of inter-cellular competition in which Gram-negative bacteria exchange polymorphic toxins using type V secretion systems. Here, we present structures of the CDI toxin from Escherichia coli NC101 in ternary complex with its cognate immunity protein and elongation factor Tu (EF-Tu). The toxin binds exclusively to domain 2 of EF-Tu, partially overlapping the site that interacts with the 3'-end of aminoacyl-tRNA (aa-tRNA). The toxin exerts a unique ribonuclease activity that cleaves the single-stranded 3'-end from tRNAs that contain guanine discriminator nucleotides. EF-Tu is required to support this tRNase activity in vitro, suggesting the toxin specifically cleaves substrate in the context of GTP·EF-Tu·aa-tRNA complexes. However, superimposition of the toxin domain onto previously solved GTP·EF-Tu·aa-tRNA structures reveals potential steric clashes with both aa-tRNA and the switch I region of EF-Tu. Further, the toxin induces conformational changes in EF-Tu, displacing a β-hairpin loop that forms a critical salt-bridge contact with the 3'-terminal adenylate of aa-tRNA. Together, these observations suggest that the toxin remodels GTP·EF-Tu·aa-tRNA complexes to free the 3'-end of aa-tRNA for entry into the nuclease active site.

Molstack-Interactive visualization tool for presentation, interpretation, and validation of macromolecules and electron density maps

Our understanding of the world of biomolecular structures is based upon the interpretation of macromolecular models, of which ∼90% are an interpretation of electron density maps. This structural information guides scientific progress and exploration in many biomedical disciplines. The Protein Data Bank's web portals have made these structures available for mass scientific consumption and greatly broaden the scope of information presented in scientific publications. The portals provide numerous quality metrics; however, the portion of the structure that is most vital for interpretation of the function may have the most difficult to interpret electron density and this ambiguity is not reflected by any single metric. The possible consequences of basing research on suboptimal models make it imperative to inspect the agreement of a model with its experimental evidence. Molstack, a web-based interactive publishing platform for structural data, allows users to present density maps and structural models by displaying a collection of maps and models, including different interpretation of one's own data, re-refinements, and corrections of existing structures. Molstack organizes the sharing and dissemination of these structural models along with their experimental evidence as an interactive session. Molstack was designed with three groups of users in mind; researchers can present the evidence of their interpretation, reviewers and readers can independently judge the experimental evidence of the authors' conclusions, and other researchers can present or even publish their new hypotheses in the context of prior results. The server is available at http://molstack.bioreproducibility.org.

Crystal Structure of Leiomodin 2 in Complex with Actin: A Structural and Functional Reexamination

Leiomodins (Lmods) are a family of actin filament nucleators related to tropomodulins (Tmods), which are pointed end-capping proteins. Whereas Tmods have alternating tropomyosin- and actin-binding sites (TMBS1, ABS1, TMBS2, ABS2), Lmods lack TMBS2 and half of ABS1, and present a C-terminal extension containing a proline-rich domain and an actin-binding Wiskott-Aldrich syndrome protein homology 2 (WH2) domain that is absent in Tmods. Most of the nucleation activity of Lmods resides within a fragment encompassing ABS2 and the C-terminal extension. This fragment recruits actin monomers into a polymerization nucleus. Here, we revise a recently reported structure of this region of Lmod2 in complex with actin and provide biochemical validation for the newly revised structure. We find that instead of two actin subunits connected by a single Lmod2 polypeptide, as reported in the original structure, the P1 unit cell contains two nearly identical copies of actin monomers, each bound to Lmod2's ABS2 and WH2 domain, with no electron density connecting these two domains. Moreover, we show that the two actin molecules in the unit cell are related to each other by a local twofold noncrystallographic symmetry axis, a conformation clearly distinct from that of actin subunits in the helical filament. We further find that a proposed actin-binding site within the missing connecting region of Lmod2, termed helix h1, does not bind actin in vitro and that the electron density assigned to it in the original structure corresponds instead to a WH2 domain with opposite backbone directionality. Polymerization assays using Lmod2 mutants of helix h1 and the WH2 domain support this conclusion. Finally, we find that deleting the C-terminal extension of Lmod1 and Lmod2 results in an approximately threefold decrease in the nucleation activity, which is only partially accounted for by the lack of the WH2 domain.

PDEδ Binding to Ras Isoforms Provides a Route to Proper Membrane Localization

To signal, Ras isoforms must be enriched at the plasma membrane (PM). It was suggested that phosphodiesterase-δ (PDEδ) can bind and shuttle some farnesylated Ras isoforms to the PM, but not all. Among these, interest focused on K-Ras4B, the most abundant oncogenic Ras isoform. To study PDEδ/Ras interactions, we modeled and simulated the PDEδ/K-Ras4B complex. We obtained structures, which were similar to two subsequently determined crystal structures. We next modeled and simulated complexes of PDEδ with the farnesylated hypervariable regions of K-Ras4A and N-Ras. Earlier data suggested that PDEδ extracts K-Ras4B and N-Ras from the PM, but surprisingly not K-Ras4A. Earlier analysis of the crystal structures advanced that the presence of large/charged residues adjacent to the farnesylated site precludes the PDEδ interaction. Here, we show that PDEδ can bind to farnesylated K-Ras4A and N-Ras like K-Ras4B, albeit not as strongly. This weaker binding, coupled with the stronger anchoring of K-Ras4A in the membrane (but not of electrostatically neutral N-Ras), can explain the observation why PDEδ is unable to effectively extract K-Ras4A. We thus propose that farnesylated Ras isoforms can bind PDEδ to fulfill the required PM enrichment, and argue that the different environments, PM versus solution, can resolve apparently puzzling Ras observations. These are novel insights that would not be expected based on the crystal structures alone, which provide an elegant rationale for previously puzzling observations of the differential effects of PDEδ on farnesylated Ras family proteins.

Databases, Repositories, and Other Data Resources in Structural Biology

Structural biology, like many other areas of modern science, produces an enormous amount of primary, derived, and "meta" data with a high demand on data storage and manipulations. Primary data come from various steps of sample preparation, diffraction experiments, and functional studies. These data are not only used to obtain tangible results, like macromolecular structural models, but also to enrich and guide our analysis and interpretation of various biomedical problems. Herein we define several categories of data resources, (a) Archives, (b) Repositories, (c) Databases, and (d) Advanced Information Systems, that can accommodate primary, derived, or reference data. Data resources may be used either as web portals or internally by structural biology software. To be useful, each resource must be maintained, curated, as well as integrated with other resources. Ideally, the system of interconnected resources should evolve toward comprehensive "hubs", or Advanced Information Systems. Such systems, encompassing the PDB and UniProt, are indispensable not only for structural biology, but for many related fields of science. The categories of data resources described herein are applicable well beyond our usual scientific endeavors.

In-crystal reaction cycle of a toluene-bound diiron hydroxylase

Electrophilic aromatic substitution is one of the most important and recognizable classes of organic chemical transformation. Enzymes create the strong electrophiles that are needed for these highly energetic reactions by using O₂, electrons, and metals or other cofactors. Although the nature of the oxidants that carry out electrophilic aromatic substitution has been deduced from many approaches, it has been difficult to determine their structures. Here we show the structure of a diiron hydroxylase intermediate formed during a reaction with toluene. Density functional theory geometry optimizations of an active site model reveal that the intermediate is an arylperoxo Fe²⁺/Fe³⁺ species with delocalized aryl radical character. The structure suggests that a carboxylate ligand of the diiron centre may trigger homolytic cleavage of the O-O bond by transferring a proton from a metal-bound water. Our work provides the spatial and electronic constraints needed to propose a comprehensive mechanism for diiron enzyme arene hydroxylation that accounts for many prior experimental results.

Diverse modes of galacto-specific carbohydrate recognition by a family 31 glycoside hydrolase from Clostridium perfringens

Clostridium perfringens is a commensal member of the human gut microbiome and an opportunistic pathogen whose genome encodes a suite of putative large, multi-modular carbohydrate-active enzymes that appears to play a role in the interaction of the bacterium with mucin-based carbohydrates. Among the most complex of these is an enzyme that contains a presumed catalytic module belonging to glycoside hydrolase family 31 (GH31). This large enzyme, which based on its possession of a GH31 module is a predicted α-glucosidase, contains a variety of non-catalytic ancillary modules, including three CBM32 modules that to date have not been characterized. NMR-based experiments demonstrated a preference of each module for galacto-configured sugars, including the ability of all three CBM32s to recognize the common mucin monosaccharide GalNAc. X-ray crystal structures of the CpGH31 CBM32s, both in apo form and bound to GalNAc, revealed the finely-tuned molecular strategies employed by these sequentially variable CBM32s in coordinating a common ligand. The data highlight that sequence similarities to previously characterized CBMs alone are insufficient for identifying the molecular mechanism of ligand binding by individual CBMs. Furthermore, the overlapping ligand binding profiles of the three CBMs provide a fail-safe mechanism for the recognition of GalNAc among the dense eukaryotic carbohydrate networks of the colonic mucosa. These findings expand our understanding of ligand targeting by large, multi-modular carbohydrate-active enzymes, and offer unique insights into of the expanding ligand-binding preferences and binding site topologies observed in CBM32s.

The integrative role of cryo electron microscopy in molecular and cellular structural biology

After gradually moving away from preparation methods prone to artefacts such as plastic embedding and negative staining for cell sections and single particles, the field of cryo electron microscopy (cryo-EM) is now heading off at unprecedented speed towards high-resolution analysis of biological objects of various sizes. This 'revolution in resolution' is happening largely thanks to new developments of new-generation cameras used for recording the images in the cryo electron microscope which have much increased sensitivity being based on complementary metal oxide semiconductor devices. Combined with advanced image processing and 3D reconstruction, the cryo-EM analysis of nucleoprotein complexes can provide unprecedented insights at molecular and atomic levels and address regulatory mechanisms in the cell. These advances reinforce the integrative role of cryo-EM in synergy with other methods such as X-ray crystallography, fluorescence imaging or focussed-ion beam milling as exemplified here by some recent studies from our laboratory on ribosomes, viruses, chromatin and nuclear receptors. Such multi-scale and multi-resolution approaches allow integrating molecular and cellular levels when applied to purified or in situ macromolecular complexes, thus illustrating the trend of the field towards cellular structural biology.

Computational Methodologies for Real-Space Structural Refinement of Large Macromolecular Complexes

The rise of the computer as a powerful tool for model building and refinement has revolutionized the field of structure determination for large biomolecular systems. Despite the wide availability of robust experimental methods capable of resolving structural details across a range of spatiotemporal resolutions, computational hybrid methods have the unique ability to integrate the diverse data from multimodal techniques such as X-ray crystallography and electron microscopy into consistent, fully atomistic structures. Here, commonly employed strategies for computational real-space structural refinement are reviewed, and their specific applications are illustrated for several large macromolecular complexes: ribosome, virus capsids, chemosensory array, and photosynthetic chromatophore. The increasingly important role of computational methods in large-scale structural refinement, along with current and future challenges, is discussed.