Hydrogen bonds with pi-acceptors in proteins: frequencies and role in stabilizing local 3D structures

Inhibition of NAMPT by PAK4 Inhibitors

The serine/threonine kinase PAK4 plays a crucial role in regulating cell proliferation, survival, migration, and invasion. Overexpression of PAK4 correlates with poor prognosis in some cancers. KPT-9274, a PAK4 inhibitor, significantly reduces the growth of triple-negative breast cancer cells and mammary tumors in mouse models, and it also inhibits the growth of several other types of cancer cells. Interestingly, although it was first identified as a PAK4 inhibitor, KPT-9274 was also found to inhibit the enzyme NAMPT (nicotinamide phosphoribosyltransferase), which is crucial for NAD (nicotinamide adenine dinucleotide) synthesis and vital for cellular energy and growth. These results made us question whether growth inhibition in response to KPT-9274 was due to PAK4 inhibition, NAMPT inhibition, or both. To address this, we tested several other PAK4 inhibitors that also inhibit cell growth, to determine whether they also inhibit NAMPT activity. Our findings confirm that multiple PAK4 inhibitors also inhibit NAMPT activity. This was assessed both in cell-free assays and in a breast cancer cell line. Molecular docking studies were also used to help us better understand the mechanism by which PAK4 inhibitors block PAK4 and NAMPT activity, and we identified specific residues on the PAK4 inhibitors that interact with NAMPT and PAK4. Our results suggest that PAK4 inhibitors may have a more complex mechanism of action than previously understood, necessitating further exploration of how they influence cancer cell growth.

Structural Basis for Antagonist Binding to Vasopressin V1b Receptor Revealed by the Molecular Dynamics Simulations

The human V1b receptor (V1bR) is primarily expressed in the corticotropic cells of the anterior pituitary where it is involved in the regulation of the hypothalamic-pituitary-adrenal (HPA) axis. The activation of V1bR induces the secretion of adrenocorticotropin hormone (ACTH) from the anterior pituitary cells which, in turn, stimulates the production of cortisol via the adrenal cortex. Clinical studies have demonstrated the chronic dysfunction of the HPA axis in patients with several psychiatric disorders. Thus, the inhibition of the V1b receptor and normalizing the HPA axis hyperactivity is a promising approach to the treatment of many stress-related disorders such as anxiety and depression. Nelivaptan is a selective V1bR antagonist that can be used for this purpose and an excellent molecule to study how antagonists interact with V1bR, especially since in recent years the experimental structures of vasopressin V2 and oxytocin receptors were solved, providing high-similarity templates for homology modeling of V1bR. Therefore, in this work, six independent molecular dynamics simulations of a V1bR-nelivaptan complex in a fully hydrated lipid bilayer, yielding a total simulation time of 6.0 μs, have been conducted. In the lowest-energy complexes obtained in this work and proposed to be the most probable structure of the V1bR-nelivaptan complex, the location of the ligand inside the receptor pocket is very similar to that of the other ligands observed in the experimental structures of the vasopressin/oxytocin receptor family. The receptor-ligand interaction has been analyzed and described, revealing the details of the molecular mechanism of this antagonist binding to V1bR and a probable contribution of L2005×40 and T2035×43 to binding selectivity.

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.

Quantification of CH and NH/π-Stacking Interactions in Cells Using Nuclear Magnetic Resonance Spectroscopy

π-Stacking, a type of noncovalent interactions involving aromatic residues, plays an important role in protein folding and function. In this work, an attempt has been made to measure CH/π and NH/π stacking interactions in a protein in Escherichia coli cells using a combined double-mutant cycle and nuclear magnetic resonance spectroscopy method. The results show that the CH/π and NH/π stacking interactions are generally weaker in cells than those in the buffer. The transient intermolecular noncovalent interactions between the protein and the complex cellular environment may compete with and thus weaken the stacking interactions in the protein. The weakening of stacking interactions can enhance the local conformational opening of proteins in E. coli cells. This is evident from the faster rates of amide hydrogen/deuterium exchange observed in cells than in the buffer, for residues that undergo local conformational opening. This study highlights the influence of the cellular environment on π-stacking and the conformational dynamics of proteins.

Robust Prediction of Relative Binding Energies for Protein-Protein Complex Mutations Using Free Energy Perturbation Calculations

Computational free energy-based methods have the potential to significantly improve throughput and decrease costs of protein design efforts. Such methods must reach a high level of reliability, accuracy, and automation to be effectively deployed in practical industrial settings in a way that impacts protein design projects. Here, we present a benchmark study for the calculation of relative changes in protein-protein binding affinity for single point mutations across a variety of systems from the literature, using free energy perturbation (FEP+) calculations. We describe a method for robust treatment of alternate protonation states for titratable amino acids, which yields improved correlation with and reduced error compared to experimental binding free energies. Following careful analysis of the largest outlier cases in our dataset, we assess limitations of the default FEP+ protocols and introduce an automated script which identifies probable outlier cases that may require additional scrutiny and calculates an empirical correction for a subset of charge-related outliers. Through a series of three additional case study systems, we discuss how Protein FEP+ can be applied to real-world protein design projects, and suggest areas of further study.

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.

Differing Effects of Nonlinearity around the Proton Acceptor on CH··O and NH··O H-Bond Strength within Proteins

The effects of deviations from nonlinearity around the carbonyl proton acceptor of an amide group are assessed by DFT quantum chemical calculations for both CH··O and NH··O H-bonds. The proton donors are the imidazole functional group of His and the indole of Trp, which are paired respectively with N-methylacetamide and acetamide. The displacement of either CH or NH group toward the carbonyl O sp2 lone pairs stabilizes the system and strengthens the H-bond. But the two donor groups differ in their response to a shift out of the amide plane. While the NH··O H-bond is weakened by this displacement, a substantial strengthening is observed when the CH donor is moved out of this plane, in one direction versus the other. This pattern is explained on the basis of simple Coulombic considerations.

Catalyst-free selective oxidation of C(sp3)-H bonds in toluene on water

The anisotropic water interfaces provide an environment to drive various chemical reactions not seen in bulk solutions. However, catalytic reactions by the aqueous interfaces are still in their infancy, with the emphasis being on the reaction rate acceleration on water. Here, we report that the oil-water interface activates and oxidizes C(sp3)-H bonds in toluene, yielding benzaldehyde with high selectivity (>99%) and conversion (>99%) under mild, catalyst-free conditions. Collision at the interface between oil-dissolved toluene and hydroxyl radicals spontaneously generated near the water-side interfaces is responsible for the unexpectedly high selectivity. Protrusion of free OH groups from interfacial water destabilizes the transition state of the OH-addition by forming π-hydrogen bonds with toluene, while the H-abstraction remains unchanged to effectively activate C(sp3)-H bonds. Moreover, the exposed free OH groups form hydrogen bonds with the produced benzaldehyde, suppressing it from being overoxidized. Our investigation shows that the oil-water interface has considerable promise for chemoselective redox reactions on water without any catalysts.

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/.

