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.

In Silico Modeling of PROTAC-Mediated Ternary Complexes: Validation and Application.

In this work, four methods are described and validated for generating in silico ensembles of PROTAC-mediated ternary complexes. Filters based on characteristics of the proposed ternary complexes are developed to identify those that resemble known crystal structures. We then show how to use these modeling techniques a priori to discriminate the PROTAC-mediated degradation behavior of a mutant protein vs its wild type, of three closely related targets, and among three different PROTAC molecules.

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.

A biomimetic pathway for vanadium-catalyzed aerobic oxidation of alcohols: evidence for a base-assisted dehydrogenation mechanism.

The first step in the catalytic oxidation of alcohols by molecular O(2), mediated by homogeneous vanadium(V) complexes [LV(V)(O)(OR)], is ligand exchange. The unusual mechanism of the subsequent intramolecular oxidation of benzyl alcoholate ligands in the 8-hydroxyquinolinato (HQ) complexes [(HQ)(2)V(V)(O)(OCH(2)C(6)H(4)-p-X)] involves intermolecular deprotonation. In the presence of triethylamine, complex 3 (X = H) reacts within an hour at room temperature to generate, quantitatively, [(HQ)(2)V(IV)(O)], benzaldehyde (0.5 equivalents), and benzyl alcohol (0.5 equivalents). The base plays a key role in the reaction: in its absence, less than 12% conversion was observed after 72 hours. The reaction is first order in both 3 and NEt(3), with activation parameters ΔH(≠)=(28±4) kJ mol(-1) and ΔS(≠)=(-169±4) J K(-1) mol(-1). A large kinetic isotope effect, 10.2±0.6, was observed when the benzylic hydrogen atoms were replaced by deuterium atoms. The effect of the para substituent of the benzyl alcoholate ligand on the reaction rate was investigated using a Hammett plot, which was constructed using σ(p). From the slope of the Hammett plot, ρ=+(1.34±0.18), a significant buildup of negative charge on the benzylic carbon atom in the transition state is inferred. These experimental findings, in combination with computational studies, support an unusual bimolecular pathway for the intramolecular redox reaction, in which the rate-limiting step is deprotonation at the benzylic position. This mechanism, that is, base-assisted dehydrogenation (BAD), represents a biomimetic pathway for transition-metal-mediated alcohol oxidations, differing from the previously identified hydride-transfer and radical pathways. It suggests a new way to enhance the activity and selectivity of vanadium catalysts in a wide range of redox reactions, through control of the outer coordination sphere.

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.

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.

A QM/QM multilayer composite methodology: The ONIOM correlation consistent composite approach (ONIOM-ccCA).

A QM/QM approach incorporating the correlation consistent composite approach (ccCA) into the ONIOM multilayer methodology has been implemented. This new multilayer composite scheme, ONIOM-ccCA, enables the accurate prediction of thermochemical properties for systems containing dozens-and potentially hundreds-of atoms. ONIOM-ccCA is used to predict the C-H bond dissociation energies of 18 anthracene and fluorene analogues, containing up to 43 atoms, to within 1.2 kcal mol(-1) of experimental values. Several density functional and basis set combinations are evaluated for use as the low level QM layer. The mean absolute deviation (1.2 kcal mol(-1)) using the most accurate ONIOM-ccCA method, ccCA:B3LYP/cc-pVTZ, is significantly lower than that of an earlier reported multilayer composite scheme [2.4 kcal mol(-1), Li, M.-J.; Liu, L.; Fu, Y.; Guo, Q.-X. J. Phys. Chem. B, 2005, 109, 13818]. Thus, through use of ONIOM-ccCA, accurate thermochemical calculations are now feasible for sizable molecular systems of chemical or biological interest.

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.

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.

Performance of density functional theory for 3d transition metal-containing complexes: utilization of the correlation consistent basis sets.

