Solution Ionic Strength Can Modulate Functional Loop Conformations in E. coli Dihydrofolate Reductase

The observation of multiple conformations of a functional loop (termed M20) in the Escherichia coli dihydrofolate reductase (ecDHFR) enzyme triggered the proposition that large-scale motions of protein structural elements contribute to enzyme catalysis. The transition of the M20 loop from a closed conformation to an occluded conformation was thought to aid the rate-limiting release of the products. However, the influence of charged species in the solution environment on the observed M20 loop conformations, independent of charged ligands bound to the enzyme, had not been considered. Molecular dynamics simulations of ecDHFR in model CaCl2 solutions of varying molar ionic strengths IM reveal a substantial free energy barrier between occluded and closed M20 loop states at IM exceeding the E. coli threshold (∼0.24 M). This barrier may facilitate crystallization of ecDHFR in the occluded state, consistent with ecDHFR structures obtained at IM exceeding 0.3 M. At lower IM (≤0.15 M), the M20 loop can explore the occluded state, but prefers an open/partially closed conformation, again consistent with ecDHFR structures. Our findings caution against using ecDHFR structures obtained at nonphysiological ionic strengths in interpreting catalytic events or in structure-based drug design.

Using protein language models for protein interaction hot spot prediction with limited data

BACKGROUND

Protein language models, inspired by the success of large language models in deciphering human language, have emerged as powerful tools for unraveling the intricate code of life inscribed within protein sequences. They have gained significant attention for their promising applications across various areas, including the sequence-based prediction of secondary and tertiary protein structure, the discovery of new functional protein sequences/folds, and the assessment of mutational impact on protein fitness. However, their utility in learning to predict protein residue properties based on scant datasets, such as protein-protein interaction (PPI)-hotspots whose mutations significantly impair PPIs, remained unclear. Here, we explore the feasibility of using protein language-learned representations as features for machine learning to predict PPI-hotspots using a dataset containing 414 experimentally confirmed PPI-hotspots and 504 PPI-nonhot spots.

RESULTS

Our findings showcase the capacity of unsupervised learning with protein language models in capturing critical functional attributes of protein residues derived from the evolutionary information encoded within amino acid sequences. We show that methods relying on protein language models can compete with methods employing sequence and structure-based features to predict PPI-hotspots from the free protein structure. We observed an optimal number of features for model precision, suggesting a balance between information and overfitting.

CONCLUSIONS

This study underscores the potential of transformer-based protein language models to extract critical knowledge from sparse datasets, exemplified here by the challenging realm of predicting PPI-hotspots. These models offer a cost-effective and time-efficient alternative to traditional experimental methods for predicting certain residue properties. However, the challenge of explaining why specific features are important for determining certain residue properties remains.

An Antibody-Drug Conjugate for Multiple Myeloma Prepared by Multi-Arm Linkers

First-line treatment of multiple myeloma, a prevalent blood cancer lacking a cure, using anti-CD38 daratumumab antibody and lenalidomide is often inadequate due to relapse and severe side effects. To enhance drug safety and efficacy, an antibody-drug conjugate, TE-1146, comprising six lenalidomide drug molecules site-specifically conjugated to a reconfigured daratumumab to deliver cytotoxic lenalidomide to tumor cells is developed. TE-1146 is prepared using the HighDAR platform, which employs i) a maleimide-containing "multi-arm linker" to conjugate multiple drug molecules creating a drug bundle, and ii) a designed peptide with a Zn2+ -binding cysteine at the C-termini of a reconfigured daratumumab for site-specific drug bundle conjugation. It is shown that TE-1146 remains intact and effectively enters CD38-expressing tumor cells, releasing lenalidomide, leading to enhanced cell-killing effects compared to lenalidomide/daratumumab alone or their combination. This reveals the remarkable potency of lenalidomide once internalized by myeloma cells. TE-1146 precisely delivers lenalidomide to target CD38-overexpressing tumor cells. In contrast, lenalidomide without daratumumab cannot easily enter cells, whereas daratumumab without lenalidomide relies on Fc-dependent effector functions to kill tumor cells.

System-Specific Parameter Optimization for Nonpolarizable and Polarizable Force Fields

The accuracy of classical force fields (FFs) has been shown to be limited for the simulation of cation-protein systems despite their importance in understanding the processes of life. Improvements can result from optimizing the parameters of classical FFs or by extending the FF formulation by terms describing charge transfer (CT) and polarization (POL) effects. In this work, we introduce our implementation of the CTPOL model in OpenMM, which extends the classical additive FF formula by adding CT and POL. Furthermore, we present an open-source parametrization tool, called FFAFFURR, that enables the (system-specific) parametrization of OPLS-AA and CTPOL models. The performance of our workflow was evaluated by its ability to reproduce quantum chemistry energies and by molecular dynamics simulations of a zinc-finger protein.

Site-specific Conjugation of 6 DOTA Chelators to a CA19-9-targeting scFv-Fc Antibody for Imaging and Therapy

Antibodies conjugated with diagnostic/therapeutic radionuclides are attractive options for inoperable cancers lacking accurate imaging methods and effective therapeutics, such as pancreatic cancer. Hence, we have produced an antibody radionuclide conjugate termed TE-1132 comprising a α-CA19-9 scFv-Fc that is site-specifically conjugated at each C-terminus to 3 DOTA chelators via a cysteine-containing peptide linker. The smaller scFv-Fc size facilitates diffusivity within solid tumors, whereas the chelator-to-antibody ratio of six enabled 177Lu-radiolabeled TE-1132 to exhibit high radioactivity up to 520 MBq/nmol. In mice bearing BxPC3 tumors, immuno-SPECT/CT imaging of [111In]In-TE-1132 and the biodistribution of [177Lu]Lu-TE-1132 showed selective tumor accumulation. Single and multiple doses of [177Lu]Lu-TE-1132 effectively inhibited the BxPC3 tumor growth and prolonged the survival of mice with no irreversible body weight loss or hematopoietic damage. The adequate pharmacokinetic parameters, prominent tumor accumulation, and efficacy with good safety in mice encourage the further investigation of theranostic TE-1132 for treating pancreatic cancer.

