General principles of binding between cell surface receptors and multi-specific ligands: A computational study

Computational simulations of bispecific T cell engagers by a multiscale model

The use of bispecific antibodies as T cell engagers can bypass the normal T cell receptor-major histocompatibility class interaction, redirect the cytotoxic activity of T cells, and lead to highly efficient tumor cell killing. However, this immunotherapy also causes significant on-target off-tumor toxicologic effects, especially when it is used to treat solid tumors. To avoid these adverse events, it is necessary to understand the fundamental mechanisms involved in the physical process of T cell engagement. We developed a multiscale computational framework to reach this goal. The framework combines simulations on the intercellular and multicellular levels. On the intercellular level, we simulated the spatial-temporal dynamics of three-body interactions among bispecific antibodies, CD3 and tumor-associated antigens (TAAs). The derived number of intercellular bonds formed between CD3 and TAAs was further transferred to the multicellular simulations as the input parameter of adhesive density between cells. Through the simulations under various molecular and cellular conditions, we were able to gain new insights into how to adopt the most appropriate strategy to maximize the drug efficacy and avoid the off-target effect. For instance, we discovered that the low antibody-binding affinity resulted in the formation of large clusters at the cell-cell interface, which could be important to control the downstream signaling pathways. We also tested different molecular architectures of the bispecific antibody and suggested the existence of an optimal length in regulating the T cell engagement. Overall, the current multiscale simulations serve as a proof-of-concept study to help in the future design of new biological therapeutics.

Cheminformatics identification of modulators of key carbohydrate-metabolizing enzymes from C. cujete for type-2 diabetes mellitus intervention

Purpose

The therapeutic use of oral hypoglycaemic agents in the management of type-2 diabetes mellitus (T2DM) is without adverse effects; thus, calls for alternative and novel candidates from natural products in medicinal plants.

Method

The study explored molecular docking and molecular dynamics (MD) simulation approaches to identify key antidiabetic metabolites from Crescentia cujete.

Results

