Cation-π Interactions and their Functional Roles in Membrane Proteins

The α/β3 complex of human voltage-gated sodium channel hNav1.7 to study mechanistic differences in presence and absence of auxiliary subunit β3

CONTEXT

In the context of structural interactomics, we generated a 3D model between α and β3 subunits for the hitherto unknown human voltage-gated sodium channel complex (hNa 1.7α/β3). We embedded our 3D model in a membrane lipid bilayer for molecular dynamics (MD) simulations of the sodium cation passage from the outer vestibule through the inner pore segment of our hNa 1.7 complex in presence and absence of auxiliary subunit β3 with remarkable changes close to electrophysiological study results. A complete passage could not be expected due to because the inactivated state of the underlying 3D template. A complete sodium ion passage would require an open state of the channel. The computed observations concerning side chain rearrangements for favorable cooperativity under evolutionary neighborhood conditions, favorable and unfavorable amino acid interactions, proline kink, loop, and helix displacements were all found in excellent keeping with the extant literature without any exception nor contradiction. Complex-stabilizing pairs of interacting amino acids with evolutionary neighborhood complementary were identified.

METHODS

The following tools were used: sequence search and alignment by FASTA and Clustal Omega; 3D model visualization and homology modeling by Vega ZZ, SPDBV, Chimera and Modeller, respectively; missing sections (loops) by Alphafold; geometry optimization prior to MD runs by GROMACS 2021.4 under the CHARMM 36 force field; local healing of bad contacts by SPDBV based on its Ramachandran plots; protein-protein docking by HDOCK 2.4; membrane insertion assisted by OPM; Berendsen V-rescaling for NVT; Parrinello-Rahman and Nose-Hoover for MPT; MD analyses by VMD and XMGRACE.

Evaluation of Orientation-Dependent Cation-π Pairwise Effects within Collagen Triple Helices

Various noncovalent interactions have been introduced to explore their impacts in folding a collagen triple helix. Among these interactions, the cation-π interaction represents one of the compelling forces stabilizing the triple helix. Still, the effects depend on the pairwise components and the orientation between the cationic and aromatic moieties. To gain more insights into this interaction within a collagen trimer, we prepared a series of collagen-mimetic peptides (CMPs) with cationic residues and aromatic residues incorporated to examine the contributions of two types of axial cation-π pairs (N → C and C → N cationic-to-aromatic pairwise) and the lateral cation-π pair. Circular dichroism (CD) measurements indicate that the N → C axial pairs have a significant stabilization effect. In contrast, the lateral and the C → N axial pairs destabilize the fold, and the lateral pairs cause the most destabilization consequences. We further designed and prepared the CMPs containing various lateral and axial cation-π pairs to investigate the coupling consequences in homotrimers and heterotrimers. From CD data, we found that the predicted differences in melting temperatures using individual cation-π pairwise contributions were comparable to the observed values for the designed homotrimers. CD and NMR measurements showed favorable cation-π interactions could effectively induce the folding of heterotrimers, in which the CMPs with more N → C axial pairs formed a more stable trimer than those containing a smaller number of N → C axial pairs. In this study, we have disclosed more valuable information about the properties of cation-π pairwise effects within a collagen triple helix, which can be considered in designing collagen-related peptides and materials.

Endowing Polycyclic Aromatic Hydrocarbons with Induced Supramolecular Chirality by Engineering Cation-π Interaction

The cation-π interaction is a crucial noncovalent binding force prevalent in many biological and chemical systems, however, are rarely reported in chiral supramolecular systems. Here, a supramolecular chiral system is reported, wherein cation-π interaction play significant roles. The interaction involves a chiral cationic palladium complex (R-Pd) co-assembling with 12 different polycyclic aromatic hydrocarbons (PAHs), resulting in the formation of supramolecular gels in mixed solvents. This co-assembly induces supramolecular chirality in the PAHs, with cation-π interaction identified as the primary driving force. The confirmation of this interaction is supported by high-resolution mass, 1H NMR, nuclear overhauser effect spectroscopy (NOESY), XRD and CD. Most importantly, with the synergy of hydrogen bonds, the cation-π interaction can be stabilized in a mixed solvent containing water, leading the co-assembly to form supramolecular gel, as well as enhanced chiral induction of PAHs. This discovery serves as a proof-of-principle example showing how the weak cation-π interaction could be a powerful guide in attaining induced chirality of PAHs.

Synthesis of Selectively 13C/2H/15N- Labeled Arginine to Probe Protein Conformation and Interaction by NMR Spectroscopy

The charged arginine side chain is unique in determining many innate properties of proteins, contributing to stability and interaction surfaces, and directing allosteric regulation and enzymatic catalysis. NMR experiments can be used to reveal these processes at the molecular level, but it often requires selective insertion of carbon-13, nitrogen-15, and deuterium at defined atomic positions. We introduce a method to endow arginine residues with defined isotope patterns, combining synthetic organic chemistry and cell-based protein overexpression. The resulting proteins feature NMR active spin systems with optimized relaxation pathways leading to simplified NMR spectra with a sensitive response to changes in the chemical environment of the nuclei observed.

Multiphasic condensates formed with mono-component of tetrapeptides via phase separation

Biomolecular condensates, formed by liquid-liquid phase separation of biomacromolecules, play crucial roles in regulating physiological events in biological systems. While multiphasic condensates have been extensively studied, those derived from a single component of short peptides have not yet been reported. Here, we report the symmetrical core-shell structural biomolecular condensates formed with a programmable tetrapeptide library via phase separation. Our findings reveal that tryptophan is essential for core-shell structure formation due to its strongest homotypical π-π interaction, enabling us to modulate the structure of condensates from core-shell to homogeneous by altering the amino acid composition. Molecular dynamics simulation combined with cryogenic focused ion beam scanning electron microscopy and cryogenic electron microscopy show that the inner core of multiphasic tetrapeptide condensates is solid-like, consisting of ordered structures. The core is enveloped by a liquid-like shell, stabilizing the core structure. Furthermore, we demonstrate control over multiphasic condensate formation through intrinsic redox reactions or post-translational modifications, facilitating the rational design of synthetic multiphasic condensates for various applications on demand.

The Cation-π Interaction in Chemistry and Biology

The cation-π interaction is an important noncovalent binding force that impacts all areas of chemistry and biology. Extensive computational and gas phase experimental studies have established the potential strength and the essential nature of the interaction. Previous reviews have emphasized studies of model systems and a variety of biological examples. This work includes discussion of those areas but emphasizes other areas that are perhaps less well appreciated. These include the novel cation-π binding ability of alkali metals in water; the application of the cation-π interaction to organic synthesis and chemical biology; cooperative behaviors of multiple cation-π interactions, including adhesive proteins from mussels and similar organisms and the formation and modulation of biomolecular condensates (phase separation); and cation-π interactions involved in recognizing DNA/RNA.

Designing Effective Antimicrobial Agents: Structural Insights into the Antibiofilm Activity of Ionic Liquids

Research concerning biofilm control is critical due to the pervasive and resilient nature of biofilms, which pose significant challenges across the industrial, environmental, and healthcare sectors. Traditional antimicrobial treatments are often ineffective against these robust structures. Here, we explore the antimicrobial properties of ionic liquids (ILs) and their efficacy in biofilm disruption. By examining the structural variations of ILs, we highlight the key role of hydrophobicity, noting that longer alkyl side chains in IL cations enhance biofilm disruption and bacterial death. However, upon reaching a certain optimal chain length─usually C12 to C14─the antimicrobial activity of ILs starts to decrease. Furthermore, the symmetry and size of anions significantly impact biofilm elimination. This Perspective addresses a critical gap in biofilm research, revealing the structure-activity relationships of ILs and providing a foundation for designing more effective biofilm-disrupting agents.

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.

