Cation-π interactions in protein-ligand binding: theory and data-mining reveal different roles for lysine and arginine

Controlling Stereoselectivity with Noncovalent Interactions in Chiral Phosphoric Acid Organocatalysis

Chiral phosphoric acids (CPAs) have emerged as highly effective Brønsted acid catalysts in an expanding range of asymmetric transformations, often through novel multifunctional substrate activation modes. Versatile and broadly appealing, these catalysts benefit from modular and tunable structures, and compatibility with additives. Given the unique types of noncovalent interactions (NCIs) that can be established between CPAs and various reactants─such as hydrogen bonding, aromatic interactions, and van der Waals forces─it is unsurprising that these catalyst systems have become a promising approach for accessing diverse chiral product outcomes. This review aims to provide an in-depth exploration of the mechanisms by which CPAs impart stereoselectivity, positioning NCIs as the central feature that connects a broad spectrum of catalytic reactions. Spanning literature from 2004 to 2024, it covers nucleophilic additions, radical transformations, and atroposelective bond formations, highlighting the applicability of CPA organocatalysis. Special emphasis is placed on the structural and mechanistic features that govern CPA-substrate interactions, as well as the tools and techniques developed to enhance our understanding of their catalytic behavior. In addition to emphasizing mechanistic details and stereocontrolling elements in individual reactions, we have carefully structured this review to provide a natural progression from these specifics to a broader, class-level perspective. Overall, these findings underscore the critical role of NCIs in CPA catalysis and their significant contributions to advancing asymmetric synthesis.

Dynamically chiral phosphonic acid-type metallo-β-lactamase inhibitors

Antibiotic resistance is a growing global health threat that risks the lives of millions. Among the resistance mechanisms, that mediated by metallo-β-lactamases is of particular concern as these bacterial enzymes dismantle most β-lactam antibiotics, which are our widest applied and cheapest to produce antibiotic agents. So far, no clinically applicable metallo-β-lactamase inhibitors are available. Aiming to adapt to structural variations, we introduce the inhibitor concept: dynamically chiral phosphonic acids. We demonstrate that they are straightforward to synthesize, penetrate bacterial membranes, inhibit the metallo-β-lactamase enzymes NDM-1, VIM-2 and GIM-1, and are non-toxic to human cells. Mimicking the transition state of β-lactam hydrolysis, they target the Zn ions of the metallo-β-lactamase active site. As a unique feature, both of their stereoisomers bind metallo-β-lactamases, which provides them unparalleled adaptability to the structural diversity of these enzymes, and may allow them to hamper bacteria's ability for resistance development.

Tetrahydrobenzimidazole TMQ0153 targets OPA1 and restores drug sensitivity in AML via ROS-induced mitochondrial metabolic reprogramming

BACKGROUND

Acute myeloid leukemia (AML) is a highly aggressive cancer with a 5-year survival rate of less than 35%. It is characterized by significant drug resistance and abnormal energy metabolism. Mitochondrial dynamics and metabolism are crucial for AML cell survival. Mitochondrial fusion protein optic atrophy (OPA)1 is upregulated in AML patients with adverse mutations and correlates with poor prognosis.

METHOD

This study investigated targeting OPA1 with TMQ0153, a tetrahydrobenzimidazole derivative, to disrupt mitochondrial metabolism and dynamics as a novel therapeutic approach to overcome treatment resistance. Effects of TMQ0153 treatment on OPA1 and mitofusin (MFN)2 protein levels, mitochondrial morphology, and function in AML cells. In this study, we examined reactive oxygen species (ROS) production, oxidative phosphorylation (OXPHOS) inhibition, mitochondrial membrane potential (MMP) depolarization, and apoptosis. Additionally, metabolic profiling was conducted to analyze changes in metabolic pathways.

RESULTS

TMQ0153 treatment significantly reduced OPA1 and mitofusin (MFN)2 protein levels and disrupted the mitochondrial morphology and function in AML cells. This increases ROS production and inhibits OXPHOS, MMP depolarization, and caspase-dependent apoptosis. Metabolic reprogramming was observed, shifting from mitochondrial respiration to glycolysis and impaired respiratory chain activity. Profiling revealed reduced overall metabolism along with changes in the glutathione (GSH)/oxidized glutathione (GSSG) and NAD⁺/NADH redox ratios. TMQ0153 treatment reduces tumor volume and weight in MV4-11 xenografts in vivo. Combination therapies with TMQ0153 and other AML drugs significantly reduced the leukemic burden and prolonged survival in NOD scid gamma (NSG) mice xenografted with U937-luc and MOLM-14-luc cells.

CONCLUSION

TMQ0153 targets mitochondrial dynamics by inhibiting OPA1, inducing metabolic reprogramming, and triggering apoptosis in AML cells. It enhances the efficacy of existing AML therapies and provides a promising combination treatment approach that exploits mitochondrial vulnerability and metabolic reprogramming to improve treatment outcomes in AML.

Structural basis of SARS-CoV-2 polymerase inhibition by nonnucleoside inhibitor HeE1-2Tyr

Targeting the RNA-dependent RNA polymerase (RdRp) of SARS-CoV-2 with small molecules is a promising therapeutic strategy against COVID-19, but potent and safe inhibitors are lacking. HeE1-2Tyr, a nonnucleoside inhibitor of Dengue virus RdRp, was also shown to inhibit SARS-CoV-2 RdRp in vitro and to have antiviral activity in cells, but the underlying mechanism remains unclear. Here, we elucidate the molecular mechanism of HeE1-2Tyr-mediated SARS-CoV-2 RdRp inhibition. Biochemical assays confirm that HeE1-2Tyr inhibits RdRp with an IC50 of 5 µM and show that it competes with RNA binding to RdRp in vitro. Structural analysis using cryo-EM reveals that a stack of three HeE1-2Tyr molecules binds to the RNA binding site of RdRp. The identification of the conserved HeE1-2Tyr binding site and its intriguing inhibition mechanism of three stacked molecules that outcompete RNA may facilitate further development of pan-corona nonnucleoside inhibitors.

Docking and structure activity relationship studies of potent and selective thiazolidinethione GSK-3 inhibitors

Glycogen synthase kinase-3 (GSK-3) plays a key role in several biochemical pathways and is an attractive target for pharmacological intervention. We prepared a series of analogs of a highly selective thiazolidinethione inhibitor of GSK-3. The structure-activity relationship indicated a precise structural requirement for potent inhibition. We used docking and bioinformatic analysis to explore the rationale for the potency and selectivity of this class of GSK-3 inhibitors. These computational studies identified residues unique to GSK-3 likely to play a role in ligand-specific induced fit interactions. Together, these studies highlight the structural stringency of specific kinase inhibition that can be achieved for GSK-3.

Unraveling the molecular grammar and the structural transitions underlying the fibrillation of a viral fibrillogenic domain

Hendra virus (HeV) is a biosafety level 4 human pathogen belonging to the Henipavirus genus within the Paramyxoviridae family. In HeV, the phosphoprotein-encoding gene also drives the synthesis of the V and W proteins that are two major players in the host innate immune response evasion. These three proteins share a common intrinsically disordered N-terminal domain (NTD) and have distinct C-terminal domains. We recently reported the ability of a short region (i.e., PNT3), located within the shared NTD, to form fibrils. We subsequently identified a PNT3 motif (EYYY) critically involved in fibrillation and deciphered the contribution of each tyrosine to the process. Herein, we combined mutational studies with various biochemical and biophysical approaches to further investigate the molecular mechanisms underlying PNT3 fibrillation. The results show that (i) lysine residues play a critical role in driving fibrillation, (ii) hydrophobic residues affect the nucleation step, and (iii) charge distribution strongly affects the fibrillation propensities. Vibrational Raman spectroscopy data further validated the role of lysine residues in promoting fibrillation and enabled documenting the formation of cross-β amyloid structures. Altogether, these results illuminate the molecular mechanisms involved in fibril formation and pave the way towards the rational design of inhibitors.