Identification of Novel GANT61 Analogs with Activity in Hedgehog Functional Assays and GLI1-Dependent Cancer Cells

Aberrant activation of hedgehog (Hh) signaling has been implicated in various cancers. Current FDA-approved inhibitors target the seven-transmembrane receptor Smoothened, but resistance to these drugs has been observed. It has been proposed that a more promising strategy to target this pathway is at the GLI1 transcription factor level. GANT61 was the first small molecule identified to directly suppress GLI-mediated activity; however, its development as a potential anti-cancer agent has been hindered by its modest activity and aqueous chemical instability. Our study aimed to identify novel GLI1 inhibitors. JChem searches identified fifty-two compounds similar to GANT61 and its active metabolite, GANT61-D. We combined high-throughput cell-based assays and molecular docking to evaluate these analogs. Five of the fifty-two GANT61 analogs inhibited activity in Hh-responsive C3H10T1/2 and Gli-reporter NIH3T3 cellular assays without cytotoxicity. Two of the GANT61 analogs, BAS 07019774 and Z27610715, reduced Gli1 mRNA expression in C3H10T1/2 cells. Treatment with BAS 07019774 significantly reduced cell viability in Hh-dependent glioblastoma and lung cancer cell lines. Molecular docking indicated that BAS 07019774 is predicted to bind to the ZF4 region of GLI1, potentially interfering with its ability to bind DNA. Our findings show promise in developing more effective and potent GLI inhibitors.

Robust prediction of relative binding energies for protein-protein complex mutations using free energy perturbation calculations

Computational free energy-based methods have the potential to significantly improve throughput and decrease costs of protein design efforts. Such methods must reach a high level of reliability, accuracy, and automation to be effectively deployed in practical industrial settings in a way that impacts protein design projects. Here, we present a benchmark study for the calculation of relative changes in protein-protein binding affinity for single point mutations across a variety of systems from the literature, using free energy perturbation (FEP+) calculations. We describe a method for robust treatment of alternate protonation states for titratable amino acids, which yields improved correlation with and reduced error compared to experimental binding free energies. Following careful analysis of the largest outlier cases in our dataset, we assess limitations of the default FEP+ protocols and introduce an automated script which identifies probable outlier cases that may require additional scrutiny and calculates an empirical correction for a subset of charge-related outliers. Through a series of three additional case study systems, we discuss how protein FEP+ can be applied to real-world protein design projects, and suggest areas of further study.

Luminescent Ruthenium-Terpyridine Complexes Coupled with Stilbene-Appended Naphthalene, Anthracene, and Pyrene Motifs Demonstrate Fluoride Ion Sensing and Reversible Trans-Cis Photoisomerization

A new family of luminescent heteroleptic Ru(II)-terpyridine complexes coupled with stilbene-appended naphthalene, anthracene, and pyrene motifs is reported. Each of the complexes features moderately intense emission at room temperature having a lifetime of 16.7 ns for naphthalene and 11.4 ns for anthracene, while a substantially elevated lifetime of 8.3 μs was observed for the pyrene derivative. All the three complexes display a reversible couple in the positive potential window due to Ru2+/Ru3+ oxidation but multiple reversible and/or quasi-reversible peaks in the negative potential domain because of the reduction of the terpyridine moieties. All the complexes selectively sense F- among the studied anions via the intermediary of different noncovalent interactions. The interaction event is monitored through absorption, emission, and 1H and 19F NMR spectroscopy. Additionally, upon utilizing the stilbene motif, reversible trans-cis isomerization of the complexes has been undertaken upon alternate treatment of visible and UV light so that the complexes can act as potential photomolecular switches. We also carried out the anion sensing characterization of the cis form of the complexes. Theoretical calculation employing density functional theory is also executed for a selective complex (naphthalene derivative) to elucidate different noncovalent interactions that are operative during the complex-fluoride interplay.

Exploring the conformational landscape, hydrogen bonding, and internal dynamics in the diallyl ether and diallyl sulfide monohydrates

The conformational spaces of the diallyl ether (DAE) and diallyl sulfide (DAS) monohydrates were explored using rotational spectroscopy from 6 to 19 GHz. Calculations at the B3LYP-D3(BJ)/aug-cc-pVTZ level suggested significant differences in their conformational behavior, with DAE-w exhibiting 22 unique conformers and DAS-w featuring three stable structures within 6 kJ mol-1. However, only transitions from the lowest energy conformer of each were experimentally observed. Spectral analysis confirmed that binding with water does not alter the conformational preference for the lowest energy structure of the monomers, but it does influence the relative stabilities of all other conformers, particularly in the case of DAE. Non-covalent interaction and quantum theory of atoms in molecules analyses showed that the observed conformer for each complex is stabilized by two intermolecular hydrogen bonds (HBs), where water primarily interacts with the central oxygen or sulfur atom of the diallyl compounds, along with secondary interactions involving the allyl groups. The nature of these interactions was further elucidated using symmetry-adapted perturbation theory, which suggests that the primary HB interaction with S in DAS is weaker and more dispersive in nature compared to the primary HB in DAE. This supports the experimental observation of a tunneling splitting exclusively in the rotational spectrum of DAS-w, as the weaker contact allows water to undergo internal motions within the complex, as shown based on calculated transition state structures for possible tunneling pathways.

Stability of Quadruple Hydrogen Bonds in an Ionic Liquid Environment

Hydrogen bonds (H-bonds) are highly sensitive to the surrounding environments owing to their dipolar nature, with polar solvents kown to significantly weaken H-bonds. Herein, the stability of the H-bonding motif ureidopyrimidinone (UPy) is investigated, embedded into a highly polar polymeric ionic liquid (PIL) consisting of pendant pyrrolidinium bis(trifluoromethylsulfonyl)imide (IL) moieties, to study the influence of such ionic environments on the UPy H-bonds. The content of the surrounding IL is changed by addition of an additional low molecular weight IL to further boost the IL content around the UPy moieties in molar ratios of UPy/IL ranging from 1/4 up to 1/113, thereby promoting the polar microenvironment around the UPy-H-bonds. Variable-temperature solid-state MAS NMR spectroscopy and FT-IR spectroscopy demonstrate that the UPy H-bonds are largely present as (UPy-) dimers, but sensitive to elevated temperatures (>70 °C). Subsequent rheology and DSC studies reveal that the ILs only solvate the polymeric chains but do not interfere with the UPy-dimer H-bonds, thus accounting for their high stability and applicability in many material systems.

A novel fully human anti-NT-ANGPTL3 antibody from phage display library exhibits potent ApoB, TG, and LDL-C lowering activities in hyperlipidemia mice