Protein Ca2+-Sites Prone to Sr2+ Substitution: Implications for Strontium Therapy

Strontium (Sr), an alkali metal with properties similar to calcium, in the form of soluble salts is used to treat osteoporosis. Despite the information accumulated on the role of Sr²⁺ as a Ca²⁺ mimetic in biology and medicine, there is no systematic study of how the outcome of the competition between the two dications depends on the physicochemical properties of (i) the metal ions, (ii) the first- and second-shell ligands, and (iii) the protein matrix. Specifically, the key features of a Ca²⁺-binding protein that enable Sr²⁺ to displace Ca²⁺ remain unclear. To address this, we studied the competition between Ca²⁺ and Sr²⁺ in protein Ca²⁺-binding sites using density functional theory combined with the polarizable continuum model. Our findings indicate that Ca²⁺-sites with multiple strong charge-donating protein ligands, including one or more bidentately bound Asp^-/Glu^- that are relatively buried and rigid are protected against Sr²⁺ attack. On the other hand, Ca²⁺-sites crowded with multiple protein ligands may be prone to Sr²⁺ displacement if they are solvent-exposed and flexible enough so that an extra backbone ligand from the outer shell can bind to Sr²⁺. In addition, solvent-exposed Ca²⁺ sites with only a few weak charge-donating ligands that can rearrange to fit the strontium's coordination requirements are susceptible to Sr²⁺ displacement. We provide the physical basis of these results and discuss potential novel protein targets of therapeutic Sr²⁺.

Molar-incisor hypomineralisation prevalence in a cohort of Australian children with type 1 diabetes

PURPOSE

Systemic diseases or drugs administered early in life may cause a disruption in amelogenesis and contribute to the qualitative defect of enamel described as molar-incisor hypomineralisation (MIH). Therefore, an increase in prevalence of MIH in children with type 1 diabetes (T1D) may be expected as this systemic disorder is commonly diagnosed in early childhood. The aim of this study was to determine the prevalence of MIH in a cohort of children with T1D and investigate diagnosis of MIH with T1D factors.

METHODS

Cross-sectional study of children with T1D recruited from paediatric diabetes clinics at the Women's and Children's Hospital (South Australia). A detailed medical history, comprehensive dental and MIH examination according to the European Academy of Paediatric Dentistry (EAPD) long form classification was collected for each child. All upper and lower first permanent molars and central incisors were scored.

RESULTS

A total number of 73 participants; 35 (47.95%) males were examined including 584 teeth. The mean age of the participants was 13.25 ± 2.58 years, with a mean age of diagnosis 7.75 ± 3.58 years, and a mean HbA1c of 8.5 ± 1.6%. 42 out of 73 children (54.8%) had enamel defects on at least one of the teeth examined. However, 19.2% met the criteria for MIH. Univariate and bivariate analyses were conducted but no significant associations were noted between MIH and risk factors including diabetes control (p > 0.1).

CONCLUSION

There was a high prevalence of enamel defects and MIH amongst children with T1D. More research is required to establish association between T1D and MIH.

A strategy for evaluating potential antiviral resistance to small molecule drugs and application to SARS-CoV-2

Alterations in viral fitness cannot be inferred from only mutagenesis studies of an isolated viral protein. To-date, no systematic analysis has been performed to identify mutations that improve virus fitness and reduce drug efficacy. We present a generic strategy to evaluate which viral mutations might diminish drug efficacy and applied it to assess how SARS-CoV-2 evolution may affect the efficacy of current approved/candidate small-molecule antivirals for Mpro, PLpro, and RdRp. For each drug target, we determined the drug-interacting virus residues from available structures and the selection pressure of the virus residues from the SARS-CoV-2 genomes. This enabled the identification of promising drug target regions and small-molecule antivirals that the virus can develop resistance. Our strategy of utilizing sequence and structural information from genomic sequence and protein structure databanks can rapidly assess the fitness of any emerging virus variants and can aid antiviral drug design for future pathogens.

Limited liver resections in the posterosuperior segments: international multicentre propensity score-matched and coarsened exact-matched analysis comparing the laparoscopic and robotic approaches

Background

Limited liver resections (LLRs) for tumours located in the posterosuperior segments of the liver are technically demanding procedures. This study compared outcomes of robotic (R) and laparoscopic (L) LLR for tumours located in the posterosuperior liver segments (IV, VII, and VIII).

Methods

This was an international multicentre retrospective analysis of patients who underwent R-LLR or L-LLR at 24 centres between 2010 and 2019. Patient demographics, perioperative parameters, and postoperative outcomes were analysed; 1 : 3 propensity score matching (PSM) and 1 : 1 coarsened exact matching (CEM) were performed.

Results

Of 1566 patients undergoing R-LLR and L-LLR, 983 met the study inclusion criteria. Before matching, 159 R-LLRs and 824 L-LLRs were included. After 1 : 3 PSM of 127 R-LLRs and 381 L-LLRs, comparison of perioperative outcomes showed that median blood loss (100 (i.q.r. 40–200) versus 200 (100–500) ml; P = 0.003), blood loss of at least 500 ml (9 (7.4 per cent) versus 94 (27.6 per cent); P < 0.001), intraoperative blood transfusion rate (4 (3.1 per cent) versus 38 (10.0 per cent); P = 0.025), rate of conversion to open surgery (1 (0.8 per cent) versus 30 (7.9 per cent); P = 0.022), median duration of Pringle manoeuvre when applied (30 (20–46) versus 40 (25–58) min; P = 0.012), and median duration of operation (175 (130–255) versus 224 (155–300); P < 0.001) were lower in the R-LLR group compared with the L-LLR group. After 1 : 1 CEM of 104 R-LLRs with 104 L-LLRs, R-LLR was similarly associated with significantly reduced blood loss and a lower rate of conversion to open surgery.