Optimizing cation-π force fields for molecular dynamics studies of competitive solvation in conjugated organosulfur polymers for lithium-sulfur batteries

Lithium-sulfur (Li/S) batteries are emerging as a next-generation energy storage technology due to their high theoretical energy density and cost-effectiveness. π-Conjugated organosulfur polymers, such as poly(4-(thiophene-3-yl)benzenethiol) (PTBT), have shown promise in overcoming challenges such as the polysulfide shuttle effect by providing a conductive framework and enabling sulfur copolymerization. In these cathodes, cation-π interactions significantly influence Li+ diffusion and storage properties in π-conjugated cathodes, but classical OPLS-AA force fields fail to capture these effects. This study employs a bottom-up approach based on density functional theory (DFT) to optimize the nonbonded interaction parameters (OPLS-AA/corr.), particularly for the Li+-π interactions with the PTBT polymer. Following prior work, we used an ion-induced dipole potential to model the cation-π interactions. The impact of the solvent on the PTBT monomers was examined by computing the potential of mean force (PMF) between PTBT monomers and Li+ ions in both explicit and implicit solvents using the Boltzmann inversion of probability distributions close to room temperature. In the implicit solvent case, the magnitude of the binding free energy decreased with increasing dielectric constant, as the dominant electrostatics scaled with the dielectric constant. In contrast, in the explicit solvent case, considering the mixtures of organic solvent DME and DOL, the binding free energy shows minimal dependence on solvent composition due to the competing interaction of TBT and Li+ with the solvent molecules. However, increasing salt concentration decreases the binding free energy due to Debye-Hückel screening effects. In general, this work suggests that the optimized parameters can be widely used in the simulation of polymers in electrolytes for the Li/S battery to enhance the representation of cation-π interactions for a fixed charge force field.

Chemically Informed Coarse-Graining of Electrostatic Forces in Charge-Rich Biomolecular Condensates

Biomolecular condensates composed of highly charged biomolecules, such as DNA, RNA, chromatin, and nucleic-acid binding proteins, are ubiquitous in the cell nucleus. The biophysical properties of these charge-rich condensates are largely regulated by electrostatic interactions. Residue-resolution coarse-grained models that describe solvent and ions implicitly are widely used to gain mechanistic insights into the biophysical properties of condensates, offering transferability, computational efficiency, and accurate predictions for multiple systems. However, their predictive accuracy diminishes for charge-rich condensates due to the implicit treatment of solvent and ions. Here, we present Mpipi-Recharged, a residue-resolution coarse-grained model that improves the description of charge effects in biomolecular condensates containing disordered proteins, multidomain proteins, and/or disordered single-stranded RNAs. Mpipi-Recharged introduces a pair-specific asymmetric Yukawa electrostatic potential, informed by atomistic simulations. We show that this asymmetric coarse-graining of electrostatic forces captures intricate effects, such as charge blockiness, stoichiometry variations in complex coacervates, and modulation of salt concentration, without requiring explicit solvation. Mpipi-Recharged provides excellent agreement with experiments in predicting the phase behavior of highly charged condensates. Overall, Mpipi-Recharged improves the computational tools available to investigate the physicochemical mechanisms regulating biomolecular condensates, enhancing the scope of computer simulations in this field.

Prediction of Single-Mutation Effects for Fluorescent Immunosensor Engineering with an End-to-End Trained Protein Language Model

A quenchbody (Q-body) is a fluorophore-labeled homogeneous immunosensor in which the fluorophore is quenched by tryptophan (Trp) residues in the vicinity of the antigen-binding paratope and dequenched in response to antigen binding. Developing Q-bodies against targets on demand remains challenging due to the large sequence space of the complementarity-determining regions (CDRs) related to antigen binding and fluorophore quenching. In this study, we pioneered a strategy using high-throughput screening and a protein language model (pLM) to predict the effects of mutations on fluorophore quenching with single amino acid resolution, thereby enhancing the performance of Q-bodies. We collected yeasts displaying nanobodies with high- and low-quenching properties for the TAMRA fluorophore from a modified large synthetic nanobody library followed by next-generation sequencing. The pretrained pLM, connected to a single-layer perceptron, was trained end-to-end on the enriched CDR sequences. The achieved quenching prediction model that focused on CDR1 + 3 performed best in the evaluation with precision-recall curves. Using this model, we predicted and validated the effective mutations in two anti-SARS-CoV-2 nanobodies, RBD1i13 and RBD10i14, which converted them into Q-bodies. For RBD1i13, three Trp mutants were predicted to have high probability scores for quenching through in silico Trp scanning. These mutants were verified via yeast surface display, and all showed enhanced quenching. For RBD10i14, mutations at four positions close to an existing Trp gave high scores through in silico saturation mutagenesis scanning. Six of eight high-score mutants, derived from two mutants at each of the four positions, exhibited deeper quenching on the yeast surface. Next, combined with the investigation of antigen binding of the mutants, we successfully achieved Q-bodies with enhanced responses. Overall, our strategy allows the prediction of fluorescence responses solely on the basis of the antibody sequence and will be essential for the rational selection and design of antibodies to achieve immunosensors with larger responses.

Antibiofilm Activities of Tritrpticin Analogs Against Pathogenic Pseudomonas aeruginosa PA01 Strains

In our previous work, we showed that short antimicrobial hexapeptides (AMPs) containing three Trp and three Arg residues had a potent antibiofilm activity against a pathogenic Gram-positive Staphylococcus aureus MRSA strain. However, the activity of these hexapeptides against a Gram-negative Pseudomonas aeruginosa PA01 strain was relatively poor. Herein, we tested the longer 13-residue synthetic AMP tritrpticin-NH2 (Tritrp) and several of its analogs as potential antibiofilm agents that can prevent biofilm formation (MBIC) and/or cause biofilm dissolution (MBEC) for two P. aeruginosa PA01 strains, one of which expressed the GFP protein. Tritrp, a porcine cathelicidin, is currently the only known naturally occurring cationic AMP that has three Trp in sequence (WWW), a feature that was found to be important in our previous study. Our results show that several Tritrp analogs were effective. In particular, analogs with Pro substitutions that had altered peptide backbone structures compared to the naturally occurring amphipathic two-turn structure showed more potent MBIC and MBEC antibiofilm activities. Selectivity of the peptides towards P. aeruginosa could be improved by introducing the non-proteinogenic amino acid 2,3-diaminopropionic acid, rather than Arg or Lys, as the positively charged residues. Using 1H NMR spectroscopy, we also reinvestigated the role of the two Pro residues in cis-trans isomerism of the peptide in aqueous solution. Overall, our results show that the WWW motif embedded in longer cationic AMPs has considerable potential to combat biofilm formation in pathogenic Gram-negative strains.

The energetic landscape of CH-π interactions in protein-carbohydrate binding