Molecular docking results identified four and/or five best compounds against each target enzyme (alpha-glucosidase, dipeptidyl peptidase-IV, aldose reductase, and protein tyrosine phosphatase-1B (PTP-1B)) implicated in diabetes. The resulting complexes (except against PTP-1B) had higher docking scores above respective standards (acarbose, Diprotin A, ranirestat). The MD simulation results revealed compounds such as benzoic acid (-48.414 kcal/mol) and phytol (-45.112 kcal/mol) as well as chlorogenic acid (-42.978 kcal/mol) and naringenin (-31.292 kcal/mol) had higher binding affinities than the standards [acarbose (-28.248 kcal/mol), ranirestat (-21.042 kcal/mol)] against alpha-glucosidase and aldose reductase, respectively while Diprotin A (-45.112 kcal/mol) and ursolic acid (-18.740 kcal/mol) presented superior binding affinities than the compounds [luteolin (-41.957 kcal/mol and naringenin (-16.518 kcal/mol)] against DPP-IV and PTP-1B respectively.

Conclusion

While isoflavone (alpha-glucosidase), xylocaine (DPP-IV), luteolin (aldose reductase,) and chlorogenic acid (PTP-1B) were affirmed as the best inhibitors of respective enzyme targets, luteolin, and chlorogenic acid may be suggested and proposed as probable candidates against T2DM and related retinopathy complication based on their structural stability, compactness and affinity for three (DPP-IV, aldose reductase, and PTP-1B) of the four targets investigated. Further studies are warranted in vitro and in vivo on the antihyperglycaemic effects of these drug candidates.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40200-023-01249-7.

Understanding the General Principles of T Cell Engagement by Multiscale Computational Simulations

The use of bispecific antibodies as T cell engagers can bypass the normal TCR-MHC interaction, redirect the cytotoxic activity of T-cells, and lead to highly efficient tumor cell killing. However, this immunotherapy also causes significant on-target off-tumor toxicologic effects, especially when they were used to treat solid tumors. In order to avoid these adverse events, it is necessary to understand the fundamental mechanisms during the physical process of T cell engagement. We developed a multiscale computational framework to reach this goal. The framework combines simulations on the intercellular and multicellular levels. On the intercellular level, we simulated the spatial-temporal dynamics of three-body interactions among bispecific antibodies, CD3 and TAA. The derived number of intercellular bonds formed between CD3 and TAA were further transferred into the multicellular simulations as the input parameter of adhesive density between cells. Through the simulations under various molecular and cellular conditions, we were able to gain new insights of how to adopt the most appropriate strategy to maximize the drug efficacy and avoid the off-target effect. For instance, we discovered that the low antibody binding affinity resulted in the formation of large clusters at the cell-cell interface, which could be important to control the downstream signaling pathways. We also tested different molecular architectures of the bispecific antibody and suggested the existence of an optimal length in regulating the T cell engagement. Overall, the current multiscale simulations serve as a prove-of-concept study to help the future design of new biological therapeutics.

SIGNIFICANCE

T-cell engagers are a class of anti-cancer drugs that can directly kill tumor cells by bringing T cells next to them. However, current treatments using T-cell engagers can cause serious side-effects. In order to reduce these effects, it is necessary to understand how T cells and tumor cells interact together through the connection of T-cell engagers. Unfortunately, this process is not well studied due to the limitations in current experimental techniques. We developed computational models on two different scales to simulate the physical process of T cell engagement. Our simulation results provide new insights into the general properties of T cell engagers. The new simulation methods can therefore serve as a useful tool to design novel antibodies for cancer immunotherapy.

Blockade of Activin Receptor IIB Protects Arthritis Pathogenesis by Non-Amplification of Activin A-ACVR2B-NOX4 Axis Pathway

Although activin receptor IIB (ACVR2B) is emerging as a novel pathogenic receptor, its ligand and assembled components (or assembly) are totally unknown in the context of osteoarthritis (OA) pathogenesis. The present results suggest that upregulation of ACVR2B and its assembly could affect osteoarthritic cartilage destruction. It is shown that the ACVR2B ligand, activin A, regulates catabolic factor expression through ACVR2B in OA development. Activin A Tg mice (Col2a1-Inhba) exhibit enhanced cartilage destruction, whereas heterozygous activin A KO mice (Inhba^+/- ) show protection from cartilage destruction. In silico analysis suggests that the Activin A-ACVR2B axis is involved in Nox4-dependent ROS production. Activin A Tg:Nox4 KO (Col2a1-Inhba:Nox4^-/- ) mice show inhibition of experimental OA pathogenesis. NOX4 directly binds to the C-terminal binding site on ACVR2B-ACVR1B and amplifies the pathogenic signal for cartilage destruction through SMAD2/3 signaling. Together, the findings reveal that the ACVR2B assembly, which comprises Activin A, ACVR2B, ACVR1B, Nox4, and AP-1-induced HIF-2α, accelerates OA development. Furthermore, it is shown that shRNA-mediated ACVR2B knockdown or trapping ligands of ACVR2B abrogate OA development by competitively disrupting the ACVR2B-Activin A interaction. These results suggest that the ACVR2B assembly is required to amplify osteoarthritic cartilage destruction and could be a potential therapeutic target in efforts to treat OA.

Dynamics and Patterning of 5-Hydroxytryptamine 2 Subtype Receptors in JC Polyomavirus Entry

The organization and dynamics of plasma membrane receptors are a critical link in virus-receptor interactions, which finetune signaling efficiency and determine cellular responses during infection. Characterizing the mechanisms responsible for the active rearrangement and clustering of receptors may aid in developing novel strategies for the therapeutic treatment of viruses. Virus-receptor interactions are poorly understood at the nanoscale, yet they present an attractive target for the design of drugs and for the illumination of viral infection and pathogenesis. This study utilizes super-resolution microscopy and related techniques, which surpass traditional microscopy resolution limitations, to provide both a spatial and temporal assessment of the interactions of human JC polyomavirus (JCPyV) with 5-hydroxytrypamine 2 receptors (5-HT2Rs) subtypes during viral entry. JCPyV causes asymptomatic kidney infection in the majority of the population and can cause fatal brain disease, and progressive multifocal leukoencephalopathy (PML), in immunocompromised individuals. Using Fluorescence Photoactivation Localization Microscopy (FPALM), the colocalization of JCPyV with 5-HT2 receptor subtypes (5-HT2A, 5-HT2B, and 5-HT2C) during viral attachment and viral entry was analyzed. JCPyV was found to significantly enhance the clustering of 5-HT2 receptors during entry. Cluster analysis of infected cells reveals changes in 5-HT2 receptor cluster attributes, and radial distribution function (RDF) analyses suggest a significant increase in the aggregation of JCPyV particles colocalized with 5-HT2 receptor clusters in JCPyV-infected samples. These findings provide novel insights into receptor patterning during viral entry and highlight improved technologies for the future development of therapies for JCPyV infection as well as therapies for diseases involving 5-HT2 receptors.

Distributed Acceleration of Adhesive Dynamics Simulations

Cell adhesion plays a critical role in processes ranging from leukocyte migration to cancer cell transport during metastasis. Adhesive cell interactions can occur over large distances in microvessel networks with cells traveling over distances much greater than the length scale of their own diameter. Therefore, biologically relevant investigations necessitate efficient modeling of large field-of-view domains, but current models are limited by simulating such geometries at the sub-micron scale required to model adhesive interactions which greatly increases the computational requirements for even small domain sizes. In this study we introduce a hybrid scheme reliant on both on-node and distributed parallelism to accelerate a fully deformable adhesive dynamics cell model. This scheme leads to performant system usage of modern supercomputers which use a many-core per-node architecture. On-node acceleration is augmented by a combination of spatial data structures and algorithmic changes to lessen the need for atomic operations. This deformable adhesive cell model accelerated with hybrid parallelization allows us to bridge the gap between high-resolution cell models which can capture the sub-micron adhesive interactions between the cell and its microenvironment, and large-scale fluid-structure interaction (FSI) models which can track cells over considerable distances. By integrating the sub-micron simulation environment into a distributed FSI simulation we enable the study of previously unfeasible research questions involving numerous adhesive cells in microvessel networks such as cancer cell transport through the microcirculation.

Expression of Recombinant CTLA-4 and PD-L1 Proteins Fused with Thioredoxin, and Determination of Their Ligand-Binding Activities

BACKGROUND

The use of chimeric proteins that selectively interact with various immune cell receptors to treat oncology patients has increased. One effective way to obtain recombinant proteins is to use the E. coli expression system. However, in eukaryotic protein production in E. coli, several difficulties arise that can be solved by fusing the target protein with thioredoxin. Thioredoxin can enhance solubility, but its large size can lead to an erroneous assessment of protein solubility, folding, and activity. The present study examined the ligand-binding activity of PD-L1, and CTLA-4 receptors fused with thioredoxin.

METHODS

The de novo synthesized genes of the extracellular domains of the PD-L1 and CTLA-4 were cloned into the pET28 and pET32 expression plasmids and used to transform E. coli BL21 cells. Purified recombinant proteins were characterized by western blotting, LC-MS/MS spectrometry, and ELISA.

RESULTS

Amino acid sequence comparisons of the recombinant proteins obtained by LC-MS/MS with the SwissProt database resulted in the highest comparison scores from 4950 to 13396. The binding efficiencies of recombinant human B7-1 Fc to rCTLA-4 and rTrx-CTLA-4 proteins in ELISA did not differ significantly. Similar results were obtained with recombinant rhesus monkey PD-1 hFc against rPD-L1 and rTrx-PD-L1.

CONCLUSION

Recombinant proteins specifically reacted with recombinant human B7-1 Fc and recombinant rhesus monkey PD-1 hFc. The fusion of thioredoxin with recombinant proteins through linkers slightly affected the activity of the extracellular domains of CTLA-4 and PD-L1.

G Protein-Coupled Receptors as Potential Intercellular Communication Mediators in Trypanosomatidae

Detection and transduction of environmental signals, constitute a prerequisite for successful parasite invasion; i.e., Leishmania transmission, survival, pathogenesis and disease manifestation and dissemination, with diverse molecules functioning as inter-cellular signaling ligands. Receptors [i.e., G protein-coupled receptors (GPCRs)] and their associated transduction mechanisms, well conserved through evolution, specialize in this function. However, canonical GPCR-related signal transduction systems have not been described in Leishmania, although orthologs, with reduced domains and function, have been identified in Trypanosomatidae. These inter-cellular communication means seem to be essential for multicellular and unicellular organism’s survival. GPCRs are flexible in their molecular architecture and may interact with the so-called receptor activity-modifying proteins (RAMPs), which modulate their function, changing GPCRs pharmacology, acting as chaperones and regulating signaling and/or trafficking in a receptor-dependent manner. In the skin, vasoactive- and neuro- peptides released in response to the noxious stimuli represented by the insect bite may trigger parasite physiological responses, for example, chemotaxis. For instance, in Leishmania (V.) braziliensis, sensory [Substance P, SP, chemoattractant] and autonomic [Vasoactive Intestinal Peptide, VIP, and Neuropeptide Y, NPY, chemorepellent] neuropeptides at physiological levels stimulate in vitro effects on parasite taxis. VIP and NPY chemotactic effects are impaired by their corresponding receptor antagonists, suggesting that the stimulated responses might be mediated by putative GPCRs (with essential conserved receptor domains); the effect of SP is blocked by [(D-Pro 2, D-Trp7,9]-Substance P (10^-6 M)] suggesting that it might be mediated by neurokinin-1 transmembrane receptors. Additionally, vasoactive molecules like Calcitonin Gene-Related Peptide [CGRP] and Adrenomedullin [AM], exert a chemorepellent effect and increase the expression of a 24 kDa band recognized in western blot analysis by (human-)-RAMP-2 antibodies. In-silico search oriented towards GPCRs-like receptors and signaling cascades detected a RAMP-2-aligned sequence corresponding to Leishmania folylpolyglutamate synthase and a RAMP-3 aligned protein, a hypothetical Leishmania protein with yet unknown function, suggesting that in Leishmania, CGRP and AM activities may be modulated by RAMP- (-2) and (-3) homologs. The possible presence of proteins and molecules potentially involved in GPCRs cascades, i.e., RAMPs, signpost conservation of ancient signaling systems associated with responses, fundamental for cell survival, (i.e., taxis and migration) and may constitute an open field for description of pharmacophores against Leishmania parasites.

Targeted Delivery of Chemotherapeutic Agents for Osteosarcoma Treatment

Since osteosarcoma (OS) is an aggressive bone cancer with unknown molecular pathways of etiology and pathophysiology, improving patient survival has long been a challenge. The conventional therapy is a complex multidisciplinary management that include radiotherapy, chemotherapy which followed by surgery and then post-operative adjuvant chemotherapy. However, they have severe side effects because the majority of the medicines used have just a minor selectivity for malignant tissue. As a result, treating tumor cells specifically without damaging healthy tissue is currently a primary goal in OS therapy. The coupling of chemotherapeutic drugs with targeting ligands is a unique therapy method for OS that, by active targeting, can overcome the aforementioned hurdles. This review focuses on advances in ligands and chemotherapeutic agents employed in targeted delivery to improve the capacity of active targeting and provide some insight into future therapeutic research for OS.

Mesenchymal stem cells transfer mitochondria to allogeneic Tregs in an HLA-dependent manner improving their immunosuppressive activity

Cell-based immunotherapies can provide safe and effective treatments for various disorders including autoimmunity, cancer, and excessive proinflammatory events in sepsis or viral infections. However, to achieve this goal there is a need for deeper understanding of mechanisms of the intercellular interactions. Regulatory T cells (Tregs) are a lymphocyte subset that maintain peripheral tolerance, whilst mesenchymal stem cells (MSCs) are multipotent nonhematopoietic progenitor cells. Despite coming from different origins, Tregs and MSCs share immunoregulatory properties that have been tested in clinical trials. Here we demonstrate how direct and indirect contact with allogenic MSCs improves Tregs’ potential for accumulation of immunosuppressive adenosine and suppression of conventional T cell proliferation, making them more potent therapeutic tools. Our results also demonstrate that direct communication between Tregs and MSCs is based on transfer of active mitochondria and fragments of plasma membrane from MSCs to Tregs, an event that is HLA-dependent and associates with HLA-C and HLA-DRB1 eplet mismatch load between Treg and MSC donors.

Regulatory T (Treg) cells and mesenchymal stem cells (MSCs) are both cell populations capable of immune tolerance induction. Here the authors show that the transfer of mitochondria from mesenchymal stem cells to allogeneic Treg cells in an HLA-dependent manner results in enhanced immunosuppressive functions of Treg cells.

Phage Display-Derived Peptide for the Specific Binding of SARS-CoV-2

Beginning from the end of 2019, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic swept all over the world and is still afflicting the whole global population. Given that the vaccine-manufacturing ability is limited and the virus can evolve quickly, vaccination alone may not be able to end the pandemic, thus developing fast and accurate diagnoses and effective therapeutics will always be unmet needs. Phage display peptide library has been used in screening antigen-specific peptides for the invention of novel mimic receptors/ligands. Here, we report that a 12-mer phage display peptide library has been screened against the SARS-CoV-2 receptor-binding domain (RBD), and five of the screened peptides show binding ability with the RBD protein by the enzyme-linked immune sorbent assay. The surface plasmon resonance assay further demonstrates that peptide no. 1 can specifically bind to SARS-CoV-2 RBD with a binding affinity constant (K d) of 5.8 μM. Transmission electron microscopy coupled with a magnetic bead assay further confirms that the screened peptide can specifically bind the inactivated SARS-CoV-2 virus. This SARS-CoV-2-specific peptide holds great promise as a new bioreceptor/ligand for the rapid and accurate detection of SARS-CoV-2.

Time dependent adhesion of cells on nanorough surfaces

Understanding and controlling the mechanisms of cell adhesion to nanomaterials is essential in tissue engineering, regenerative medicine, the development of experimental models for the study of neurodegenerative diseases. Nonetheless, despite the great many of studies that have examined how cells interact with nanoscale surfaces, little is known about the temporal dimension of the process of adhesion. In a previous work, Decuzzi and Ferrari, by examining how the energy of a cell changes while binding to a nanoscale surface, determined a criterion to decide whether nanoroughness can either enhance or retard cell adhesion. While accurate, however their model template disregards the time variable. Here, starting from the work of Decuzzi and Ferrari, we have developed a mathematical model based on chemotaxis that describes how cells adhere to a nanorough surface over time. Relaxing the originating constraint of a fixed density of ligand molecules expressed by the cell membrane, we show that the strength of adhesion depends on time and that, for certain values of the model parameters, a cell can arrive to establish a stable adhesion to a substrate even if the process of binding is initially energetically unfavourable. We show that, for a cell-membrane stiffness of 10kPa, an initial density of receptors of 500bonds/μm², a specific and non-specific energy density of adhesion of 10^-5J/m² and 10^-7J/m², and roughness in the low nanometer range, cell adhesion forces can be completely activated from few seconds to some tens of minutes from the initial contact with the surface.

Stimulus-specific responses in innate immunity: Multilayered regulatory circuits

Immune sentinel cells initiate immune responses to pathogens and tissue injury and are capable of producing highly stimulus-specific responses. Insight into the mechanisms underlying such specificity has come from the identification of regulatory factors and biochemical pathways, as well as the definition of signaling circuits that enable combinatorial and temporal coding of information. Here, we review the multi-layered molecular mechanisms that underlie stimulus-specific gene expression in macrophages. We categorize components of inflammatory and anti-pathogenic signaling pathways into five layers of regulatory control and discuss unifying mechanisms determining signaling characteristics at each layer. In this context, we review mechanisms that enable combinatorial and temporal encoding of information, identify recurring regulatory motifs and principles, and present strategies for integrating experimental and computational approaches toward the understanding of signaling specificity in innate immunity.

A DNA Molecular Robot that Autonomously Walks on the Cell Membrane to Drive Cell Motility

Synthetic molecular robots can execute sophisticated molecular tasks at nanometer resolution. However, a molecular robot capable of controlling cellular behavior remains unexplored. Herein, we report a self-propelled DNA robot operating on the cell membrane to control the migration of a cell. Driven by DNAzyme catalytic activity, the DNA robot could autonomously and stepwise move on the membrane-floating cell-surface receptors in a stochastic manner and simultaneously trigger the receptor-dimerization to activate downstream signaling for cell motility. The cell membrane-associated continuous motion and operation of a DNA robot allowed for the ultrasensitive regulation of MET/AKT signaling and cytoskeleton remodeling to enhance cell migration. Finally, we designed distinct conditional DNA robots to orthogonally manipulate the cell migration in a coculture of mixed cell populations. We have developed a novel strategy to engineer a cell-driving molecular robot, representing a promising avenue for precise cell manipulation with nanoscale resolution.

Estrogen Receptor β Contributes to Both Hypertension and Hypothalamic Plasticity in a Mouse Model of Peri-Menopause

Hypertension susceptibility in women increases at the transition to menopause, termed perimenopause, a state characterized by erratic estrogen fluctuation and extended hormone cycles. Elucidating the role of estrogen signaling in the emergence of hypertension during perimenopause has been hindered by animal models that are confounded by abrupt estrogen cessation or effects of aging. In the present study, accelerated ovarian failure (AOF) in estrogen receptor β (ERβ) reporter mice was induced by 4-vinylcyclohexene diepoxide in young mice to model early-stage ovarian failure (peri-AOF) characteristic of peri-menopause. It was found that administering ERβ agonists suppressed elevated blood pressure in a model of neurogenic hypertension induced by angiotensin II (AngII) in peri-AOF, but not in age-matched male mice. It was also found that ERβ agonist administration in peri-AOF females, but not males, suppressed the heightened NMDAR signaling and reactive oxygen production in ERβ neurons in the hypothalamic paraventricular nucleus (PVN), a critical neural regulator of blood pressure. It was further shown that deleting ERβ in the PVN of gonadally intact females produced a phenotype marked by a sensitivity to AngII hypertension. These results suggest that ERβ signaling in the PVN plays an important role in blood pressure regulation in female mice and contributes to hypertension susceptibility in females at an early stage of ovarian failure comparable to human perimenopause.

Mechanistic dissection of spatial organization in NF-κB signaling pathways by hybrid simulations

The nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) is one of the most important transcription factors involved in the regulation of inflammatory signaling pathways. Inappropriate activation of these pathways has been linked to autoimmunity and cancers. Emerging experimental evidences have been showing the existence of elaborate spatial organizations for various molecular components in the pathways. One example is the scaffold protein tumor necrosis factor receptor associated factor (TRAF). While most TRAF proteins form trimeric quaternary structure through their coiled-coil regions, the N-terminal region of some members in the family can further be dimerized. This dimerization of TRAF trimers can drive them into higher-order clusters as a response to receptor stimulation, which functions as a spatial platform to mediate the downstream poly-ubiquitination. However, the molecular mechanism underlying the TRAF protein clustering and its functional impacts are not well-understood. In this article, we developed a hybrid simulation method to tackle this problem. The assembly of TRAF-based signaling platform at the membrane-proximal region is modeled with spatial resolution, while the dynamics of downstream signaling network, including the negative feedbacks through various signaling inhibitors, is simulated as stochastic chemical reactions. These two algorithms are further synchronized under a multiscale simulation framework. Using this computational model, we illustrated that the formation of TRAF signaling platform can trigger an oscillatory NF-κB response. We further demonstrated that the temporal patterns of downstream signal oscillations are closely regulated by the spatial factors of TRAF clustering, such as the geometry and energy of dimerization between TRAF trimers. In general, our study sheds light on the basic mechanism of NF-κB signaling pathway and highlights the functional importance of spatial regulation within the pathway. The simulation framework also showcases its potential of application to other signaling pathways in cells.