Conclusion

Based on a matched analysis of well selected patients, both robotic and laparoscopic access could be undertaken safely with good outcomes for tumours in the posterosuperior liver segments.

Hepatic transcriptome profiling according to growth rate reveals acclimation in metabolic regulatory mechanisms to cyclic heat stress in broiler chickens

Climate change has numerous effects on poultry that result in welfare concerns and economic losses in agricultural industries. However, the mechanisms underlying the acclimation to heat stress in poultry have not been comprehensively defined. Therefore, identifying associated patterns of gene regulation and understanding the molecular mechanisms of acclimation to a warmer environment will provide insights into the acclimation system of broiler chickens. We profiled differentially expressed genes (DEGs) associated with differences in growth performance under heat stress conditions in the liver tissues of broilers based on RNA sequencing data. The DEGs were identified by comparison to the gene expression levels of broilers exhibiting average growth at 28 d of age (D28A) and D36A relative to those at D21A. In D36A, 507 and 312 DEGs were up- and downregulated, respectively, whereas 400 and 156 DEGs were up- and downregulated in D28A, respectively. Pathway enrichment analysis further revealed that “fatty acid degradation” and “heat shock protein expression” were upregulated in broilers exhibiting a higher growth and weight, whereas “cell cycle arrest” and “amino acid metabolism” were downregulated. Transcriptome profiling revealed that the acclimatized group supplied fat and energy from the liver to tissues through the breakdown of fatty acids. Furthermore, homeostasis was maintained via heat shock proteins and antioxidant enzymes. The characterized candidate genes and mechanisms associated with the response to heat stress might serve as a foundation for improving the ability of broilers to acclimatize under heat stress conditions.

Efficient Strategy to Design Protease Inhibitors: Application to Enterovirus 71 2A Protease

One strategy to counter viruses that persistently cause outbreaks is to design molecules that can specifically inhibit an essential multifunctional viral protease. Herein, we present such a strategy using well-established methods to first identify a region present only in viral (but not human) proteases and find peptides that can bind specifically to this "unique" region by maximizing the protease-peptide binding free energy iteratively using single-point mutations starting with the substrate peptide. We applied this strategy to discover pseudosubstrate peptide inhibitors for the multifunctional 2A protease of enterovirus 71 (EV71), a key causative pathogen for hand-foot-and-mouth disease affecting young children, along with coxsackievirus A16. Four peptide candidates predicted to bind EV71 2A protease more tightly than the natural substrate were experimentally validated and found to inhibit protease activity. Furthermore, the crystal structure of the best pseudosubstrate peptide bound to the EV71 2A protease was determined to provide a molecular basis for the observed inhibition. Since the 2A proteases of EV71 and coxsackievirus A16 share nearly identical sequences and structures, our pseudosubstrate peptide inhibitor may prove useful in inhibiting the two key pathogens of hand-foot-and-mouth disease.

Factors allowing small monovalent Li+ to displace Ca2+ in proteins

Because Li⁺ and Ca²⁺ differ in both charge and size, the possibility that monovalent Li⁺ could dislodge the bulkier, divalent Ca²⁺ in Ca²⁺ proteins had not been considered. However, our recent density functional theory/continuum dielectric calculations predicted that Li⁺ could displace the native Ca²⁺ from the C2 domain of cytosolic PKCα/γ. This would reduce electrostatic interactions between the Li⁺-bound C2 domain and the membrane, consistent with experimental studies showing that Li⁺ can inhibit the translocation of cytoplasmic PKC to membranes. Besides the trinuclear Ca²⁺-site in the PKCα/γ C2 domain, it is not known whether other Ca²⁺-sites in human proteins may be susceptible to Li⁺ substitution. Furthermore, it is unclear what factors determine the outcome of the competition between divalent Ca²⁺ and monovalent Li⁺. Here we show that the net charge of residues in the first and second coordination shell is a key determinant of the selectivity for divalent Ca²⁺ over monovalent Li⁺ in proteins: neutral/anionic Ca²⁺-carboxylate sites are protected against Li⁺ attack. They are further protected by outer-shell Asp^-/Glu^- and the protein matrix rigidifying the Ca²⁺-site or limiting water entry. In contrast, buried, cationic Ca²⁺-sites surrounded by Arg⁺/Lys⁺, which are found in the C2 domains of PKCα/γ, as well as certain synaptotagmins, are prone to Li⁺ attack.

Allosteric coupling between transmembrane segment 4 and the selectivity filter of TALK1 potassium channels regulates their gating by extracellular pH

Opening of two-pore domain K+ channels (K2Ps) is regulated by various external cues, such as pH, membrane tension, or temperature, which allosterically modulate the selectivity filter (SF) gate. However, how these cues cause conformational changes in the SF of some K2P channels remains unclear. Herein, we investigate the mechanisms by which extracellular pH affects gating in an alkaline-activated K2P channel, TALK1, using electrophysiology and molecular dynamics (MD) simulations. We show that R233, located at the N-terminal end of transmembrane segment 4, is the primary pHo sensor. This residue distally regulates the orientation of the carbonyl group at the S1 potassium-binding site through an interacting network composed of residues on transmembrane segment 4, the pore helix domain 1, and the SF. Moreover, in the presence of divalent cations, we found the acidic pH-activated R233E mutant recapitulates the network interactions of protonated R233. Intriguingly, our data further suggested stochastic coupling between R233 and the SF gate, which can be described by an allosteric gating model. We propose that this allosteric model could predict the hybrid pH sensitivity in heterodimeric channels with alkaline-activated and acidic-activated K2P subunits.

From quantum-derived principles underlying cysteine reactivity to combating the COVID-19 pandemic