CH-π interactions between carbohydrates and aromatic amino acids play an essential role in biological systems that span all domains of life. Quantifying the strength and importance of these CH-π interactions is challenging because these interactions involve several atoms and can exist in many distinct orientations. To identify an orientational landscape of CH-π interactions, we constructed a dataset of close contacts formed between β-d-galactose residues and the aromatic amino acids, tryptophan, tyrosine, and phenylalanine, across crystallographic structures deposited in the Protein Data Bank. We carried out quantum mechanical calculations to quantify their interaction strengths. The data indicate that tryptophan-containing CH-π interactions have more favorable interaction energies than those formed by tyrosine or phenylalanine. The energetic differences between these amino acids are caused by the aromatic ring system electronics and size. We use individual distance and angle features to train random forest models to successfully predict the first-principles computed energetics of CH-π interactions. Using insights from our models, we define a tradeoff in CH-π interaction strength arising from the proximity of galactose carbons 1 and 2 versus carbons 4 and 6 to the aromatic amino acid. Our work demonstrates that a feature of CH-π stacking interactions is that numerous orientations allow for highly favorable interaction strengths.

Mesoscale properties of biomolecular condensates emerging from protein chain dynamics

Biomolecular condensates form by phase separation of biological polymers and have important functions in the cell - functions that are inherently connected to their physical properties. A remarkable aspect of such condensates is that their viscoelastic properties can vary by orders of magnitude, but it has remained unclear how these pronounced differences are rooted in the nanoscale dynamics at the molecular level. Here we investigate a series of condensates formed by complex coacervation that span about two orders of magnitude in molecular dynamics, diffusivity, and viscosity. We find that the nanoscale chain dynamics on the nano- to microsecond timescale can be accurately related to both translational diffusion and mesoscale condensate viscosity by analytical relations from polymer physics. Atomistic simulations reveal that the observed differences in friction - a key quantity underlying these relations - are caused by differences in inter-residue contact lifetimes, leading to the vastly different dynamics among the condensates. The rapid exchange of inter-residue contacts we observe may be a general mechanism for preventing dynamic arrest in compartments densely packed with polyelectrolytes, such as the cell nucleus.

Lowering the affinity of single-chain monovalent BBB shuttle scFc-scFv8D3 prolongs its half-life and increases brain concentration

Monoclonal antibody therapeutics is a massively growing field. Progress in providing monoclonal antibody therapeutics to treat brain disorders is complicated, due to the impermeability of the blood-brain barrier (BBB) to large macromolecular structures. To date, the most successful approach for delivering antibody therapeutics to the brain is by targeting the transferrin receptor (TfR) using anti-TfR BBB shuttles, with the 8D3 antibody being one of the most extensively studied in the field. The strategy of fine-tuning TfR binding affinity has shown promise, with previous results showing an improved brain delivery of bivalent 8D3-BBB constructs. In the current study, a fine-tuning TfR affinity strategy has been employed to improve single-chain variable fragment (scFv) 8D3 (scFv8D3) affinity mutants. Initially, in silico protein-protein docking analysis was performed to identify amino acids (AAs) likely to contribute to 8D3s TfR binding affinity. Mutating the identified AAs resulted in decreased TfR binding affinity, increased blood half-life and increased brain concentration. As monovalent BBB shuttles are seemingly superior for delivering antibodies at therapeutically relevant doses, our findings and approach may be relevant for optimizing brain delivery.

Extracellular Peptide-Ligand Dimerization Actuator Receptor Design for Reversible and Spatially Dosed 3D Cell-Material Communication

Transmembrane receptors that endow mammalian cells with the ability to sense and respond to biomaterial-bound ligands will prove instrumental in bridging the fields of synthetic biology and biomaterials. Materials formed with thiol-norbornene chemistry are amenable to thiol-peptide patterning, and this study reports the rational design of synthetic receptors that reversibly activate cellular responses based on peptide-ligand recognition. This transmembrane receptor platform, termed Extracellular Peptide-ligand Dimerization Actuator (EPDA), consists of stimulatory or inhibitory receptor pairs that come together upon extracellular peptide dimer binding with corresponding monobody receptors. Intracellularly, Stimulatory EPDAs phosphorylate a substrate that merges two protein halves, whereas Inhibitory EPDAs revert split proteins back to their unmerged, inactive state via substrate dephosphorylation. To identify ligand-receptor pairs, over 2000 candidate monobodies were built in silico using PETEI, a novel computational algorithm we developed. The top 30 monobodies based on predicted peptide binding affinity were tested experimentally, and monobodies that induced the highest change in protein merging (green fluorescent protein, GFP) were incorporated in the final EPDA receptor design. In soluble form, stimulatory peptides induce intracellular GFP merging in a time- and concentration-dependent manner, and varying levels of green fluorescence were observed based on stimulatory and inhibitory peptide-ligand dosing. EPDA-programmed cells encapsulated in thiol-norbornene hydrogels patterned with stimulatory and inhibitory domains exhibited 3D activation or deactivation based on their location within peptide-patterned hydrogels. EPDA receptors can recognize a myriad of peptide-ligands bound to 3D materials, can reversibly induce cellular responses beyond fluorescence, and are widely applicable in biological research and regenerative medicine.

Structural basis for membrane association and catalysis by phosphatidylserine synthase in Escherichia coli

Phosphatidylserine synthase (PssA) is essential in the biosynthesis of phosphatidylethanolamine, a major phospholipid of bacterial membranes. A peripheral membrane protein PssA can associate with the cellular membrane in its active state or exist in the cytosol in an inactive form. The membrane-bound enzyme acts on cytidine diphosphate diacylglycerol (CDP-DG) to form cytidine monophosphate and a covalent intermediate, which is subsequently targeted by serine to produce phosphatidylserine. Here, we present two crystal structures of Escherichia coli PssA, one complexed with CDP-DG and the other without. The lipid-bound structure mimics the Michaelis complex before the formation of a covalent intermediate, revealing key determinants for substrate recognition and catalysis. Notably, membrane-free PssA is in a monomer-dimer equilibrium, with only the monomer capable of associating with the membrane, suggesting a regulatory mechanism for phospholipid biosynthesis dependent on the oligomerization state of the enzyme.

A recombinant L-threonine aldolase with high catalytic efficiency for the asymmetric synthesis of L-threo-phenylserine and L-threo-4-fluorophenylserine

OBJECTIVES

To develop robust variants of L-threonine aldolases (L-TAs), potent catalysts for synthesizing asymmetric β-hydroxy-α-amino acids, it is necessary to identify critical residues beyond the known active site residues.

RESULTS

Through virtual screening, a neglected residue Asn305, was identified as critical for catalytic efficiency. Subsequent site-saturation mutagenesis led to a potent variant N305R which exhibited excellent conversions of 88%conv (87%de) and 80%conv (94%de) for the synthesis of L-threo-phenylserine and L-threo-4-fluorophenylserine respectively. This variant not only outperformed the template enzyme, but also represented a promising L-TA for synthesizing the two β-hydroxy-α-amino acids. It was suggested that Arg305 of the variant N305R generated strong cation-arene interaction and electrostatic force with the intermediates, leading to strengthened binding, enhanced L-threo favored orientation and wider entrance.

CONCLUSIONS

Our work not only provided an excellent variant N305R, but also suggested the crucial function of a neglected residue Asn305, which offered valuable experiences for other L-TA studies.

Molecular Determinants for Guanine Binding in GTP-Binding Proteins: A Data Mining and Quantum Chemical Study