Sticky Polyelectrolyte Shield for Enhancing Biological Half-Life of Growth Factors

Delivery of secretomes, which includes growth factors, cytokines, and mRNA, is critical in regenerative medicine for cell-to-cell communication. However, the harsh in vivo environment presents significant challenges for secretome delivery. Proteolytic enzymes shorten secretomes' half-lives, and secretomes tend to rapidly diffuse at defect sites. Therefore, a delivery system that ensures prolonged retention and enhanced therapeutic efficacy of secretomes is required. In this study, a Coating Optimized Drug Delivery Enhancement (COD2E) system, a coacervate composed of dopamine functionalized fucoidan and poly-l-lysine, was fabricated for secretome delivery. The dopamine modification significantly enhanced adhesive strength (>7-fold) compared to that of the neat coacervates, which enabled rapid (5 min) and uniform coating ability on collagen sponges. The COD2E system was able to encapsulate fibroblast growth factor (FGF2) and prolong the half-life of FGF2. Notably, its efficacy, demonstrated through a single application of FGF2 encapsulated COD2E on collagen sponge, in a wound model demonstrated a successful tissue repair. The COD2E system is an effective growth factor delivery vehicle since it can protect growth factors, has an antioxidant ability, adheres on various material surfaces, and is hemocompatible.

Antibiofilm activities of lactoferricin-related Trp- and Arg-rich antimicrobial hexapeptides against pathogenic Staphylococcus aureus and Pseudomonas aeruginosa strains

Recently, several antimicrobial peptides (AMPs), varying in length from 12 to 37 residues, have been shown to act as antibiofilm agents. Here, we report a study of 23 hexapeptides modeled after four different Trp- and Arg-rich AMPs, including the RRWQWR-NH2 peptide, derived from bovine lactoferrin. They were tested against the pathogenic Gram-negative Pseudomonas aeruginosa PAO1 strain and a Gram-positive Staphylococcus aureus MRSA strain. Both strains were engineered to express the green fluorescent protein (GFP) protein, and fluorescence detection was used to measure the ability of the peptides to prevent biofilm formation (minimum biofilm inhibitory concentration (MBIC)) or to cause the breakdown of established biofilms (minimum biofilm eradication concentration (MBEC)). Similar antibiofilm activities were obtained with the standard crystal violet dye assay. Most Trp- and Arg-rich hexapeptides displayed a potent antibiofilm activity against the Gram-positive S. aureus MRSA strain. In particular, hexapeptides with 3 Arg and 3 Trp were very effective, especially when they contained the three Trp in sequence. Somewhat unexpectedly, the antimicrobial (MIC) values correlated with the MBIC and MBEC values, which has not been seen for several other AMP/antibiofilm peptides. Our results demonstrate that short Trp- and Arg-rich peptides merit further studies as antibiofilm agents that could be deployed to address part of the antimicrobial resistance problem.

pH-Dependent Assembly and Stability of Toll-Like Receptor 3/dsRNA Signaling Complex: Insights from Constant pH Molecular Dynamics and Metadynamics Simulations

The pH-dependent assembly of Toll-like receptors (TLRs), which triggers a threshold-like response, is a key principle in immune signaling. While crystallography has revealed the intricate structure of these assembly complexes, the mechanisms underlying their pH dependency remain unclear. Herein, constant pH simulations and metadynamics are employed to investigate the pH-dependent assembly and stability of the TLR3/dsRNA signaling complex. The findings demonstrate that system pH regulates complex assembly and stability by modulating the protonation and charge states of histidines. Histidines in TLR3 act as pH-dependent, positively charged binding sites that capture negatively charged dsRNA. Additionally, these histidines form a [H682⁺]-[E626⁻] dipole, facilitating the assembly of two TLR3 molecules into an antisymmetric dimer through dipole-dipole interactions. Surprisingly, TLR3 can shift the pKa values of key histidines from their model pKa of 6.5, increasing protonation likelihood and enhancing ligand binding. Notably, the aromatic residue Phe84, located within the dsRNA binding site [His39⁺-His60⁺-Phe84-His108⁺], alters the pKa of His60 through cation-π interactions with its protonated state. This study offers new insights into the molecular mechanisms underlying pH-dependent immune signaling via higher-order assemblies and suggests potential applications for histidine in self-assembling biomaterials.

Appendix

Since the first Chapter dealt with the well-known charge-charge interactions familiar to biologists, this concluding Chapter introduces some key electrical forces, probably much less familiar, perhaps even unknown. LLPS (liquid liquid phase separation) which we have seen involved in so much of cell biology depends on multivalent, π-π and cation-π electrical forces. How these arise is dealt with here and may be especially useful as an aide memoir to return to when such issues arise within the bulk of the text.

Insights into the adsorption of emerging organic contaminant by low-cost readily separable modified jute fiber

A high-efficiency biosorbent based on the low-priced jute fiber was developed, characterized, and applied to remove the emerging organic contaminant diclofenac from aqueous solutions. Jute fiber was treated by NaOH (named AJF) followed by grafting different amounts of trimethyl[3-(trimethoxysilyl) propyl] ammonium chloride (named AJF-TTSAC). The composition, morphology, porosity, and adsorption features of the neat and modified jute fiber were evaluated and compared. The surface of neat JF was smooth, nonporous, and free of cracks. NaOH treatment increased the fibrillation, created cracks and grooves, and increased the oxygen content, total pore volume, and surface area. In comparison to AJF, grafting TTSAC filled in the crevices, grooves, and spaces between fibrillates, and decreased the total pore volume and surface area. The adsorption of diclofenac by the neat and modified JF occurred at highly acidic pHo and peaked at pHo 3. Among the neat and modified JF, AJF-TTSAC5 was the most efficient followed by AJF. The efficiency of AJF and AJF-TTSAC5 was highest using 1.00 g/L, at 35 °C and was not affected by the presence of NaCl. The Elovich, pseudo-first-order, and pseudo-second-order models described the adsorption kinetic satisfactorily with the marginal advantage of Elovich for AJF and pseudo-second-order for AJF-TTSAC5. The isotherm study exposed the multilayer and physisorption nature of the adsorption of diclofenac onto AJF and AJF-TTSAC5. The Langmuir monolayer saturation capacity of AJF-TTSAC5 was 37.43 mg/g which revealed its great potential relative to other adsorbents in the literature. The AJF and AJF-TTSAC5 were easily regenerated using distilled water and kept good performance for 5 repetitive cycles.

Tunable Coassembly of Octaarginine with Thiazolyl Benzenesulfonamides Exerts Variable Antibacterial Activity

The cationic peptide octaarginine (R8) is a prominent cell-penetrating peptide and has been extensively researched as a carrier of diverse cell-destined cargo. In this work, we describe the coassembly of R8 with small molecule thiazolyl benzenesulfonamide (TBS) derivatives. Physical complexation of R8 with three TBS derivatives across a range of weight ratios results in the formation of a distinctive set of nano- and microstructures. A detailed structural characterization of the R8:TBS-derivative coassemblies has been performed by a combination of FTIR, XRD, SEM, and DSC. The major functional groups that facilitate coassembly include sulfonamide SO2 and NH groups of the TBS derivatives, and the guanidinium of R8, via a combination of cation-π and hydrogen-bonding interactions. The R8:4F-TBS coassembly displays singular topological features compared to R8:4Br-TBS and R8:4CH3-TBS complexes. These differences are attributed to the changes in the preferred orientation of the guanidino groups of R8 with respect to the π-surface of TBS derivatives. The modulation of forces of interaction across the R8:TBS-derivative coassemblies aligns with their respective thermal stabilities. The single-crystal structure of bare 4F-TBS has been subjected to Hirshfeld and 2D fingerprinting analysis and indicates notable variations from the crystal packing of the R8:4F-TBS coassembly. The structural differences among the R8:TBS-derivative coassemblies correlate with distinctive profiles of antibacterial activity in each case. The coassembled structures exert a variable extent of bacterial membrane disruption and damage based on the unique disposition of R8 and the potency of small molecule in each case. The aqueous suspension of R8:4F-TBS displays significant outer membrane disruption and bacterial killing compared with the other complexes. This work successfully demonstrates the hitherto unreported potential for coassembly of cell-penetrating peptides with other entities. The coassembly of R8 with small molecules highlights an attractive strategy for tuning the functional properties of each component.

