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.

Ligand Binding to the Collagen VI Receptor Triggers a Talin-to-RhoA Switch that Regulates Receptor Endocytosis

Capillary morphogenesis gene 2 (CMG2/ANTXR2) is a cell surface receptor for both collagen VI and anthrax toxin. Biallelic loss-of-function mutations in CMG2 lead to a severe condition, hyaline fibromatosis syndrome (HFS). We have here dissected a network of dynamic interactions between CMG2 and various actin interactors and regulators, describing a different behavior from other extracellular matrix receptors. CMG2 binds talin, and thereby the actin cytoskeleton, only in its ligand-free state. Extracellular ligand binding leads to src-dependent talin release and recruitment of the actin cytoskeleton regulator RhoA and its effectors. These sequential interactions of CMG2 are necessary for the control of oriented cell division during fish development. Finally, we demonstrate that effective switching between talin and RhoA binding is required for the intracellular degradation of collagen VI in human fibroblasts, which explains why HFS mutations in the cytoskeleton-binding domain lead to dysregulation of extracellular matrix homeostasis.

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.

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.

An experiment-informed signal transduction model for the role of the Staphylococcus aureus MecR1 protein in β-lactam resistance

The treatment of hospital- and community-associated infections by methicillin-resistant Staphylococcus aureus (MRSA) is a perpetual challenge. This Gram-positive bacterium is resistant specifically to β-lactam antibiotics, and generally to many other antibacterial agents. Its resistance mechanisms to β-lactam antibiotics are activated only when the bacterium encounters a β-lactam. This activation is regulated by the transmembrane sensor/signal transducer proteins BlaR1 and MecR1. Neither the transmembrane/metalloprotease domain, nor the complete MecR1 and BlaR1 proteins, are isolatable for mechanistic study. Here we propose a model for full-length MecR1 based on homology modeling, residue coevolution data, a new extensive experimental mapping of transmembrane topology, partial structures, molecular simulations, and available NMR data. Our model defines the metalloprotease domain as a hydrophilic transmembrane chamber effectively sealed by the apo-sensor domain. It proposes that the amphipathic helices inserted into the gluzincin domain constitute the route for transmission of the β-lactam-binding event in the extracellular sensor domain, to the intracellular and membrane-embedded zinc-containing active site. From here, we discuss possible routes for subsequent activation of proteolytic action. This study provides the first coherent model of the structure of MecR1, opening routes for future functional investigations on how β-lactam binding culminates in the proteolytic degradation of MecI.

Modelling structures in cryo-EM maps

Recent advances in structure determination of sub-cellular structures using cryo-electron microscopy and tomography have enabled us to understand their architecture in a more detailed manner and gain insight into their function. The choice of approach to use for atomic model building, fitting, refinement and validation in the 3D map resulting from these experiments depends primarily on the resolution of the map and the prior information on the corresponding model. Here, we survey some of such methods and approaches and highlight their uses in specific recent examples.

KAP1 is an antiparallel dimer with a functional asymmetry

This study reveals the architecture of human KAP1 by integrating molecular modeling with small-angle X-ray scattering and single-molecule experiments. KAP1 dimers feature a structural asymmetry at the C-terminal domains that has functional implications for recruitment of HP1.

KAP1 (KRAB domain–associated protein 1) plays a fundamental role in regulating gene expression in mammalian cells by recruiting different transcription factors and altering the chromatin state. In doing so, KAP1 acts both as a platform for macromolecular interactions and as an E3 small ubiquitin modifier ligase. This work sheds light on the overall organization of the full-length protein combining solution scattering data, integrative modeling, and single-molecule experiments. We show that KAP1 is an elongated antiparallel dimer with an asymmetry at the C-terminal domains. This conformation is consistent with the finding that the Really Interesting New Gene (RING) domain contributes to KAP1 auto-SUMOylation. Importantly, this intrinsic asymmetry has key functional implications for the KAP1 network of interactions, as the heterochromatin protein 1 (HP1) occupies only one of the two putative HP1 binding sites on the KAP1 dimer, resulting in an unexpected stoichiometry, even in the context of chromatin fibers.

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.