Metallic nanoparticles as drug delivery system for the treatment of cancer

INTRODUCTION

The targeted delivery of anticancer agents to tumor is a major challenge because most of the drugs show off-target effect resulting in nonspecific cell death. Multifunctionalized metallic nanoparticles (NPs) are explored as new carrier system in the era of cancer therapeutics. Researchers investigated the potential of metallic NPs to target tumor cells by active and passive mechanisms, thereby reducing off-target effects of anticancer agents. Moreover, photocatalytic activity of upconversion nanoparticles (UCNPs) and the enhanced permeation and retention (EPR) effect have also gained wide potential in cancer treatment. Recent advancement in the field of nanotechnology highlights their potency for cancer therapy.

AREAS COVERED

This review summarizes the types of gold and silver metallic NPs with targeting mechanisms and their potentiality in cancer therapy.

EXPERT OPINION

Recent advances in the field of nanotechnology for cancer therapy offer high specificity and targeting efficiency. Targeting tumor cells through mechanistic pathways using metallic NPs for the disruption/alteration of molecular profile and survival rate of the tumor cells has led to an effective approach for cancer therapeutics. This alteration in the survival rate of the tumor cells might decrease the proliferation thereby resulting in more efficient management in the treatment of cancer.

A computational study of co-inhibitory immune complex assembly at the interface between T cells and antigen presenting cells