The COVID‐19 pandemic poses a challenge in coming up with quick and effective means to counter its cause, the SARS‐CoV‐2. Here, we show how the key factors governing cysteine reactivity in proteins derived from combined quantum mechanical/continuum calculations led to a novel multi‐targeting strategy against SARS‐CoV‐2, in contrast to developing potent drugs/vaccines against a single viral target such as the spike protein. Specifically, they led to the discovery of reactive cysteines in evolutionary conserved Zn²⁺‐sites in several SARS‐CoV‐2 proteins that are crucial for viral polypeptide proteolysis as well as viral RNA synthesis, proofreading, and modification. These conserved, reactive cysteines, both free and Zn²⁺‐bound, can be targeted using the same Zn‐ejector drug (disulfiram/ebselen), which enables the use of broad‐spectrum anti‐virals that would otherwise be removed by the virus's proofreading mechanism. Our strategy of targeting multiple, conserved viral proteins that operate at different stages of the virus life cycle using a Zn‐ejector drug combined with other broad‐spectrum anti‐viral drug(s) could enhance the barrier to drug resistance and antiviral effects, as compared to each drug alone. Since these functionally important nonstructural proteins containing reactive cysteines are highly conserved among coronaviruses, our proposed strategy has the potential to tackle future coronaviruses.

This article is categorized under:

Structure and Mechanism > Reaction Mechanisms and Catalysis

Structure and Mechanism > Computational Biochemistry and Biophysics

Electronic Structure Theory > Density Functional Theory

Probiotic characterization of Lactiplantibacillus plantarum HOM3204 and its restoration effect on antibiotic-induced dysbiosis in mice

The purpose of this study was to evaluate the probiotic characteristics of Lactiplantibacillus plantarum HOM3204 isolated from homemade pickled cabbage and to examine its restoration effect on antibiotic-induced dysbiosis in mice. Lact. plantarum HOM3204 tolerated simulated gastric and intestinal juices with a 99.38% survival rate. It also showed strong adhesion ability (3.45%) to Caco-2 cells and excellent antimicrobial activity against foodborne pathogens in vitro. For safety (antibiotic susceptibility) of this strain, it was susceptible to all the tested seven antibiotics. Lact. plantarum HOM3204 had good stability during storage, especially in cold and frozen conditions. Furthermore, Lact. plantarum HOM3204 significantly restored the gut microbiota composition by increasing the abundance of Lactobacilli and Bifidobacteria and decreasing Enterococci, and improved antioxidative function by raising the concentrations of glutathione peroxidase (GSH-Px) and superoxide dismutase (SOD) in serum of antibiotic-induced dysbiosis in mice. These results suggest that Lact. plantarum HOM3204 could be a potential probiotic as a functional food ingredient.

PPI-HotspotDB: Database of Protein-Protein Interaction Hot Spots

Single-point mutations of certain residues (so-called hot spots) impair/disrupt protein-protein interactions (PPIs), leading to pathogenesis and drug resistance. Conventionally, a PPI-hot spot is identified when its replacement decreased the binding free energy significantly, generally by ≥2 kcal/mol. The relatively few mutations with such a significant binding free energy drop limited the number of distinct PPI-hot spots. By defining PPI-hot spots based on mutations that have been manually curated in UniProtKB to significantly impair/disrupt PPIs in addition to binding free energy changes, we have greatly expanded the number of distinct PPI-hot spots by an order of magnitude. These experimentally determined PPI-hot spots along with available structures have been collected in a database called PPI-Hotspot^DB. We have applied the PPI-Hotspot^DB to create a nonredundant benchmark, PPI-Hotspot+PDB^BM, for assessing methods to predict PPI-hot spots using the free structure as input. PPI-Hotspot^DB will benefit the design of mutagenesis experiments and development of PPI-hot spot prediction methods. The database and benchmark are freely available at https://ppihotspot.limlab.dnsalias.org.

1104 Plasma Circulating Cell-Free Messenger RNA as a Potential Biomarker of Melanoma

Introduction

Blood-borne cell-free nucleic acids are increasingly emerging as significant non-invasive adjuncts to current methods of disease status evaluation in cancer patients. In this study, we sought to examine whether significant differences exist in the plasma transcriptomic profile of advanced melanoma patients with a high disease burden compared to patients with therapeutic response.

Method

Plasma circulating cell-free messenger ribonucleic acid (ccfmRNA) was extracted from twenty patients with stage IV melanoma receiving immunotherapy. Pathway focused gene expression analysis was performed. Patients were assessed with paired blood sample collection and CT scan assessments at baseline and at 3 months follow up. Response Evaluation Criteria in Solid Tumours (RECIST 1.1) was used for tumour burden estimation.

Results

In stage IV melanoma patients, CCL5, GZMB and MYD88 genes were significantly over-expressed (p < 0.05 versus healthy controls). In patients with high disease burden or progressive disease, CCL18, CCR1, CCR4, CD274, CSF2, and GBP1 genes were significantly over-expressed (p < 0.05 versus patients with therapeutic response). Finally, in stage IV melanoma patients with brain metastases, CCL18, CCR1, CCR4, CD274, CSF2, EGF, GBP1, and PTGS2 genes were significantly over-expressed (p < 0.05, versus patients without brain metastasis).

Conclusions

Significant differences were observed in the plasma transcriptomic profile between the various melanoma patient groups, and we postulate that these differences may be exploited to identify novel therapeutic targets or biomarkers relevant to melanoma. CCL4 and CCL5 are prognostic in melanoma, both genes had significantly higher expressions in low disease burden patients compared to patients with a high disease burden.

1095 The Utility of Plasma Circulating Cell-Free Messenger RNA as A Biomarker of Glioma: A Pilot Study

Introduction

Research into the potential utility of plasma-derived circulating-cell-free nucleic acids as non-invasive adjuncts to radiological imaging has been occasioned by the invasive nature of brain tumour biopsy. Circulating-cell-free messenger RNAs are short fragments of RNA present in blood. The objective of this study was to determine whether significant differences exist in the plasma transcriptomic profile of glioma patients relative to differences in their tumour characteristics, and also whether any observed differences were representative of synchronously obtained glioma samples and TCGA glioma derived RNA.