Abstract P250: Unveiling Binding Pocket Structure Of Mas Receptor And Its Interaction With Angiotensin-(1-7)

Substantial amount of preclinical data points out Mas receptor as a promising pharmacological target to treat cardiovascular diseases. Here, we aimed to characterize Mas binding pocket and its interaction with Angiotensin-(Ang)-1-7, in order to provide insights for potential structure-based drug discovery. Homology modeling based on AT1 and AT2 crystal structures, combined with peptide docking (Schrodinger TM ) and molecular dynamics simulations (Gromacs) were employed to generate hypothesis about crucial interactions in protein-ligand complex. Proposed model was further validated in vitro using HEK cells transiently transfected with wild-type Mas or mutant Mas (Arg245Ala). Mas-receptor activation was assessed by FOXO1 dephosphorylation at Ser256 after 5 minutes incubation with Ang-1-7 (10 -7 M). Model of Mas/Ang-(1-7) complex revealed electrostatic interaction between Ang-(1-7) C-terminus and Mas Arg245. Interestingly, while Ang-1-7 incubation reduced FOXO1 phosphorylation by 35 % in HEK-Mas wild type (1.000 vs. 0.654 A.U., p<0.01, n=8), this effect was abolished in HEK-Mas mutant (1.000 vs. 0.990 A.U., n.s. , n=5). In sum, Arg245 of Mas is essential for the receptor activation, which is in accordance with A-779 antagonism (C-terminus modification of Ang-(1-7)), endorsing our proposed structure model. Further elucidation of the ligand-receptor interactions may drive rational drug discovery and lead to development of improved strategies for Mas targeting.

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.

Functional assays for the assessment of the pathogenicity of variants of GOSR2, an ER-to-Golgi SNARE involved in progressive myoclonus epilepsies

Progressive myoclonus epilepsies (PMEs) are inherited disorders characterized by myoclonus, generalized tonic-clonic seizures, and ataxia. One of the genes that is associated with PME is the ER-to-Golgi Q_b-SNARE GOSR2, which forms a SNARE complex with syntaxin-5, Bet1 and Sec22b. Most PME patients are homo­zygous for a p.Gly144Trp mutation and develop similar clinical presentations. Recently, a patient who was compound heterozygous for p.Gly144Trp and a previously unseen p.Lys164del mutation was identified. Because this patient presented with a milder disease phenotype, we hypothesized that the p.Lys164del mutation may be less severe compared to p.Gly144Trp. To characterize the effect of the p.Gly144Trp and p.Lys164del mutations, both of which are present in the SNARE motif of GOSR2, we examined the corresponding mutations in the yeast ortholog Bos1. Yeasts expressing the orthologous mutants in Bos1 showed impaired growth, suggesting a partial loss of function, which was more severe for the Bos1 p.Gly176Trp mutation. Using anisotropy and gel filtration, we report that Bos1 p.Gly176Trp and p.Arg196del are capable of complex formation, but with partly reduced activity. Molecular dynamics (MD) simulations showed that the hydrophobic core, which triggers SNARE complex formation, is compromised due to the glycine-to-tryptophan substitution in both GOSR2 and Bos1. In contrast, the deletion of residue p.Lys164 (or p.Arg196del in Bos1) interferes with the formation of hydrogen bonds between GOSR2 and syntaxin-5. Despite these perturbations, all SNARE complexes stayed intact during longer simulations. Thus, our data suggest that the milder course of disease in compound heterozygous PME is due to less severe impairment of the SNARE function.

Summary: Mutations in the Q_b-SNARE GOSR2 cause progressive myoclonus epilepsies. The authors report the effect of two mutations on SNARE function to investigate their correlation with progression and severity of disease.

Conserved Lipid and Small-Molecule Modulation of COQ8 Reveals Regulation of the Ancient Kinase-like UbiB Family