The activation and differentiation of T-cells are mainly directly by their co-regulatory receptors. T lymphocyte-associated protein-4 (CTLA-4) and programed cell death-1 (PD-1) are two of the most important co-regulatory receptors. Binding of PD-1 and CTLA-4 with their corresponding ligands programed cell death-ligand 1 (PD-L1) and B7 on the antigen presenting cells (APC) activates two central co-inhibitory signaling pathways to suppress T cell functions. Interestingly, recent experiments have identified a new cis-interaction between PD-L1 and B7, suggesting that a crosstalk exists between two co-inhibitory receptors and the two pairs of ligand-receptor complexes can undergo dynamic oligomerization. Inspired by these experimental evidences, we developed a coarse-grained model to characterize the assembling of an immune complex consisting of CLTA-4, B7, PD-L1 and PD-1. These four proteins and their interactions form a small network motif. The temporal dynamics and spatial pattern formation of this network was simulated by a diffusion-reaction algorithm. Our simulation method incorporates the membrane confinement of cell surface proteins and geometric arrangement of different binding interfaces between these proteins. A wide range of binding constants was tested for the interactions involved in the network. Interestingly, we show that the CTLA-4/B7 ligand-receptor complexes can first form linear oligomers, while these oligomers further align together into two-dimensional clusters. Similar phenomenon has also been observed in other systems of cell surface proteins. Our test results further indicate that both co-inhibitory signaling pathways activated by B7 and PD-L1 can be down-regulated by the new cis-interaction between these two ligands, consistent with previous experimental evidences. Finally, the simulations also suggest that the dynamic and the spatial properties of the immune complex assembly are highly determined by the energetics of molecular interactions in the network. Our study, therefore, brings new insights to the co-regulatory mechanisms of T cell activation.