Dyslipidemia is characterized by elevated plasma levels of low-density lipoprotein cholesterol (LDL-C), triglycerides (TG), and TG-rich lipoprotein (TGRLs) in circulation, and is closely associated with the incidence and development of cardiovascular disease. Angiopoietin-like protein 3 (ANGPTL3) deficiency has been identified as a cause of familial combined hypolipidemia in humans, which allows it to be an important therapeutic target for reducing plasma lipids. Here, we report the discovery and characterization of a novel fully human antibody F1519-D95aA against N-terminal ANGPTL3 (NT-ANGPTL3), which potently inhibits NT-ANGPTL3 with a KD as low as 9.21 nM. In hyperlipidemic mice, F1519-D95aA shows higher apolipoprotein B (ApoB) and TG-lowering, and similar LDL-C reducing activity as compared to positive control Evinacumab (56.50% vs 26.01% decrease in serum ApoB levels, 30.84% vs 25.28% decrease in serum TG levels, 23.32% vs 22.52% decrease in serum LDLC levels, relative to vehicle group). Molecular docking and binding energy calculations reveal that the F1519-D95aA-ANGPTL3 complex (10 hydrogen bonds, -65.51 kcal/mol) is more stable than the Evinacumab-ANGPTL3 complex (4 hydrogen bonds, -63.76 kcal/mol). Importantly, F1519-D95aA binds to ANGPTL3 with different residues in ANGPTL3 from Evinacumab, suggesting that F1519-D95aA may be useful for the treatment of patients resistant to Evinacumab. In conclusion, F1519-D95aA is a novel fully human anti-NT-ANGPTL3 antibody with potent plasma ApoB, TG, and LDL-C lowering activities, which can potentially serve as a therapeutic agent for hyperlipidemia and relevant cardiovascular diseases.

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.

An In Silico Study Investigating Camptothecin-Analog Interaction with Human Protein Tyrosine Phosphatase, SHP2 (PTPN11)

The human PTPN11 gene encodes for the src tyrosine phosphatase protein (SHP2) is now gaining much attention in many disorders, particularly its oncogenic involvement in many types of cancer. Efforts in developing molecules targeting SHP2 with high efficacy are the future of drug discovery and chemotherapy. However, the interaction of a new camptothecin analog with the catalytic domain of SHP2 protein remains unknown. Therefore, this study aims to provide in silico rationale for the recognition and binding of FL118 and irinotecan with the catalytic domain of human protein tyrosine phosphatase-SHP2 (PTPc-SH2-SHP2, chain A). The docking interaction of the human SHP2 protein's catalytic domain as well as Y279C and R465G mutants with FL118 and irinotecan ligands were calculated and analyzed using the Autodock 4.2 programme, setting the docking grid to target the protein's active site. The camptothecin analog FL118 had the best lowest negative affinity energies with PTPc-SHP2 wildtype and SHP2-Y279C mutant model (-7.54 Kcal/mol and -6.94 Kcal/mol, respectively). Moreover, the protein-ligand complexes revealed several hydrogen bond interactions reflecting the degree of stability that each structure possesses, with the FL118-SHP2-wildtype forming the most stable complex among the structures. In addition, the FL118-SHP2 wildtype complex was validated for RMSD, RMSF, hydrogen bonds, and salt bridges. This revealed that the complex generated became stable over time. This in silico rationale identifies the novel FL118 camptothecin analog as a potent selective inhibitor of PTPc-SH2 domain of SHP2 protein, paving way for further in vitro investigations into the interactions and binding activity of analogs with SHP2 for potential therapeutic applications in PTPN11-associated disorders.

Subfunctionalization of a monolignol to a phytoalexin glucosyltransferase is accompanied by substrate inhibition

Uridine diphosphate-dependent glycosyltransferases (UGTs) mediate the glycosylation of plant metabolites, thereby altering their physicochemical properties and bioactivities. Plants possess numerous UGT genes, with the encoded enzymes often glycosylating multiple substrates and some exhibiting substrate inhibition kinetics, but the biological function and molecular basis of these phenomena are not fully understood. The promiscuous monolignol/phytoalexin glycosyltransferase NbUGT72AY1 exhibits substrate inhibition (K_i) at 4 μM scopoletin, whereas the highly homologous monolignol StUGT72AY2 is inhibited at 190 μM. We therefore used hydrogen/deuterium exchange mass spectrometry and structure-based mutational analyses of both proteins and introduced NbUGT72AY1 residues into StUGT72AY2 and vice versa to study promiscuity and substrate inhibition of UGTs. A single F87I and chimeric mutant of NbUGT72AY1 showed significantly reduced scopoletin substrate inhibition, whereas its monolignol glycosylation activity was almost unaffected. Reverse mutations in StUGT72AY2 resulted in increased scopoletin glycosylation, leading to enhanced promiscuity, which was accompanied by substrate inhibition. Studies of 3D structures identified open and closed UGT conformers, allowing visualization of the dynamics of conformational changes that occur during catalysis. Previously postulated substrate access tunnels likely serve as drainage channels. The results suggest a two-site model in which the second substrate molecule binds near the catalytic site and blocks product release. Mutational studies showed that minor changes in amino acid sequence can enhance the promiscuity of the enzyme and add new capabilities such as substrate inhibition without affecting existing functions. The proposed subfunctionalization mechanism of expanded promiscuity may play a role in enzyme evolution and highlights the importance of promiscuous enzymes in providing new functions.

Probing a Hydrogen-π Interaction Involving a Trapped Water Molecule in the Solid State

The detection and characterization of trapped water molecules in chemical entities and biomacromolecules remains a challenging task for solid materials. We herein present proton-detected solid-state Nuclear Magnetic Resonance (NMR) experiments at 100 kHz magic-angle spinning and at high static magnetic-field strengths (28.2 T) enabling the detection of a single water molecule fixed in the calix[4]arene cavity of a lanthanide complex by a combination of three types of non-covalent interactions. The water proton resonances are detected at a chemical-shift value close to zero ppm, which we further confirm by quantum-chemical calculations. Density Functional Theory calculations pinpoint to the sensitivity of the proton chemical-shift value for hydrogen-π interactions. Our study highlights how proton-detected solid-state NMR is turning into the method-of-choice in probing weak non-covalent interactions driving a whole branch of molecular-recognition events in chemistry and biology.

Entropically-Driven Co-assembly of l-Histidine and l-Phenylalanine to Form Supramolecular Materials