Method

Blood samples were collected from twenty-nine patients prior to tumour resection. Plasma-ccfmRNA and glioma derived RNA were extracted and profiled.

Results

BCL2L1, CXCL5, GZMB, HLA-A, HLA-C, IRF1, MYD88, TGFB1, TLR2, and TP53 genes were significantly over-expressed in glioma (high-grade-glioma-HGG and low-grade-glioma-LGG) patients (p < 0.05, versus control). BCL2L1, GZMB and HLA-A genes were significantly over-expressed in HGG patients (p < 0.05, versus LGG patients). There was positive correlation between the magnitude of fold change of differentially expressed genes in plasma and glioma derived RNA (Spearman r = 0.6344, n = 14, p = 0.017), and with the mean FPKM of TCGA glioma derived RNA samples (Spearman r = 0.4614, n = 19, p = 0.047). There was positive correlation between glioma radiographic tumour burden and the magnitude of fold change of CSF3 gene (r = 0.9813, n = 20, p < 0.001).

Conclusions

We identified significant differential expression of genes involved in cancer inflammation and immunity among patients with different glioma grades, and we identified positive correlation between the plasma transcriptomic profile and tumour samples, and with TCGA glioma derived RNA.

Calcium in Signaling: Its Specificity and Vulnerabilities toward Biogenic and Abiogenic Metal Ions

Divalent calcium ion (Ca²⁺) plays an indispensable role as a second messenger in a myriad of signal transduction processes. Of utmost importance for the faultless functioning of calcium-modulated signaling proteins is their binding selectivity of the native metal cation over rival biogenic/abiogenic metal ion contenders in the intra/extracellular fluids. In this Perspective, we summarize recent findings on the competition between the cognate Ca²⁺ and other biogenic or abiogenic divalent cations for binding to Ca²⁺-signaling proteins or organic cofactors. We describe the competition between the two most abundant intracellular biogenic metal ions (Mg²⁺ and Ca²⁺) for Ca²⁺-binding sites in signaling proteins, followed by the rivalry between native Ca²⁺ and "therapeutic" Li⁺ as well as "toxic" Pb²⁺. We delineate the key factors governing the rivalry between the native and non-native cations in proteins and highlight key implications for the biological performance of the respective proteins/organic cofactors.

Trinuclear Calcium Site in the C2 Domain of PKCα/γ Is Prone to Lithium Attack

Lithium (Li⁺) is the first-line therapy for bipolar disorder and a candidate drug for various diseases such as amyotrophic lateral sclerosis, multiple sclerosis, and stroke. Despite being the captivating subject of many studies, the mechanism of lithium's therapeutic action remains unclear. To date, it has been shown that Li⁺ competes with Mg²⁺ and Na⁺ to normalize the activity of inositol and neurotransmitter-related signaling proteins, respectively. Furthermore, Li⁺ may co-bind with Mg²⁺-loaded adenosine or guanosine triphosphate to alter the complex's susceptibility to hydrolysis and mediate cellular signaling. Bipolar disorder patients exhibit abnormally high cytosolic Ca²⁺ levels and protein kinase C (PKC) hyperactivity that can be downregulated by long-term Li⁺ treatment. However, the possibility that monovalent Li⁺ could displace the bulkier divalent Ca²⁺ and inhibit PKC activity has not been considered. Here, using density functional theory calculations combined with continuum dielectric methods, we show that Li⁺ may displace the native dication from the positively charged trinuclear site in the C2 domain of cytosolic PKCα/γ. This would affect the membrane-docking ability of cytosolic PKCα/γ and reduce the abnormally high membrane-associated active PKCα/γ levels, thus downregulating the PKC hyperactivity found in bipolar patients.

Influence of solution ionic strength on the stabilities of M20 loop conformations in apo E. coli dihydrofolate reductase

Interactions among ions and their specific interactions with macromolecular solutes are known to play a central role in biomolecular stability. However, similar effects in the conformational stability of protein loops that play functional roles, such as binding ligands, proteins, and DNA/RNA molecules, remain relatively unexplored. A well-characterized enzyme that has such a functional loop is Escherichia coli dihydrofolate reductase (ecDHFR), whose so-called M20 loop has been observed in three ordered conformations in crystal structures. To explore how solution ionic strengths may affect the M20 loop conformation, we proposed a reaction coordinate that could quantitatively describe the loop conformation and used it to classify the loop conformations in representative ecDHFR x-ray structures crystallized in varying ionic strengths. The Protein Data Bank survey indicates that at ionic strengths (I) below the intracellular ion concentration-derived ionic strength in E. coli (I ≤ 0.237M), the ecDHFR M20 loop tends to adopt open/closed conformations, and rarely an occluded loop state, but when I is >0.237M, the loop tends to adopt closed/occluded conformations. Distance-dependent electrostatic potentials around the most mobile M20 loop region from molecular dynamics simulations of ecDHFR in equilibrated CaCl₂ solutions of varying ionic strengths show that high ionic strengths (I = 0.75/1.5M) can preferentially stabilize the loop in closed/occluded conformations. These results nicely correlate with conformations derived from ecDHFR structures crystallized in varying ionic strengths. Altogether, our results suggest caution in linking M20 loop conformations derived from crystal structures solved at ionic strengths beyond that tolerated by E. coli to the ecDHFR function.

Correction: Multi-targeting of functional cysteines in multiple conserved SARS-CoV-2 domains by clinically safe Zn-ejectors

[This corrects the article DOI: 10.1039/D0SC02646H.].

Synergistic Inhibition of SARS-CoV-2 Replication Using Disulfiram/Ebselen and Remdesivir