Genetic Code Expansion for Mechanistic Studies in Ion Channels: An (Un)natural Union of Chemistry and Biology

Ion channels play central roles in biology and human health by catalyzing the transmembrane flow of electrical charge. These proteins are ideal targets for genetic code expansion (GCE) methods because it is feasible to measure ion channel activity from miniscule amounts of protein and to analyze the resulting data via rigorous, established biophysical methods. In an ideal scenario, the encoding of synthetic, noncanonical amino acids via GCE allows the experimenter to ask questions inaccessible to traditional methods. For this reason, GCE has been successfully applied to a variety of ligand- and voltage-gated channels wherein extensive structural, functional, and pharmacological data exist. Here, we provide a comprehensive summary of GCE as applied to ion channels. We begin with an overview of the methods used to encode noncanonical amino acids in channels and then describe mechanistic studies wherein GCE was used for photochemistry (cross-linking; caged amino acids) and atomic mutagenesis (isosteric manipulation of charge and aromaticity; backbone mutation). Lastly, we cover recent advances in the encoding of fluorescent amino acids for the real-time study of protein conformational dynamics.

Structural Insights and Catalytic Mechanism of 3-Hydroxybutyryl-CoA Dehydrogenase from Faecalibacterium Prausnitzii A2-165

Atopic dermatitis (AD) is characterized by a T-helper cell type 2 (Th2) inflammatory response leading to skin damage with erythema and edema. Comparative fecal sample analysis has uncovered a strong correlation between AD and Faecalibacterium prausnitzii strain A2-165, specifically associated with butyrate production. Therefore, understanding the functional mechanisms of crucial enzymes in the butyrate pathway, such as 3-hydroxybutyryl-CoA dehydrogenase of A2-165 (A2HBD), is imperative. Here, we have successfully elucidated the three-dimensional structure of A2HBD in complex with acetoacetyl-CoA and NAD+ at a resolution of 2.2Å using the PAL-11C beamline (third generation). Additionally, X-ray data of A2HBD in complex with acetoacetyl-CoA at a resolution of 1.9 Å were collected at PAL-XFEL (fourth generation) utilizing Serial Femtosecond Crystallography (SFX). The monomeric structure of A2HBD consists of two domains, N-terminal and C-terminal, with cofactor binding occurring at the N-terminal domain, while the C-terminal domain facilitates dimerization. Our findings elucidate the binding mode of NAD+ to A2HBD. Upon acetoacetyl-CoA binding, the crystal structure revealed a significant conformational change in the Clamp-roof domain (root-mean-square deviation of 2.202 Å). Notably, residue R143 plays a critical role in capturing the adenine phosphate ring, underlining its significance in substrate recognition and catalytic activity. The binding mode of acetoacetyl-CoA was also clarified, indicating its lower stability compared to NAD+. Furthermore, the conformational change of hydrophobic residues near the catalytic cavity upon substrate binding resulted in cavity shrinkage from an open to closed conformation. This study confirms the conformational changes of catalytic triads involved in the catalytic reaction and presents a proposed mechanism for substrate reduction based on structural observations.

The Interaction Of Homodimer Styrylcyanine Dyes With Sodium Dodecyl Sulfate And Triton X-100

Solubilization of the styrylcyanine dye Sbt ((E)-2-(4-(dimethylamino)styryl)-3-methylbenzo[d]thiazol-3-ium iodide) and its homodimers Dbt-5 and Dbt-10 in aqueous solution of sodium dodecyl sulfate and Triton X-100 has been studied by steady state and picosecond time-resolved fluorescence spectroscopy. At low concentration of sodium dodecyl sulfate in solution, between Sbt, Dbt-5 dyes molecules and surfactant ion pairs are formed followed by the formation non-luminescent H-aggregates. The nature of the interaction between molecules of dyes and surfactants has been revealed. The binding constants Ks of the dyes to the surfactants, free energy changes (ΔG0), the number of dye molecules (n) included in a single micelle and photophysical parameters have been determined. The degree of solubilization of dyes in micellar solution of Triton X-100 is higher as compared to sodium dodecyl sulfate and depends on the molecular weight and size of both dye molecules and micelles.

Aromatic Residues in Proteins: Re-Evaluating the Geometry and Energetics of π-π, Cation-π, and CH-π Interactions

Aromatic residues can participate in various biomolecular interactions, such as π-π, cation-π, and CH-π interactions, which are essential for protein structure and function. Here, we re-evaluate the geometry and energetics of these interactions using quantum mechanical (QM) calculations, focusing on pairwise interactions involving the aromatic amino acids Phe, Tyr, and Trp and the cationic amino acids Arg and Lys. Our findings reveal that π-π interactions, while energetically favorable, are less abundant in structured proteins than commonly assumed and are often overshadowed by previously underappreciated, yet prevalent, CH-π interactions. Cation-π interactions, particularly those involving Arg, show strong binding energies and a specific geometric preference toward stacked conformations, despite the global QM minimum, suggesting that a rather perpendicular T-shape conformation should be more favorable. Our results support a more nuanced understanding of protein stabilization via interactions involving aromatic residues. On the one hand, our results challenge the traditional emphasis on π-π interactions in structured proteins by showing that CH-π and cation-π interactions contribute significantly to their structure. On the other hand, π-π interactions appear to be key stabilizers in solvated regions and thus may be particularly important to the stabilization of intrinsically disordered proteins.

Ion-combination specific effects driving the enzymatic activity of halophilic alcohol dehydrogenase 2 from Haloferax volcanii in aqueous ionic liquid solvent mixtures

Biocatalysis in ionic liquids enables novel routes for bioprocessing. Enzymes derived from extremophiles promise greater stability and activity under ionic liquid (IL) influence. Here, we probe the enzyme alcohol dehydrogenase 2 from the halophilic archaeon Haloferax volcanii in thirteen different ion combinations for relative activity and analyse the results against molecular dynamics (MD) simulations of the same IL systems. We probe the ionic liquid property space based on ion polarizability and molecular electrostatic potential. Using the radial distribution functions, survival probabilities and spatial distribution functions of ions, we show that cooperative ion-ion interactions determine ion-protein interactions, and specifically, strong ion-ion interactions equate to higher enzymatic activity if neither of the ions interact strongly with the protein surface. We further demonstrate a tendency for cations interacting with the protein surface to be least detrimental to enzymatic activity if they show a low polarizability when combined with small hydrophilic anions. We also find that the IL ion influence is not mitigated by the surplus of negatively charged residues of the halophilic enzyme. This is shown by free energy landscape analysis in root mean square deviation and distance variation plots of active site gating residues (Trp43 and His273) demonstrating no protection of specific structural elements relevant to preserving enzymatic activity. On the other hand, we observe a general effect across all IL systems that a tight binding of water at acidic residues is preferentially interrupted at these residues through the increased presence of potassium ions. Overall, this study demonstrates a co-ion interaction dependent influence on allosteric surface residues controlling the active/inactive conformation of halophilic alcohol dehydrogenase 2 and the necessity to engineer ionic liquid systems for enzymes that rely on the integrity of functional surface residues regardless of their halophilicity or thermophilicity for use in bioprocessing.