GTP-binding proteins are essential molecular switches that regulate a wide range of cellular processes. Their function relies on the specific recognition and binding of guanine within their binding pockets. This study aims to elucidate the molecular determinants underlying this recognition. A large-scale data mining of the Protein Data Bank yielded 298 GTP-binding protein complexes, which provided a structural foundation for a systematic analysis of the intermolecular interactions that are responsible for the molecular recognition of guanine in proteins. It was found that multiple modes of non-bonded interactions including hydrogen bonding, cation-π interactions, and π-π stacking interactions are employed by GTP-binding proteins for binding. Subsequently, the strengths of non-bonded interaction energies between guanine and its surrounding protein residues were quantified by means of the double-hybrid DFT method B2PLYP-D3/cc-pVDZ. Hydrogen bonds, particularly those involving the N2 and O6 atoms of guanine, confer specificity to guanine recognition. Cation-π interactions between the guanine ring and basic residues (Lys and Arg) provide significant electrostatic stabilization. π-π stacking interactions with aromatic residues (Phe, Tyr, and Trp) further contribute to the overall binding affinity. This synergistic interplay of multiple interaction modes enables GTP-binding proteins to achieve high specificity and stability in guanine recognition, ultimately underpinning their crucial roles in cellular signaling and regulation. Notably, the NKXD motif, while historically considered crucial for guanine binding in GTP-binding proteins, is not universally required. Our study revealed significant variability in hydrogen bonding patterns, with many proteins lacking the NKXD motif but still effectively binding guanine through alternative arrangements of interacting residues.

Glycerol 3-Phosphate Dehydrogenase Catalyzed Hydride Transfer: Enzyme Activation by Cofactor Pieces

Glycerol 3-phosphate dehydrogenase catalyzes reversible hydride transfer from glycerol 3-phosphate (G3P) to NAD+ to form dihydroxyacetone phosphate; from the truncated substrate ethylene glycol to NAD+ in a reaction activated by the phosphite dianion substrate fragment; and from G3P to the truncated nicotinamide riboside cofactor in a reaction activated by adenosine 5'-diphosphate, adenosine 5'-monophosphate, and ribose 5-phosphate cofactor fragments. The sum of the stabilization of the transition state for GPDH-catalyzed hydride transfer reactions of the whole substrates by the phosphodianion fragment of G3P and the ADP fragment of NAD+ is 25 kcal/mol. Fourteen kcal/mol of this transition state stabilization is recovered as phosphite dianion and AMP activation of the reactions of the substrate and cofactor fragments. X-ray crystal structures for unliganded GPDH, for a binary GPDH·NAD+ complex, and for a nonproductive ternary GPDH·NAD+·DHAP complex show that the ligand binding energy is utilized to drive an extensive protein conformational change that creates a caged complex for these ligands. The phosphite dianion and AMP fragments are proposed to activate GPDH for the catalysis of hydride transfer by stabilization of this active caged complex. The closure of a conserved loop [292-LNGQKL-297] during substrate binding stabilizes the G3P and NAD+ complexes by interactions, respectively, with the Q295 and K296 loop side chains. The appearance and apparent conservation of two side chains that interact with the hydride donor and acceptor to stabilize the active closed enzyme are proposed to represent a significant improvement in the catalytic performance of GPDH.

In Silico Assessment for Risk of Possible Human Transmission of FCoV-23

Since the pandemic in 2019, coronaviruses (CoVs) have been a great concern for public health burden. The fact that CoVs can infect all animals including domestic ones and livestock points to a future pandemic even though interaction between human and wildlife animals is restricted. Moreover, interspecies transmission abilities of CoVs by mutations make them drastically risky not only for humans but also for animal health. Recently, a new CoV outbreak in cats in Cyprus, the so-called FCoV-23, has been realized. In addition to worries over animal health, any possible transmission to humans is now controversial. However, there have been limited characterization studies on FCoV-23. Thus, we aimed to assess the possible transmission of FCoV-23 to humans using in silico prediction tools. Accordingly, we first checked the binding affinities of receptor binding domain (RBD) of FCoV-23 against feline target protein and its human homolog. Next, we randomly and rationally created mutations on the RBD sequence and evaluated the binding affinities using protein docking tools. Our results underlined that multiple mutations at the same time were needed for increased binding affinity towards human target protein, demonstrating that the probability of transmission to humans was extremely low when mutation rates were regarded. Still, further molecular studies are required to comprehensively conclude the possible transmission risk.

Computational insights into intrinsically disordered regions in protein-nucleic acid complexes

We study transitions in intrinsically disordered regions (IDRs) upon complex formation, utilizing X-ray-solved structural dataset of protein-DNA and protein-RNA complexes, along with their available unbound protein forms. The identified IDRs are categorized into three classes: Disordered-to-Ordered (D-O), Disordered-to-Partial Ordered (D-PO) and Disordered-to-Disordered (D-D) after comparing them in unbound and complex forms. In the D-O class, IDRs form secondary structures like coils, helices, and strands upon binding to nucleic acids. Though a majority of these IDRs are present at the surface of the complexes, a significant number of IDRs are also observed at the interfaces and are involved in polar interactions. The hydrogen bonds made by the interface IDRs (B_IDRs) with phosphates and bases of nucleic acids are comparatively more than those formed with sugars. B_IDRs form more H-bonds with the ribose in protein-RNA than with the deoxyribose in protein-DNA. Among the B_IDRs, Arg and Lys prefer to interact with the major and minor grooves of DNA and RNA, respectively. Ser, however, prefers the minor groove in both the nucleic acids. Interestingly, we report 61 and 48 IDRs in 31 protein-DNA and 22 protein-RNA complexes, respectively, suggesting nucleic acid binding to proteins may also result in ordered-to-disordered transitions.

The [(bathophenanthroline)3:Fe2+] complex as an aromatic non-polymeric medium for purification of human lactoferrin

We describe a non-chromatographic, ligand-free platform for the efficient purification of recombinant human lactoferrin (LF). The platform consists of a [metal:chelator] complex precipitate in the presence of osmotically active polyethylene glycol 6000 (PEG-6000). Purification is achieved in three stages. Following formation of the complex, LF is captured under neutral conditions by the aggregated complexes (Step I), a washing step follows (Step II) and then, (Step III) LF is extracted in pure form with 100 mM tribasic Na citrate buffer (pH 7). Of the four complexes investigated, [bathophenanthroline (batho)3:Fe2+] was determined to be the most efficient. LF is recovered with high yield (∼90%, by densitometry) and purity (≥97%, by SDS polyacrylamide gel electrophoresis (SDS-PAGE)) from an artificial contamination background comprising E. coli lysate proteins. Purified LF is demonstrated to be monomeric by dynamic light scattering (DLS); to preserve its native secondary structure by circular dichroism (CD) spectroscopy; and, as apo-LF, to efficiently inhibit bacterial growth. Process yield is not affected by a 45-fold increase in LF concentration from 0.2 to 9 mg/mL. We provide evidence that protein capture relies on [cation:π] interactions between the lysine and arginine residues of LF with the fully aromatic [(batho)3:Fe2+] complexes. The use of [metal:chelator] complex aggregates is demonstrated to provide an economical and efficient avenue for LF purification.

Supramolecular Complex Surfactants and Structured Liquids Enabled by Cation-π and Charge-Transfer Interactions

Cation-π and charge-transfer (CT) interactions are pervasive with significant implications in the fields of chemistry, materials science, and biology. However, much less is known about the construction of interfacial assemblies based on the two interactions. Here, by combining cation-π and CT interactions between an acceptor molecule, dicationic naphthalenediimide, and an aromatic donor, pyrene-terminated poly-l-lactic acid, we report the generation of supramolecular complex surfactants (SCSs) in situ at the toluene-water interface. The utilization of SCSs as building blocks enables the fabrication of interfacial assemblies including 2D films, emulsions, and structured liquids. By modification of the redox state of the acceptor molecules under chemical stimulus, the association/assembly and dissociation/disassembly of SCSs can be precisely regulated, imparting intriguing redox-responsive properties to the resulting assemblies.

A suicidal and extensively disordered luciferase with a bright luminescence