The performance of 44 density functionals used in conjunction with the correlation consistent basis sets (cc-pVnZ where n = T and Q) has been assessed for the gas-phase enthalpies of formation at 298.15 K of 3d transition metal (TM) containing systems. Nineteen molecules were examined: ScS, VO, VO(2), Cr(CO)(6), MnS, MnCl(2), Mn(CO)(5)Cl, FeCl(3), Fe(CO)(5), CoH(CO)(4), NiCl(2), Ni(CO)(4), CuH, CuF, CuCl, ZnH, ZnO, ZnCl, and Zn(CH(3))(2). Of the functionals examined, the functionals that resulted in the smallest mean absolute deviation (MAD, in parentheses, kcal mol(-1)) from experiment were B97-1 (6.9), PBE1KCIS (8.1), TPSS1KCIS (9.6), B97-2 (9.7), and B98 (10.7). All five of these functionals include some degree of Hartree-Fock (HF) exchange. The impact of increasing the basis set from cc-pVTZ to cc-pVQZ was found to be slight for the generalized gradient approximation (GGA) and meta-GGA (MGGA) functionals studied, indicating basis set saturation at the triple-zeta level. By contrast, for most of the generalized gradient exchange (GGE), hybrid GGA (HGGA), and hybrid meta-GGA (HMGGA) functionals considered, improvements in the average MAD of 2-3 kcal mol(-1) were seen upon progressing to a quadruple-zeta level basis set. Overall, it was found that the functionals that include Hartree-Fock exchange performed best overall, but those with greater than 40% HF exchange exhibit significantly poor performance for the prediction of enthalpies of formation for 3d TM complexes. Carbonyl-containing complexes, a mainstay in organometallic TM chemistry, are demonstrated to be exceedingly difficult to describe accurately with all but 2 of the 44 functionals considered. The most accurate functional, for both CO-containing and CO-free compounds, is B97-1/cc-pVQZ, which is shown to be capable of yielding results within 1 kcal mol(-1) of high-level ab initio composite methodologies.

Actin polymerization controls cilia-mediated signaling.

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.

Molecular Control of Atypical Protein Kinase C: Tipping the Balance between Self-Renewal and Differentiation.

Complex organisms are faced with the challenge of generating and maintaining diverse cell types, ranging from simple epithelia to neurons and motile immune cells [1-3]. To meet this challenge, a complex set of regulatory pathways controls nearly every aspect of cell growth and function, including genetic and epigenetic programming, cytoskeleton dynamics, and protein trafficking. The far reach of cell fate specification pathways makes it particularly catastrophic when they malfunction, both during development and for tissue homeostasis in adult organisms. Furthermore, the therapeutic promise of stem cells derives from their ability to deftly navigate the multitude of pathways that control cell fate [4]. How the molecular components making up these pathways function to specify cell fate is beginning to become clear. Work from diverse systems suggests that the atypical Protein Kinase C (aPKC) is a key regulator of cell fate decisions in metazoans [5-7]. Here, we examine some of the diverse physiological outcomes of aPKC's function in differentiation, along with the molecular pathways that control aPKC and those that are responsive to changes in its catalytic activity.

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.

The importance of secondary structure in determining CO2-protein binding patterns.

One potential means to decrease the level of atmospheric carbon dioxide is through the utilization of protein-CO(2) interactions. A recent bioinformatics analysis [Cundari TR et al. (2009) J Chem Inf Model 49:2111-2115] of these interactions revealed a marked disparity in CO(2) affinity between α-helices and ß-sheets. In order to understand this difference, a series of molecular dynamics simulations was performed on polypeptide model systems. Numerous factors that may influence CO(2) affinity were systematically investigated, including the specific location of the amino acids within the secondary structural elements (SSEs), the partial charges on CO(2), chemical modifications made to the protein backbone, the inclusion of singly, doubly, and many functionalized residues, and the effect of solvent water. The differing abilities of the secondary structure types to participate in hydrogen bonding along the backbone were identified as being a crucial influence on CO(2) affinity; the lesser role of polypeptide-CO(2) electrostatic interactions was also noted. The effect of incorporating functionalized amino acid side chains, such as those possessed by Arg and His, on the affinity differs between the two structure types, and also strongly depends on the number included and the distance between them. The inclusion of explicit water molecules was found to attenuate all interactions, but did not change the overall trends in CO(2) affinity. Collectively, these results highlight the role of the backbone atoms in binding the CO(2) ligand, which will have important implications for efforts to ameliorate atmospheric carbon dioxide through the use of natural, designed, and modified proteins.

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.

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.

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.