Cooperative regulation of C1-domain membrane recruitment polarizes atypical protein kinase C

Recruitment of the Par complex protein atypical protein kinase C (aPKC) to a specific membrane domain is a key step in the polarization of animal cells. While numerous proteins and phospholipids interact with aPKC, how these interactions cooperate to control its membrane recruitment has been unknown. Here, we identify aPKC's C1 domain as a phospholipid interaction module that targets aPKC to the membrane of Drosophila neural stem cells (NSCs). The isolated C1 binds the NSC membrane in an unpolarized manner during interphase and mitosis and is uniquely sufficient among aPKC domains for targeting. Other domains, including the catalytic module and those that bind the upstream regulators Par-6 and Bazooka, restrict C1's membrane targeting activity-spatially and temporally-to the apical NSC membrane during mitosis. Our results suggest that aPKC polarity results from cooperative activation of autoinhibited C1-mediated membrane binding activity.

Theoretical Investigation of Repurposed Drugs Potentially Capable of Binding to the Catalytic Site and the Secondary Binding Pocket of Subunit A of Ricin

Recently, we reported a library of 82 compounds, selected from different databanks through virtual screening and docking studies, and pointed to 6 among them as potential repurposed dual binders to both the catalytic site and the secondary binding pockets of subunit A of ricin (RTA). Here, we report additional molecular modeling studies of an extended list of compounds from the original library. Rounds of flexible docking followed by molecular dynamics simulations and further rounds of MM-PBSA calculations using a more robust protocol, enabled a better investigation of the interactions of these compounds inside RTA, the elucidation of their dynamical behaviors, and updating the list of the most important residues for the ligand binding. Four compounds were pointed as potential repurposed ricin inhibitors that are worth being experimentally investigated.

Single-cell analysis of human basal cell carcinoma reveals novel regulators of tumor growth and the tumor microenvironment

How basal cell carcinoma (BCC) interacts with its tumor microenvironment to promote growth is unclear. We use singe-cell RNA sequencing to define the human BCC ecosystem and discriminate between normal and malignant epithelial cells. We identify spatial biomarkers of tumors and their surrounding stroma that reinforce the heterogeneity of each tissue type. Combining pseudotime, RNA velocity-PAGA, cellular entropy, and regulon analysis in stromal cells reveals a cancer-specific rewiring of fibroblasts, where STAT1, TGF-β, and inflammatory signals induce a noncanonical WNT5A program that maintains the stromal inflammatory state. Cell-cell communication modeling suggests that tumors respond to the sudden burst of fibroblast-specific inflammatory signaling pathways by producing heat shock proteins, whose expression we validated in situ. Last, dose-dependent treatment with an HSP70 inhibitor suppresses in vitro vismodegib-resistant BCC cell growth, Hedgehog signaling, and in vivo tumor growth in a BCC mouse model, validating HSP70's essential role in tumor growth and reinforcing the critical nature of tumor microenvironment cross-talk in BCC progression.

Improved Accuracy for Modeling PROTAC-Mediated Ternary Complex Formation and Targeted Protein Degradation via New In Silico Methodologies

Extending upon our previous publication [Drummond, M.; J. Chem. Inf. Model. 2019, 59, 1634-1644], two additional computational methods are presented to model PROTAC-mediated ternary complex structures, which are then used to predict the efficacy of any accompanying protein degradation. Method 4B, an extension to one of our previous approaches, incorporates a clustering procedure uniquely suited for considering ternary complexes. Method 4B yields the highest proportion to date of crystal-like poses in modeled ternary complex ensembles, nearing 100% in two cases and always giving a hit rate of at least 10%. Techniques to further improve this performance for particularly troublesome cases are suggested and validated. This demonstrated ability to reliably reproduce known crystallographic ternary complex structures is further established through modeling of a newly released crystal structure. Moreover, for the far more common scenario where the structure of the ternary complex intermediate is unknown, the methods detailed in this work nonetheless consistently yield results that reliably follow experimental protein degradation trends, as established through seven retrospective case studies. These various case studies cover challenging yet common modeling situations, such as when the precise orientation of the PROTAC binding moiety in one (or both) of the protein pockets has not been experimentally established. Successful results are presented for one PROTAC targeting many proteins, for different PROTACs targeting the same protein, and even for degradation effected by an E3 ligase that has not been structurally characterized in a ternary complex. Overall, the computational modeling approaches detailed in this work should greatly facilitate PROTAC screening and design efforts, so that the many advantages of a PROTAC-based degradation approach can be effectively utilized both rapidly and at reduced cost.