Human COQ8A (ADCK3) and Saccharomyces cerevisiae Coq8p (collectively COQ8) are UbiB family proteins essential for mitochondrial coenzyme Q (CoQ) biosynthesis. However, the biochemical activity of COQ8 and its direct role in CoQ production remain unclear, in part due to lack of known endogenous regulators of COQ8 function and of effective small molecules for probing its activity in vivo. Here, we demonstrate that COQ8 possesses evolutionarily conserved ATPase activity that is activated by binding to membranes containing cardiolipin and by phenolic compounds that resemble CoQ pathway intermediates. We further create an analog-sensitive version of Coq8p and reveal that acute chemical inhibition of its endogenous activity in yeast is sufficient to cause respiratory deficiency concomitant with CoQ depletion. Collectively, this work defines lipid and small-molecule modulators of an ancient family of atypical kinase-like proteins and establishes a chemical genetic system for further exploring the mechanistic role of COQ8 in CoQ biosynthesis.

Recruitment of the amyloid precursor protein by γ-secretase at the synaptic plasma membrane

Γ-secretase is a membrane-embedded protease that cleaves single transmembrane helical domains of various integral membrane proteins. The amyloid precursor protein (APP) is an important substrate due to its pathological relevance to Alzheimer's disease. The mechanism of the cleavage of APP by γ-secretase that leads to accumulation of Alzheimer's disease causing amyloid-β (Aβ) is still unknown. Coarse-grained molecular dynamics simulations in this study reveal initial lipids raft formation near the catalytic site of γ-secretase as well as changes in dynamic behavior of γ-secretase once interacting with APP. The results suggest a precursor of the APP binding mode and hint at conformational changes of γ-secretase in the nicastrin (NCT) domain upon APP binding.

Structural, physicochemical and dynamic features conserved within the aerolysin pore-forming toxin family

Aerolysin is the founding member of a major class of β-pore-forming toxins (β-PFTs) found throughout all kingdoms of life. PFTs are cytotoxic proteins produced as soluble monomers, which oligomerize at the membrane of target host cells forming pores that may lead to osmotic lysis and cell death. Besides their role in microbial infection, they have become interesting for their potential as biotechnological sensors and delivery systems. Using an approach that integrates bioinformatics with molecular modeling and simulation, we looked for conserved features across this large toxin family. The cell surface-binding domains present high variability within the family to provide membrane receptor specificity. On the contrary, the novel concentric double β-barrel structure found in aerolysin is highly conserved in terms of sequence, structure and conformational dynamics, which likely contribute to preserve a common transition mechanism from the prepore to the mature pore within the family.Our results point to the key role of several amino acids in the conformational changes needed for oligomerization and further pore formation, such as Y221, W227, P248, Q263 and L277, which we propose are involved in the release of the stem loop and the two adjacent β-strands to form the transmembrane β-barrel.

Molecular Effects of Concentrated Solutes on Protein Hydration, Dynamics, and Electrostatics

Most studies of protein structure and function are performed in dilute conditions, but proteins typically experience high solute concentrations in their physiological scenarios and biotechnological applications. High solute concentrations have well-known effects on coarse protein traits like stability, diffusion, and shape, but likely also perturb other traits through finer effects pertinent at the residue and atomic levels. Here, NMR and molecular dynamics investigations on ubiquitin disclose variable interactions with concentrated solutes that lead to localized perturbations of the protein's surface, hydration, electrostatics, and dynamics, all dependent on solute size and chemical properties. Most strikingly, small polar uncharged molecules are sticky on the protein surface, whereas charged small molecules are not, but the latter still perturb the internal protein electrostatics as they diffuse nearby. Meanwhile, interactions with macromolecular crowders are favored mainly through hydrophobic, but not through polar, surface patches. All the tested small solutes strongly slow down water exchange at the protein surface, whereas macromolecular crowders do not exert such strong perturbation. Finally, molecular dynamics simulations predict that unspecific interactions slow down microsecond- to millisecond-timescale protein dynamics despite having only mild effects on pico- to nanosecond fluctuations as corroborated by NMR. We discuss our results in the light of recent advances in understanding proteins inside living cells, focusing on the physical chemistry of quinary structure and cellular organization, and we reinforce the idea that proteins should be studied in native-like media to achieve a faithful description of their function.