Molecular self- and co-assembly allow the formation of diverse and well-defined supramolecular structures with notable physical properties. Among the associating molecules, amino acids are especially attractive due to their inherent biocompatibility and simplicity. The biologically active enantiomer of l-histidine (l-His) plays structural and functional roles in proteins but does not self-assemble to form discrete nanostructures. In order to expand the structural space to include l-His-containing materials, we explored the co-assembly of l-His with all aromatic amino acids, including phenylalanine (Phe), tyrosine (Tyr), and tryptophan (Trp), all in both enantiomeric forms. In contrast to pristine l-His, the combination of this building block with all aromatic amino acids resulted in distinct morphologies including fibers, rods, and flake-like structures. Electrospray ionization mass spectrometry (ESI-MS) indicated the formation of supramolecular co-assemblies in all six combinations, but time-of-flight secondary-ion mass spectrometry (ToF-SIMS) indicated the best seamless co-assembly occurs between l-His and l-Phe while in the other cases, different degrees of phase separation could be observed. Indeed, isothermal titration calorimetry (ITC) suggested the highest affinity between l-His and l-Phe where the formation of co-assembled structures was driven by entropy. In accordance, among all the combinations, the co-assembly of l-His and l-Phe produced single crystals. The structure revealed the formation of a 3D network with nanocavities stabilized by hydrogen bonding between -N (l-His) and -NH (l-Phe). Taken together, using the co-assembly approach we expanded the field of amino acid nanomaterials and showed the ability to obtain discrete supramolecular nanostructures containing l-His based on its specific interactions with l-Phe.

Insights into the interactions between cellulose and biological molecules

Understanding the interactions between carbohydrate polymer molecules and biomolecules is of primary significance for its application. In this paper, the interaction between cellulose and biomolecules was studied using density functional theory method, in which cellobiose, nucleobases, and aromatic amino acids were employed as the structural models of cellulose, DNA, and protein, respectively. Quantitative molecular surface electrostatic potential (ESP) results well represented how cellulose perceived by organism during the recognition. The structural and energetic studies of cellulose with biomolecules complexes show that weak interactions, such as hydrogen bonding interaction, vdW interaction, and pi-H interaction, play an important role in stabilizing these complexes. Through systematic wavefunction analysis, including reduced density gradient (RDG) and natural bond orbital (NBO) methods, the nature of these weak interactions was revealed and further graphically visualized. In-depth understanding of the interaction between cellobiose with biological model molecules may shed lights on the application of carbohydrate polymer-based materials in biological fields.

Cation-π Interactions and Their Role in Assembling Collagen Triple Helices

Cation-π interactions play a significant role in the stabilization of globular proteins. However, their role in collagen triple helices is less well understood and they have rarely been used in de novo designed collagen mimetic systems. In this study, we analyze the stabilizing and destabilizing effects in pairwise amino acid interactions between cationic and aromatic residues in both axial and lateral sequential relationships. Thermal unfolding experiments demonstrated that only axial pairs are stabilizing, while the lateral pairs are uniformly destabilizing. Molecular dynamics simulations show that pairs with an axial relationship can achieve a near-ideal interaction distance, but pairs in a lateral relationship do not. Arginine-π systems were found to be more stabilizing than lysine-π and histidine-π. Arginine-π interactions were then studied in more chemically diverse ABC-type heterotrimeric helices, where arginine-tyrosine pairs were found to form the best helix. This work helps elucidate the role of cation-π interactions in triple helices and illustrates their utility in designing collagen mimetic peptides.

Structural Achievability of an NH-π Interaction between Gln and Phe in a Crystal Structure of a Collagen-like Peptide

NH-π interactions between polar and aromatic residues are well distributed in proteins whose stabilizing effects have been investigated in globular and fibrous proteins. In order to gain structural insights into side chain NH-π interactions, we solved a crystal structure of a collagen-like peptide containing Gln-Phe pairs. The Gln-Phe NH-π interactions were further characterized by quantum calculations, molecular simulations, and structural bioinformatics. The analyses indicated that the NH-π interactions are robust under various solvent conditions, can be distributed either on the protein surface or in its hydrophobic core and can form at a wide range of distances between residues. This study suggested that NH-π interactions can play a versatile role in protein design, including engineering hydrophobic cores, solvent accessible surfaces, and protein-protein interfaces.

Understanding the anomaly of cis-trans isomerism in Pro-His sequence

The anomalous absence of cisPro stabilizing C^αH^α_Xaa···π_Aro interactions at Xaa-Pro-Aro exclusively when Aro is His, is understood by NMR structural analyses of model peptides, as due to i → i backbone-side chain C₆ H-bond that forms uniquely when Aro is His, which significantly decreases its χ₁-g^- population essential for C^αH^α_Xaa···π_Aro formation.

Pressure-Induced Hydrogen Transfer in 2-Butyne via a Double CH···π Aromatic Transition State

Hydrogen transfer (H-transfer) is an important elementary reaction in chemistry and bioscience. It is often facilitated by the hydrogen bonds between the H-donor and acceptor. Here, at room temperature and high pressure, we found that solid 2-butyne experienced a concerted two-in-two-out intermolecular CH···π H-transfer, which initiated the subsequent polymerization. Such double H-transfer goes through an aromatic Hückel six-membered ring intermediate state via intermolecular CH···π interactions enhanced by external pressure. Our work shows that H-transfer can occur via the CH···π route in appropriate conformations under high pressure, which gives important insights into the H-transfer in solid-state hydrocarbons.

Donor-Acceptor Pyridinium Salts for Photo-Induced Electron-Transfer-Driven Modification of Tryptophan in Peptides, Proteins, and Proteomes Using Visible Light

Tryptophan (Trp) plays a variety of critical functional roles in protein biochemistry; however, owing to its low natural frequency and poor nucleophilicity, the design of effective methods for both single protein bioconjugation at Trp as well as for in situ chemoproteomic profiling remains a challenge. Here, we report a method for covalent Trp modification that is suitable for both scenarios by invoking photo-induced electron transfer (PET) as a means of driving efficient reactivity. We have engineered biaryl N-carbamoyl pyridinium salts that possess a donor-acceptor relationship that enables optical triggering with visible light whilst simultaneously attenuating the probe's photo-oxidation potential in order to prevent photodegradation. This probe was assayed against a small bank of eight peptides and proteins, where it was found that micromolar concentrations of the probe and short irradiation times (10-60 min) with violet light enabled efficient reactivity toward surface exposed Trp residues. The carbamate transferring group can be used to transfer useful functional groups to proteins including affinity tags and click handles. DFT calculations and other mechanistic analyses reveal correlations between excited state lifetimes, relative fluorescence quantum yields, and chemical reactivity. Biotinylated and azide-functionalized pyridinium salts were used for Trp profiling in HEK293T lysates and in situ in HEK293T cells using 440 nm LED irradiation. Peptide-level enrichment from live cell labeling experiments identified 290 Trp modifications, with 82% selectivity for Trp modification over other π-amino acids, demonstrating the ability of this method to identify and quantify reactive Trp residues from live cells.

Visualizing protein breathing motions associated with aromatic ring flipping