Single cell transcriptomics of human epidermis identifies basal stem cell transition states

How stem cells give rise to epidermis is unclear despite the crucial role the epidermis plays in barrier and appendage formation. Here we use single cell-RNA sequencing to interrogate basal stem cell heterogeneity of human interfollicular epidermis and find four spatially distinct stem cell populations at the top and bottom of rete ridges and transitional positions between the basal and suprabasal epidermal layers. Cell-cell communication modeling suggests that basal cell populations serve as crucial signaling hubs to maintain epidermal communication. Combining pseudotime, RNA velocity, and cellular entropy analyses point to a hierarchical differentiation lineage supporting multi-stem cell interfollicular epidermal homeostasis models and suggest that transitional basal stem cells are stable states essential for proper stratification. Finally, alterations in differentially expressed transitional basal stem cell genes result in severe thinning of human skin equivalents, validating their essential role in epidermal homeostasis and reinforcing the critical nature of basal stem cell heterogeneity.

Actin polymerization controls cilia-mediated signaling

Actin polymerization is important to generate primary cilia. Drummond et al. show that upstream actin regulators are necessary for this process by controlling aPKC and Src kinase activity to promote Hedgehog signaling and restrict primary cilia.

Primary cilia are polarized organelles that allow detection of extracellular signals such as Hedgehog (Hh). How the cytoskeleton supporting the cilium generates and maintains a structure that finely tunes cellular response remains unclear. Here, we find that regulation of actin polymerization controls primary cilia and Hh signaling. Disrupting actin polymerization, or knockdown of N-WASp/Arp3, increases ciliation frequency, axoneme length, and Hh signaling. Cdc42, a potent actin regulator, recruits both atypical protein pinase C iota/lambda (aPKC) and Missing-in-Metastasis (MIM) to the basal body to maintain actin polymerization and restrict axoneme length. Transcriptome analysis implicates the Src pathway as a major aPKC effector. aPKC promotes whereas MIM antagonizes Src activity to maintain proper levels of primary cilia, actin polymerization, and Hh signaling. Hh pathway activation requires Smoothened-, Gli-, and Gli1-specific activation by aPKC. Surprisingly, longer axonemes can amplify Hh signaling, except when aPKC is disrupted, reinforcing the importance of the Cdc42–aPKC–Gli axis in actin-dependent regulation of primary cilia signaling.

sm"FISH"ing for Hedgehog

Patched (Ptch) receptors are critical negative regulators of Hedgehog signaling, where Ptch1 loss causes basal cell carcinoma and Ptch1;Ptch2 loss disrupts skin and hair follicle development. Adolphe et al. use single molecule fluorescent in situ hybridization to show quantitatively that Ptch receptors create a Hedgehog signaling gradient that may specify hair follicle development.

Nature of protein-CO2 interactions as elucidated via molecular dynamics