Signal Sensing and Transduction by Histidine Kinases as Unveiled through Studies on a Temperature Sensor

Histidine kinases (HK) are the sensory proteins of two-component systems, responsible for a large fraction of bacterial responses to stimuli and environmental changes. Prototypical HKs are membrane-bound proteins that phosphorylate cognate response regulator proteins in the cytoplasm upon signal detection in the membrane or periplasm. HKs stand as potential drug targets but also constitute fascinating systems for studying proteins at work, specifically regarding the chemistry and mechanics of signal detection, transduction through the membrane, and regulation of catalytic outputs. In this Account, we focus on Bacillus subtilis DesK, a membrane-bound HK part of a two-component system that maintains appropriate membrane fluidity at low growth temperatures. Unlike most HKs, DesK has no extracytoplasmic signal-sensing domains; instead, sensing is carried out by 10 transmembrane helices (coming from two protomers) arranged in an unknown structure. The fifth transmembrane helix from each protomer connects, without any of the intermediate domains found in other HKs, into the dimerization and histidine phosphotransfer (DHp) domain located in the cytoplasm, which is followed by the ATP-binding domains (ABD). Throughout the years, genetic, biochemical, structural, and computational studies on wild-type, mutant, and truncated versions of DesK allowed us to dissect several aspects of DesK's functioning, pushing forward a more general understanding of its own structure/function relationships as well as those of other HKs. We have shown that the sensing mechanism is rooted in temperature-dependent membrane properties, most likely a combination of thickness, fluidity, and water permeability, and we have proposed possible mechanisms by which DesK senses these properties and transduces the signals. X-ray structures and computational models have revealed structural features of TM and cytoplasmic regions in DesK's kinase- and phosphatase-competent states. Biochemical and genetic experiments and molecular simulations further showed that reversible formation of a two-helix coiled coil in the fifth TM segment and the N-terminus of the cytoplasmic domain is essential for the sensing and signal transduction mechanisms. Together with other structural and functional works, the emerging picture suggests that diverse HKs possess distinct sensing and transduction mechanisms but share as rather general features (i) a symmetric phosphatase state and an asymmetric kinase state and (ii) similar functional outputs on the conserved DHp and ABD domains, achieved through different mechanisms that depend on the nature of the initial signal. We here advance (iii) an important role for TM prolines in transducing the initial signals to the cytoplasmic coiled coils, based on simulations of DesK's TM helices and our previous work on a related HK, PhoQ. Lastly, evidence for DesK, PhoQ, BvgS, and DctB HKs shows that (iv) overall catalytic output is tuned by a delicate balance between hydration potentials, coiled coil stability, and exposure of hydrophobic surface patches at their cytoplasmic coiled coils and at the N-terminal and C-terminal sides of their TM helices. This balance is so delicate that small perturbations, either physiological signals or induced by mutations, lead to large remodeling of the underlying conformational landscape achieving clear-cut changes in catalytic output, mirroring the required response speed of these systems for proper biological function.

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.

Protein post-translational modifications: In silico prediction tools and molecular modeling

Post-translational modifications (PTMs) occur in almost all proteins and play an important role in numerous biological processes by significantly affecting proteins' structure and dynamics. Several computational approaches have been developed to study PTMs (e.g., phosphorylation, sumoylation or palmitoylation) showing the importance of these techniques in predicting modified sites that can be further investigated with experimental approaches. In this review, we summarize some of the available online platforms and their contribution in the study of PTMs. Moreover, we discuss the emerging capabilities of molecular modeling and simulation that are able to complement these bioinformatics methods, providing deeper molecular insights into the biological function of post-translational modified proteins.

Disentangling constraints using viability evolution principles in integrative modeling of macromolecular assemblies

Predicting the structure of large molecular assemblies remains a challenging task in structural biology when using integrative modeling approaches. One of the main issues stems from the treatment of heterogeneous experimental data used to predict the architecture of native complexes. We propose a new method, applied here for the first time to a set of symmetrical complexes, based on evolutionary computation that treats every available experimental input independently, bypassing the need to balance weight components assigned to aggregated fitness functions during optimization.