Structure and functionality of a multimeric human COQ7:COQ9 complex.

Coenzyme Q (CoQ) is a redox-active lipid essential for core metabolic pathways and antioxidant defense. CoQ is synthesized upon the mitochondrial inner membrane by an ill-defined "complex Q" metabolon. Here, we present structure-function analyses of a lipid-, substrate-, and NADH-bound complex comprising two complex Q subunits: the hydroxylase COQ7 and the lipid-binding protein COQ9. We reveal that COQ7 adopts a ferritin-like fold with a hydrophobic channel whose substrate-binding capacity is enhanced by COQ9. Using molecular dynamics, we further show that two COQ7:COQ9 heterodimers form a curved tetramer that deforms the membrane, potentially opening a pathway for the CoQ intermediates to translocate from the bilayer to the proteins' lipid-binding sites. Two such tetramers assemble into a soluble octamer with a pseudo-bilayer of lipids captured within. Together, these observations indicate that COQ7 and COQ9 cooperate to access hydrophobic precursors within the membrane and coordinate subsequent synthesis steps toward producing CoQ.

Structural mechanism of mitochondrial membrane remodelling by human OPA1.

Distinct morphologies of the mitochondrial network support divergent metabolic and regulatory processes that determine cell function and fate1-3. The mechanochemical GTPase optic atrophy 1 (OPA1) influences the architecture of cristae and catalyses the fusion of the mitochondrial inner membrane4,5. Despite its fundamental importance, the molecular mechanisms by which OPA1 modulates mitochondrial morphology are unclear. Here, using a combination of cellular and structural analyses, we illuminate the molecular mechanisms that are key to OPA1-dependent membrane remodelling and fusion. Human OPA1 embeds itself into cardiolipin-containing membranes through a lipid-binding paddle domain. A conserved loop within the paddle domain inserts deeply into the bilayer, further stabilizing the interactions with cardiolipin-enriched membranes. OPA1 dimerization through the paddle domain promotes the helical assembly of a flexible OPA1 lattice on the membrane, which drives mitochondrial fusion in cells. Moreover, the membrane-bending OPA1 oligomer undergoes conformational changes that pull the membrane-inserting loop out of the outer leaflet and contribute to the mechanics of membrane remodelling. Our findings provide a structural framework for understanding how human OPA1 shapes mitochondrial morphology and show us how human disease mutations compromise OPA1 functions.

An Isoprene Lipid-Binding Protein Promotes Eukaryotic Coenzyme Q Biosynthesis.

The biosynthesis of coenzyme Q presents a paradigm for how cells surmount hydrophobic barriers in lipid biology. In eukaryotes, CoQ precursors-among nature's most hydrophobic molecules-must somehow be presented to a series of enzymes peripherally associated with the mitochondrial inner membrane. Here, we reveal that this process relies on custom lipid-binding properties of COQ9. We show that COQ9 repurposes the bacterial TetR fold to bind aromatic isoprenes with high specificity, including CoQ intermediates that likely reside entirely within the bilayer. We reveal a process by which COQ9 associates with cardiolipin-rich membranes and warps the membrane surface to access this cargo. Finally, we identify a molecular interface between COQ9 and the hydroxylase COQ7, motivating a model whereby COQ9 presents intermediates directly to CoQ enzymes. Overall, our results provide a mechanism for how a lipid-binding protein might access, select, and deliver specific cargo from a membrane to promote biosynthesis.

Bicc1 ribonucleoprotein complexes specifying organ laterality are licensed by ANKS6-induced structural remodeling of associated ANKS3.

Organ laterality of vertebrates is specified by accelerated asymmetric decay of Dand5 mRNA mediated by Bicaudal-C1 (Bicc1) on the left side, but whether binding of this or any other mRNA to Bicc1 can be regulated is unknown. Here, we found that a CRISPR-engineered truncation in ankyrin and sterile alpha motif (SAM)-containing 3 (ANKS3) leads to symmetric mRNA decay mediated by the Bicc1-interacting Dand5 3' UTR. AlphaFold structure predictions of protein complexes and their biochemical validation by in vitro reconstitution reveal a novel interaction of the C-terminal coiled coil domain of ANKS3 with Bicc1 that inhibits binding of target mRNAs, depending on the conformation of ANKS3 and its regulation by ANKS6. The dual regulation of RNA binding by mutually opposing structured protein domains in this multivalent protein network emerges as a novel mechanism linking associated laterality defects and possibly other ciliopathies to perturbed dynamics in Bicc1 ribonucleoparticle (RNP) formation.