Rising global temperatures require innovative measures to reduce atmospheric concentrations of CO(2). The most successful carbon capture technology on Earth is the enzymatic capture of CO(2) and its sequestration in the form of glucose. Efforts to improve upon or mimic this naturally occurring process will require a rich understanding of protein-CO(2) interactions. Toward that end, extensive all-atom molecular dynamics (MD) simulations were performed on the CO(2)-utilizing enzyme phosphoenolpyruvate carboxykinase (PEPCK). Preliminary simulations were performed using implicit and explicit solvent models, which yielded similar results: arginine, lysine, tyrosine, and asparagine enhance the ability of a protein to bind carbon dioxide. Extensive explicit solvent simulations were performed for both wild-type PEPCK and five single-point PEPCK mutants, revealing three prevalent channels by which CO(2) enters (or exits) the active site cleft, as well as a fourth channel (observed only once), the existence of which can be rationalized in terms of the position of a key Arg residue. The strongest CO(2) binding sites in these simulations consist of appropriately positioned hydrogen bond donors and acceptors. Interactions between CO(2) and both Mn(2+) and Mg(2+) present in PEPCK are minimal due to the stable protein- and solvent-based coordination environments of these cations. His 232, suggested by X-ray crystallography as being a potential important CO(2) binding site, is indeed found to be particularly "CO(2)-philic" in these simulations. Finally, a recent mechanism, proposed on the basis of X-ray crystallography, for PEPCK active site lid closure is discussed in light of the MD trajectories. Overall, the results of this work will prove useful not only to scientists investigating PEPCK, but also to groups seeking to develop an environmentally benign, protein-based carbon capture, sequestration, and utilization system.

Carbon Dioxide Migration Pathways in Proteins

Some of the most important biological processes, such as carbon fixation, are dependent on protein-gas interactions. The motion of CO2 through the enzyme phosphoenolpyruvate carboxykinase was investigated using extensive all-atom molecular dynamics simulations. Three discrete migration pathways were located, suggesting the protein directs the movement of CO2. The chemical nature of these pathways is discussed, as are their biotechnological ramifications.

Aspartate 458 of human glutathione synthetase is important for cooperativity and active site structure

Human glutathione synthetase (hGS) catalyzes the second ATP-dependent step in the biosynthesis of glutathione (GSH) and is negatively cooperative to the γ-glutamyl substrate. The hGS active site is composed of three highly conserved catalytic loops, notably the alanine rich A-loop. Experimental and computational investigations of the impact of mutation of Asp458 are reported, and thus the role of this A-loop residue on hGS structure, activity, negativity cooperativity and stability is defined. Several Asp458 hGS mutants (D458A, D458N and D458R) were constructed using site-directed mutagenesis and their activities determined (10%, 15% and 7% of wild-type hGS, respectively). The Michaelis-Menten constant (K(m)) was determined for all three substrates (glycine, GAB and ATP): glycine K(m) increased by 30-115-fold, GAB K(m) decreased by 8-17-fold, and the ATP K(m) was unchanged. All Asp458 mutants display a change in cooperativity from negative cooperativity to non-cooperative. All mutants show similar stability as compared to wild-type hGS, as determined by differential scanning calorimetry. The findings indicate that Asp458 is essential for hGS catalysis and that it impacts the allostery of hGS.

Identification of novel proteins in Neospora caninum using an organelle purification and monoclonal antibody approach

Neospora caninum is an important veterinary pathogen that causes abortion in cattle and neuromuscular disease in dogs. Neospora has also generated substantial interest because it is an extremely close relative of the human pathogen Toxoplasma gondii, yet does not appear to infect humans. While for Toxoplasma there are a wide array of molecular tools and reagents available for experimental investigation, relatively few reagents exist for Neospora. To investigate the unique biological features of this parasite and exploit the recent sequencing of its genome, we have used an organelle isolation and monoclonal antibody approach to identify novel organellar proteins and develop a wide array of probes for subcellular localization. We raised a panel of forty-six monoclonal antibodies that detect proteins from the rhoptries, micronemes, dense granules, inner membrane complex, apicoplast, mitochondrion and parasite surface. A subset of the proteins was identified by immunoprecipitation and mass spectrometry and reveal that we have identified and localized many of the key proteins involved in invasion and host interaction in Neospora. In addition, we identified novel secretory proteins not previously studied in any apicomplexan parasite. Thus, this organellar monoclonal antibody approach not only greatly enhances the tools available for Neospora cell biology, but also identifies novel components of the unique biological characteristics of this important veterinary pathogen.

Drosophila type II neuroblast lineages keep Prospero levels low to generate large clones that contribute to the adult brain central complex