The SARS-CoV-2 replication and transcription complex (RTC) comprising nonstructural protein (nsp) 2-16 plays crucial roles in viral replication, reducing the efficacy of broad-spectrum nucleoside analog drugs such as remdesivir and evading innate immune responses. Most studies target a specific viral component of the RTC such as the main protease or the RNA-dependent RNA polymerase. In contrast, our strategy is to target multiple conserved domains of the RTC to prevent SARS-CoV-2 genome replication and to create a high barrier to viral resistance and/or evasion of antiviral drugs. We show that the clinically safe Zn-ejector drugs disulfiram and ebselen can target conserved Zn2+ sites in SARS-CoV-2 nsp13 and nsp14 and inhibit nsp13 ATPase and nsp14 exoribonuclease activities. As the SARS-CoV-2 nsp14 domain targeted by disulfiram/ebselen is involved in RNA fidelity control, our strategy allows coupling of the Zn-ejector drug with a broad-spectrum nucleoside analog that would otherwise be excised by the nsp14 proofreading domain. As proof-of-concept, we show that disulfiram/ebselen, when combined with remdesivir, can synergistically inhibit SARS-CoV-2 replication in Vero E6 cells. We present a mechanism of action and the advantages of our multitargeting strategy, which can be applied to any type of coronavirus with conserved Zn2+ sites.

Metal Affinity/Selectivity of Monophosphate-Containing Signaling/Lipid Molecules

Monophosphate, an essential component of nucleic acids, as well as cell membranes and signaling molecules, is often bound to metal cations. Despite the biological importance of monophosphate-containing cell-signaling or lipid molecules, their propensity to bind the two most abundant cellular dications, Mg²⁺ and Ca²⁺, in a particular mode (inner/outer shell, mono/bidentate) is not well understood. Whether they prefer binding to Mg²⁺ than to Ca²⁺ and if they can outcompete the carboxylates of excitatory Asp/Glu and inhibitory gamma-aminobutyric acid (GABA) neurotransmitters in binding to Mg²⁺/Ca²⁺ remain unclear. To address these questions, we modeled cyclic adenosine/guanosine monophosphate (cAMP/cGMP), nucleoside 2',3'-cyclic phosphate, phosphatidylinositol (PI), phosphatidylserine (PS), and phosphatidylethanolamine (PEA) and determined their most stable metal-binding modes, including those of Asp/Glu and GABA, as well as their selectivity for Mg²⁺/Ca²⁺ using density functional theory combined with the polarizable continuum model. The results obtained, which are consistent with the available experimental findings, reveal that the structurally and functionally diverse monophosphate-containing ligands studied prefer monodentate coordination of Mg²⁺ because of the greater strain encountered upon bidentate coordination, whereas the larger Ca²⁺ imposes less strain upon bidentate binding and has reduced/no preference for monodentate coordination. We further show that in a low-dielectric environment, negatively charged monophosphate-containing ligands favor the better charge-accepting dication, that is, Mg²⁺ rather than Ca²⁺. By promoting Mg²⁺ over Ca²⁺ binding, signaling monophosphates (cAMP/cGMP) do not entrap cellular Ca²⁺ and interfere with signal transduction processes employing Ca²⁺ as a second messenger. In regions with high glutamate cytoplasmic concentration, glutamate may sequester Mg²⁺ bound to isolated five-/six-membered ring phosphates, PI, or neutral PEA, but not anionic phospholipids constituting the inner leaflet of the cell membrane.

A case study exploring the 'real world' process of 'naturalizing' school playgrounds

Background: The aim of this study was to examine the lived experiences of teachers and administrators at two schools (one elementary Kindergarten-Grade 5 and one middle school, grades 6-8) following the 'naturalizing' of a portion of their schoolyard. Methods: A qualitative case study design was used. Focus groups and interviews allowed participants to express their rationale, identify implementation facilitators and barriers and critical processes/steps for realizing their naturalization goals. Researcher questions explored (1) factors which led to naturalizing the schoolyard, (2) key planning and process steps (3) challenges and successes experienced, and, (4) lessons learned. Findings: Six themes were identified as rationale for naturalizing including: re-integration of outdoor play into children's daily lives, pro-active mental health, opportunity for inclusion, nature as a positive space for learning, it fit with our values and principles, we had a need. Implementation was facilitated by having champions and support from multiple levels, adopting a phased approach, resources and having the buy-in to persevere over time. Significant challenges and 'how to' steps were also identified. Conclusions: Naturalizing schoolyards does not happen without a supportive context, significant time and effort. Implementation steps and issues were those commonly found in other school-based health promotion efforts.

Sensitive and Specific Cadmium Biosensor Developed by Reconfiguring Metal Transport and Leveraging Natural Gene Repositories

Whole-cell biosensors are useful for monitoring heavy metal toxicity in public health and ecosystems, but their development has been hindered by intrinsic trade-offs between sensitivity and specificity. Here, we demonstrated an effective engineering solution by building a sensitive, specific, and high-response biosensor for carcinogenic cadmium ions. We genetically programmed the metal transport system of Escherichia coli to enrich intracellular cadmium ions and deprive interfering metal species. We then selected 16 cadmium-sensing transcription factors from the GenBank database and tested their reactivity to 14 metal ions in the engineered E. coli using the expression of the green fluorescent protein as the readout. The resulting cadmium biosensor was highly specific and showed a detection limit of 3 nM, a linear increase in fluorescent intensities from 0 to 200 nM, and a maximal 777-fold signal change. Using this whole-cell biosensor, a smartphone, and low-tech equipment, we developed a simple assay capable of measuring cadmium ions at the same concentration range in irrigation water and human urine. This method is user-friendly and cost-effective, making it affordable to screen large amounts of samples for cadmium toxicity in agriculture and medicine. Moreover, our work highlights natural gene repositories as a treasure chest for bioengineering.

The novel Mechanical Ventilator Milano for the COVID-19 pandemic