Nucleo-cytosolic Shuttling of ARGONAUTE1 Prompts a Revised Model of the Plant MicroRNA Pathway.

Unlike in metazoans, plant microRNAs (miRNAs) undergo stepwise nuclear maturation before engaging cytosolic, sequence-complementary transcripts in association with the silencing effector protein ARGONAUTE1 (AGO1). Since their discovery, how and under which form plant miRNAs translocate to the cytosol has remained unclear, as has their sub-cellular AGO1 loading site(s). Here, we show that the N termini of all plant AGO1s contain a nuclear-localization (NLS) and nuclear-export signal (NES) that, in Arabidopsis thaliana (At), enables AtAGO1 nucleo-cytosolic shuttling in a Leptomycin-B-inhibited manner, diagnostic of CRM1(EXPO1)/NES-dependent nuclear export. Nuclear-only AtAGO1 contains the same 2'O-methylated miRNA cohorts as its nucleo-cytosolic counterpart, but it preferentially interacts with the miRNA loading chaperone HSP90. Furthermore, mature miRNA translocation and miRNA-mediated silencing both require AtAGO1 nucleo-cytosolic shuttling. These findings lead us to propose a substantially revised view of the plant miRNA pathway in which miRNAs are matured, methylated, loaded into AGO1 in the nucleus, and exported to the cytosol as AGO1:miRNA complexes in a CRM1(EXPO1)/NES-dependent manner.

Cerebellar Ataxia and Coenzyme Q Deficiency through Loss of Unorthodox Kinase Activity.

The UbiB protein kinase-like (PKL) family is widespread, comprising one-quarter of microbial PKLs and five human homologs, yet its biochemical activities remain obscure. COQ8A (ADCK3) is a mammalian UbiB protein associated with ubiquinone (CoQ) biosynthesis and an ataxia (ARCA2) through unclear means. We show that mice lacking COQ8A develop a slowly progressive cerebellar ataxia linked to Purkinje cell dysfunction and mild exercise intolerance, recapitulating ARCA2. Interspecies biochemical analyses show that COQ8A and yeast Coq8p specifically stabilize a CoQ biosynthesis complex through unorthodox PKL functions. Although COQ8 was predicted to be a protein kinase, we demonstrate that it lacks canonical protein kinase activity in trans. Instead, COQ8 has ATPase activity and interacts with lipid CoQ intermediates, functions that are likely conserved across all domains of life. Collectively, our results lend insight into the molecular activities of the ancient UbiB family and elucidate the biochemical underpinnings of a human disease.

State-of-the-art web services for de novo protein structure prediction.

Residue coevolution estimations coupled to machine learning methods are revolutionizing the ability of protein structure prediction approaches to model proteins that lack clear homologous templates in the Protein Data Bank (PDB). This has been patent in the last round of the Critical Assessment of Structure Prediction (CASP), which presented several very good models for the hardest targets. Unfortunately, literature reporting on these advances often lacks digests tailored to lay end users; moreover, some of the top-ranking predictors do not provide webservers that can be used by nonexperts. How can then end users benefit from these advances and correctly interpret the predicted models? Here we review the web resources that biologists can use today to take advantage of these state-of-the-art methods in their research, including not only the best de novo modeling servers but also datasets of models precomputed by experts for structurally uncharacterized protein families. We highlight their features, advantages and pitfalls for predicting structures of proteins without clear templates. We present a broad number of applications that span from driving forward biochemical investigations that lack experimental structures to actually assisting experimental structure determination in X-ray diffraction, cryo-EM and other forms of integrative modeling. We also discuss issues that must be considered by users yet still require further developments, such as global and residue-wise model quality estimates and sources of residue coevolution other than monomeric tertiary structure.