Tissue homeostasis depends on the ability of stem cells to properly regulate self-renewal versus differentiation. Drosophila neural stem cells (neuroblasts) are a model system to study self-renewal and differentiation. Recent work has identified two types of larval neuroblasts that have different self-renewal/differentiation properties. Type I neuroblasts bud off a series of small basal daughter cells (ganglion mother cells) that each generate two neurons. Type II neuroblasts bud off small basal daughter cells called intermediate progenitors (INPs), with each INP generating 6 to 12 neurons. Type I neuroblasts and INPs have nuclear Asense and cytoplasmic Prospero, whereas type II neuroblasts lack both these transcription factors. Here we test whether Prospero distinguishes type I/II neuroblast identity or proliferation profile, using several newly characterized Gal4 lines. We misexpress prospero using the 19H09-Gal4 line (expressed in type II neuroblasts but no adjacent type I neuroblasts) or 9D11-Gal4 line (expressed in INPs but not type II neuroblasts). We find that differential prospero expression does not distinguish type I and type II neuroblast identities, but Prospero regulates proliferation in both type I and type II neuroblast lineages. In addition, we use 9D11 lineage tracing to show that type II lineages generate both small-field and large-field neurons within the adult central complex, a brain region required for locomotion, flight, and visual pattern memory.

CO2-formatics: how do proteins bind carbon dioxide?

The rising atmospheric concentration of CO(2) has motivated researchers to seek routes for improved utilization, increased mitigation, and enhanced sequestration of this greenhouse gas. Through a combination of bioinformatics, molecular modeling, and first-principles quantum mechanics the binding of carbon dioxide to proteins is analyzed. It is concluded that acid/base interactions are the principal chemical force by which CO(2) is bound inside proteins. With respect to regular secondary structural elements, beta-sheets show a marked preference for CO(2) binding compared to alpha-helices. The data also support the inference that while either or both oxygens of CO(2) are generally tightly bound in the protein environment, the carbon is much less "sequestered." First principles and more approximate modeling techniques are assessed for quantifying CO(2) binding thermodynamics.

Use of Drug Discovery Tools in Rational Organometallic Catalyst Design

A computational procedure is detailed where techniques common in the drug discovery process-2D- and 3D-quantitative structure-activity relationships (QSAR)-are applied to rationalize the catalytic activity of a synthetically flexible, Ti-N=P ethylene polymerization catalyst system. Once models relating molecular properties to catalyst activity are built with the two QSAR approaches, two database mining approaches are used to select a small number of ligands from a larger database that are likely to produce catalysts with high activity when grafted onto the Ti-N=P framework. The software employed throughout this work is freely available, is easy to use, and was applied in a "black box" approach to highlight areas where the drug discovery tools, designed to address organic molecules, have difficulty in addressing issues arising from the presence of a metal atom. In general, 3D-QSAR offers an efficient way to screen new potential ligands and separate those likely to lead to poor catalysts from those that are likely to contribute to highly active catalysts. The results for 2D-QSAR appear to be quantitatively unreliable, likely due to the presence of a metal atom; nonetheless, there is evidence that qualitative predictions from different models may be reliable. Pitfalls in the database mining techniques are identified, none of which are insurmountable. The lessons learned about the potential uses and drawbacks of the techniques described herein are readily applicable to other catalyst frameworks, thereby enabling a rational approach to catalyst improvement and design.

Structure and stability of small boron and boron oxide clusters

To rationally design and explore a potential energy source based on the highly exothermic oxidation of boron, density functional theory (DFT) was used to characterize small boron clusters with 0-3 oxygen atoms and a total of up to ten atoms. The structures, vibrational frequencies, and stabilities were calculated for each of these clusters. A quantum molecular dynamics procedure was used to locate the global minimum for each species, which proved to be crucial given the unintuitive structure of many of the most stable isomers. Additionally, due to the plane-wave, periodic DFT code used in this study a straightforward comparison of these clusters to the bulk boron and B2O3 structures was possible despite the great structural and energetic differences between the two forms. Through evaluation of previous computational and experimental work, the relevant low-energy structures of all but one of the pure boron clusters can be assigned with great certainty. Nearly all of the boron oxide clusters are described here for the first time, but there are strong indications that the DFT procedure chosen is particularly well suited for the task. Insight into the trends in boron and boron oxide cluster stabilities, as well as the ultimate limits of growth for each, are also provided. The work reported herein provides crucial information towards understanding the oxidation of boron at a molecular level.