The activation of a T cell can be regulated by the receptors on its surface, such as CTLA-4 and PD-1. People used to think that these two receptors inhibit T cell activation through distinct pathways. However, recent experiments discovered that the ligands of these two receptors, B7 and PD-L1, can interact with each other on the same surface of antigen presenting cells. Here we utilized computational simulations to investigate functional roles of this newly discovered interaction in T cell coregulation. The specific environment of interface between T cell and antigen presenting cell has been taken into account of our model. Ligand and receptors randomly diffuse within this interface area. They further involve in different types of interactions, with each other from the same side or the opposite side of cell surface. Using this method, we found ligands and receptors can not only form complexes, but also aggregate into large-scale clusters. We also demonstrated that the engagement between B7 and PD-L1 can reduce the interactions with their corresponding receptors. This study, therefore, offers new insights to our understanding of signal regulation in T cells.

Understanding the impacts of cellular environments on ligand binding of membrane receptors by computational simulations

Binding of cell surface receptors with their extracellular ligands initiates various intracellular signaling pathways. However, our understanding of the cellular functions of these receptors is very limited due to the fact that in vivo binding between ligands and receptors has only been successfully measured in a very small number of cases. In living cells, receptors are anchored on surfaces of the plasma membrane, which undergoes thermal undulations. Moreover, it has been observed in various systems that receptors can be organized into oligomers prior to ligand binding. It is not well understood how these cellular factors play roles in regulating the dynamics of ligand-receptor interactions. Here, we tackled these problems by using a coarse-grained kinetic Monte Carlo simulation method. Using this method, we demonstrated that the membrane undulations cause a negative effect on ligand-receptor interactions. We further found that the preassembly of membrane receptors on the cell surface can not only accelerate the kinetics of ligand binding but also reduce the noises during the process. In general, our study highlights the importance of membrane environments in regulating the function of membrane receptors in cells. The simulation method can be potentially applied to specific receptor systems involved in cell signaling.