Palmitoylated acyl protein thioesterase APT2 deforms membranes to extract substrate acyl chains.

Many biochemical reactions require controlled recruitment of proteins to membranes. This is largely regulated by posttranslational modifications. A frequent one is S-acylation, which consists of the addition of acyl chains and can be reversed by poorly understood acyl protein thioesterases (APTs). Using a panel of computational and experimental approaches, we dissect the mode of action of the major cellular thioesterase APT2 (LYPLA2). We show that soluble APT2 is vulnerable to proteasomal degradation, from which membrane binding protects it. Interaction with membranes requires three consecutive steps: electrostatic attraction, insertion of a hydrophobic loop and S-acylation by the palmitoyltransferases ZDHHC3 or ZDHHC7. Once bound, APT2 is predicted to deform the lipid bilayer to extract the acyl chain bound to its substrate and capture it in a hydrophobic pocket to allow hydrolysis. This molecular understanding of APT2 paves the way to understand the dynamics of APT2-mediated deacylation of substrates throughout the endomembrane system.

CLoNe: automated clustering based on local density neighborhoods for application to biomolecular structural ensembles.

MOTIVATION: Proteins are intrinsically dynamic entities. Flexibility sampling methods, such as molecular dynamics or those arising from integrative modeling strategies, are now commonplace and enable the study of molecular conformational landscapes in many contexts. Resulting structural ensembles increase in size as technological and algorithmic advancements take place, making their analysis increasingly demanding. In this regard, cluster analysis remains a go-to approach for their classification. However, many state-of-the-art algorithms are restricted to specific cluster properties. Combined with tedious parameter fine-tuning, cluster analysis of protein structural ensembles suffers from the lack of a generally applicable and easy to use clustering scheme. RESULTS: We present CLoNe, an original Python-based clustering scheme that builds on the Density Peaks algorithm of Rodriguez and Laio. CLoNe relies on a probabilistic analysis of local density distributions derived from nearest neighbors to find relevant clusters regardless of cluster shape, size, distribution and amount. We show its capabilities on many toy datasets with properties otherwise dividing state-of-the-art approaches and improves on the original algorithm in key aspects. Applied to structural ensembles, CLoNe was able to extract meaningful conformations from membrane binding events and ligand-binding pocket opening as well as identify dominant dimerization motifs or inter-domain organization. CLoNe additionally saves clusters as individual trajectories for further analysis and provides scripts for automated use with molecular visualization software. AVAILABILITY AND IMPLEMENTATION: www.epfl.ch/labs/lbm/resources, github.com/LBM-EPFL/CLoNe. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Visualization, Interactive Handling and Simulation of Molecules in Commodity Augmented Reality in Web Browsers Using moleculARweb's Virtual Modeling Kits.