Spiers memorial lecture. The quantum chemistry of d- and f-element complexes: from an approximate existence to functional happiness

The field of modern quantum inorganic chemistry is just over 50 years old, dating back to 1951, when quantitative LCAO molecular orbital theory was developed and ferrocene was discovered. This Lecture provides a survey of the development of the field through about 1980, which has led to its current state. The application of modern quantum chemical techniques are illustrated via two disparate examples from the authors' research group. First, the recent discovery of uranium-noble gas bonds is discussed including the synergy between the theoretical and experimental investigations of this phenomenon. New theoretical results using coupled-cluster [CCSD(T)] methodology is contrasted to the original scalar-relativistic density functional theory results. Second, new applications of time-dependent density functional theory to the rich photochemistry of a dinuclear organometallic complex, (eta5-C5H5)2Fe2(mu-CO)2(CO)2, are discussed.

Human and murine osteocalcin gene expression: conserved tissue restricted expression and divergent responses to 1,25-dihydroxyvitamin D3 in vivo

Human and murine osteocalcin genes demonstrate similar cell-specific expression patterns despite significant differences in gene locus organization and sequence variations in cis-acting regulatory elements. To investigate whether differences in these regulatory regions result in an altered response to 1,25-dihydroxyvitamin D3 [1,25-(OH)2D3] in vivo, we compared the response of the endogenous mouse osteocalcin gene to a bacterial reporter gene directed by flanking regions of the human osteocalcin gene in transgenic mice. Transgene expression colocalized with endogenous osteocalcin expression in serial sections, being detected in osteoblasts, osteocytes and hypertrophic chondrocytes. In calvarial cell culture lysates from transgenic and nontransgenic mice, the endogenous mouse osteocalcin gene did not respond to 1,25-(OH)2D3 treatment. Despite this, transgene activity was significantly increased in the same cells. Similarly, Northern blots of total cellular RNA and in situ hybridization studies of transgenic animals demonstrated a maximal increase in transgene expression at 6 h after 1,25-(OH)2D3 injection (23.6+/-3.6-fold) with a return to levels equivalent to uninjected animals by 24 h (1.2+/-0.1-fold). This increase in transgene expression was also observed at 6 h after 1,25-(OH)2D3 treatment in animals on a low calcium diet (25.2+/-7.7-fold) as well as in transgenic mice fed a vitamin D-deficient diet containing strontium chloride to block endogenous 1,25-(OH)2D3 production (7.5+/-0.9-fold). In contrast to the increased transgene expression levels, neither endogenous mouse osteocalcin mRNA levels nor serum osteocalcin levels were significantly altered after 1,25-(OH)2D3 injection in transgenic or nontransgenic mice, regardless of dietary manipulations, supporting evidence for different mechanisms regulating the response of human and mouse osteocalcin genes to 1,25-(OH)2D3. Although the cis- and trans-acting mechanisms directing cell-specific gene expression appear to be conserved in the mouse and human osteocalcin genes, responsiveness to 1,25-(OH)2D3 is not. The mouse osteocalcin genes do not respond to 1,25-(OH)2D3 treatment, but the human osteocalcin-directed transgene is markedly upregulated under the same conditions and in the same cells. The divergent responses of these homologous genes to 1,25-(OH)2D3 are therefore likely to be due to differences in mouse and human osteocalcin-regulatory sequences rather than to variation in the complement of trans-acting factors present in mouse osteoblastic cells. Increased understanding of these murine-human differences in osteocalcin regulation may shed light on the function of osteocalcin and its regulation by vitamin D in bone physiology.