Model Mechanism for Lipid Uptake by the Human STARD2/PC-TP Phosphatidylcholine Transfer Protein

The human StAR-related lipid transfer domain protein 2 (STARD2), also known as phosphatidylcholine (PC) transfer protein, is a single-domain lipid transfer protein thought to transfer PC lipids between intracellular membranes. We performed extensive μs-long molecular dynamics simulations of STARD2 of its apo and holo forms in the presence or absence of complex lipid bilayers. The simulations in water reveal ligand-dependent conformational changes. In the 2 μs-long simulations of apo STARD2 in the presence of a lipid bilayer, we observed spontaneous reproducible PC lipid uptake into the protein hydrophobic cavity. We propose that the lipid extraction mechanism involves one to two metastable states stabilized by choline-tyrosine or choline-tryptophane cation-π interactions. Using free energy perturbation, we evaluate that PC-tyrosine cation-π interactions contribute 1.8 and 2.5 kcal/mol to the affinity of a PC-STARD2 metastable state, thus potentially providing a significant decrease of the energy barrier required for lipid desorption.

Histidine in Proteins: pH-Dependent Interplay between π-π, Cation-π, and CH-π Interactions

Histidine (His) stands out as the most versatile natural amino acid due to its side chain's facile propensity to protonate at physiological pH, leading to a transition from aromatic to cationic characteristics and thereby enabling diverse biomolecular interactions. In this study, our objective was to quantify the energetics and geometries of pairwise interactions involving His at varying pH levels. Through quantum chemical calculations, we discovered that His exhibits robust participation in both π-π and cation-π interactions, underscoring its ability to adopt a π or cationic nature, akin to other common residues. Of particular note, we found that the affinity of protonated His for aromatic residues (via cation-π interactions) is greater than the affinity of neutral His for either cationic residues (also via cation-π interactions) or aromatic residues (via π-π interactions). Furthermore, His frequently engages in CH-π interactions, and notably, depending on its protonation state, we found that some instances of hydrogen bonding by His exhibit greater stability than is typical for interamino acid hydrogen bonds. The strength of the pH-dependent pairwise energies of His with aromatic residues is supported by the abundance of pairwise interactions with His of low and high predicted pKa values. Overall, our findings illustrate the contribution of His interactions to protein stability and its potential involvement in conformational changes despite its relatively low abundance in proteins.

Relationships between conformational and vibrational features of tryptophan characteristic Raman markers

Tryptophan (Trp) residue provides characteristic vibrational markers to the middle wavenumber spectral region of the Raman spectra recorded from peptides and proteins. In this report, we were particularly interested in eight Trp Raman markers, referred to as Wi (i = 1,…,8). All responsible for pronounced Raman lines, these markers originate from indole moiety, a bicyclic conjugated segment involved in the Trp structure. Numerous investigations have previously attempted to relate the variations observed in the spectral features of these markers to the environmental changes of Trp residues. To emphasize the most important points we can mention (i) the variations in the Raman profile of W4 (∼1360 cm-1) and W5 (∼1340 cm-1), frequently observed as a doublet with variable intensity ratio. These two markers were thought to result from a Fermi-resonance effect between certain planar and nonplanar modes; (ii) the changes observed in the wavenumbers and relative intensities of W4, W7 (∼880 cm-1) and W8 (∼760 cm-1) were supposed to be related to the accessibility of Trp to surrounding water molecules; and (iii) the wavenumber fluctuations of W3 (∼1550 cm-1), taken as a Trp side chain orientational marker. However, some ambiguities still exist regarding the interpretation of these markers, needing further clarification. Herein, upon a joint experimental and theoretical analysis based on a multiconformational approach, attention was paid to the relationships between structural and vibrational features of three indole-containing compounds with increasing structural complexity, i.e., skatole (3-methylindole), tryptophan, and tripeptide Gly-Trp-Gly. This study clearly shows that the existing assignments given to certain Trp Raman markers should be reconsidered, especially those based on the Fermi-resonance origin of W4-W5 (∼1360-1340 cm-1) doublet, as well as the purely environmental dependence of W7 and W8 markers.

Exploring the Therapeutic Potential of 5-Fluorouracil-Loaded Calcium Carbonate Nanoparticles Combined with Natural Compound Thymoquinone for Colon Cancer Treatment

Given the need for novel and effective therapies for colon cancer, this study aimed to investigate the effects of 5-fluorouracil-loaded calcium carbonate nanoparticles (5FU-CaCO3np) combined with thymoquinone (TQ) against colon cancer. A shaking incubator and a high-speed homogenizer were used to prepare the optimal 5FU-CaCO3np, with characterizations of physicochemical properties, in vitro drug release profile, and biocompatibility. In vitro experiments and molecular docking were employed to evaluate the therapeutic potential of the combination for colon cancer treatment. Study results revealed that 5FU-CaCO3np with a size of approximately 130 nm was synthesized using the high-speed homogenizer. Its favorable biocompatibility, pH sensitivity, and sustained release properties facilitated reduced toxic side effects of 5-FU on NIH3T3 normal cells and enhanced inhibitory effects on CT26 colon cancer cells. The combination of 5FU-CaCO3np (1.875 μM) and TQ (30 μM) showed significantly superior anti-colon cancer effects to 5FU-CaCO3np alone in terms of cell proliferation and migration inhibition, cell apoptosis induction, and spheroid growth suppression in CT26 cells (p < 0.05), with strong interactions between the drugs and targets (E-cadherin, Bcl-2, PCNA, and MMP-2). These results provide evidence for 5FU-CaCO3np as a novel regimen against colon cancer. Combining 5FU-CaCO3np and TQ may offer a new perspective for colon cancer therapy.

RING 4.0: faster residue interaction networks with novel interaction types across over 35,000 different chemical structures

Residue interaction networks (RINs) are a valuable approach for representing contacts in protein structures. RINs have been widely used in various research areas, including the analysis of mutation effects, domain-domain communication, catalytic activity, and molecular dynamics simulations. The RING server is a powerful tool to calculate non-covalent molecular interactions based on geometrical parameters, providing high-quality and reliable results. Here, we introduce RING 4.0, which includes significant enhancements for identifying both covalent and non-covalent bonds in protein structures. It now encompasses seven different interaction types, with the addition of π-hydrogen, halogen bonds and metal ion coordination sites. The definitions of all available bond types have also been refined and RING can now process the complete PDB chemical component dictionary (over 35000 different molecules) which provides atom names and covalent connectivity information for all known ligands. Optimization of the software has improved execution time by an order of magnitude. The RING web server has been redesigned to provide a more engaging and interactive user experience, incorporating new visualization tools. Users can now visualize all types of interactions simultaneously in the structure viewer and network component. The web server, including extensive help and tutorials, is available from URL: https://ring.biocomputingup.it/.

Deciphering the etiology of undiagnosed ocular anomalies along with systemic alterations in pediatric patients through whole exome sequencing