Gaussia luciferase (GLuc) is one of the most luminescent luciferases known and is widely used as a reporter in biochemistry and cell biology. During catalysis, GLuc undergoes inactivation by irreversible covalent modification. The mechanism by which GLuc generates luminescence and how it becomes inactivated are however not known. Here, we show that GLuc unlike other enzymes has an extensively disordered structure with a minimal hydrophobic core and no apparent binding pocket for the main substrate, coelenterazine. From an alanine scan, we identified two Arg residues required for light production. These residues separated with an average of about 22 Å and a major structural rearrangement is required if they are to interact with the substrate simultaneously. We furthermore show that in addition to coelenterazine, GLuc also can oxidize furimazine, however, in this case without production of light. Both substrates result in the formation of adducts with the enzyme, which eventually leads to enzyme inactivation. Our results demonstrate that a rigid protein structure and substrate-binding site are no prerequisites for high enzymatic activity and specificity. In addition to the increased understanding of enzymes in general, the findings will facilitate future improvement of GLuc as a reporter luciferase.

Expanding the molecular language of protein liquid-liquid phase separation

Understanding the relationship between a polypeptide sequence and its phase separation has important implications for analysing cellular function, treating disease and designing novel biomaterials. Several sequence features have been identified as drivers for protein liquid-liquid phase separation (LLPS), schematized as a 'molecular grammar' for LLPS. Here we further probe how sequence modulates phase separation and the material properties of the resulting condensates, targeting sequence features previously overlooked in the literature. We generate sequence variants of a repeat polypeptide with either no charged residues, high net charge, no glycine residues or devoid of aromatic or arginine residues. All but one of 12 variants exhibited LLPS, albeit to different extents, despite substantial differences in composition. Furthermore, we find that all the condensates formed behaved like viscous fluids, despite large differences in their viscosities. Our results support the model of multiple interactions between diverse residue pairs-not just a handful of residues-working in tandem to drive the phase separation and dynamics of condensates.

Cellular and Biophysical Applications of Genetic Code Expansion

Despite their diverse functions, proteins are inherently constructed from a limited set of building blocks. These compositional constraints pose significant challenges to protein research and its practical applications. Strategically manipulating the cellular protein synthesis system to incorporate novel building blocks has emerged as a critical approach for overcoming these constraints in protein research and application. In the past two decades, the field of genetic code expansion (GCE) has achieved significant advancements, enabling the integration of numerous novel functionalities into proteins across a variety of organisms. This technological evolution has paved the way for the extensive application of genetic code expansion across multiple domains, including protein imaging, the introduction of probes for protein research, analysis of protein-protein interactions, spatiotemporal control of protein function, exploration of proteome changes induced by external stimuli, and the synthesis of proteins endowed with novel functions. In this comprehensive Review, we aim to provide an overview of cellular and biophysical applications that have employed GCE technology over the past two decades.

Aromatic and arginine content drives multiphasic condensation of protein-RNA mixtures

Multiphasic architectures are found ubiquitously in biomolecular condensates and are thought to have important implications for the organization of multiple chemical reactions within the same compartment. Many of these multiphasic condensates contain RNA in addition to proteins. Here, we investigate the importance of different interactions in multiphasic condensates comprising two different proteins and RNA using computer simulations with a residue-resolution coarse-grained model of proteins and RNA. We find that in multilayered condensates containing RNA in both phases, protein-RNA interactions dominate, with aromatic residues and arginine forming the key stabilizing interactions. The total aromatic and arginine content of the two proteins must be appreciably different for distinct phases to form, and we show that this difference increases as the system is driven toward greater multiphasicity. Using the trends observed in the different interaction energies of this system, we demonstrate that we can also construct multilayered condensates with RNA preferentially concentrated in one phase. The "rules" identified can thus enable the design of synthetic multiphasic condensates to facilitate further study of their organization and function.

One touch is all it takes: the supramolecular interaction between ubiquitin and lanthanide complexes revisited by paramagnetic NMR and molecular dynamics

The supramolecular interaction between lanthanide complexes and proteins is at the heart of numerous chemical and biological studies. Some of these complexes have demonstrated remarkable interaction properties with proteins or peptides in solution and in the crystalline state. Here we have used the paramagnetism of lanthanide ions to characterize the affinity of two lanthanide complexes for ubiquitin. As the interaction process is dynamic, the acquired NMR data only reflect the time average of the different steps. We have used molecular dynamics (MD) simulations to get a deeper insight into the detailed interaction scenario at the microsecond scale. This NMR/MD approach enabled us to establish that the tris-dipicolinate complex interacts specifically with arginines and lysines, while the crystallophore explores the protein surface through weak interactions with carboxylates. These observations shed new light on the dynamic interaction properties of these complexes, which will ultimately enable us to propose a crystallization mechanism.

Intramolecular Cation-π Interactions Organize Bowl-Shaped, Luminescent Molecular Containers

Molecules with nonplanar architectures are highly desirable due to their unique topological structures and functions. We report here the synthesis of two molecular containers (1 ⋅ 3Br- and 1 ⋅ 3Cl-), which utilize intramolecular cation-π interactions to enforce macrocylic arrangements and exhibit high binding affinity and luminescent properties. Remarkably, the geometry of the cation-π interaction can be flexibly tailored to achieve a precise ring arrangement, irrespective of the angle of the noncovalent bonds. Additionally, the C-H⋅⋅⋅Br- hydrogen bonds within the container are also conducive to stabilizing the bowl-shaped conformation. These bowl-shaped conformations were confirmed both in solution through NMR spectroscopy and in the solid state by X-ray studies. 1 ⋅ 3Br- shows high binding affinity and selectivity: F->Cl-, through C-H⋅⋅⋅X- (X=F, Cl) hydrogen bonds. Additionally, these containers exhibited blue fluorescence in solution and yellow room-temperature phosphorescence (RTP) in the solid state. Our findings illustrate the utility of cation-π interactions in designing functional molecules.

Reexamining the diverse functions of arginine in biochemistry

Arginine in a free-state and as part of peptides and proteins shows distinct tendency to form clusters. In free-form, it has been found useful in cryoprotection, as a drug excipient for both solid and liquid formulations, as an aggregation suppressor, and an eluent in protein chromatography. In many cases, the mechanisms by which arginine acts in all these applications is either debatable or at least continues to attract interest. It is quite possible that arginine clusters may be involved in many such applications. Furthermore, it is possible that such clusters are likely to behave as intrinsically disordered polypeptides. These considerations may help in understanding the roles of arginine in diverse applications and may even lead to better strategies for using arginine in different situations.

Bioisoteres for carboxylic acids: From ionized isosteres to novel unionized replacements

Carboxylic acids are key pharmacophoric elements in many molecules. They can be seen as a problem by some, due to perceived permeability challenges, potential for high plasma protein binding and the risk of forming reactive metabolites due to acyl-glucuronidation. By others they are viewed more favorably as they can decrease lipophilicity by adding an ionizable center which can be beneficial for solubility, and can add enthalpic interactions with the target protein. However, there are many instances where the replacement of a carboxylic acid with a bioisosteric group is required. This has led to the development of a number of ionizable groups which sufficiently mimic the carboxylic acid functionality whilst improving, for example, the metabolic profile of the molecule in question. An alternative strategy involves replacement of the carboxylate by neutral functional groups. This review initially details carefully selected examples whereby tetrazoles, acyl sulfonamides or isoxazolols have been beneficially utilized as carboxylic acid bioisosteres altering physicohemical properties, interactions with the target and metabolism and/or pharmacokinetics, before delving further into the binding mode of carboxylic acid derivatives with their target proteins. This analysis highlights new ways to consider the replacement of carboxylic acids by neutral bioisosteric groups which either rely on hydrogen bonds or cation-π interactions. It should serve as a useful guide for scientists working in drug discovery.