This paper presents the Mechanical Ventilator Milano (MVM), a novel intensive therapy mechanical ventilator designed for rapid, large-scale, low-cost production for the COVID-19 pandemic. Free of moving mechanical parts and requiring only a source of compressed oxygen and medical air to operate, the MVM is designed to support the long-term invasive ventilation often required for COVID-19 patients and operates in pressure-regulated ventilation modes, which minimize the risk of furthering lung trauma. The MVM was extensively tested against ISO standards in the laboratory using a breathing simulator, with good agreement between input and measured breathing parameters and performing correctly in response to fault conditions and stability tests. The MVM has obtained Emergency Use Authorization by U.S. Food and Drug Administration (FDA) for use in healthcare settings during the COVID-19 pandemic and Health Canada Medical Device Authorization for Importation or Sale, under Interim Order for Use in Relation to COVID-19. Following these certifications, mass production is ongoing and distribution is under way in several countries. The MVM was designed, tested, prepared for certification, and mass produced in the space of a few months by a unique collaboration of respiratory healthcare professionals and experimental physicists, working with industrial partners, and is an excellent ventilator candidate for this pandemic anywhere in the world.

Cardiovascular phenotyping of the first mouse model of Sarcoidosis

Introduction

Sarcoidosis is a potentially life-threatening, inflammatory, granulomatous disease that affects multiple organs including the heart. Heretofore, its unknown etiology had hindered the creation of experimental models and the understanding of the molecular mechanisms of pathogenesis behind it.

Purpose

To extensively phenotype the heart of the first mouse model of sarcoidosis created through deletion of the tuberous sclerosis 2 (Tsc2) gene in the CD11c-positive macrophage population.

Methods

Tsc2 fl/fl CD11c Cre+ (Tsc2-KO; n=7) and Tsc2 fl/fl CD11c Cre- (Tsc2-WT; n=7) mice were subjected to echocardiography at 25 weeks of age (woa) to assess myocardial dimensions and function. Hearts of 13 and 25woa animals were subjected to histological and immunological stains to assess tissue changes, subtype inflammatory infiltrates and examine the localization of key proteins shown to be re-distributed in patients.

Results