Aromatic residues cluster in the core of folded proteins, where they stabilize the structure through multiple interactions. Nuclear magnetic resonance (NMR) studies in the 1970s showed that aromatic side chains can undergo ring flips-that is, 180° rotations-despite their role in maintaining the protein fold1-3. It was suggested that large-scale 'breathing' motions of the surrounding protein environment would be necessary to accommodate these ring flipping events1. However, the structural details of these motions have remained unclear. Here we uncover the structural rearrangements that accompany ring flipping of a buried tyrosine residue in an SH3 domain. Using NMR, we show that the tyrosine side chain flips to a low-populated, minor state and, through a proteome-wide sequence analysis, we design mutants that stabilize this state, which allows us to capture its high-resolution structure by X-ray crystallography. A void volume is generated around the tyrosine ring during the structural transition between the major and minor state, and this allows fast flipping to take place. Our results provide structural insights into the protein breathing motions that are associated with ring flipping. More generally, our study has implications for protein design and structure prediction by showing how the local protein environment influences amino acid side chain conformations and vice versa.

Conformational preferences in a series of α-hydroxy ketone derivatives: interplay of conformational energies and lattice cohesive energies

A detailed structural analysis of a series of α-hydroxy ketone derivatives in terms of lattice energies and conformational energies reveals the competing factors of inter- and intramolecular interactions that direct conformational preferences.

N-Substituted Valiolamine Derivatives as Potent Inhibitors of Endoplasmic Reticulum α-Glucosidases I and II with Antiviral Activity

Most enveloped viruses rely on the host cell endoplasmic reticulum (ER) quality control (QC) machinery for proper folding of glycoproteins. The key ER α-glucosidases (α-Glu) I and II of the ERQC machinery are attractive targets for developing broad-spectrum antivirals. Iminosugars based on deoxynojirimycin have been extensively studied as ER α-glucosidase inhibitors; however, other glycomimetic compounds are less established. Accordingly, we synthesized a series of N-substituted derivatives of valiolamine, the iminosugar scaffold of type 2 diabetes drug voglibose. To understand the basis for up to 100,000-fold improved inhibitory potency, we determined high-resolution crystal structures of mouse ER α-GluII in complex with valiolamine and 10 derivatives. The structures revealed extensive interactions with all four α-GluII subsites. We further showed that N-substituted valiolamines were active against dengue virus and SARS-CoV-2 in vitro. This study introduces valiolamine-based inhibitors of the ERQC machinery as candidates for developing potential broad-spectrum therapeutics against the existing and emerging viruses.

Beyond the Molecule: Intermolecular Forces from Gas Liquefaction to X-H⋅⋅⋅π Hydrogen Bonds

Interest toward molecule-molecule interactions developed during the first half of the 19th century with studies on gas liquefaction. In 1869, Andrews carried out the first accurate study on the effects of temperature and pressure on the behaviour of real gases, and in 1873, van der Waals formulated an equation capable of accounting for critical phenomena of each individual gas The nature of the intermolecular forces responsible for aggregation of gases was investigated in the early 20th century by Keesom and Debye (electrostatic interactions between dipoles) and by London (interaction between instantaneous dipoles, studied by quantum mechanics). The hydrogen bond theory, a particular case of dipolar interaction, originated in the 1920s in Berkeley in the institute directed by G. N. Lewis. Later it was made clear that hydrogen bonding is responsible for the aggregation of most biological systems, from proteins to nucleic acids. This Perspective describes the most significant steps through which the science of intermolecular interactions has progressed over the last two centuries, revisiting classical experiments and theoretical formulations, which led to the complex state of art of today.

Platinum Indolylphosphine Fluorido and Polyfluorido Complexes: An Interplay between Cyclometallation, Fluoride Migration, and Hydrogen Bonding

The reaction of [PtCl2 (COD)] (COD=1,5-cyclooctadiene) with diisopropyl-2-(3-methyl)indolylphosphine (iPr2 P(C9 H8 N)) led to the formation of the platinum(ii) chlorido complexes, cis-[PtCl2 {iPr2 P(C9 H8 N)}2 ] (1) and trans-[PtCl2 {iPr2 P(C9 H8 N)}2 ] (2). The cis-complex 1 reacted with NEt3 yielding the complex cis-[PtCl{κ2 -(P,N)-iPr2 P(C9 H7 N)}{iPr2 P(C9 H8 N)}] (3) bearing a cyclometalated κ2 -(P,N)-phosphine ligand, while the isomer 2 with a trans-configuration did not show any reactivity towards NEt3 . Treatment of 1 or 3 with (CH3 )4 NF (TMAF) resulted in the formation of the twofold cyclometalated complex cis-[Pt{κ2 -(P,N)-iPr2 P(C9 H7 N)}2 ] (4). The molecular structures of the complexes 1-4 were determined by single-crystal X-ray diffraction. The fluorido complex cis-[PtF{κ2 -(P,N)-iPr2 P(C9 H7 N)}{iPr2 P(C9 H8 N)}] ⋅ (HF)4 (5 ⋅ (HF)4 ) was formed when complex 4 was treated with different hydrogen fluoride sources. The Pt(ii) fluorido complex 5 ⋅ (HF)4 exhibits intramolecular hydrogen bonding in its outer coordination sphere between the fluorido ligand and the NH group of the 3-methylindolyl moiety. In contrast to its chlorido analogue 3, complex 5 ⋅ (HF)4 reacted with CO or the ynamide 1-(2-phenylethynyl)-2-pyrrolidinone to yield the complexes trans-[Pt(CO){κ2 -(P,C)-iPr2 P(C9 H7 NCO)}{iPr2 P(C9 H8 N)}][F(HF)4 ] (7) and a complex, which we suggest to be cis-[Pt{C=C(Ph)OCN(C3 H6 )}{κ2 -(P,N)-iPr2 P(C9 H7 N)}{iPr2 P(C9 H8 N)}][F(HF)4 ] (9), respectively. The structure of 9 was assigned on the basis of DFT calculations as well as NMR and IR data. Hydrogen bonding of HF and NH to fluoride was proven to be crucial for the existence of 7 and 9.

A Non-Covalent Dimer Formation of Quaternary Ammonium Cation with Unusual Charge Neutralization in Electrospray-Ionization Mass Spectrometry

Specific and nonspecific non-covalent molecular association of biomolecules is characteristic for electrospray-ionization mass spectrometry analysis of biomolecules. Understanding the interaction between two associated molecules is of significance not only from the biological point of view but also gas phase analysis by mass spectrometry. Here we reported a formation of non-covalent dimer of quaternary ammonium denatonium cation with +1 charge detected in the positive ion mode electrospray ionization mass spectrometry analysis of denatonium benzoate. Hydrogen deuterium exchange of amide and carbon-bonded hydrogens revealed that charge neutralization of one denatonium cation is the consequence of amide hydrogen dissociation. DFT (Density Functional Theory) calculations proved high thermodynamic stable of formed dimer stabilized by the short and strong N..H-N hydrogen bond. The signal intensity of the peak characterizing non-covalent dimer is low intensity and does not depend on the sample concentration. Additionally, dimer observation was found to be instrument-dependent. The current investigation is the first experimental and theoretical study on the quaternary ammonium ions dimer. Thus the present study has great significance for understanding the structures of the biomolecules as well as materials.