Situating the phosphonated calixarene-cytochrome C association by molecular dynamics simulations

Protein-calixarenes binding plays an increasingly central role in many applications, spanning from molecular recognition to drug delivery strategies and protein inhibition. These ligands obey a specific bio-supramolecular chemistry, which can be revealed by computational approaches, such as molecular dynamics simulations. In this paper, we rely on all-atom, explicit-solvent molecular dynamics simulations to capture the electrostatically driven association of a phosphonated calix-[4]-arene with cytochome-C, which critically relies on surface-exposed paired lysines. Beyond two binding sites identified in direct agreement with the x-ray structure, the association has a larger structural impact on the protein dynamics. Then, our simulations allow a direct comparison to analogous calixarenes, namely, sulfonato, similarly reported as "molecular glue." Our work can contribute to a robust in silico predictive tool to assess binding sites for any given protein of interest for crystallization, with the specificity of a macromolecular cage whose endo/exo orientation plays a role in the binding.

Omniligase-1-Mediated Phage-Peptide Library Modification and Insulin Engineering

Chemical and enzymatic modifications of peptide-displayed libraries have been successfully employed to expand the phage display library. However, the requirement of specific epitopes and scaffolds has limited the scope of protein engineering using phage display. In this study, we present a novel approach utilizing omniligase-1-mediated selective and specific ligation on the phage pIII protein, offering a high conversion rate and compatibility with commercially available phage libraries. We applied this method to perform high-throughput engineering of insulin analogues with randomized B chain C-terminal regions. Insulin analogues with different B chain C-terminal segments were selected and exhibited biological activity equivalent to that of human insulin. Molecular dynamics studies of insulin analogues revealed a novel interaction between the insulin B27 residue and insulin receptor L1 domain. In summary, our findings highlight the potential of omniligase-1-mediated phage display in the development and screening of disulfide-rich peptides and proteins. This approach holds promise for the creation of novel insulin analogues with enhanced therapeutic properties and exhibits potential for the development of other therapeutic compounds.

Biological importance of arginine: A comprehensive review of the roles in structure, disorder, and functionality of peptides and proteins

Arginine shows Jekyll and Hyde behavior in several respects. It participates in protein folding via ionic and H-bonds and cation-pi interactions; the charge and hydrophobicity of its side chain make it a disorder-promoting amino acid. Its methylation in histones; RNA binding proteins; chaperones regulates several cellular processes. The arginine-centric modifications are important in oncogenesis and as biomarkers in several cardiovascular diseases. The cross-links involving arginine in collagen and cornea are involved in pathogenesis of tissues but have also been useful in tissue engineering and wound-dressing materials. Arginine is a part of active site of several enzymes such as GTPases, peroxidases, and sulfotransferases. Its metabolic importance is obvious as it is involved in production of urea, NO, ornithine and citrulline. It can form unusual functional structures such as molecular tweezers in vitro and sprockets which engage DNA chains as part of histones in vivo. It has been used in design of cell-penetrating peptides as drugs. Arginine has been used as an excipient in both solid and injectable drug formulations; its role in suppressing opalescence due to liquid-liquid phase separation is particularly very promising. It has been known as a suppressor of protein aggregation during protein refolding. It has proved its usefulness in protein bioseparation processes like ion-exchange, hydrophobic and affinity chromatographies. Arginine is an amino acid, whose importance in biological sciences and biotechnology continues to grow in diverse ways.

Adenosine Triphosphate Mediates Phase Separation of Disordered Basic Proteins by Bridging Intermolecular Interaction Networks

Adenosine triphosphate (ATP) is an abundant molecule with crucial cellular roles as the energy currency and a building block of nucleic acids and for protein phosphorylation. Here we show that ATP mediates the phase separation of basic intrinsically disordered proteins (bIDPs). In the resulting condensates, ATP is highly concentrated (apparent partition coefficients up to 7700) and serves as bridges between bIDP chains. These liquid-like droplets have some of the lowest interfacial tension (∼25 pN/μm) but high zero-shear viscosities (1-15 Pa s) due to the bridged protein networks, and yet their fusion has some of the highest speeds (∼1 μm/ms). The rapid fusion manifests extreme shear thinning, where the apparent viscosity is lower than zero-shear viscosity by over 100-fold, made possible by fast reformation of the ATP bridges. At still higher concentrations, ATP does not dissolve bIDP droplets but results in aggregates and fibrils.

Design principles for cyclin K molecular glue degraders

Molecular glue degraders are an effective therapeutic modality, but their design principles are not well understood. Recently, several unexpectedly diverse compounds were reported to deplete cyclin K by linking CDK12-cyclin K to the DDB1-CUL4-RBX1 E3 ligase. Here, to investigate how chemically dissimilar small molecules trigger cyclin K degradation, we evaluated 91 candidate degraders in structural, biophysical and cellular studies and reveal all compounds acquire glue activity via simultaneous CDK12 binding and engagement of DDB1 interfacial residues, in particular Arg928. While we identify multiple published kinase inhibitors as cryptic degraders, we also show that these glues do not require pronounced inhibitory properties for activity and that the relative degree of CDK12 inhibition versus cyclin K degradation is tuneable. We further demonstrate cyclin K degraders have transcriptional signatures distinct from CDK12 inhibitors, thereby offering unique therapeutic opportunities. The systematic structure-activity relationship analysis presented herein provides a conceptual framework for rational molecular glue design.

Effects of arginine on coenzyme-Q10 micelle uptake for mitochondria-targeted nanotherapy in phenylketonuria

Phenylketonuria (PKU) is a rare inherited metabolic disease characterized by phenylalanine hydroxylase enzyme deficiency. In PKU patients, coenzyme Q10 (CoQ10) levels were found low. Therefore, we focused on the modification of CoQ10 to load the micelles and increase entry of micelles into the cell and mitochondria, and it is taking a part in ATP turnover. Micelles had produced by comparing two different production methods (thin-film layer and direct-dissolution), and characterization studies were performed (zeta potential, size, and encapsulation efficiency). Then, L-arginine (LARG) and poly-arginine (PARG) were incorporated with the micelles for subsequential release and PKU cell studies. The effects of these components on intracellular uptake and their use in the cellular cycle were analyzed by ELISA, Western blot, membrane potential measurement, and flow cytometry methods. In addition, both effects of LARG and PARG micelles on pharmacokinetics at the cellular level and their cell binding rate were determined. The thin-film method was found superior in micelle preparation. PARG/LARG-modified micelles showed sustained release. In the cellular and mitochondrial uptake of CoQ10, CoQ10-micelle + PARG > CoQ10-micelle + LARG > CoQ10-micelle > CoQ10 was found. This increased localization caused lowering of oxygen consumption rates, but maintaining mitochondrial membrane potential. The study results had showed that besides micelle formulation, PARG and LARG are effective in cellular and mitochondrial targeting.

A Transient π-π or Cation-π Interaction between Degron and Degrader Dual Residues: A Key Step for the Substrate Recognition and Discrimination in the Processive Degradation of SulA by ClpYQ (HslUV) Protease in Escherichia coli