Regulatory role of insulin-like growth factor-binding proteins in odontogenic mineralization in rats

Much information is currently available for molecules in early odontogenesis, but there is limited knowledge regarding terminal cytodifferentiation of ameloblasts and odontoblasts for the determination of normal crown morphology. The present differential display PCR (DD-PCR) revealed that insulin-like growth factor-binding protein 5 (IGFBP5) was differentially expressed in molar tooth germs between the cap (before crown mineralization) and root formation (after crown mineralization) stages. Real-time PCR confirmed that the expression levels of IGFBP1-4 were not significantly changed but those of IGFBP5-7 were upregulated in a time-dependent manner. Immunoreactivities for IGFBP5-7 were hardly seen in molar germs at the cap/early bell stage and protective-stage ameloblasts at the root formation stage. However, the reactivity was strong in odontoblasts and maturation-stage ameloblasts, which are morphologically and functionally characterized by wide intercellular space and active enamel matrix mineralization. The localization of each IGFBP was temporospatial. IGFBP5 was localized in the nuclei of fully differentiated odontoblasts and ameloblasts, while IGFBP6 was localized in the apical cytoplasm of ameloblasts and odontoblasts with dentinal tubules, and IGFBP7 was mainly found in the whole cytoplasm of odontoblasts and the intercellular space of ameloblasts. IGFBP silencing using specific siRNAs upregulated representative genes for dentinogenesis and amelogenesis, such as DMP1 and amelogenin, respectively, and augmented the differentiation media-induced mineralization, which was confirmed by alizarin red s and alkaline phosphatase staining. These results suggest that IGFBP5-7 may play independent and redundant regulatory roles in late-stage odontogenesis by modulating the functional differentiation of ameloblasts and odontoblasts.

Jennings C, Parekh SH. Effect of ambient temperature on respiratory tract cells exposed to SARS-CoV-2 viral mimicking nanospheres-An experimental study

The novel coronavirus caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has reached more than 160 countries and has been declared a pandemic. SARS-CoV-2 infects host cells by binding to the angiotensin-converting enzyme 2 (ACE-2) surface receptor via the spike (S) receptor-binding protein (RBD) on the virus envelope. Global data on a similar infectious disease spread by SARS-CoV-1 in 2002 indicated improved stability of the virus at lower temperatures facilitating its high transmission in the community during colder months (December-February). Seasonal viral transmissions are strongly modulated by temperatures, which can impact viral trafficking into host cells; however, an experimental study of temperature-dependent activity of SARS-CoV-2 is still lacking. We mimicked SARS-CoV-2 with polymer beads coated with the SARS-CoV-2 S protein to study the effect of seasonal temperatures on the binding of virus-mimicking nanospheres to lung epithelia. The presence of the S protein RBD on nanosphere surfaces led to binding by Calu-3 airway epithelial cells via the ACE-2 receptor. Calu-3 and control fibroblast cells with S-RBD-coated nanospheres were incubated at 33 and 37 °C to mimic temperature fluctuations in the host respiratory tract, and we found no temperature dependence in contrast to nonspecific binding of bovine serum ablumin-coated nanospheres. Moreover, the ambient temperature changes from 4 to 40 °C had no effect on S-RBD-ACE-2 ligand-receptor binding and minimal effect on the S-RBD protein structure (up to 40 °C), though protein denaturing occurred at 51 °C. Our results suggest that ambient temperatures from 4 to 40 °C have little effect on the SARS-CoV-2-ACE-2 interaction in agreement with the infection data currently reported.

Understanding the Targeting Mechanisms of Multi-Specific Biologics in Immunotherapy with Multiscale Modeling

Immunotherapeutics are frequently associated with adverse side effects due to the elicitation of global immune modulation. To lower the risk of these side effects, recombinant DNA technology is employed to enhance the selectivity of cell targeting by genetically fusing different biomolecules, yielding new species referred to as multi-specific biologics. The design of new multi-specific biologics is a central challenge for the realization of new immunotherapies. To understand the molecular determinants responsible for regulating the binding between multi-specific biologics and surface-bound membrane receptors, we developed a multiscale computational framework that integrates various simulation approaches covering different timescales and spatial resolutions. Our model system of multi-specific biologics contains two natural ligands of immune receptors, which are covalently tethered by a peptide linker. Using this method, a number of interesting features of multi-specific biologics were identified. Our study therefore provides an important strategy to design the next-generation biologics for immunotherapy.

•

Two proteins are connected by different linkers as a model of bispecific biologics

•

Conformational dynamics of biologics are captured by microsecond MD simulations

•

Coarse-grained simulations are used to test binding between biologics and receptors

•

Biologics with long and flexible linkers are more efficient in targeting receptors

Human recombinant lysosomal β-Hexosaminidases produced in Pichia pastoris efficiently reduced lipid accumulation in Tay-Sachs fibroblasts