Unexpectedly rigid short peptide foldamers in which NH-π and CH-π interactions are preserved in solution

NH-π and CH-π interactions, due to their weak character, are not easily identified in solution. We report a group of isolable short peptides with stable folds, in which NH-π and CH-π main chain-side chain interactions can be detected in solution by means of NMR and ATR-IR spectroscopy.

Ultrafast vibrational dynamics of the tyrosine ring mode and its application to enkephalin insertion into phospholipid membranes as probed by two-dimensional infrared spectroscopy

Enkephalins are small opioid peptides whose binding conformations are catalyzed by phospholipid membranes. Binding to opioid receptors is determined by the orientation of tyrosine and phenylalanine side chains. In this work, we investigate the effects of different charged phospholipid headgroups on the insertion of the tyrosine side chain into a lipid bilayer using a combination of 2D IR spectroscopy, anharmonic DFT calculations, and third order response function modeling. The insertion is probed by using the ∼1515 cm^-1 tyrosine ring breathing mode, which we found exhibits rich vibrational dynamics on the picosecond timescale. These dynamics include rapid intramolecular vibrational energy redistribution (IVR), where some of the energy ends up in a dark state that shows up as an anharmonically shifted combination band. The waiting-time dependent 2D IR spectra also show an unusual line shape distortion that affects the extraction of the frequency-frequency correlation function (FFCF), which is the dynamic observable of interest that reflects the tyrosine side chain's insertion into the lipid bilayer. We proposed three models to account for this distortion: a hot-state exchange model, a local environment dependent IVR model, and a coherence transfer model. A qualitative analysis of these models suggests that the local environment dependent IVR rate best explains the line shape distortion, while the coherence transfer model best reproduced the effects on the FFCF. Even with these complex dynamics, we found that the tyrosine ring mode's FFCF is qualitatively correlated with the degree of insertion expected from the different phospholipid headgroups.

Controlling Amphiphilic Polymer Folding beyond the Primary Structure with Protein-Mimetic Di(Phenylalanine)

While methods for polymer synthesis have proliferated, their functionality pales in comparison to natural biopolymers-strategies are limited for building the intricate network of noncovalent interactions necessary to elicit complex, protein-like functions. Using a bioinspired di(phenylalanine) acrylamide (FF) monomer, we explored the impact of various noncovalent interactions in generating ordered assembled structures. Amphiphilic copolymers were synthesized that exhibit β-sheet-like local structure upon collapsing into single-chain assemblies in aqueous environments. Systematic analysis of a series of amphiphilic copolymers illustrated that the global collapse is primarily driven by hydrophobic forces. Hydrogen-bonding and aromatic interactions stabilize local structure, as β-sheet-like interactions were identified via circular dichroism and thioflavin T fluorescence. Similar analysis of phenylalanine (F) and alanine-phenylalanine acrylamide (AF) copolymers found that distancing the aromatic residue from the polymer backbone is sufficient to induce β-sheet-like local structure akin to the FF copolymers; however, the interactions between AF subunits are less stable than those formed by FF. Further, hydrogen-bond donating hydrophilic monomers disrupt internal structure formed by FF within collapsed assemblies. Collectively, these results illuminate design principles for the facile incorporation of multiple facets of protein-mimetic, higher-order structure within folded synthetic polymers.

In silico investigation of potential inhibitors to main protease and spike protein of SARS-CoV-2 in propolis

Docking analysis of propolis's natural compound was successfully performed against SARS-CoV-2 main protease (Mpro) and spike protein subunit 2 (S2). Initially, the propolis's protein was screened using chromatography analysis and successfully identified 22 compounds in the propolis. Four compounds were further investigated, i.e., neoblavaisoflavone, methylophiopogonone A, 3'-Methoxydaidzin, and genistin. The binding affinity of 3'-Methoxydaidzin was -7.7 kcal/mol, which is similar to nelfinavir (control), while the others were -7.6 kcal/mol. However, we found the key residue of Glu A:166 in the methylophiopogonone A and genistin, even though the predicted binding energy slightly higher than nelfinavir. In contrast, the predicted binding affinity of neoblavaisoflavone, methylophiopogonone A, 3'-Methoxydaidzin, and genistin against S2 were -8.1, -8.2, -8.3, and -8.3 kcal/mol, respectively, which is far below of the control (pravastatin, -7.3 kcal/mol). Instead of conventional hydrogen bonding, the π bonding influenced the binding affinity against S2. The results reveal that this is the first report about methylophiopogonone A, 3'-Methoxydaidzin, and genistin as candidates for anti-viral agents. Those compounds can then be further explored and used as a parent backbone molecule to develop a new supplementation for preventing SARS-CoV-2 infections during COVID-19 outbreaks.

Design and synthesis of polyacrylic acid/deoxycholic acid-modified chitosan copolymer and a close inspection of human growth hormone-copolymer interactions: An experimental and computational study

Despite efforts to achieve a long-acting formulation for human growth hormone (hGH), daily injections are still prescribed for children with growth hormone deficiency. To grapple with the issue, acquiring a deep knowledge of the protein and understanding its interaction mechanism with the carrier can be beneficial. Herein, we designed and synthesized a novel chitosan-based copolymer and investigated its interaction with hGH using a combination of experimental and computational strategies. To construct the amphiphilic triblock copolymers (CDP), we grafted deoxycholic acid (DCA) and polyacrylic acid (PAA) onto the chitosan chains, and Fourier-transform infrared (FTIR) analysis confirmed the proper formation of CDP. Circular dichroism (CD) demonstrated the preservation of the secondary structure of hGH interacting with CDP, and, further, fluorescence spectroscopy proved the stability of the tertiary structure of the protein. Applying molecular dynamics simulation (MD), we examined the dynamics and integrity of hGH in the presence of the copolymer and compared its behavior with the protein in aquatic environments. Additionally, energy and contact analysis illustrated that the residues involved in the interaction were located predominantly in the connecting loops, and van der Waals (vdW) and electrostatic interactions were the main driving forces of the polymer-protein complex formation. This research's main aim was to trace the protein-polymer interaction's mechanism. We anticipate that the utility of the copolymer can address the challenges of fabricating a new sustained-release delivery platform for therapeutic proteins.

PEGylation Increases the Strength of a Nearby NH-π Hydrogen Bond in the WW Domain

Here we show that an NH-π interaction between a highly conserved Asn and a nearby Trp stabilizes the WW domain of the human protein Pin1. The strength of this NH-π interaction depends on the structure of the arene, with NH-π interactions involving Trp or naphthylalanine being substantially more stabilizing than those involving Tyr or Phe. Calculations suggest arene size and polarizability are key structural determinants of NH-π interaction strength. Methylation or PEGylation of the Asn side-chain amide nitrogen each strengthens the associated NH-π interaction, though likely for different reasons. We hypothesize that methylation introduces steric clashes that destabilize conformations in which the NH-π interaction is not possible, whereas PEGylation strengthens the NH-π interaction via localized desolvation of the protein surface.