Inherited and developmental eye diseases are quite diverse and numerous, and determining their genetic cause is challenging due to their high allelic and locus heterogeneity. New molecular approaches, such as whole exome sequencing (WES), have proven to be powerful molecular tools for addressing these cases. The present study used WES to identify the genetic etiology in ten unrelated Mexican pediatric patients with complex ocular anomalies and other systemic alterations of unknown etiology. The WES approach allowed us to identify five clinically relevant variants in the GZF1, NFIX, TRRAP, FGFR2 and PAX2 genes associated with Larsen, Malan, developmental delay with or without dysmorphic facies and autism, LADD1 and papillorenal syndromes. Mutations located in GZF1 and NFIX were classified as pathogenic, those in TRRAP and FGFR2 were classified as likely pathogenic variants, and those in PAX2 were classified as variants of unknown significance. Protein modeling of the two missense FGFR2 p.(Arg210Gln) and PAX2 p.(Met3Thr) variants showed that these changes could induce potential structural alterations in important functional regions of the proteins. Notably, four out of the five variants were not previously reported, except for the TRRAP gene. Consequently, WES enabled the identification of the genetic cause in 40% of the cases reported. All the syndromes reported herein are very rare, with phenotypes that may overlap with other genetic entities.

Hydrated cation-π interactions of π-electrons with hydrated Mg2+ and Ca2+ cations

Hydrated cation-π interactions at liquid-solid interfaces between hydrated cations and aromatic ring structures of carbon-based materials are pivotal in many material, biological, and chemical processes, and water serves as a crucial mediator in these interactions. However, a full understanding of the hydrated cation-π interactions between hydrated alkaline earth cations and aromatic ring structures, such as graphene remains elusive. Here, we present a molecular picture of hydrated cation-π interactions for Mg2+ and Ca2+ by using the density functional theory methods. Theoretical results show that the graphene sheet can distort the hydration shell of the hydrated Ca2+ to interact with Ca2+ directly, which is water-cation-π interactions. In contrast, the hydration shell of the hydrated Mg2+ is quite stable and the graphene sheet interacts with Mg2+ indirectly, mediated by water molecules, which is the cation-water-π interactions. These results lead to the anomalous order of adsorption energies for these alkaline earth cations, with hydrated Mg2+-π < hydrated Ca2+-π when the number of water molecules is large (n ≥ 6), contrary to the order observed for cation-π interactions in the absence of water molecules (n = 0). The behavior of hydrated alkaline earth cations adsorbed on a graphene surface is mainly attributed to the competition between the cation-π interactions and hydration effects. These findings provide valuable details of the structures and the adsorption energy of hydrated alkaline earth cations adsorbed onto the graphene surface.

Efficient overexpression and purification of severe acute respiratory syndrome coronavirus 2 nucleocapsid proteins in Escherichia coli

The fundamental biology of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) nucleocapsid protein (Ncap), its use in diagnostic assays and its potential application as a vaccine component have received considerable attention since the outbreak of the Covid19 pandemic in late 2019. Here we report the scalable expression and purification of soluble, immunologically active, SARS-CoV-2 Ncap in Escherichia coli. Codon-optimised synthetic genes encoding the original Ncap sequence and four common variants with an N-terminal 6His affinity tag (sequence MHHHHHHG) were cloned into an inducible expression vector carrying a regulated bacteriophage T5 synthetic promoter controlled by lac operator binding sites. The constructs were used to express Ncap proteins and protocols developed which allow efficient production of purified Ncap with yields of over 200 mg per litre of culture media. These proteins were deployed in ELISA assays to allow comparison of their responses to human sera. Our results suggest that there was no detectable difference between the 6His-tagged and untagged original Ncap proteins but there may be a slight loss of sensitivity of sera to other Ncap isolates.

Complex biophysical changes and reduced neuronal firing in an SCN8A variant associated with developmental delay and epilepsy

Mutations in the SCN8A gene, encoding the voltage-gated sodium channel NaV1.6, are associated with a range of neurodevelopmental syndromes. The p.(Gly1625Arg) (G1625R) mutation was identified in a patient diagnosed with developmental epileptic encephalopathy (DEE). While most of the characterized DEE-associated SCN8A mutations were shown to cause a gain-of-channel function, we show that the G1625R variant, positioned within the S4 segment of domain IV, results in complex effects. Voltage-clamp analyses of NaV1.6G1625R demonstrated a mixture of gain- and loss-of-function properties, including reduced current amplitudes, increased time constant of fast voltage-dependent inactivation, a depolarizing shift in the voltage dependence of activation and inactivation, and increased channel availability with high-frequency repeated depolarization. Current-clamp analyses in transfected cultured neurons revealed that these biophysical properties caused a marked reduction in the number of action potentials when firing was driven by the transfected mutant NaV1.6. Accordingly, computational modeling of mature cortical neurons demonstrated a mild decrease in neuronal firing when mimicking the patients' heterozygous SCN8A expression. Structural modeling of NaV1.6G1625R suggested the formation of a cation-π interaction between R1625 and F1588 within domain IV. Double-mutant cycle analysis revealed that this interaction affects the voltage dependence of inactivation in NaV1.6G1625R. Together, our studies demonstrate that the G1625R variant leads to a complex combination of gain and loss of function biophysical changes that result in an overall mild reduction in neuronal firing, related to the perturbed interaction network within the voltage sensor domain, necessitating personalized multi-tiered analysis for SCN8A mutations for optimal treatment selection.

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.

Preparation and Preliminary Structure-Activity Relationship Studies of Schwarzinicine A Analogs as Vasorelaxant Agents

Schwarzinicines A-D, a series of alkaloids recently discovered from Ficus schwarzii, exhibit pronounced vasorelaxant activity in rat isolated aorta. Building on this finding, a concise synthesis of schwarzinicines A and B has been reported, allowing further investigations into their biological properties. Herein, a preliminary exploration of the chemical space surrounding the structure of schwarzinicine A (1) was carried out aiming to identify structural features that are essential for vasorelaxant activity. A total of 57 analogs were synthesized and tested for vasorelaxant activity in rat isolated aorta. Both efficacy (Emax) and potency (EC50) of these analogs were compared. In addition to identifying structural features that are required for activity or associated with potency enhancement effect, four analogs showed significant potency improvements of up to 40.2-fold when compared to 1. Molecular dynamics simulation of a tetrameric 44-bound transient receptor potential canonical-6 (TRPC6) protein indicated that 44 could potentially form important interactions with the residues Glu509, Asp530, Lys748, Arg758, and Tyr521. These results may serve as a foundation for guiding further structural optimization of the schwarzinicine A scaffold, aiming to discover even more potent analogs.

Regulation of inositol 5-phosphatase activity by the C2 domain of SHIP1 and SHIP2

SHIP1, an inositol 5-phosphatase, plays a central role in cellular signaling. As such, it has been implicated in many conditions. Exploiting SHIP1 as a drug target will require structural knowledge and the design of selective small molecules. We have determined apo, and magnesium and phosphate-bound structures of the phosphatase and C2 domains of SHIP1. The C2 domains of SHIP1 and the related SHIP2 modulate the activity of the phosphatase domain. To understand the mechanism, we performed activity assays, hydrogen-deuterium exchange mass spectrometry, and molecular dynamics on SHIP1 and SHIP2. Our findings demonstrate that the influence of the C2 domain is more pronounced for SHIP2 than SHIP1. We determined 91 structures of SHIP1 with fragments bound, with some near the interface between the two domains. We performed a mass spectrometry screen and determined four structures with covalent fragments. These structures could act as starting points for the development of potent, selective probes.

Computer-assisted discovery and evaluation of potential ribosomal protein S6 kinase beta 2 inhibitors