GM2 gangliosidosis, Tay-Sachs and Sandhoff diseases, are lysosomal storage disorders characterized by the lysosomal accumulation of GM2 gangliosides. This accumulation is due to deficiency in the activity of the β-hexosaminidases Hex-A or Hex-B, which are dimeric hydrolases formed by αβ or ββ subunits, respectively. These disorders show similar clinical manifestations that range from mild systemic symptoms to neurological damage and premature death. There is still no effective therapy for GM2 gangliosidoses, but some therapeutic alternatives, as enzyme replacement therapy, have being evaluated. Previously, we reported the production of active human recombinant β-hexosaminidases (rhHex-A and rhHex-B) in the methylotrophic yeast Pichia pastoris. In this study, we evaluated in vitro the cellular uptake, intracellular delivery to lysosome, and reduction of stored substrates. Both enzymes were taken-up via endocytic pathway mediated by mannose and mannose-6-phosphate receptors and delivered to lysosomes. Noteworthy, rhHex-A diminished the levels of stored lipids and lysosome mass in fibroblasts from Tay-Sachs patients. Overall, these results confirm the potential of P. pastoris as host to produce recombinant β-hexosaminidases intended to be used in the treatment of GM2 gangliosidosis.

A Multiscale and Comparative Model for Receptor Binding of 2019 Novel Coronavirus and the Implication of its Life Cycle in Host Cells

The respiratory syndrome caused by a new type of coronavirus has been emerging from China and caused more than one million death globally since December 2019. This new virus, called severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) uses the same receptor called Angiotensin-converting enzyme 2 (ACE2) to attack humans as the coronavirus that caused the severe acute respiratory syndrome (SARS) seventeen years ago. Both viruses recognize ACE2 through the spike proteins (S-protein) on their surfaces. It was found that the S-protein from the SARS coronavirus (SARS-CoV) bind stronger to ACE2 than SARS-CoV-2. However, function of a bio-system is often under kinetic, rather than thermodynamic, control. To address this issue, we constructed a structural model for complex formed between ACE2 and the S-protein from SARS-CoV-2, so that the rate of their association can be estimated and compared with the binding of S-protein from SARS-CoV by a multiscale simulation method. Our simulation results suggest that the association of new virus to the receptor is slower than SARS, which is consistent with the experimental data obtained very recently. We further integrated this difference of association rate between virus and receptor into a mathematical model which describes the life cycle of virus in host cells and its interplay with the innate immune system. Interestingly, we found that the slower association between virus and receptor can result in longer incubation period, while still maintaining a relatively higher level of viral concentration in human body. Our computational study therefore provides, from the molecular level, one possible explanation that this new pandemic by far spread much faster than SARS.

C1188D mutation abolishes specific recognition between MLL1-CXXC domain and CpG site by inducing conformational switch of flexible N-terminal

Mixed lineage leukemia protein (MLL1 protein) recognizes the CpG site via its CXXC domain and is frequently associated with leukemia. The specific recognition is abolished by C1188D mutation, which also prevents MLL-related leukemia. In this paper, multiple molecular dynamic (MD) simulations were performed to investigate the mechanism of recognition and influences of C1188D mutation. Started from fully dissociated DNA and MLL1-CXXC domain, remarkably, the center of mass (COM) of MLL1-CXXC domain quickly concentrates on the vicinity of the CpG site in all 53 short MD simulations. Extended simulations of the wild type showed that the native complex formed in 500 ns among 4 of 53 simulations. In contrast, the C1188D mutant COM distributed broadly around the DNA and the native complex was not observed in any of the extended simulations. Simulations on the apo MLL1-CXXC domain further suggest that the wild type protein remained predominantly in an open form that closely resembles its structure in the native complex whereas C1188D mutant formed predominantly compact structures in which the N- terminal bends to D1188. This conformational switch hinders the formation of encounter complex, thus abolishes the recognition. Our study also provides clues to the study mechanism of recognition, by the CXXC domain from proteins like DNA methyltransferase and ten-eleven translocation enzymes.

Surface engineering of nanoparticles with ligands for targeted delivery to osteosarcoma

Osteosarcoma is one of the most common malignant bone tumors which affect adolescents. Neoadjuvant chemotherapy followed by operation has become recommended for osteosarcoma treatment. Whereas, the effects of conventional chemotherapy are unsatisfactory because of multidrug resistance, fast clearance rate, nontargeted delivery, side effects and so on. Accordingly, Nanoparticle-mediated targeted drug delivery system (NTDDS) is recommended to be a novel treatment strategy for osteosarcoma. NTDDS can overcome the above obstacles by enhanced permeability and retention effect and active targeting. The active targeting of the delivery system is mainly based on ligands. In this study, we investigate and summarize the most common ligands used in the latest NTDDS for osteosarcoma. It might provide new insights into nanomedicine for osteosarcoma treatment.

A computational model for understanding the oligomerization mechanisms of TNF receptor superfamily

By recognizing members in the tumor necrosis factor (TNF) receptor superfamily, TNF ligand proteins function as extracellular cytokines to activate various signaling pathways involved in inflammation, proliferation, and apoptosis. Most ligands in TNF superfamily are trimeric and can simultaneously bind to three receptors on cell surfaces. It has been experimentally observed that the formation of these molecular complexes further triggers the oligomerization of TNF receptors, which in turn regulate the intracellular signaling processes by providing transient compartmentalization in the membrane proximal regions of cytoplasm. In order to decode the molecular mechanisms of oligomerization in TNF receptor superfamily, we developed a new computational method that can physically simulate the spatial-temporal process of binding between TNF ligands and their receptors. The simulations show that the TNF receptors can be organized into hexagonal oligomers. The formation of this spatial pattern is highly dependent not only on the molecular properties such as the affinities of trans and cis binding, but also on the cellular factors such as the concentration of TNF ligands in the extracellular area or the density of TNF receptors on cell surfaces. Moreover, our model suggests that if TNF receptors are pre-organized into dimers before ligand binding, these lateral interactions between receptor monomers can play a positive role in stabilizing the ligand-receptor interactions, as well as in regulating the kinetics of receptor oligomerization. Altogether, this method throws lights on the mechanisms of TNF ligand-receptor interactions in cellular environments.