D936Y and Other Mutations in the Fusion Core of the SARS-CoV-2 Spike Protein Heptad Repeat 1: Frequency, Geographical Distribution, and Structural Effect

The crown of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is constituted by its spike (S) glycoprotein. S protein mediates the SARS-CoV-2 entry into the host cells. The “fusion core” of the heptad repeat 1 (HR1) on S plays a crucial role in the virus infectivity, as it is part of a key membrane fusion architecture. While SARS-CoV-2 was becoming a global threat, scientists have been accumulating data on the virus at an impressive pace, both in terms of genomic sequences and of three-dimensional structures. On 15 February 2021, from the SARS-CoV-2 genomic sequences in the GISAID resource, we collected 415,673 complete S protein sequences and identified all the mutations occurring in the HR1 fusion core. This is a 21-residue segment, which, in the post-fusion conformation of the protein, gives many strong interactions with the heptad repeat 2, bringing viral and cellular membranes in proximity for fusion. We investigated the frequency and structural effect of novel mutations accumulated over time in such a crucial region for the virus infectivity. Three mutations were quite frequent, occurring in over 0.1% of the total sequences. These were S929T, D936Y, and S949F, all in the N-terminal half of the HR1 fusion core segment and particularly spread in Europe and USA. The most frequent of them, D936Y, was present in 17% of sequences from Finland and 12% of sequences from Sweden. In the post-fusion conformation of the unmutated S protein, D936 is involved in an inter-monomer salt bridge with R1185. We investigated the effect of the D936Y mutation on the pre-fusion and post-fusion state of the protein by using molecular dynamics, showing how it especially affects the latter one.

Biocomputational Screening of Natural Compounds against Acetylcholinesterase

Alzheimer's disease (AD) is the most common form of dementia and is characterized by irreversible and progressive neurodegeneration. Cholinergic dysfunction has been reported in AD, and several cholinesterase inhibitors, including natural compounds and synthetic analogs, have been developed to treat the disease. However, there is currently no treatment for AD, as most drug-like compounds have failed in clinical trials. Acetylcholinesterase (AChE) is the target of most drugs used commercially to treat AD. This work focused on screening natural compounds obtained from the ZINC database (224, 205 compounds) against AChE to identify those possibly capable of enabling the management of AD. Indirubin and dehydroevodiamine were the best potential AChE inhibitors with free binding energies of -10.03 and -9.00 kcal/mol, respectively. The key residue (His447) of the active site of AChE was found to participate in complex interactions with these two molecules. Six H-bonds were involved in the 'indirubin-AChE' interaction and three H-bonds in the 'dehydroevodiamine-AChE' interaction. These compounds were predicted to cross the blood-brain barrier (BBB) and to exhibit high levels of intestinal absorption. Furthermore, 'indirubin-AChE' and 'dehydroevodiamine-AChE' complexes were found to be stable, as determined by root mean square deviation (RMSD) during a 50 ns molecular dynamics simulation study. Based on the free binding energies and stabilities obtained by simulation studies, we recommend that experimental studies be undertaken on indirubin and dehydroevodiamine with a view towards their potential use as treatments for AD.

Charge-Free, Stabilizing Amide-π Interactions Can Be Used to Control Collagen Triple-Helix Self-Assembly

There is a noted lack of understood, controllable interactions for directing the organization of collagen triple helices. While the field has had success using charge-pair interactions and cation-π interactions in helix design, these alone are not adequate for achieving the degree of specificity desirable for these supramolecular structures. Furthermore, because of the reliance on electrostatic interactions, designed heterotrimeric systems have been heavily charged, a property undesirable in some applications. Amide-π interactions are a comparatively understudied class of charge-free interactions, which could potentially be harnessed for triple-helix design. Herein, we propose, validate, and utilize pairwise amino acid amide-π interactions in collagen triple-helix design. Glutamine-phenylalanine pairs, when arranged in an axial geometry, are found to exhibit a moderately stabilizing effect, while in the lateral geometry, this pair is destabilizing. Together this allows glutamine-phenylalanine pairs to effectively set the register of triple helices. In contrast, interactions between asparagine and phenylalanine appear to have little effect on triple-helical stability. After deconvoluting the contributions of these amino acids to triple-helix stability, we demonstrate these new glutamine-phenylalanine interactions in the successful design of a heterotrimeric triple helix. The results of all of these analyses are used to update our collagen triple-helix thermal stability prediction algorithm, Scoring function for Collagen Emulating Peptides' Temperature of Transition (SCEPTTr).

Accurate prediction for adsorption rate of peptides with high ACE-inhibitory activity from sericin hydrolysate on thiophene hypercross-linked polymer using CoMSIA in 3D-QSAR model

Efficient screening of angiotensin converting enzyme inhibitory (ACE-I) peptides from agricultural or edible sources attract increasing attention. However, their purification process from the complex natural system is still semi-empirical or even uncontrollable, which has seriously reduced their screening efficiency. Herein, inspired by the prediction of ACE-I activity, 3D-QSAR was proposed to predict the adsorption performance of peptides from sericin hydrolysate (SH) having high ACE-I activity on porous hypercross-linked polymers according to their molecular structures. Thiophene hypercross-linked polymer (T-HCP) possessing better screening capacity for ACE-I peptides was chosen as our research object in this work. The sequence and relative adsorption rate of 101 peptides in SH were analyzed by LC-MS and was used as a database to construct the relationship of peptide's chemical structure and adsorption performance on T-HCP by Comparative molecular similarity indices analysis (CoMSIA) from 3D-QSAR. Optimum CoMSIA revealed that enhanced interaction of hydrophobicity and H-bond between T-HCP and the peptide was conducive to increase the adsorption performance of di- to hexa- peptides. Based on these relationships, the adsorption capability of 24 designed peptides with distinguished hydrophobic and H-bond fields was predicted on T-HCP by using optimum CoMSIA and the results of half of these were verified, which showed high consistency with their predicted adsorption rate. Interestingly, these peptides having higher adsorption capacities on T-HCP also possessed higher ACE-I activity. This can be attributed to the high concentration of aromatic surface with π-π interaction and weak-polar CSC group with H-bond interaction on T-HCP material, which is ideal for the selective adsorption of peptides with higher ACE-I activity from SH. This study provides important theoretical guidance for the industrial screening of bio-functional peptides from complex protein mixtures.

Novel diagnostic and therapeutic techniques reveal changed metabolic profiles in recurrent focal segmental glomerulosclerosis