The Escherichia coli ATP-dependent ClpYQ protease constitutes ClpY ATPase/unfoldase and ClpQ peptidase. The Tyr91st residue within the central pore-I site of ClpY-hexamer is important for unfolding and translocating substrates into the catalytic site of ClpQ. We have identified the degron site (GFIMRP147th) of SulA, a cell-division inhibitor recognized by ClpYQ and that the Phe143rd residue in degron site is necessary for SulA native folded structure. However, the functional association of this degron site with the ClpYQ degrader is unknown. Here, we investigated the molecular insights into substrate recognition and discrimination by the ClpYQ protease. We found that the point mutants ClpYY91FQ, ClpYY91HQ, and ClpYY91WQ, carrying a ring structure at the 91st residue of ClpY, efficiently degraded their natural substrates, evidenced by the suppressed bacterial methyl-methane-sulfonate (MMS) sensitivity, the reduced β-galactosidase activity of cpsB::lacZ, and the lowest amounts of MBP-SulA in both in vivo and in vitro degradation analyses. Alternatively, mimicking the wild-type SulA, SulAF143H, SulAF143K and SulAF143W, harboring a ring structure or a cation side-group in 143rd residue of SulA, were efficiently degraded by ClpYQ in the bacterial cells, also revealing shorter half-lives at 41 °C and higher binding affinities towards ClpY in pull-down assays. Finally, ClpYY91FQ and ClpYY91HQ, were capable of effectively degrading SulAF143H and SulAF143K, highlighting a correspondingly functional interaction between the SulA 143rd and ClpY 91st residues. According to the interchangeable substituted amino acids, our results uniquely indicate that a transient π-π or cation-π interaction between the SulA 143rd and ClpY 91st residues could be aptly gripped between the degron site of substrates and the pore site of proteases (degraders) for substrate recognition and discrimination of the processive degradation.

Identification of allosteric inhibitor against AKT1 through structure-based virtual screening

AKT (serine/threonine protein kinase) is a potential therapeutic target for many types of cancer as it plays a vital role in cancer progression. Many AKT inhibitors are already in practice under single and combinatorial therapy. However, most of these inhibitors are orthosteric / pan-AKT that are non-selective and non-specific to AKT kinase and their isoforms. Hence, researchers are searching for novel allosteric inhibitors that bind in the interface between pH and kinase domain. In this study, we performed structure-based virtual screening from the afroDB (a diverse natural compounds library) to find the potential inhibitor targeting the AKT1. These compounds were filtered through Lipinski, ADMET properties, combined with a molecular docking approach to obtain the 8 best compounds. Then we performed molecular dynamics simulation for apoprotein, AKT1 with 8 complexes, and AKT1 with the positive control (Miransertib). Molecular docking and simulation analysis revealed that Bianthracene III (hit 1), 10-acetonyl Knipholonecyclooxanthrone (hit 2), Abyssinoflavanone VII (hit 5) and 8-c-p-hydroxybenzyldiosmetin (hit 6) had a better binding affinity, stability, and compactness than the reference compound. Notably, hit 1, hit 2 and hit 5 had molecular features required for allosteric inhibition.

1,2,3-Triazoles in Biomolecular Crystallography: A Geometrical Data-Mining Approach

The 1,2,3-triazole scaffold has become very attractive to identify new chemical entities in drug discovery projects. Despite the widespread use of click chemistry to synthesize numerous 123Ts, there are few drugs on the market that incorporate this scaffold as a substructure. To investigate the true potential of 123Ts in protein-ligand interactions, we examined the noncovalent interactions between the 1,2,3-triazole ring and amino acids in protein-ligand cocrystals using a geometrical approach. For this purpose, we constructed a nonredundant database of 220 PDB IDs from available 123T-protein cocrystal structures. Subsequently, using the Protein Ligand Interaction Profiler web platform (PLIP), we determined whether 1,2,3-triazoles primarily act as linkers or if they can be considered interactive scaffolds. We then manually analyzed the geometrical descriptors from 333 interactions between 1,4-disubstituted 123T rings and amino acid residues in proteins. This study demonstrates that 1,2,3-triazoles exhibit diverse preferred interactions with amino acids, which contribute to protein-ligand binding.

Electrostatic Interactions Dictate Bile Salt Hydrolase Substrate Preference

The human intestines are colonized by trillions of microbes, comprising the gut microbiota, which produce diverse small molecule metabolites and modify host metabolites, such as bile acids, that regulate host physiology. Biosynthesized in the liver, bile acids are conjugated with glycine or taurine and secreted into the intestines, where gut microbial bile salt hydrolases (BSHs) deconjugate the amino acid to produce unconjugated bile acids that serve as precursors for secondary bile acid metabolites. Among these include a recently discovered class of microbially conjugated bile acids (MCBAs), wherein alternative amino acids are conjugated onto bile acids. To elucidate the metabolic potential of MCBAs, we performed detailed kinetic studies to investigate the preference of BSHs for host-conjugated bile acids and MCBAs. We identified a BSH that exhibits positive cooperativity uniquely for MCBAs containing an aromatic side chain. Further molecular modeling and phylogenetic analyses indicated that the BSH preference for aromatic MCBAs is due to a substrate-specific cation-π interaction and is predicted to be widespread among human gut microbial BSHs.

Dynamic extracellular vestibule of human SERT: Unveiling druggable potential with high-affinity allosteric inhibitors

The serotonin transporter (SERT) tightly regulates synaptic serotonin levels and has been the primary target of antidepressants. Binding of inhibitors to the allosteric site of human SERT (hSERT) impedes the dissociation of antidepressants bound at the central site and may enhance the efficacy of such antidepressants to potentially reduce their dosage and side effects. Here, we report the identification of a series of high-affinity allosteric inhibitors of hSERT in a unique scaffold, with the lead compound, Lu AF88273 (3-(1-(2-(1H-indol-3-yl)ethyl)piperidin-4-yl)-6-chloro-1H-indole), having 2.1 nM allosteric potency in inhibiting imipramine dissociation. In addition, we find that Lu AF88273 also inhibits serotonin transport in a noncompetitive manner. The binding pose of Lu AF88273 in the allosteric site of hSERT is determined with extensive molecular dynamics simulations and rigorous absolute binding free energy perturbation (FEP) calculations, which show that a part of the compound occupies a dynamically formed small cavity. The predicted binding location and pose are validated by site-directed mutagenesis and can explain much of the structure-activity relationship of these inhibitors using the relative binding FEP calculations. Together, our findings provide a promising lead compound and the structural basis for the development of allosteric drugs targeting hSERT. Further, they demonstrate that the divergent allosteric sites of neurotransmitter transporters can be selectively targeted.

The Prenucleation Equilibrium of the Parathyroid Hormone Determines the Critical Aggregation Concentration and Amyloid Fibril Nucleation

Nucleation and growth of amyloid fibrils were found to only occur in supersaturated solutions above a critical concentration (ccrit ). The biophysical meaning of ccrit remained mostly obscure, since typical low values of ccrit in the sub-μM range hamper investigations of potential oligomeric states and their structure. Here, we investigate the parathyroid hormone PTH84 as an example of a functional amyloid fibril forming peptide with a comparably high ccrit of 67±21 μM. We describe a complex concentration dependent prenucleation ensemble of oligomers of different sizes and secondary structure compositions and highlight the occurrence of a trimer and tetramer at ccrit as possible precursors for primary fibril nucleation. Furthermore, the soluble state found in equilibrium with fibrils adopts to the prenucleation state present at ccrit . Our study sheds light onto early events of amyloid formation directly related to the critical concentration and underlines oligomer formation as a key feature of fibril nucleation. Our results contribute to a deeper understanding of the determinants of supersaturated peptide solutions. In the current study we present a biophysical approach to investigate ccrit of amyloid fibril formation of PTH84 in terms of secondary structure, cluster size and residue resolved intermolecular interactions during oligomer formation. Throughout the investigated range of concentrations (1 μM to 500 μM) we found different states of oligomerization with varying ability to contribute to primary fibril nucleation and with a concentration dependent equilibrium. In this context, we identified the previously described ccrit of PTH84 to mark a minimum concentration for the formation of homo-trimers/tetramers. These investigations allowed us to characterize molecular interactions of various oligomeric states that are further converted into elongation competent fibril nuclei during the lag phase of a functional amyloid forming peptide.