S6K2 is an important protein in mTOR signaling pathway and cancer. To identify potential S6K2 inhibitors for mTOR pathway treatment, a virtual screening of 1,575,957 active molecules was performed using PLANET, AutoDock GPU, and AutoDock Vina, with their classification abilities compared. The MM/PB(GB)SA method was used to identify four compounds with the strongest binding energies. These compounds were further investigated using molecular dynamics (MD) simulations to understand the properties of the S6K2/ligand complex. Due to a lack of available 3D structures of S6K2, OmegaFold served as a reliable 3D predictive model with higher evaluation scores in SAVES v6.0 than AlphaFold, AlphaFold2, and RoseTTAFold2. The 150 ns MD simulation revealed that the S6K2 structure in aqueous solvation experienced compression during conformational relaxation and encountered potential energy traps of about 19.6 kJ mol-1. The virtual screening results indicated that Lys75 and Lys99 in S6K2 are key binding sites in the binding cavity. Additionally, MD simulations revealed that the ligands remained attached to the activation cavity of S6K2. Among the compounds, compound 1 induced restrictive dissociation of S6K2 in the presence of a flexible region, compound 8 achieved strong stability through hydrogen bonding with Lys99, compound 9 caused S6K2 tightening, and the binding of compound 16 was heavily influenced by hydrophobic interactions. This study suggests that these four potential inhibitors with different mechanisms of action could provide potential therapeutic options.

Analysing the effect caused by increasing the molecular volume in M1-AChR receptor agonists and antagonists: a structural and computational study

M1 muscarinic acetylcholine receptor (M1-AChR), a member of the G protein-coupled receptors (GPCR) family, plays a crucial role in learning and memory, making it an important drug target for Alzheimer's disease (AD) and schizophrenia. M1-AChR activation and deactivation have shown modifying effects in AD and PD preclinical models, respectively. However, understanding the pharmacology associated with M1-AChR activation or deactivation is complex, because of the low selectivity among muscarinic subtypes, hampering their therapeutic applications. In this regard, we constructed two quantitative structure-activity relationship (QSAR) models, one for M1-AChR agonists (total and partial), and the other for the antagonists. The binding mode of 59 structurally different compounds, including agonists and antagonists with experimental binding affinity values (pKi), were analyzed employing computational molecular docking over different structures of M1-AChR. Furthermore, we considered the interaction energy (Einter), the number of rotatable bonds (NRB), and lipophilicity (ilogP) for the construction of the QSAR model for agonists (R2 = 89.64, QLMO2 = 78, and Qext2 = 79.1). For the QSAR model of antagonists (R2 = 88.44, QLMO2 = 82, and Qext2 = 78.1) we considered the Einter, the fraction of sp3 carbons fCsp3, and lipophilicity (MlogP). Our results suggest that the ligand volume is a determinant to establish its biological activity (agonist or antagonist), causing changes in binding energy, and determining the affinity for M1-AChR.

Dynamics of activation in the voltage-sensing domain of Ciona intestinalis phosphatase Ci-VSP

The Ciona intestinalis voltage-sensing phosphatase (Ci-VSP) is a membrane protein containing a voltage-sensing domain (VSD) that is homologous to VSDs from voltage-gated ion channels responsible for cellular excitability. Previously published crystal structures of Ci-VSD in putative resting and active conformations suggested a helical-screw voltage sensing mechanism in which the S4 helix translocates and rotates to enable exchange of salt-bridge partners, but the microscopic details of the transition between the resting and active conformations remained unknown. Here, by combining extensive molecular dynamics simulations with a recently developed computational framework based on dynamical operators, we elucidate the microscopic mechanism of the resting-active transition at physiological membrane potential. Sparse regression reveals a small set of coordinates that distinguish intermediates that are hidden from electrophysiological measurements. The intermediates arise from a noncanonical helical-screw mechanism in which translocation, rotation, and side-chain movement of the S4 helix are only loosely coupled. These results provide insights into existing experimental and computational findings on voltage sensing and suggest ways of further probing its mechanism.

Probing the functional hotspots inside protein hydrophobic pockets by in situ photochemical trifluoromethylation and mass spectrometry

Due to the complex high-order structures and interactions of proteins within an aqueous solution, a majority of chemical functionalizations happen on the hydrophilic sites of protein external surfaces which are naturally exposed to the solution. However, the hydrophobic pockets inside proteins are crucial for ligand binding and function as catalytic centers and transporting tunnels. Herein, we describe a reagent pre-organization and in situ photochemical trifluoromethylation strategy to profile the functional sites inside the hydrophobic pockets of native proteins. Unbiased mass spectrometry profiling was applied for the characterization of trifluoromethylated sites with high sensitivity. Native proteins including myoglobin, trypsin, haloalkane dehalogenase, and human serum albumin have been engaged in this mild photochemical process and substantial hydrophobic site-specific and structure-selective trifluoromethylation substitutes are obtained without significant interference to their bioactivity and structures. Sodium triflinate is the only reagent required to functionalize the unprotected proteins with wide pH-range tolerance and high biocompatibility. This "in-pocket" activation model provides a general strategy to modify the potential binding pockets and gain essential structural insights into the functional hotspots inside protein hydrophobic pockets.

Trimethyllysine Reader Proteins Exhibit Widespread Charge-Agnostic Binding via Different Mechanisms to Cationic and Neutral Ligands

In the last 40 years, cation-π interactions have become part of the lexicon of noncovalent forces that drive protein binding. Indeed, tetraalkylammoniums are universally bound by aromatic cages in proteins, suggesting that cation-π interactions are a privileged mechanism for binding these ligands. A prominent example is the recognition of histone trimethyllysine (Kme3) by the conserved aromatic cage of reader proteins, dictating gene expression. However, two proteins have recently been suggested as possible exceptions to the conventional understanding of tetraalkylammonium recognition. To broadly interrogate the role of cation-π interactions in protein binding interactions, we report the first large-scale comparative evaluation of reader proteins for a neutral Kme3 isostere, experimental and computational mechanistic studies, and structural analysis. We find unexpected widespread binding of readers to a neutral isostere with the first examples of readers that bind the neutral isostere more tightly than Kme3. We find that no single factor dictates the charge selectivity, demonstrating the challenge of predicting such interactions. Further, readers that bind both cationic and neutral ligands differ in mechanism: binding Kme3 via cation-π interactions and the neutral isostere through the hydrophobic effect in the same aromatic cage. This discovery explains apparently contradictory results in previous studies, challenges traditional understanding of molecular recognition of tetraalkylammoniums by aromatic cages in myriad protein-ligand interactions, and establishes a new framework for selective inhibitor design by exploiting differences in charge dependence.

Saturation Transfer Difference NMR and Molecular Docking Interaction Study of Aralkyl-Thiodigalactosides as Potential Inhibitors of the Human-Galectin-3 Protein

Human Galectin-3 (hGal-3) is a protein that selectively binds to β-galactosides and holds diverse roles in both normal and pathological circumstances. Therefore, targeting hGal-3 has become a vibrant area of research in the pharmaceutical chemistry. As a step towards the development of novel hGal-3 inhibitors, we synthesized and investigated derivatives of thiodigalactoside (TDG) modified with different aromatic substituents. Specifically, we describe a high-yielding synthetic route of thiodigalactoside (TDG); an optimized procedure for the synthesis of the novel 3,3'-di-O-(quinoline-2-yl)methyl)-TDG and three other known, symmetric 3,3'-di-O-TDG derivatives ((naphthalene-2yl)methyl, benzyl, (7-methoxy-2H-1-benzopyran-2-on-4-yl)methyl). In the present study, using competition Saturation Transfer Difference (STD) NMR spectroscopy, we determined the dissociation constant (Kd) of the former three TDG derivatives produced to characterize the strength of the interaction with the target protein (hGal-3). Based on the Kd values determined, the (naphthalen-2-yl)methyl, the (quinolin-2-yl)methyl and the benzyl derivatives bind to hGal-3 94, 30 and 24 times more strongly than TDG. Then, we studied the binding modes of the derivatives in silico by molecular docking calculations. Docking poses similar to the canonical binding modes of well-known hGal-3 inhibitors have been found. However, additional binding forces, cation-π interactions between the arginine residues in the binding pocket of the protein and the aromatic groups of the ligands, have been established as significant features. Our results offer a molecular-level understanding of the varying affinities observed among the synthesized thiodigalactoside derivatives, which can be a key aspect in the future development of more effective ligands of hGal-3.