Idiopathic forms of Focal Segmental Glomerulosclerosis (FSGS) are caused by circulating permeability factors, which can lead to early recurrence of FSGS and kidney failure after kidney transplantation. In the past three decades, many research endeavors were undertaken to identify these unknown factors. Even though some potential candidates have been recently discussed in the literature, "the" actual factor remains elusive. Therefore, there is an increased demand in FSGS research for the use of novel technologies that allow us to study FSGS from a yet unexplored angle. Here, we report the successful treatment of recurrent FSGS in a patient after living-related kidney transplantation by removal of circulating factors with CytoSorb apheresis. Interestingly, the classical published circulating factors were all in normal range in this patient but early disease recurrence in the transplant kidney and immediate response to CytoSorb apheresis were still suggestive for pathogenic circulating factors. To proof the functional effects of the patient's serum on podocytes and the glomerular filtration barrier we used a podocyte cell culture model and a proteinuria model in zebrafish to detect pathogenic effects on the podocytes actin cytoskeleton inducing a functional phenotype and podocyte effacement. We then performed Raman spectroscopy in the < 50 kDa serum fraction, on cultured podocytes treated with the FSGS serum and in kidney biopsies of the same patient at the time of transplantation and at the time of disease recurrence. The analysis revealed changes in podocyte metabolome induced by the FSGS serum as well as in focal glomerular and parietal epithelial cell regions in the FSGS biopsy. Several altered Raman spectra were identified in the fractionated serum and metabolome analysis by mass spectrometry detected lipid profiles in the FSGS serum, which were supported by disturbances in the Raman spectra. Our novel innovative analysis reveals changed lipid metabolome profiles associated with idiopathic FSGS that might reflect a new subtype of the disease.

Investigation of the Substrate-Binding Site of a Prostaglandin E Synthase in Bombyx mori

Prostaglandin E synthase (PGES) catalyzes the conversion of prostaglandin H₂ to prostaglandin E₂ in the presence of glutathione (GSH) in mammals. Amid the limited knowledge on prostaglandin and its related enzymes in insects, we recently identified PGES from the silkworm Bombyx mori (bmPGES) and determined its crystal structure complexed with GSH. In the current study, we investigated the substrate-binding site of bmPGES by site-directed mutagenesis and X-ray crystallography. We found that the residues Tyr107, Val155, Met159, and Glu203 are located in the catalytic pockets of bmPGES, and mutagenesis of each residue reduced the bmPGES activity. Our results suggest that these four residues contribute to the catalytic activity of bmPGES. Overall, this structure-function study holds implications in controlling pests by designing rational and efficient pesticides.

When are two hydrogen bonds better than one? Accurate first-principles models explain the balance of hydrogen bond donors and acceptors found in proteins

Hydrogen bonds (HBs) play an essential role in the structure and catalytic action of enzymes, but a complete understanding of HBs in proteins challenges the resolution of modern structural (i.e., X-ray diffraction) techniques and mandates computationally demanding electronic structure methods from correlated wavefunction theory for predictive accuracy. Numerous amino acid sidechains contain functional groups (e.g., hydroxyls in Ser/Thr or Tyr and amides in Asn/Gln) that can act as either HB acceptors or donors (HBA/HBD) and even form simultaneous, ambifunctional HB interactions. To understand the relative energetic benefit of each interaction, we characterize the potential energy surfaces of representative model systems with accurate coupled cluster theory calculations. To reveal the relationship of these energetics to the balance of these interactions in proteins, we curate a set of 4000 HBs, of which >500 are ambifunctional HBs, in high-resolution protein structures. We show that our model systems accurately predict the favored HB structural properties. Differences are apparent in HBA/HBD preference for aromatic Tyr versus aliphatic Ser/Thr hydroxyls because Tyr forms significantly stronger O–H⋯O HBs than N–H⋯O HBs in contrast to comparable strengths of the two for Ser/Thr. Despite this residue-specific distinction, all models of residue pairs indicate an energetic benefit for simultaneous HBA and HBD interactions in an ambifunctional HB. Although the stabilization is less than the additive maximum due both to geometric constraints and many-body electronic effects, a wide range of ambifunctional HB geometries are more favorable than any single HB interaction.

Correlated wavefunction theory predicts and high-resolution crystal structure analysis confirms the important, stabilizing effect of simultaneous hydrogen bond donor and acceptor interactions in proteins.

Joint neutron/X-ray crystal structure of a mechanistically relevant complex of perdeuterated urate oxidase and simulations provide insight into the hydration step of catalysis

Cofactor-independent urate oxidase (UOX) is an ∼137 kDa tetrameric enzyme essential for uric acid (UA) catabolism in many organisms. UA is first oxidized by O2 to de-hydro-isourate (DHU) via a peroxo intermediate. DHU then undergoes hydration to 5-hy-droxy-isourate (5HIU). At different stages of the reaction both catalytic O2 and water occupy the 'peroxo hole' above the organic substrate. Here, high-resolution neutron/X-ray crystallographic analysis at room temperature has been integrated with molecular dynamics simulations to investigate the hydration step of the reaction. The joint neutron/X-ray structure of perdeuterated Aspergillus flavus UOX in complex with its 8-azaxanthine (8AZA) inhibitor shows that the catalytic water molecule (W1) is present in the peroxo hole as neutral H2O, oriented at 45° with respect to the ligand. It is stabilized by Thr57 and Asn254 on different UOX protomers as well as by an O-H⋯π interaction with 8AZA. The active site Lys10-Thr57 dyad features a charged Lys10-NH3 + side chain engaged in a strong hydrogen bond with Thr57OG1, while the Thr57OG1-HG1 bond is rotationally dynamic and oriented toward the π system of the ligand, on average. Our analysis offers support for a mechanism in which W1 performs a nucleophilic attack on DHUC5 with Thr57HG1 central to a Lys10-assisted proton-relay system. Room-temperature crystallography and simulations also reveal conformational heterogeneity for Asn254 that modulates W1 stability in the peroxo hole. This is proposed to be an active mechanism to facilitate W1/O2 exchange during catalysis.

Identification of Persuasive Antiviral Natural Compounds for COVID-19 by Targeting Endoribonuclease NSP15: A Structural-Bioinformatics Approach

SARS-CoV-2 is a positive-stranded RNA virus that bundles its genomic material as messenger-sense RNA in infectious virions and replicates these genomes through RNA intermediates. Several virus-encoded nonstructural proteins play a key role during the viral life cycle. Endoribonuclease NSP15 is vital for the replication and life cycle of the virus, and is thus considered a compelling druggable target. Here, we performed a combination of multiscoring virtual screening and molecular docking of a library of 1624 natural compounds (Nuclei of Bioassays, Ecophysiology and Biosynthesis of Natural Products (NuBBE) database) on the active sites of NSP15 (PDB:6VWW). After sequential high-throughput screening by LibDock and GOLD, docking optimization by CDOCKER, and final scoring by calculating binding energies, top-ranked compounds NuBBE-1970 and NuBBE-242 were further investigated via an indepth molecular-docking and molecular-dynamics simulation of 60 ns, which revealed that the binding of these two compounds with active site residues of NSP15 was sufficiently strong and stable. The findings strongly suggest that further optimization and clinical investigations of these potent compounds may lead to effective SARS-CoV-2 treatment.