Electrostatic Interactions Dictate Bile Salt Hydrolase Substrate Preference

The human intestines are colonized by trillions of microbes, comprising the gut microbiota, which produce diverse small molecule metabolites and modify host metabolites, such as bile acids, that regulate host physiology. Biosynthesized in the liver, bile acids are conjugated with glycine or taurine and secreted into the intestines, where gut microbial bile salt hydrolases (BSHs) deconjugate the amino acid to produce unconjugated bile acids that serve as precursors for secondary bile acid metabolites. Among these include a recently discovered class of microbially-conjugated bile acids (MCBAs), wherein alternative amino acids are conjugated onto bile acids. To elucidate the metabolic potential of MCBAs, we performed detailed kinetic studies to investigate the preference of BSHs for host- and microbially-conjugated bile acids. We identified a BSH that exhibits positive cooperativity uniquely for MCBAs containing an aromatic sidechain. Further molecular modeling and phylogenetic analyses indicated that BSH preference for aromatic MCBAs is due to a substrate-specific cation-π interaction and is predicted to be widespread among human gut microbial BSHs.

Host Membranes as Drivers of Virus Evolution

The molecular mechanisms controlling the adaptation of viruses to host cells are generally poorly documented. An essential issue to resolve is whether host membranes, and especially lipid rafts, which are usually considered passive gateways for many enveloped viruses, also encode informational guidelines that could determine virus evolution. Due to their enrichment in gangliosides which confer an electronegative surface potential, lipid rafts impose a first control level favoring the selection of viruses with enhanced cationic areas, as illustrated by SARS-CoV-2 variants. Ganglioside clusters attract viral particles in a dynamic electrostatic funnel, the more cationic viruses of a viral population winning the race. However, electrostatic forces account for only a small part of the energy of raft-virus interaction, which depends mainly on the ability of viruses to form a network of hydrogen bonds with raft gangliosides. This fine tuning of virus-ganglioside interactions, which is essential to stabilize the virus on the host membrane, generates a second level of selection pressure driven by a typical induced-fit mechanism. Gangliosides play an active role in this process, wrapping around the virus spikes through a dynamic quicksand-like mechanism. Viruses are thus in an endless race for access to lipid rafts, and they are bound to evolve perpetually, combining speed (electrostatic potential) and precision (fine tuning of amino acids) under the selective pressure of the immune system. Deciphering the host membrane guidelines controlling virus evolution mechanisms may open new avenues for the design of innovative antivirals.

A Route to Potent, Selective, and Biased Salvinorin Chemical Space

The salvinorins serve as templates for next generation analgesics, antipruritics, and dissociative hallucinogens via selective and potent agonism of the kappa-opioid receptor (KOR). In contrast to most opioids, the salvinorins lack basic amines and bind with high affinity and selectivity via complex polyoxygenated scaffolds that have frustrated deep-seated modification by synthesis. Here we describe a short asymmetric synthesis that relies on a sterically confined organocatalyst to dissociate acidity from reactivity and effect Robinson annulation of an unactivated nucleophile/unstable electrophile pair. Combined with a cobalt-catalyzed polarized diene-alkyne cycloaddition, the route allows divergent access to a focused library of salvinorins. We appraise the synthesis by its generation of multiple analogs that exceed the potency, selectivity, stability, and functional bias of salvinorin A itself.

ATP Mediates Phase Separation of Disordered Basic Proteins by Bridging Intermolecular Interaction Networks

ATP is an abundant molecule with crucial cellular roles as the energy currency and a building block of nucleic acids and for protein phosphorylation. Here we show that ATP mediates the phase separation of basic intrinsically disordered proteins (bIDPs). In the resulting condensates, ATP is highly concentrated (apparent partition coefficients at 200-5000) and serves as bridges between bIDP chains. These liquid-like droplets have some of the lowest interfacial tension (~25 pN/μm) but high zero-shear viscosities (1-15 Pa s) due to the bridged protein networks, and yet their fusion has some of the highest speeds (~1 μm/ms). The rapid fusion manifests extreme shear thinning, where the apparent viscosity is lower than zero-shear viscosity by over 100-fold, made possible by fast reformation of the ATP bridges. At still higher concentrations, ATP does not dissolve bIDP droplets but results in aggregates and fibrils.

Design, synthesis, and mechanistic study of 2-piperazineone-bearing peptidomimetics as novel HIV capsid modulators

HIV-1 capsid (CA) is an attractive target for its indispensable roles in the viral life cycle. We report the design, synthesis, and mechanistic study of a novel series of 2-piperazineone peptidomimetics as HIV capsid modulators by mimicking the structure of host factors binding to CA. F-Id-3o was the most potent compound from the synthesized series, with an anti-HIV-1 EC50 value of 6.0 μM. However, this series of compounds showed a preference for HIV-2 inhibitory activity, in which Id-3o revealed an EC50 value of 2.5 μM (anti-HIV-2 potency), an improvement over PF74. Interestingly, F-Id-3o did bind HIV-1 CA monomers and hexamers with comparable affinity, unlike PF74, consequently showing antiviral activity in the early and late stages of the HIV-1 lifecycle. Molecular dynamics simulations shed light on F-Id-3o and Id-3o binding modes within the HIV-1/2 CA protein and provide a possible explanation for the increased anti-HIV-2 potency. Metabolic stability assays in human plasma and human liver microsomes indicated that although F-Id-3o has enhanced metabolic stability over PF74, further optimization is necessary. Moreover, we utilized computational prediction of drug-like properties and metabolic stability of F-Id-3o and PF74, which correlated well with experimentally derived metabolic stability, providing an efficient computational pipeline for future preselection based on metabolic stability prediction. Overall, the 2-piperazineone-bearing peptidomimetics are a promising new chemotype in the CA modulators class with considerable optimization potential.

PIQLE: protein-protein interface quality estimation by deep graph learning of multimeric interaction geometries

Motivation

Accurate modeling of protein-protein interaction interface is essential for high-quality protein complex structure prediction. Existing approaches for estimating the quality of a predicted protein complex structural model utilize only the physicochemical properties or energetic contributions of the interacting atoms, ignoring evolutionarily information or inter-atomic multimeric geometries, including interaction distance and orientations.

Results

Here, we present PIQLE, a deep graph learning method for protein-protein interface quality estimation. PIQLE leverages multimeric interaction geometries and evolutionarily information along with sequence- and structure-derived features to estimate the quality of individual interactions between the interfacial residues using a multi-head graph attention network and then probabilistically combines the estimated quality for scoring the overall interface. Experimental results show that PIQLE consistently outperforms existing state-of-the-art methods including DProQA, TRScore, GNN-DOVE and DOVE on multiple independent test datasets across a wide range of evaluation metrics. Our ablation study and comparison with the self-assessment module of AlphaFold-Multimer repurposed for protein complex scoring reveal that the performance gains are connected to the effectiveness of the multi-head graph attention network in leveraging multimeric interaction geometries and evolutionary information along with other sequence- and structure-derived features adopted in PIQLE.

Availability and implementation

An open-source software implementation of PIQLE is freely available at https://github.com/Bhattacharya-Lab/PIQLE.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.