Expression and purification of fluorinated proteins from mammalian suspension culture

The site-specific encoding of noncanonical amino acids allows for the introduction of rationalized chemistry into a target protein. Of the methods that enable this technology, evolved tRNA and synthetase pairs offer the potential for expanded protein production and purification. Such an approach combines the versatility of solid-phase peptide synthesis with the scalable features of recombinant protein production. We describe the large scale production and purification of eukaryotic proteins bearing fluorinated phenylalanine in mammalian suspension cell preparations. Downstream applications of this approach include scalable recombinant protein preparation for ligand binding assays with small molecules and ligands, protein structure determination, and protein stability assays.

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.

Insights into the recognition mechanism in the UBR box of UBR4 for its specific substrates

The N-end rule pathway is a proteolytic system involving the destabilization of N-terminal amino acids, known as N-degrons, which are recognized by N-recognins. Dysregulation of the N-end rule pathway results in the accumulation of undesired proteins, causing various diseases. The E3 ligases of the UBR subfamily recognize and degrade N-degrons through the ubiquitin-proteasome system. Herein, we investigated UBR4, which has a distinct mechanism for recognizing type-2 N-degrons. Structural analysis revealed that the UBR box of UBR4 differs from other UBR boxes in the N-degron binding sites. It recognizes type-2 N-terminal amino acids containing an aromatic ring and type-1 N-terminal arginine through two phenylalanines on its hydrophobic surface. We also characterized the binding mechanism for the second ligand residue. This is the report on the structural basis underlying the recognition of type-2 N-degrons by the UBR box with implications for understanding the N-end rule pathway.

In Silico Molecular Docking and Dynamics Simulation Analysis of Potential Histone Lysine Methyl Transferase Inhibitors for Managing β-Thalassemia

A decreased hemoglobin synthesis is contemplated as a pathological indication of β-thalassemia. Recent studies show that EPZ035544 from Epizyme could induce fetal hemoglobin (HbF) levels due to its proven capability to inhibit euchromatin histone lysine methyl transferase (EHMT2). Therefore, the development of EHMT2 inhibitors is considered promising in managing β-thalassemia. Our strategy to find novel compounds that are EHMT2 inhibitors relies on the virtual screening of ligands that have a structural similarity to N2-[4-methoxy-3-(2,3,4,7-tetrahydro-1H-azepin-5-yl) phenyl]-N4,6-dimethyl-pyrimidine-2,4-diamine (F80) using the PubChem database. In silico docking studies using Autodock Vina were employed to screen a library of 985 compounds and evaluate their binding ability with EHMT2. The selection of hit compounds was based on the docking score and mode of interaction with the protein. The top two ranked compounds were selected for further investigations, including pharmacokinetic properties analysis and molecular dynamics simulations (MDS). Based on the obtained docking score and interaction analysis, N-(4-methoxy-3-methylphenyl)-4,6-diphenylpyrimidin-2-amine (TP1) and 2-N-[4-methoxy-3-(5-methoxy-3H-indol-2-yl)phenyl]-4-N,6-dimethylpyrimidine-2,4-diamine (TP2) were found to be promising candidates, and TP1 exhibited better stability in the MDS study compared to TP2. In summary, our approach helps identify potential EHMT2 inhibitors, and further validation using in vitro and in vivo experiments could certainly enable this molecule to be used as a therapeutic drug in managing β-thalassemia disease.

Amine-recognizing domain in diverse receptors from bacteria and archaea evolved from the universal amino acid sensor

Bacteria possess various receptors that sense different signals and transmit information to enable an optimal adaptation to the environment. A major limitation in microbiology is the lack of information on the signal molecules that activate receptors. Signals recognized by sensor domains are poorly reflected in overall sequence identity, and therefore, the identification of signals from the amino acid sequence of the sensor alone presents a challenge. Biogenic amines are of great physiological importance for microorganisms and humans. They serve as substrates for aerobic and anaerobic growth and play a role of neurotransmitters and osmoprotectants. Here, we report the identification of a sequence motif that is specific for amine-sensing sensor domains that belong to the Cache superfamily of the most abundant extracellular sensors in prokaryotes. We identified approximately 13,000 sensor histidine kinases, chemoreceptors, receptors involved in second messenger homeostasis and Ser/Thr phosphatases from 8,000 bacterial and archaeal species that contain the amine-recognizing motif. The screening of compound libraries and microcalorimetric titrations of selected sensor domains confirmed their ability to specifically bind biogenic amines. Mutants in the amine-binding motif or domains that contain a single mismatch in the binding motif had either no or a largely reduced affinity for amines. We demonstrate that the amine-recognizing domain originated from the universal amino acid-sensing Cache domain, thus providing insight into receptor evolution. Our approach enables precise "wet"-lab experiments to define the function of regulatory systems and therefore holds a strong promise to enable the identification of signals stimulating numerous receptors.

Neurobiochemical characteristics of arginine-rich peptides explain their potential therapeutic efficacy in neurodegenerative diseases

Neurodegenerative diseases, including Alzheimer̕ s disease (AD), Parkinson̕ s disease (PD), Huntington̕ s disease (HD), and Amyotrophic Lateral Sclerosis (ALS) require special attention to find new potential treatment methods. This review aims to summarize the current knowledge of the relationship between the biochemical properties of arginine-rich peptides (ARPs) and their neuroprotective effects to deal with the harmful effects of risk factors. It seems that ARPs have portrayed a promising and fantastic landscape for treating neurodegeneration-associated disorders. With multimodal mechanisms of action, ARPs play various unprecedented roles, including as the novel delivery platforms for entering the central nervous system (CNS), the potent antagonists for calcium influx, the invader molecules for targeting mitochondria, and the protein stabilizers. Interestingly, these peptides inhibit the proteolytic enzymes and block protein aggregation to induce pro-survival signaling pathways. ARPs also serve as the scavengers of toxic molecules and the reducers of oxidative stress agents. They also have anti-inflammatory, antimicrobial, and anti-cancer properties. Moreover, by providing an efficient nucleic acid delivery system, ARPs can play an essential role in developing various fields, including gene vaccines, gene therapy, gene editing, and imaging. ARP agents and ARP/cargo therapeutics can be raised as an emergent class of neurotherapeutics for neurodegeneration. Part of the aim of this review is to present recent advances in treating neurodegenerative diseases using ARPs as an emerging and powerful therapeutic tool. The applications and progress of ARPs-based nucleic acid delivery systems have also been discussed to highlight their usefulness as a broad-acting class of drugs.

AQcalc: A web server that identifies weak molecular interactions in protein structures