At 13 woa, Tsc2-KO animals show inflammatory infiltrates; subtyped mainly as macrophages as well as evidence of myocyte destruction. At 25 woa, the number of inflammatory cells is significantly higher and there is heavy fibrotic replacement primarily in the septum and trabeculae. Older animals also show giant cells and non-necrotizing granulomas. The hearts show heterogeneous gap junction remodeling known to constitute an arrhythmogenic substrate and lack of immunoreactive signal for the desmosomal protein plakoglobin from the cell-cell junctions just as described in patients. The left ventricular ejection fraction and LV morphology was not significantly different between the two groups (EF: 64±4% in Tsc2-KO vs 64±2% in Tsc2-WT; LV end-systolic diameter: 4.51±0.54 mm in Tsc2-KO vs 4.59±0.29 mm in Tsc2-WT). However, there was a strong trend towards increasing filling pressure (E/e'ratio; 14.24±4.01 vs 12.15±2.54) and mean pulmonary pressure (21±6 vs 18±3 mmHg) in Tsc2-KO mice compared to controls suggesting diastolic dysfunction.

Conclusion

Hearts of the Tsc2 fl/fl CD11c Cre+ animals show a phenotype highly reminiscent of cardiac sarcoidosis in patients. We anticipate that this model will be very useful in deciphering molecular mechanisms of pathogenesis as well as testing much-needed mechanism-based therapies.

Funding Acknowledgement

Type of funding source: Foundation. Main funding source(s): British Heart Foundation - PG/18/27/33616

Multi-targeting of functional cysteines in multiple conserved SARS-CoV-2 domains by clinically safe Zn-ejectors

We present a near-term treatment strategy to tackle pandemic outbreaks of coronaviruses with no specific drugs/vaccines by combining evolutionary and physical principles to identify conserved viral domains containing druggable Zn-sites that can be targeted by clinically safe Zn-ejecting compounds. By applying this strategy to SARS-CoV-2 polyprotein-1ab, we predicted multiple labile Zn-sites in papain-like cysteine protease (PLpro), nsp10 transcription factor, and nsp13 helicase. These are attractive drug targets because they are highly conserved among coronaviruses and play vital structural/catalytic roles in viral proteins indispensable for virus replication. We show that five Zn-ejectors can release Zn2+ from PLpro and nsp10, and clinically-safe disulfiram and ebselen can not only covalently bind to the Zn-bound cysteines in both proteins, but also inhibit PLpro protease. We propose combining disulfiram/ebselen with broad-spectrum antivirals/drugs to target different conserved domains acting at various stages of the virus life cycle to synergistically inhibit SARS-CoV-2 replication and reduce the emergence of drug resistance.

How the Local Environment of Functional Sites Regulates Protein Function

Proteins form complex biological machineries whose functions in the cell are highly regulated at both the cellular and molecular levels. Cellular regulation of protein functions involves differential gene expressions, post-translation modifications, and signaling cascades. Molecular regulation, on the other hand, involves tuning an optimal local protein environment for the functional site. Precisely how a protein achieves such an optimal environment around a given functional site is not well understood. Herein, by surveying the literature, we first summarize the various reported strategies used by certain proteins to ensure their correct functioning. We then formulate three key physicochemical factors for regulating a protein's functional site, namely, (i) its immediate interactions, (ii) its solvent accessibility, and (iii) its conformational flexibility. We illustrate how these factors are applied to regulate the functions of free/metal-bound Cys and Zn sites in proteins.

Sensitivity of Functional Loop Conformations on Long-Range Electrostatics: Implications for M20 Loop Dynamics in E. coli Dihydrofolate Reductase

In E. coli dihydrofolate reductase, unusual conformational motions of a functional M20 loop that interacts with substrate and coenzyme have been construed as evidence for dynamical effects in enzyme catalysis. By computing this loop's conformational free energies in the apoenzyme, we show that it is sensitive to the treatment of long-range electrostatic interactions and the solvation box size in modeling/simulations. These results provide important guidelines for computing reaction/binding free energy profiles of proteins with functional loops.

Gene Silencing Mechanisms Revealed by Dynamics of Guide, Target, and Duplex Binding to Argonaute

Argonaute (Ago) protein plays a central role in silencing gene expression by binding a "guide" strand to the base-pair with a complementary mRNA and degrading the mRNA. The current understanding of how Ago-guide and Ago-guide-mRNA complexes assemble is based mainly on static crystal structures; the associated kinetic pathways remain unknown/unclear. By simulating the successive binding of guide/target strand to Thermus thermophilus Ago (TtAgo) and computing the respective free energy landscapes, we directly visualize how TtAgo silencing complexes form and function. We show that the guide binding rate depends on its initial loading position onto TtAgo. Subsequent target recognition beyond the scissile 10-11 nucleotides must overcome a substantial energy barrier for TtAgo's nucleotide-binding groove to expand widely. This work reveals novel roles for the core TtAgo domains and shows how the kinetic barriers that must be overcome for critical structural changes to occur lead to target repression/cleavage.

Why Cellular Di/Triphosphates Preferably Bind Mg2+ and Not Ca2

Di/triphosphates perform a multitude of essential tasks, being important components of many vital organic cofactors such as adenosine/guanosine di/triphosphate (ADP/GDP, ATP/GTP), flavin adenine dinucleotide, and nicotinamide adenine dinucleotide and its phosphate derivative. They are generally bound to cations inside cells, in particular Mg²⁺ in the case of ATP/GTP. Yet how their metal-binding modes depend on the number, charge, and solvent exposure of the polyphosphate group and how Mg²⁺and Ca²⁺ dications that coexist in cellular fluids compete for di/triphosphates in biological systems remain elusive. Using density functional theory calculations combined with a polarizable continuum model, we have determined the relative free energies and stabilities of the different binding modes of di- and triphosphate groups to Mg²⁺ and Ca²⁺. We show that the thermodynamic outcome of the competition between Mg²⁺ and Ca²⁺ for cellular di/triphosphates depends mainly on the oligomericity/charge and metal-binding mode of the phosphate ligand as well as the solvent exposure of the binding site. Increasing the charge and thus denticity of the phosphate ligand from bi- to tridentate in a buried binding pocket enhances the affinity of the host system for the stronger charge acceptor, Mg²⁺. The cellular di/triphosphates's intrinsic properties and the protein matrix allowing them to bind a dication bi/tridentately, along with the higher cytosolic concentration of Mg²⁺ compared to Ca²⁺, enables Mg²⁺ to outcompete Ca²⁺ in binding to these highly charged anions. This suggests an explanation for why nature has chosen Mg²⁺ but not Ca²⁺ to perform most of the essential tasks associated with biological triphosphates.

Free and Bound Therapeutic Lithium in Brain Signaling

Lithium, a first-line therapy for bipolar disorder, is effective in preventing suicide and new depressive/manic episodes. Yet, how this beguilingly simple monocation with only two electrons could yield such profound therapeutic effects remains unclear. An in-depth understanding of lithium's mechanisms of actions would help one to develop better treatments limiting its adverse side effects and repurpose lithium for treating traumatic brain injury and chronic neurodegenerative diseases. In this Account, we begin with a comparison of the physicochemical properties of Li⁺ and its key native rivals, Na⁺ and Mg²⁺, to provide physical grounds for their competition in protein binding sites. Next, we review the abnormal signaling pathways and proteins found in bipolar patients, who generally have abnormally high intracellular Na⁺ and Ca²⁺ concentrations, high G-protein levels, and hyperactive phosphatidylinositol signaling and glycogen synthase kinase-3β (GSK3β) activity. We briefly summarize experimental findings on how lithium, at therapeutic doses, modulates these abnormal signaling pathways and proteins. Following this survey, we address the following aspects of lithium's therapeutic actions: (1) Can Li⁺ displace Na⁺ from the allosteric Na⁺-binding sites in neurotransmitter transporters and G-protein coupled receptors (GPCRs); if so, how would this affect the host protein's function? (2) Why are certain Mg²⁺-dependent enzymes targeted by Li⁺? (3) How does Li⁺ binding to Mg²⁺-bound ATP/GTP (denoted as NTP) in solution affect the cofactor's conformation and subsequent recognition by the host protein? (4) How do NTP-Mg-Li complexes modulate the properties of the respective cellular receptors and signal-transducing proteins? We show that Li⁺ may displace Na⁺ from allosteric Na⁺-binding sites in certain GPCRs and stabilize inactive conformations, preventing these receptors from relaying signal to the respective G-proteins. It may also displace Mg²⁺ in enzymes containing highly cationic Mg²⁺-binding sites such as GSK3β, but not in enzymes containing Mg²⁺-binding sites with low or zero charge. We further show that Li⁺ binding to Mg²⁺-NTP in water does not alter the NTP conformation, which is locked by all three phosphates binding to Mg²⁺. However, bound lithium in the form of [NTP-Mg-Li]^2- dianions can activate or inhibit the host protein depending on the NTP-binding pocket's shape, which determines the metal-binding mode: The ATP-binding pocket's shape in the P2X receptor is complementary to the native ATP-Mg solution conformation and nicely fits [ATP-Mg-Li]^2-. However, since the ATP βγ phosphates bind Li⁺, bimetallic [ATP-Mg-Li]^2- may be more resistant to hydrolysis than the native cofactor, enabling ATP to reside longer in the binding site and elicit a prolonged P2X response. In contrast, the elongated GTP-binding pockets in G-proteins allow only two GTP phosphates to bind Mg²⁺, so the GTP conformation is no longer "triply-locked". Consequently, Li⁺ binding to GTP-Mg can significantly alter the native cofactor's structure, lowering the activated G-protein level, thus attenuating hyperactive G-protein-mediated signaling in bipolar patients. In summary, we have presented a larger "connected" picture of lithium's diverse effects based on its competition as a free monocation with native cations or as a phosphate-bound polyanionic complex modulating the host protein function.