A Multiscale Model for the Self-Assembly of Coat Proteins in Bacteriophage MS2

The self-assembly of viral capsids is an essential step to the formation of infectious viruses. Elucidating the kinetic mechanisms of how a capsid or virus-like particle assembles could advance our knowledge about the viral lifecycle, as well as the general principles in self-assembly of biomaterials. However, current understanding of capsid assembly remains incomplete for many viruses due to the fact that the transient intermediates along the assembling pathways are experimentally difficult to be detected. In this paper, we constructed a new multiscale computational framework to simulate the self-assembly of virus-like particles. We applied our method to the coat proteins of bacteriophage MS2 as a specific model system. This virus-like particle of bacteriophage MS2 has a unique feature that its 90 sequence-identical dimers can be classified into two structurally various groups: one is the symmetric CC dimer, and the other is the asymmetric AB dimer. The homotypic interactions between AB dimers result in a 5-fold symmetric contact, while the heterotypic interactions between AB and CC dimers result in 6-fold symmetric contact. We found that the assembly can be described as a physical process of phase transition that is regulated by various factors such as concentration and specific stoichiometry between AB and CC dimers. Our simulations also demonstrate that heterotypic and homotypic interfaces play distinctive roles in modulating the assembling kinetics. The interaction between AB and CC dimers is much more dynamic than that between two AB dimers. We therefore suggest that the alternate growth of viral capsid through the heterotypic dimer interactions dominates the assembling pathways. This is, to the best of our knowledge, the first multiscale model to simulate the assembling process of coat proteins in bacteriophage MS2. The generality of this approach opens the door to its further applications in assembly of other viral capsids, virus-like particles, and novel drug delivery systems.

Etiologic Role of Kinases in the Progression of Human Cancers and Its Targeting Strategies

Cancer is one of the dominant causes of death worldwide while lifelong prognosis is still inauspicious. The maturation of the cancer is seen as a process of transformation of a healthy cell into a tumor-sensitive cell, which is held entirely at the cellular, molecular, and genetic levels of the organism. Tyrosine kinases can play a major, etiologic role in the inception of malignancy and devote to the uncontrolled proliferation of cancerous cells and the progression of a tumor as well as the development of metastatic disease. Angiogenesis and oncogene activation are the major event in cell proliferation. The growth of a tumor and metastasis are fully depending on angiogenesis and lymphangiogenesis triggered by chemical signals from tumor cells in a phase of rapid growth. Tyrosine kinase inhibitors are compounds that inhibit tyrosine kinases and effective in targeting angiogenesis and blocking the signaling pathways of oncogenes. Small molecule tyrosine kinase inhibitors like afatinib, erlotinib, crizotinib, gefitinib, and cetuximab are shown to a selective cut off tactic toward the constitutive activation of an oncogene in tumor cells, and thus contemplated as promising therapeutic approaches for the diagnosis of cancer and malignancies.

A Model for the Binding of Fluorescently Labeled Anti-Human CD4 Monoclonal Antibodies to CD4 Receptors on Human Lymphocytes

The CD4 glycoprotein is a component of the T cell receptor complex which plays an important role in the human immune response. This manuscript describes the measurement and modeling of the binding of fluorescently labeled anti-human CD4 monoclonal antibodies (mAb; SK3 clone) to CD4 receptors on the surface of human peripheral blood mononuclear cells (PBMC). CD4 mAb fluorescein isothiocyanate (FITC) and CD4 mAb allophycoerythrin (APC) conjugates were obtained from commercial sources. Four binding conditions were performed, each with the same PBMC sample and different CD4 mAb conjugate. Each binding condition consisted of the PBMC sample incubated for 30 min in labeling solutions containing progressively larger concentrations of the CD4 mAb-label conjugate. After the incubation period, the cells were re-suspended in PBS-based buffer and analyzed using a flow cytometer to measure the mean fluorescence intensity (MFI) of the labeled cell populations. A model was developed to estimate the equilibrium concentration of bound CD4 mAb-label conjugates to CD4 receptors on PBMC. A set of parameters was obtained from the best fit of the model to the measured MFI data and the known number of CD4 receptors on PBMC surface. Divalent and monovalent binding had to be invoked for the APC and FITC CD4 mAb conjugates, respectively. This suggests that the mAb binding depends on the size of the label, which has significant implications for quantitative flow cytometry. The study supports the National Institute of Standards and Technology program to develop quantitative flow cytometry measurements.

Integrating Structural Information to Study the Dynamics of Protein-Protein Interactions in Cells

The information of how two proteins interact is embedded in the atomic details of their binding interfaces. These interactions, spatial-temporally coordinating each other as a network in a variable cytoplasmic environment, dominate almost all biological functions. A feasible and reliable computational model is highly demanded to realistically simulate these cellular processes and unravel the complexities beneath them. We therefore present a multiscale framework that integrates simulations on two different scales. The higher-resolution model incorporates structural information of proteins and energetics of their binding, while the lower-resolution model uses a highly simplified representation of proteins to capture the long-time-scale dynamics of a system with multiple proteins. Through a systematic benchmark test and two practical applications of biomolecular systems with specific cellular functions, we demonstrated that this method could be a powerful approach to understand molecular mechanisms of dynamic interactions between biomolecules and their functional impacts with high computational efficiency.