Weak molecular interactions play an important role in protein structure and function. Computational tools that identify weak molecular interactions are, therefore, valuable for the study of proteins. Here, we present AQcalc, a web server (https://aqcalcbiocomputing.com/) that can be used to identify anion-quadrupole (AQ) interactions, which are weak interactions involving aromatic residue (Trp, Tyr, and Phe) ring edges and anions (Asp, Glu, and phosphate ion) both within proteins and at their interfaces (protein-protein, protein-nucleic acids, and protein-lipid bilayer). AQcalc identifies AQ interactions as well as clusters involving AQ, cation-π, and salt bridges, among others. Utilizing AQcalc we analyzed weak interactions in protein models, even in the absence of experimental structures, to understand the contributions of weak interactions to deleterious structural changes, including those associated with oncogenic and germline disease variants. We identified several deleterious variants with disrupted AQ interactions (comparable in frequency to cation-π disruptions). Amyloid fibrils utilize AQ to bury anions at frequencies that far exceed those observed for globular proteins. AQ interactions were detected three and five times more frequently than the hydrogen-bonded AQ (HBAQ) in fibril structures and protein-lipid bilayer interfaces, respectively. By contrast, AQ and HBAQ interactions were detected with similar frequencies in globular proteins. Collectively, these findings suggest AQcalc will be effective in facilitating fine structural analysis. As other web utilities designed to identify protein residue interaction networks do not report AQ interactions, wide use of AQcalc will enrich our understanding of residue interaction networks and facilitate hypothesis testing by identifying and experimentally characterizing these comparably weak but important interactions.

Xanomeline displays concomitant orthosteric and allosteric binding modes at the M4 mAChR

The M4 muscarinic acetylcholine receptor (M4 mAChR) has emerged as a drug target of high therapeutic interest due to its expression in regions of the brain involved in the regulation of psychosis, cognition, and addiction. The mAChR agonist, xanomeline, has provided significant improvement in the Positive and Negative Symptom Scale (PANSS) scores in a Phase II clinical trial for the treatment of patients suffering from schizophrenia. Here we report the active state cryo-EM structure of xanomeline bound to the human M4 mAChR in complex with the heterotrimeric Gi1 transducer protein. Unexpectedly, two molecules of xanomeline were found to concomitantly bind to the monomeric M4 mAChR, with one molecule bound in the orthosteric (acetylcholine-binding) site and a second molecule in an extracellular vestibular allosteric site. Molecular dynamic simulations supports the structural findings, and pharmacological validation confirmed that xanomeline acts as a dual orthosteric and allosteric ligand at the human M4 mAChR. These findings provide a basis for further understanding xanomeline's complex pharmacology and highlight the myriad of ways through which clinically relevant ligands can bind to and regulate GPCRs.

Stem Cell Bioengineering with Bioportides: Inhibition of Planarian Head Regeneration with Peptide Mimetics of Eyes Absent Proteins

Djeya1 (RKLAFRYRRIKELYNSYR) is a very effective cell penetrating peptide (CPP) that mimics the α5 helix of the highly conserved Eya domain (ED) of eyes absent (Eya) proteins. The objective of this study was to bioengineer analogues of Djeya1 that, following effective translocation into planarian tissues, would reduce the ability of neoblasts (totipotent stem cells) and their progeny to regenerate the anterior pole in decapitated S. mediterranea. As a strategy to increase the propensity for helix formation, molecular bioengineering of Djeya1 was achieved by the mono-substitution of the helicogenic aminoisobutyric acid (Aib) at three species-variable sites: 10, 13, and 16. CD analyses indicated that Djeya1 is highly helical, and that Aib-substitution had subtle influences upon the secondary structures of bioengineered analogues. Aib-substituted Djeya1 analogues are highly efficient CPPs, devoid of influence upon cell viability or proliferation. All three peptides increase the migration of PC-3 cells, a prostate cancer line that expresses high concentrations of Eya. Two peptides, [Aib13]Djeya1 and [Aib16]Djeya1, are bioportides which delay planarian head regeneration. As neoblasts are the only cell population capable of division in planaria, these data indicate that bioportide technologies could be utilised to directly manipulate other stem cells in situ, thus negating any requirement for genetic manipulation.

EMC chaperone-CaV structure reveals an ion channel assembly intermediate

Voltage-gated ion channels (VGICs) comprise multiple structural units, the assembly of which is required for function1,2. Structural understanding of how VGIC subunits assemble and whether chaperone proteins are required is lacking. High-voltage-activated calcium channels (CaVs)3,4 are paradigmatic multisubunit VGICs whose function and trafficking are powerfully shaped by interactions between pore-forming CaV1 or CaV2 CaVα1 (ref. 3), and the auxiliary CaVβ5 and CaVα2δ subunits6,7. Here we present cryo-electron microscopy structures of human brain and cardiac CaV1.2 bound with CaVβ3 to a chaperone-the endoplasmic reticulum membrane protein complex (EMC)8,9-and of the assembled CaV1.2-CaVβ3-CaVα2δ-1 channel. These structures provide a view of an EMC-client complex and define EMC sites-the transmembrane (TM) and cytoplasmic (Cyto) docks; interaction between these sites and the client channel causes partial extraction of a pore subunit and splays open the CaVα2δ-interaction site. The structures identify the CaVα2δ-binding site for gabapentinoid anti-pain and anti-anxiety drugs6, show that EMC and CaVα2δ interactions with the channel are mutually exclusive, and indicate that EMC-to-CaVα2δ hand-off involves a divalent ion-dependent step and CaV1.2 element ordering. Disruption of the EMC-CaV complex compromises CaV function, suggesting that the EMC functions as a channel holdase that facilitates channel assembly. Together, the structures reveal a CaV assembly intermediate and EMC client-binding sites that could have wide-ranging implications for the biogenesis of VGICs and other membrane proteins.

Dissecting the Self-assembly Dynamics of Imperfect Repeats in α-Synuclein

The pathological aggregation of α-synuclein (αS) into amyloid fibrils is the hallmark of Parkinson's disease (PD). The self-assembly and membrane interactions of αS are mainly governed by the seven imperfect 11-residue repeats of the XKTKEGVXXXX motif around residues 1-95. However, the particular role of each repeat in αS fibrillization remains unclear. To answer this question, we studied the aggregation dynamics of each repeat with up to 10 peptides in silico by conducting multiple independent micro-second atomistic discrete molecular dynamics simulations. Our simulations revealed that only repeats R3 and R6 readily self-assembled into β-sheet-rich oligomers, while the other repeats remained as unstructured monomers with weak self-assembly and β-sheet propensities. The self-assembly process of R3 featured frequent conformational changes with β-sheet formation mainly in the non-conserved hydrophobic tail, whereas R6 spontaneously self-assembled into extended and stable cross-β structures. These results of seven repeats are consistent with their structures and organization in recently solved αS fibrils. As the primary amyloidogenic core, R6 was buried inside the central cross-β core of all αS fibrils, attracting the hydrophobic tails of adjacent R4, R5, and R7 repeats forming β-sheets around R6 in the core. Further away from R6 in the sequence but with a moderate amyloid aggregation propensity, the R3 tail could serve as a secondary amyloidogenic core and form independent β-sheets in the fibril. Overall, our results demonstrate the critical role of R3 and R6 repeats in αS amyloid aggregation and suggest their potential as targets for the peptide-based and small-molecule amyloid inhibitors.

The anaphylatoxin C5a: Structure, function, signaling, physiology, disease, and therapeutics

The complement system is one of the oldest known tightly regulated host defense systems evolved for efficiently functioning cell-based immune systems and antibodies. Essentially, the complement system acts as a pivot between the innate and adaptive arms of the immune system. The complement system collectively represents a cocktail of ∼50 cell-bound/soluble glycoproteins directly involved in controlling infection and inflammation. Activation of the complement cascade generates complement fragments like C3a, C4a, and C5a as anaphylatoxins. C5a is the most potent proinflammatory anaphylatoxin, which is involved in inflammatory signaling in a myriad of tissues. This review provides a comprehensive overview of human C5a in the context of its structure and signaling under several pathophysiological conditions, including the current and future therapeutic applications targeting C5a.