moleculARweb (https://molecularweb.epfl.ch) began as a website for education and outreach in chemistry and structural biology through augmented reality (AR) content that runs in the web browsers of regular devices like smartphones, tablets, and computers. Here we present two evolutions of moleculARweb's Virtual Modeling Kits (VMK), tools where users can build and view molecules, and explore their mechanics, in 3D AR by handling the molecules in full 3D with custom-printed cube markers (VMK 2.0) or by moving around a simulated scene with mouse or touch gestures (VMK 3.0). Upon simulation the molecules experience visually realistic torsions, clashes, and hydrogen-bonding interactions that the user can manually switch on and off to explore their effects. Moreover, by manually tuning a fictitious temperature the users can accelerate conformational transitions or 'freeze' specific conformations for careful inspection in 3D. Even some phase transitions and separations can be simulated. We here showcase these and other features of the new VMKs connecting them to possible specific applications to teaching and self-learning of concepts from general, organic, biological and physical chemistry; and in assisting with small tasks in molecular modelling for research. Last, in a short discussion section we overview what future developments are needed for the 'dream tool' for the future of chemistry education and work.

Specific Protein-Membrane Interactions Promote Packaging of Metallo-ß-Lactamases into Outer Membrane Vesicles.

Outer membrane vesicles (OMVs) act as carriers of bacterial products such as plasmids and resistance determinants, including metallo-ß-lactamases. The lipidated, membrane-anchored metallo-ß-lactamase NDM-1 can be detected in Gram-negative OMVs. The soluble domain of NDM-1 also forms electrostatic interactions with the membrane. Here, we show that these interactions promote its packaging into OMVs produced by Escherichia coli. We report that favorable electrostatic protein-membrane interactions are also at work in the soluble enzyme IMP-1 while being absent in VIM-2. These interactions correlate with an enhanced incorporation of IMP-1 compared to VIM-2 into OMVs. Disruption of these interactions in NDM-1 and IMP-1 impairs their inclusion into vesicles, confirming their role in defining the protein cargo in OMVs. These results also indicate that packaging of metallo-ß-lactamases into vesicles in their active form is a common phenomenon that involves cargo selection based on specific molecular interactions.

Salt Enhances the Thermostability of Enteroviruses by Stabilizing Capsid Protein Interfaces.

Enteroviruses are common agents of infectious disease that are spread by the fecal-oral route. They are readily inactivated by mild heat, which causes the viral capsid to disintegrate or undergo conformational change. While beneficial for the thermal treatment of food or water, this heat sensitivity poses challenges for the stability of enterovirus vaccines. The thermostability of an enterovirus can be modulated by the composition of the suspending matrix, though the effects of the matrix on virus stability are not understood. Here, we determined the thermostability of four enterovirus strains in solutions with various concentrations of NaCl and different pH values. The experimental findings were combined with molecular modeling of the protein interaction forces at the pentamer and the protomer interfaces of the viral capsids. While pH only had a modest effect on thermostability, increasing NaCl concentrations raised the breakpoint temperatures of all viruses tested by up to 20°C. This breakpoint shift could be explained by an enhancement of the van der Waals attraction forces at the two protein interfaces. In comparison, the (net repulsive) electrostatic interactions were less affected by NaCl. Depending on the interface considered, the breakpoint temperature shifted by 7.5 or 5.6°C per 100-kcal/(mol·Å) increase in protein interaction force.IMPORTANCE The genus Enterovirus encompasses important contaminants of water and food (e.g., coxsackieviruses), as well as viruses of acute public health concern (e.g., poliovirus). Depending on the properties of the surrounding matrix, enteroviruses exhibit different sensitivities to heat, which in turn influences their persistence in the environment, during food treatment, and during vaccine storage. Here, we determined the effect of NaCl and pH on the heat stability of different enteroviruses and related the observed effects to changes in protein interaction forces in the viral capsid. We demonstrate that NaCl renders enteroviruses thermotolerant and that this effect stems from an increase in van der Waals forces at different protein subunits in the viral capsid. This work sheds light on the mechanism by which salt enhances virus stability.

Active Site-Induced Evolutionary Constraints Follow Fold Polarity Principles in Soluble Globular Enzymes.

A recent analysis of evolutionary rates in >500 globular soluble enzymes revealed pervasive conservation gradients toward catalytic residues. By looking at amino acid preference profiles rather than evolutionary rates in the same data set, we quantified the effects of active sites on site-specific constraints for physicochemical traits. We found that conservation gradients respond to constraints for polarity, hydrophobicity, flexibility, rigidity and structure in ways consistent with fold polarity principles; while sites far from active sites seem to experience no physicochemical constraint, rather being highly variable and favoring amino acids of low metabolic cost. Globally, our results highlight that amino acid variation contains finer information about protein structure than usually regarded in evolutionary models, and that this information is retrievable automatically with simple fits. We propose that analyses of the kind presented here incorporated into models of protein evolution should allow for better description of the physical chemistry that underlies molecular evolution.

Transcriptional regulatory logic of the diurnal cycle in the mouse liver.

Many organisms exhibit temporal rhythms in gene expression that propel diurnal cycles in physiology. In the liver of mammals, these rhythms are controlled by transcription-translation feedback loops of the core circadian clock and by feeding-fasting cycles. To better understand the regulatory interplay between the circadian clock and feeding rhythms, we mapped DNase I hypersensitive sites (DHSs) in the mouse liver during a diurnal cycle. The intensity of DNase I cleavages cycled at a substantial fraction of all DHSs, suggesting that DHSs harbor regulatory elements that control rhythmic transcription. Using chromatin immunoprecipitation followed by DNA sequencing (ChIP-seq), we found that hypersensitivity cycled in phase with RNA polymerase II (Pol II) loading and H3K27ac histone marks. We then combined the DHSs with temporal Pol II profiles in wild-type (WT) and Bmal1-/- livers to computationally identify transcription factors through which the core clock and feeding-fasting cycles control diurnal rhythms in transcription. While a similar number of mRNAs accumulated rhythmically in Bmal1-/- compared to WT livers, the amplitudes in Bmal1-/- were generally lower. The residual rhythms in Bmal1-/- reflected transcriptional regulators mediating feeding-fasting responses as well as responses to rhythmic systemic signals. Finally, the analysis of DNase I cuts at nucleotide resolution showed dynamically changing footprints consistent with dynamic binding of CLOCK:BMAL1 complexes. Structural modeling suggested that these footprints are driven by a transient heterotetramer binding configuration at peak activity. Together, our temporal DNase I mappings allowed us to decipher the global regulation of diurnal transcription rhythms in the mouse liver.

Will Cryo-Electron Microscopy Shift the Current Paradigm in Protein Structure Prediction?

Protein dynamics is undoubtedly a pervasive ingredient in all biological functions. However, structural biology has been strongly driven by a static-centered view of protein architecture. We argue that the recent advances of cryo-electron microscopy (EM) have the potential to more broadly explore the conformational landscapes of protein complexes and therefore will enhance our ability to predict the diverse conformations of tertiary and quaternary protein structures that are functionally relevant in physiological conditions.

A further leap of improvement in tertiary structure prediction in CASP13 prompts new routes for future assessments.

We present our assessment of tertiary structure predictions for hard targets in Critical Assessment of Structure Prediction round 13 (CASP13). The analysis includes (a) assignment and discussion of best models through scores-aided visual inspection of models for each evaluation unit (EU); (b) ranking of predictors resulting from this evaluation and from global scores; and (c) evaluation of progress, state of the art, and current limitations of protein structure prediction. We witness a sizable improvement in tertiary structure prediction building on the progress observed from CASP11 to CASP12, with (a) top models reaching backbone RMSD <3 å for several EUs of size <150 residues, contributed by many groups; (b) at least one model that roughly captures global topology for all EUs, probably unprecedented in this track of CASP; and (c) even quite good models for full, unsplit targets. Better structure predictions are brought about mainly by improved residue-residue contact predictions, and since this CASP also by distance predictions, achieved through state-of-the-art machine learning methods which also progressed to work with slightly shallower alignments compared to CASP12. As we reach a new realm of tertiary structure prediction quality, new directions are proposed and explored for future CASPs: (a) dropping splitting into EUs, (b) rethinking difficulty metrics probably in terms of contact and distance predictions, (c) assessing also side chains for models of high backbone accuracy, and (d) assessing residue-wise and possibly residue-residue quality estimates.

Immobilization of the N-terminal helix stabilizes prefusion paramyxovirus fusion proteins.

Parainfluenza virus 5 (PIV5) is an enveloped, single-stranded, negative-sense RNA virus of the Paramyxoviridae family. PIV5 fusion and entry are mediated by the coordinated action of the receptor-binding protein, hemagglutinin-neuraminidase (HN), and the fusion protein (F). Upon triggering by HN, F undergoes an irreversible ATP- and pH-independent conformational change, going down an energy gradient from a metastable prefusion state to a highly stable postfusion state. Previous studies have highlighted key conformational changes in the F-protein refolding pathway, but a detailed understanding of prefusion F-protein metastability remains elusive. Here, using two previously described F-protein mutations (S443D or P22L), we examine the capacity to modulate PIV5 F stability and the mechanisms by which these point mutants act. The S443D mutation destabilizes prefusion F proteins by disrupting a hydrogen bond network at the base of the F-protein globular head. The introduction of a P22L mutation robustly rescues destabilized F proteins through a local hydrophobic interaction between the N-terminal helix and a hydrophobic pocket. Prefusion stabilization conferred by a P22L-homologous mutation is demonstrated in the F protein of Newcastle disease virus, a paramyxovirus of a different genus, suggesting a conserved stabilizing structural element within the paramyxovirus family. Taken together, the available data suggest that movement of the N-terminal helix is a necessary early step for paramyxovirus F-protein refolding and presents a novel target for structure-based drug design.

Assessment of data-assisted prediction by inclusion of crosslinking/mass-spectrometry and small angle X-ray scattering data in the 12th Critical Assessment of protein Structure Prediction experiment.

Integrative modeling approaches attempt to combine experiments and computation to derive structure-function relationships in complex molecular assemblies. Despite their importance for the advancement of life sciences, benchmarking of existing methodologies is rather poor. The 12th round of the Critical Assessment of protein Structure Prediction (CASP) offered a unique niche to benchmark data and methods from two kinds of experiments often used in integrative modeling, namely residue-residue contacts obtained through crosslinking/mass-spectrometry (CLMS), and small-angle X-ray scattering (SAXS) experiments. Upon assessment of the models submitted by predictors for 3 targets assisted by CLMS data and 11 targets by SAXS data, we observed no significant improvement when compared to the best data-blind models, although most predictors did improve relative to their own data-blind predictions. Only for target Tx892 of the CLMS-assisted category and for target Ts947 of the SAXS-assisted category, there was a net, albeit mild, improvement relative to the best data-blind predictions. We discuss here possible reasons for the relatively poor success, which point rather to inconsistencies in the data sources rather than in the methods, to which a few groups were less sensitive. We conclude with suggestions that could improve the potential of data integration in future CASP rounds in terms of experimental data production, methods development, data management and prediction assessment.

Assessment of hard target modeling in CASP12 reveals an emerging role of alignment-based contact prediction methods.

We present our assessment of CASP12 modeling efforts for targets with no obvious templates of high sequence/structure similarity in the PDB, that is for evaluation units of the free modeling (FM) and free modeling/template-based modeling (FM/TBM) categories. Models were clustered and ranked using the Global Distance Test-Total Score and 5 additional metrics developed in previous CASP rounds, producing short lists of models that were subject to visual inspection in comparison to the target structures. The whole procedure was implemented as a web app that facilitates model selection and visual inspection, and could become useful to facilitate and standardize future assessments. We describe cases of (1) targets with remarkably good predictions, (2) targets whose models captured some global shape and topology features, and (3) targets for which models fail to capture even coarse features. We note that despite this CASP being among the most challenging ones, a measurable improvement of the top predictions is apparent, that we attribute to the emergence of accurate contact prediction methods and the increased number of available sequences. We also briefly discuss current limitations in tertiary structure prediction exemplified by CASP12 targets. Overall, the Baker, Zhang, and Lee manual groups and servers were identified as the top global performing groups.

Definition and classification of evaluation units for tertiary structure prediction in CASP12 facilitated through semi-automated metrics.

For assessment purposes, CASP targets are split into evaluation units. We herein present the official definition of CASP12 evaluation units (EUs) and their classification into difficulty categories. Each target can be evaluated as one EU (the whole target) or/and several EUs (separate structural domains or groups of structural domains). The specific scenario for a target split is determined by the domain organization of available templates, the difference in server performance on separate domains versus combination of the domains, and visual inspection. In the end, 71 targets were split into 96 EUs. Classification of the EUs into difficulty categories was done semi-automatically with the assistance of metrics provided by the Prediction Center. These metrics account for sequence and structural similarities of the EUs to potential structural templates from the Protein Data Bank, and for the baseline performance of automated server predictions. The metrics readily separate the 96 EUs into 38 EUs that should be straightforward for template-based modeling (TBM) and 39 that are expected to be hard for homology modeling and are thus left for free modeling (FM). The remaining 19 borderline evaluation units were dubbed FM/TBM, and were inspected case by case. The article also overviews structural and evolutionary features of selected targets relevant to our accompanying article presenting the assessment of FM and FM/TBM predictions, and overviews structural features of the hardest evaluation units from the FM category. We finally suggest improvements for the EU definition and classification procedures.