Salt bridges: geometrically specific, designable interactions

Understanding the Specific Implications of Amino Acids in the Antibody Development

As the demand for immunotherapy to treat and manage cancers, infectious diseases and other disorders grows, a comprehensive understanding of amino acids and their intricate role in antibody engineering has become a prime requirement. Naturally produced antibodies may not have the most suitable amino acids at the complementarity determining regions (CDR) and framework regions, for therapeutic purposes. Therefore, to enhance the binding affinity and therapeutic properties of an antibody, the specific impact of certain amino acids on the antibody's architecture must be thoroughly studied. In antibody engineering, it is crucial to identify the key amino acid residues that significantly contribute to improving antibody properties. Therapeutic antibodies with higher binding affinity and improved functionality can be achieved through modifications or substitutions with highly suitable amino acid residues. Here, we have indicated the frequency of amino acids and their association with the binding free energy in CDRs. The review also analyzes the experimental outcome of two studies that reveal the frequency of amino acids in CDRs and provides their significant correlation between the outcomes. Additionally, it discusses the various bond interactions within the antibody structure and antigen binding. A detailed understanding of these amino acid properties should assist in the analysis of antibody sequences and structures needed for designing and enhancing the overall performance of therapeutic antibodies.

Development of transgenic Paulownia trees expressing antimicrobial thionin genes for enhanced resistance to fungal infections using chitosan nanoparticles

There is an increasing focus on genetically altering Paulownia trees to enhance their resistance against fungal infections, given their rapid growth and quality wood production. The aim of this research was to establish a technique for incorporating two antimicrobial thionin genes, namely thionin-60 (thio-60) and thionin-63 (thio-63), into Paulownia tomentosa and Paulownia hybrid 9501 through the utilization of chitosan nanoparticles. The outcomes revealed the successful gene transfer into Paulownia trees utilizing chitosan nanoparticles. The effectiveness of thionin proteins against plant pathogens Fusarium and Aspergillus was examined, with a specific focus on Fusarium equiseti due to limited available data. In non-transgenic Paulownia species, the leaf weight inhibition percentage varied from 25 to 36 %, whereas in transgenic species, it ranged from 22 to 7 %. In general, Paulownia species expressing thio-60 displayed increased resistance to F. equiseti, while those expressing thio-63 exhibited heightened resistance to A. niger infection. The thionin proteins displayed a strong affinity for the phospholipid bilayer of the fungal cell membrane, demonstrating their capability to disrupt its structure. The transgenic plants created through this technique showed increased resistance to fungal infections. Thionin-60 demonstrated superior antifungal properties in comparison to thio-63, being more effective at disturbing the fungal cell membrane. These findings indicate that thio-60 holds potential as a novel antifungal agent and presents a promising approach for enhancing the antimicrobial traits of genetically modified Paulownia trees.

Balancing the Affinity and Tumor Cell Binding of a Two-in-One Antibody Simultaneously Targeting EGFR and PD-L1

The optimization of the affinity of monoclonal antibodies is crucial for the development of drug candidates, as it can impact the efficacy of the drug and, thus, the dose and dosing regimen, limit adverse effects, and reduce therapy costs. Here, we present the affinity maturation of an EGFR×PD-L1 Two-in-One antibody for EGFR binding utilizing site-directed mutagenesis and yeast surface display. The isolated antibody variants target EGFR with a 60-fold-improved affinity due to the replacement of a single amino acid in the CDR3 region of the light chain. The binding properties of the Two-in-One variants were confirmed using various methods, including BLI measurements, real-time antigen binding measurements on surfaces with a mixture of both recombinant proteins and cellular binding experiments using flow cytometry as well as real-time interaction cytometry. An AlphaFold-based model predicted that the amino acid exchange of tyrosine to glutamic acid enables the formation of a salt bridge to an arginine at EGFR position 165. This easily adaptable approach provides a strategy for the affinity maturation of bispecific antibodies with respect to the binding of one of the two antigens.

Studying the ssDNA loaded adeno-associated virus aggregation using coarse-grained molecular dynamics simulations

The aggregation of adeno-associated viral (AAV) capsids in an aqueous environment was investigated via coarse-grained molecular dynamics (CG-MD) simulations. The primary driving force and mechanism of the aggregation were investigated with or without single-strand DNA (ssDNA) loaded at various process temperatures. Capsid aggregation appeared to involve multiple residue interactions (i.e., hydrophobic, polar and charged residues) leading to complex protein aggregation. In addition, two aggregation mechanisms (i.e., the fivefold face-to-face contact and the edge-to-edge contact) were identified from this study. The ssDNA with its asymmetric structure could be the reason for destabilizing protein subunits and enhancing the interaction between the charged residues, and further result in the non-reversible face-to-face contact. At higher temperature, the capsid structure was found to be unstable with the significant size expansion of the loaded ssDNA which could be attributed to reduced number of intramolecular hydrogen bonds, the increased conformational deviations of protein subunits and the higher residue fluctuations. The CG-MD model was further validated with previous experimental and simulation data, including the full capsid size measurement and the capsid internal pressure. Thus, a good understanding of AAV capsid aggregation, instability and the role of ssDNA were revealed by applying the developed computational model.

Thermophilic β-mannanases from bacteria: production, resources, structural features and bioengineering strategies

β-mannanases are pivotal enzymes that cleave the mannan backbone to release short chain mannooligosaccharides, which have tremendous biotechnological applications including food/feed, prebiotics and biofuel production. Due to the high temperature conditions in many industrial applications, thermophilic mannanases seem to have great potential to overcome the thermal impediments. Thus, structural analysis of thermostable β-mannanases is extremely important, as it could open up new avenues for genetic engineering, and protein engineering of these enzymes with enhanced properties and catalytic efficiencies. Under this scope, the present review provides a state-of-the-art discussion on the thermophilic β-mannanases from bacterial origin, their production, engineering and structural characterization. It covers broad insights into various molecular biology techniques such as gene mutagenesis, heterologous gene expression, and protein engineering, that are employed to improve the catalytic efficiency and thermostability of bacterial mannanases for potential industrial applications. Further, the bottlenecks associated with mannanase production and process optimization are also discussed. Finally, future research related to bioengineering of mannanases with novel protein expression systems for commercial applications are also elaborated.

Systematic enhancement of protein crystallization efficiency by bulk lysine-to-arginine (KR) substitution

Structural genomics consortia established that protein crystallization is the primary obstacle to structure determination using x-ray crystallography. We previously demonstrated that crystallization propensity is systematically related to primary sequence, and we subsequently performed computational analyses showing that arginine is the most overrepresented amino acid in crystal-packing interfaces in the Protein Data Bank. Given the similar physicochemical characteristics of arginine and lysine, we hypothesized that multiple lysine-to-arginine (KR) substitutions should improve crystallization. To test this hypothesis, we developed software that ranks lysine sites in a target protein based on the redundancy-corrected KR substitution frequency in homologs. This software can be run interactively on the worldwide web at https://www.pxengineering.org/. We demonstrate that three unrelated single-domain proteins can tolerate 5-11 KR substitutions with at most minor destabilization, and, for two of these three proteins, the construct with the largest number of KR substitutions exhibits significantly enhanced crystallization propensity. This approach rapidly produced a 1.9 Å crystal structure of a human protein domain refractory to crystallization with its native sequence. Structures from Bulk KR-substituted domains show the engineered arginine residues frequently make hydrogen-bonds across crystal-packing interfaces. We thus demonstrate that Bulk KR substitution represents a rational and efficient method for probabilistic engineering of protein surface properties to improve crystallization.

Determining the nanostructure and main axis of gly-his-gly fibrils using the amide I' bands in FTIR, VCD, and Raman spectra

The zwitterionic tripeptide glycyl-histidine-glycine (GHG) has been shown to self-assemble into visible crystalline fibrils that form a gel-supporting network with a very high storage modulus. Here we elaborate on the theory and experimental setup behind our novel approach employed to determining the main fibril axis for these gel-forming fibrils by simulating the amide I band profile for infrared absorption (IR), vibrational circular dichroism (VCD), and visible Raman scattering. We also highlight that combining these three vibrational spectroscopies can help in validating structures that are solved using powder x-ray diffraction analysis (PXRD). The PXRD analysis yielded a GHG fibril unit cell with P21 symmetry containing two peptide monomers and two water molecules. The monomers adopt a conformation reminiscent of the distorted polyproline II conformation obtained for tri-lysine in aqueous solution. Stabilization occurs primarily through peptide-peptide intermolecular hydrogen bond interactions, while the role of water in peptide hydration is minimal. The comparison of simulated and experimental amide I' band profiles suggests that the xz plane of the crystal unit cell is being predominantly probed in the experimental IR and VCD spectra, with the x axis of the unit cell pointing in the direction of the main fibril axis. The monomer peptide in the unit cell interacts with six adjacent peptides forming hydrophobic channels by edge-to-face and parallel-displaced ππstacking in the y direction. These cores are further stabilized by a plethora of intermolecular interactions in the x and z directions. Our result suggests that the hydrophobic xz-surfaces would be a good target for the adsorption of hydrophobic drugs.

Molecular Dynamics and Machine Learning reveal distinguishing mechanisms of Competitive Ligands to perturb α,β-Tubulin

The mechanisms of action of ligands competing for the Colchicine Binding Site (CBS) of the α,β-Tubulin are non-standard compared to the commonly witnessed ligand-induced inhibition of proteins. This is because their potencies are not solely judged by the binding affinity itself, but also by their capacity to bias the conformational states of the dimer. Regarding the latter requirement, it is observed that ligands competing for the same pocket that binds colchicine exhibit divergence in potential clinical outcomes. Molecular dynamics-based ∼5.2 µs sampling of α,β-Tubulin complexed with four different ligands has revealed that each ligand has its customized way of influencing the complex. Primarily, it is the proportion of twisting and/or bending characteristic of modes of the intrinsic dynamics which is revealed to be 'fundamental' to tune this variation in the mechanism. The milder influence of 'bending' makes a ligand (TUB092), better classifiable under the group of vascular disrupting agents (VDAs), which are phenotypically effective on cytoskeletons; whereas a stronger impact of 'bending' makes the classical ligand Colchicine (COL) a better Anti-Mitotic agent (AMA). Two other ligands BAL27862 (2RR) and Nocodazole (NZO) fall in the intermediate zone as they fail to explicitly induce bending modes. Random Forest Classification method and K-means Clustering is applied to reveal the efficiency of Machine Learning methods in classifying the Tubulin conformations according to their ligand-specific perturbations and to highlight the significance of specific amino acid residues, mostly positioned in the α-β and β-β interfaces involved in the mechanism. These key residues responsible to yield discriminative actions of the ligands are likely to be highly useful in future endeavours to design more precise inhibitors.

Biophysical chemistry behind sickle cell anemia and the mechanism of voxelotor action

Sickle cell anemia disease has been a great challenge to the world in the present situation. It occurs only due to the polymerization of sickle hemoglobin (HbS) having Pro-Val-Glu typed mutation, while the polymerization does not occur in normal hemoglobin (HbA) having Pro-Glu-Glu peptides. It is also well confirmed that the oxygenated HbS (OHbS) does not participate in the polymerization, while the deoxygenated HbS (dHbS) does, which causes the shape of red blood cells sickled. After polymerization, the blood has a low oxygen affinity. Keeping this fact into consideration, only those drugs are being synthesized that stabilize the OHbS structure so that the polymerization of HbS can be stopped. The literature data showed no systematic description of the changes occurring during the OHbS conversion to dHbS before polymerization. Hence, an innovative reasonable study between HbA and HbS, when they convert into their deoxygenated forms, was done computationally. In this evaluation, physiochemical parameters in HbA/HbS before and after deoxygenation were studied and compared deeply. The computationally collected data was used to understand the abnormal behaviour of dHbS arising due to the replacement of Glu6 with Val6. Consequently, during the presented computational study, the changes occurring in HbS were found opposite/abnormal as compared to HbA after the deoxygenation of both. The mechanism of Voxelotor (GBT-440) action to stop the HbS polymerization was also explained with the help of computationally collected data. Besides, a comparative study between GBT-440 and another suggested drug was also done to know their antisickling strength. Additionally, the effect of pH, CO, CO2, and 2,3-diphosphoglycerate (2,3-DPG) on HbS structure was also studied computationally.

Rational design engineering of a more thermostable Sulfurihydrogenibium yellowstonense carbonic anhydrase for potential application in carbon dioxide capture technologies

Implementing hyperthermostable carbonic anhydrases into CO2 capture and storage technologies in order to increase the rate of CO2 absorption from the industrial flue gases is of great importance from technical and economical points of view. The present study employed a combination of in silico tools to further improve thermostability of a known thermostable carbonic anhydrase from Sulfurihydrogenibium yellowstonense. Experimental results showed that our rationally engineered K100G mutant not only retained the overall structure and catalytic efficiency but also showed a 3 °C increase in the melting temperature and a two-fold improvement in the enzyme half-life at 85 °C. Based on the molecular dynamics simulation results, rearrangement of salt bridges and hydrogen interactions network causes a reduction in local flexibility of the K100G variant. In conclusion, our study demonstrated that thermostability can be improved through imposing local structural rigidity by engineering a single-point mutation on the surface of the enzyme.

The tubulin database: Linking mutations, modifications, ligands and local interactions

Microtubules are polymeric filaments, constructed of α-β tubulin heterodimers that underlie critical subcellular structures in eukaryotic organisms. Four homologous proteins (γ-, δ-, ε- and ζ-tubulin) additionally contribute to specialized microtubule functions. Although there is an immense volume of publicly available data pertaining to tubulins, it is difficult to assimilate all potentially relevant information across diverse organisms, isotypes, and categories of data. We previously assembled an extensive web-based catalogue of published missense mutations to tubulins with >1,500 entries that each document a specific substitution to a discrete tubulin, the species where the mutation was described and the associated phenotype with hyperlinks to the amino acid sequence and citation(s) for research. This report describes a significant update and expansion of our online resource (TubulinDB.bio.uci.edu) to nearly 18,000 entries. It now encompasses a cross-referenced catalog of post-translational modifications (PTMs) to tubulin drawn from public datasets, primary literature, and predictive algorithms. In addition, tubulin protein structures were used to define local interactions with bound ligands (GTP, GDP and diverse microtubule-targeting agents) and amino acids at the intradimer interface, within the microtubule lattice and with associated proteins. To effectively cross-reference these datasets, we established a universal tubulin numbering system to map entries into a common framework that accommodates specific insertions and deletions to tubulins. Indexing and cross-referencing permitted us to discern previously unappreciated patterns. We describe previously unlinked observations of loss of PTM sites in the context of cancer cells and tubulinopathies. Similarly, we expanded the set of clinical substitutions that may compromise MAP or microtubule-motor interactions by collecting tubulin missense mutations that alter amino acids at the interface with dynein and doublecortin. By expanding the database as a curated resource, we hope to relate model organism data to clinical findings of pathogenic tubulin variants. Ultimately, we aim to aid researchers in hypothesis generation and design of studies to dissect tubulin function.

Conformational features and interaction mechanisms of VH H antibodies with β-hairpin CDR3: A case of Nb8-HigB2 interaction

The β-hairpin conformation is regarded as an important basic motif to form and regulate protein-protein interactions. Single-domain VH H antibodies are potential therapeutic and diagnostic tools, and the third complementarity-determining regions of the heavy chains (CDR3s) of these antibodies are critical for antigen recognition. Although the sequences and conformations of the CDR3s are diverse, CDR3s sometimes adopt β-hairpin conformations. However, characteristic features and interaction mechanisms of β-hairpin CDR3s remain to be fully elucidated. In this study, we investigated the molecular recognition of the anti-HigB2 VH H antibody Nb8, which has a CDR3 that forms a β-hairpin conformation. The interaction was analyzed by evaluation of alanine-scanning mutants, molecular dynamics simulations, and hydrogen/deuterium exchange mass spectrometry. These experiments demonstrated that positions 93 and 94 (Chothia numbering) in framework region 3, which is right outside CDR3 by definition, play pivotal roles in maintaining structural stability and binding properties of Nb8. These findings will facilitate the design and optimization of single-domain antibodies.

Multi-epitope vaccine candidates based on mycobacterial membrane protein large (MmpL) proteins against Mycobacterium ulcerans

Buruli ulcer (BU) is a neglected tropical disease. It is caused by the bacterium Mycobacterium ulcerans and is characterized by skin lesions. Several studies were performed testing the Bacillus Calmette-Guérin (BCG) vaccine in human and animal models and M. ulcerans-specific vaccines in animal models. However, there are currently no clinically accepted vaccines to prevent M. ulcerans infection. The aim of this study was to identify T-cell and B-cell epitopes from the mycobacterial membrane protein large (MmpL) proteins of M. ulcerans. These epitopes were analysed for properties including antigenicity, immunogenicity, non-allergenicity, non-toxicity, population coverage and the potential to induce cytokines. The final 8 CD8+, 12 CD4+ T-cell and 5 B-cell epitopes were antigenic, non-allergenic and non-toxic. The estimated global population coverage of the CD8+ and CD4+ epitopes was 97.71%. These epitopes were used to construct five multi-epitope vaccine constructs with different adjuvants and linker combinations. The constructs underwent further structural analyses and refinement. The constructs were then docked with Toll-like receptors. Three of the successfully docked complexes were structurally analysed. Two of the docked complexes successfully underwent molecular dynamics simulations (MDS) and post-MDS analysis. The complexes generated were found to be stable. However, experimental validation of the complexes is required.

Engineering Alkaline-Stable Barley Stripe Mosaic Virus-Like Particles for Efficient Surface Modification

Viruses and virus-like particles are powerful templates for materials synthesis because of their capacity for precise protein engineering and diverse surface functionalization. We recently developed a recombinant bacterial expression system for the production of barley stripe mosaic virus-like particles (BSMV VLPs). However, the applicability of this biotemplate was limited by low stability in alkaline conditions and a lack of chemical handles for ligand attachment. Here, we identify and validate novel residues in the BSMV Caspar carboxylate clusters that mediate virion disassembly through repulsive interactions at high pH. Point mutations of these residues to create attractive interactions that increase rod length ~2 fold, with an average rod length of 91 nm under alkaline conditions. To enable diverse chemical surface functionalization, we also introduce reactive lysine residues at the C-terminus of BSMV coat protein, which is presented on the VLP surface. Chemical conjugation reactions with this lysine proceed more quickly under alkaline conditions. Thus, our alkaline-stable VLP mutants are more suitable for rapid surface functionalization of long nanorods. This work validates novel residues involved in BSMV VLP assembly and demonstrates the feasibility of chemical functionalization of BSMV VLPs for the first time, enabling novel biomedical and chemical applications.

Conservation and divergence of the G-interfaces of Drosophila melanogaster septins

Septins possess a conserved guanine nucleotide-binding (G) domain that participates in the stabilization of organized hetero-oligomeric complexes which assemble into filaments, rings and network-like structures. The fruit fly, Drosophila melanogaster, has five such septin genes encoding Sep1, Sep2, Sep4, Sep5 and Pnut. Here, we report the crystal structure of the heterodimer formed between the G-domains of Sep1 and Sep2, the first from an insect to be described to date. A G-interface stabilizes the dimer (in agreement with the expected arrangement for the Drosophila hexameric particle) and this bears significant resemblance to its human counterparts, even down to the level of individual amino acid interactions. On the other hand, a model for the G-interface formed between the two copies of Pnut which occupy the centre of the hexamer, shows important structural differences, including the loss of a highly favourable bifurcated salt-bridge network. Whereas wild-type Pnut purifies as a monomer, the reintroduction of the salt-bridge network results in stabilizing the dimeric interface in solution as shown by size exclusion chromatography and thermal stability measurements. Adaptive steered molecular dynamics reveals an unzipping mechanism for dimer dissociation which initiates at a point of electrostatic repulsion within the switch II region. Overall, the data contribute to a better understanding of the molecular interactions involved in septin assembly/disassembly.

Molecular Structural Analysis of Porcine CMAH-Native Ligand Complex and High Throughput Virtual Screening to Identify Novel Inhibitors

Porcine meat is the most consumed red meat worldwide. Pigs are also vital tools in biological and medical research. However, xenoreactivity between porcine's N-glycolylneuraminic acid (Neu5Gc) and human anti-Neu5Gc antibodies poses a significant challenge. On the one hand, dietary Neu5Gc intake has been connected to particular human disorders. On the other hand, some pathogens connected to pig diseases have a preference for Neu5Gc. The Cytidine monophospho-N-acetylneuraminic acid hydroxylase (CMAH) catalyses the conversion of N-acetylneuraminic acid (Neu5Ac) to Neu5Gc. In this study, we predicted the tertiary structure of CMAH, performed molecular docking, and analysed the protein-native ligand complex. We performed a virtual screening from a drug library of 5M compounds and selected the two top inhibitors with Vina scores of -9.9 kcal/mol for inhibitor 1 and -9.4 kcal/mol for inhibitor 2. We further analysed their pharmacokinetic and pharmacophoric properties. We conducted stability analyses of the complexes with molecular dynamic simulations of 200 ns and binding free energy calculations. The overall analyses revealed the inhibitors' stable binding, which was further validated by the MMGBSA studies. In conclusion, this result may pave the way for future studies to determine how to inhibit CMAH activities. Further in vitro studies can provide in-depth insight into these compounds' therapeutic potential.

Ins and Outs of Rocker Switch Mechanism in Major Facilitator Superfamily of Transporters

The major facilitator superfamily (MFS) of transporters consists of three classes of membrane transporters: symporters, uniporters, and antiporters. Despite such diverse functions, MFS transporters are believed to undergo similar conformational changes within their distinct transport cycles, known as the rocker-switch mechanism. While the similarities between conformational changes are noteworthy, the differences are also important since they could potentially explain the distinct functions of symporters, uniporters, and antiporters of the MFS superfamily. We reviewed a variety of experimental and computational structural data on a select number of antiporters, symporters, and uniporters from the MFS family to compare the similarities and differences of the conformational dynamics of three different classes of transporters.

Enhanced stability of a disaggregated Aβ fibril on removal of ligand inhibits refibrillation: An all atom Molecular Dynamics simulation study

The extraneuronally deposited senile plaques, composed of neurotoxic aggregates of Aβ fibril, define Alzheimer's disease (AD). Natural compounds have been tested for their destabilization potential on Aβ fibril, thereby curing AD. However, the resultant destabilized Aβ fibril, needs to be checked for its irreversibility to the native organized state after removal of the ligand. Herein, we assessed the stability of a destabilized fibril after the ligand (ellagic acid represented as REF) is removed from the complex. The study has been conducted via Molecular Dynamics (MD) simulation of 1 μs for both Aβ-Water (control) and Aβ-REF″ (test or REF removed) system. The increased value of RMSD, R_g, SASA, lower β-sheet content and reduced number of H-bonds explains enhanced destabilization observed in Aβ-REF″ system. The increased inter-chain distance demonstrates breaking of the residual contacts, testifying the drift of terminal chains from the pentamer. The increased SASA along with the ∆G_ps(polar solvation energy) accounts for the reduced interaction amongst residues, and more with solvent molecules, governing irreversibility to native state. The higher Gibb's free energy of the misaligned structure of Aβ-REF″ ensures irreversibility to the organized structure due to its inability to cross such high energy barrier. The observed stability of the disaggregated structure, despite ligand elimination, establishes the effectiveness of the destabilization technique as a promising therapeutic approach towards treating AD.

The protein family of pyruvate:quinone oxidoreductases: Amino acid sequence conservation and taxonomic distribution

Pyruvate:quinone oxidoreductases (PQOs) catalyse the oxidative decarboxylation of pyruvate to acetate and concomitant reduction of quinone to quinol with the release of CO₂. They are thiamine pyrophosphate (TPP) and flavin-adenine dinucleotide (FAD) containing enzymes, which interact with the membrane in a monotopic way. PQOs are considered as part of alternatives to most recognized pyruvate catabolizing pathways, and little is known about their taxonomic distribution and structural/functional relationship. In this bioinformatics work we tackled these gaps in PQO knowledge. We used the KEGG database to identify PQO coding genes, performed a multiple sequence analysis which allowed us to study the amino acid conservation on these enzymes, and looked at their possible cellular function. We observed that PQOS are enzymes exclusively present in prokaryotes with most of the sequences identified in bacteria. Regarding the amino acid sequence conservation, we found that 75 amino acid residues (out of 570, on average) have a conservation over 90 %, and that the most conserved regions in the protein are observed around the TPP and FAD binding sites. We systematized the presence of conserved features involved in Mg²⁺, TPP and FAD binding, as well as residues directly linked to the catalytic mechanism. We also established the presence of a new motif named "HEH lock", possibly involved in the dimerization process. The results here obtained for the PQO protein family contribute to a better understanding of the biochemistry of these respiratory enzymes.

Structure-guided stabilization improves the ability of the HIV-1 gp41 hydrophobic pocket to elicit neutralizing antibodies

The hydrophobic pocket found in the N-heptad repeat (NHR) region of HIV-1 gp41 is a highly conserved epitope that is the target of various HIV-1-neutralizing monoclonal antibodies. Although the high conservation of the pocket makes it an attractive vaccine candidate, it has been challenging to elicit potent anti-NHR antibodies via immunization. Here, we solved a high-resolution structure of the NHR mimetic IQN17, and, consistent with previous ligand-bound gp41 pocket structures, we observed remarkable conformational plasticity of the pocket. The high malleability of this pocket led us to test whether we could improve the immunogenicity of the gp41 pocket by stabilizing its conformation. We show that the addition of five amino acids at the C terminus of IQN17, to generate IQN22, introduces a stabilizing salt bridge at the base of the peptide that rigidifies the pocket. Mice immunized with IQN22 elicited higher avidity antibodies against the gp41 pocket and a more potent, albeit still weak, neutralizing response against HIV-1 compared with IQN17. Stabilized epitope-focused immunogens could serve as the basis for future HIV-1 fusion-inhibiting vaccines.

Integrated ATAC-Seq and RNA-Seq Data Analysis to Reveal OsbZIP14 Function in Rice in Response to Heat Stress

Transcription factors (TFs) play critical roles in mediating the plant response to various abiotic stresses, particularly heat stress. Plants respond to elevated temperatures by modulating the expression of genes involved in diverse metabolic pathways, a regulatory process primarily governed by multiple TFs in a networked configuration. Many TFs, such as WRKY, MYB, NAC, bZIP, zinc finger protein, AP2/ERF, DREB, ERF, bHLH, and brassinosteroids, are associated with heat shock factor (Hsf) families, and are involved in heat stress tolerance. These TFs hold the potential to control multiple genes, which makes them ideal targets for enhancing the heat stress tolerance of crop plants. Despite their immense importance, only a small number of heat-stress-responsive TFs have been identified in rice. The molecular mechanisms underpinning the role of TFs in rice adaptation to heat stress still need to be researched. This study identified three TF genes, including OsbZIP14, OsMYB2, and OsHSF7, by integrating transcriptomic and epigenetic sequencing data analysis of rice in response to heat stress. Through comprehensive bioinformatics analysis, we demonstrated that OsbZIP14, one of the key heat-responsive TF genes, contained a basic-leucine zipper domain and primarily functioned as a nuclear TF with transcriptional activation capability. By knocking out the OsbZIP14 gene in the rice cultivar Zhonghua 11, we observed that the knockout mutant OsbZIP14 exhibited dwarfism with reduced tiller during the grain-filling stage. Under high-temperature treatment, it was also demonstrated that in the OsbZIP14 mutant, the expression of the OsbZIP58 gene, a key regulator of rice seed storage protein (SSP) accumulation, was upregulated. Furthermore, bimolecular fluorescence complementation (BiFC) experiments uncovered a direct interaction between OsbZIP14 and OsbZIP58. Our results suggested that OsbZIP14 acts as a key TF gene through the concerted action of OsbZIP58 and OsbZIP14 during rice filling under heat stress. These findings provide good candidate genes for genetic improvement of rice but also offer valuable scientific insights into the mechanism of heat tolerance stress in rice.

Ion-Selective Microporous Polymer Membranes with Hydrogen-Bond and Salt-Bridge Networks for Aqueous Organic Redox Flow Batteries

Redox flow batteries (RFBs) have great potential for long-duration grid-scale energy storage. Ion-conducting membranes are a crucial component in RFBs, allowing charge-carrying ions to transport while preventing the cross-mixing of redox couples. Commercial Nafion membranes are widely used in RFBs, but their unsatisfactory ionic and molecular selectivity, as well as high costs, limit the performance and the widespread deployment of this technology. To extend the longevity and reduce the cost of RFB systems, inexpensive ion-selective membranes that concurrently deliver low ionic resistance and high selectivity toward redox-active species are highly desired. Here, high-performance RFB membranes are fabricated from blends of carboxylate- and amidoxime-functionalized polymers of intrinsic microporosity, which exploit the beneficial properties of both polymers. The enthalpy-driven formation of cohesive interchain interactions, including hydrogen bonds and salt bridges, facilitates the microscopic miscibility of the blends, while ionizable functional groups within the sub-nanometer pores allow optimization of membrane ion-transport functions. The resulting microporous membranes demonstrate fast cation conduction with low crossover of redox-active molecular species, enabling improved power ratings and reduced capacity fade in aqueous RFBs using anthraquinone and ferrocyanide as redox couples.

Stability of the Aβ42 Peptide in Mixed Solutions of Denaturants and Proline

Amyloid β-peptide (Aβ) is responsible for the neuronal damage and death of a patient with Alzheimer's disease (AD). Aβ42 oligomeric forms are dominant neurotoxins and are related to neurodegeneration. Their different forms are related to various pathological conditions in the brain. We investigated Aβ42 peptides in different environments of proline, urea, and GdmCl solutions (in pure and mixed binary forms) through atomistic molecular dynamics simulations. Preferential exclusion from the protein surface and facile formation of a large number of weak molecular interactions are the driving forces for the osmolyte's action. We have focused on these interactions between peptide monomers and pure/mixed osmolytes and denaturants. Urea, as usual, denatures the peptide strongly compared to the GdmCl by accumulation around the peptide. GdmCl shows lesser build-up around protein in contrast to urea but is involved in destabilizing the salt bridge formation of Asp23 and Lys28. Proline as an osmolyte protects the peptide from aggregation when mixed with urea and GdmCl solutions. In mixed solutions of two denaturants and osmolyte plus denaturant, the peptide shows enhanced stability as compared to pure denaturant urea solution. The enhanced stability of peptides in proline may be attributed to its exclusion from the peptide surface and favoring salt bridge formation.

Immunoinformatics Study: Multi-Epitope Based Vaccine Design from SARS-CoV-2 Spike Glycoprotein

The coronavirus disease 2019 outbreak has become a huge challenge in the human sector for the past two years. The coronavirus is capable of mutating at a higher rate than other viruses. Thus, an approach for creating an effective vaccine is still needed to induce antibodies against multiple variants with lower side effects. Currently, there is a lack of research on designing a multiepitope of the COVID-19 spike protein for the Indonesian population with comprehensive immunoinformatic analysis. Therefore, this study aimed to design a multiepitope-based vaccine for the Indonesian population using an immunoinformatic approach. This study was conducted using the SARS-CoV-2 spike glycoprotein sequences from Indonesia that were retrieved from the GISAID database. Three SARS-CoV-2 sequences, with IDs of EIJK-61453, UGM0002, and B.1.1.7 were selected. The CD8+ cytotoxic T-cell lymphocyte (CTL) epitope, CD4+ helper T lymphocyte (HTL) epitope, B-cell epitope, and IFN-γ production were predicted. After modeling the vaccines, molecular docking, molecular dynamics, in silico immune simulations, and plasmid vector design were performed. The designed vaccine is antigenic, non-allergenic, non-toxic, capable of inducing IFN-γ with a population reach of 86.29% in Indonesia, and has good stability during molecular dynamics and immune simulation. Hence, this vaccine model is recommended to be investigated for further study.

Structure of Polymer-Grafted Nanocellulose in the Colloidal Dispersion System

Clarifying the primary structure of nanomaterials is invaluable to understand how the nanostructures lead to macroscopic material functions. Nanocellulose is attracting attention as a sustainable building block in materials science. The surface of nanocellulose is often chemically modified by polymer grafting to tune the material properties, such as the viscoelastic properties in rheology modifiers and the reinforcement effect in composites. However, the structure, such as molecular conformation of the grafted polymer and the twist of the core nanocellulose, is not well understood. Here, we investigated the structure of polymer-grafted nanocellulose in the colloidal dispersion system by combining small-angle X-ray scattering measurement and all-atom molecular dynamics simulation. We demonstrated formation of the polymer brush layer on the nanocellulose surface in solvents, which explains the excellent colloidal stability. We also found that twisting of the nanocellulose in the core is suppressed by the existence of the polymer brush layer.

Destabilization of Aβ fibrils by omega-3 polyunsaturated fatty acids: a molecular dynamics study

The senile plaques of neurotoxic aggregates of Aβ protein, deposited extraneuronally, mark the pathological hallmark of Alzheimer's disease (AD). The natural compounds such as omega-3 (ω-3) polyunsaturated fatty acids (PUFAs), which can access blood-brain barrier, are believed to be potential disruptors of preformed Aβ fibrils to cure AD with unknown mechanism. Herein, we present the destabilization potential of three ω-3 PUFAs, viz. Eicosapentaenoic acid (EPA), Docosahexaenoic acid (HXA), and α-linolenic acid (LNL) by molecular dynamics simulation. After an initial testing of 300 ns, EPA and HXA have been considered further for extended production run time, 500 ns. The increased value of root mean square deviation (RMSD), radius of gyration, and solvent-accessible surface area (SASA), the reduced number of H-bonds and β-sheet content, and disruption of salt bridges and hydrophobic contacts establish the binding of these ligands to Aβ fibril leading to destabilization. The polar head was found to interact with positively charged lysine (K28) residue in the fibril. However, the hydrophobicity of the long aliphatic tail competes with the intrinsic hydrophobic interactions of Aβ fibril. This amphiphilic nature of EPA and HXA led to the breaking of inherent hydrophobic contacts and formation of new bonds between the tail of PUFA and hydrophobic residues of Aβ fibril, leading to the destabilization of fibril. The Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) results explain the binding of EPA and HXA to Aβ fibril by interacting with different residues. The destabilization potential of EPA and HXA establishes them as promising drug leads to cure AD, and encourages prospecting of other fatty acids for therapeutic intervention in AD.Communicated by Ramaswamy H. Sarma.

Understanding the Role of Intramolecular Ion-Pair Interactions in Conformational Stability Using an Ab Initio Thermodynamic Cycle

Intramolecular ion-pair interactions yield shape and functionality to many molecules. With proper orientation, these interactions overcome steric factors and are responsible for the compact structures of several peptides. In this study, we present a thermodynamic cycle based on isoelectronic and alchemical mutation to estimate the intramolecular ion-pair interaction energy. We determine these energies for 26 benchmark molecules with common ion-pair combinations and compare them with results obtained using intramolecular symmetry-adapted perturbation theory. For systems with long linkers, the ion-pair energies evaluated using both approaches deviate by less than 2.5% in the vacuum phase. The thermodynamic cycle based on density functional theory facilitates calculations of salt-bridge interactions in model tripeptides with continuum/microsolvation modeling and four large peptides: 1EJG (crambin), 1BDK (bradykinin), 1L2Y (a mini-protein with a tryptophan cage), and 1SCO (a toxin from the scorpion venom).

Acyl-CoA-dependent and acyl-CoA-independent avocado acyltransferases positively influence oleic acid content in nonseed triacylglycerols

In higher plants, acyl-CoA:diacylglycerol acyltransferase (DGAT) and phospholipid:diacylglycerol acyltransferase (PDAT) catalyze the terminal step of triacylglycerol (TAG) synthesis in acyl-CoA-dependent and -independent pathways, respectively. Avocado (Persea americana) mesocarp, a nonseed tissue, accumulates significant amounts of TAG (~70% by dry weight) that is rich in heart-healthy oleic acid (18:1). The oil accumulation stages of avocado mesocarp development coincide with high expression levels for type-1 DGAT (DGAT1) and PDAT1, although type-2 DGAT (DGAT2) expression remains low. The strong preference for oleic acid demonstrated by the avocado mesocarp TAG biosynthetic machinery represents lucrative biotechnological opportunities, yet functional characterization of these three acyltransferases has not been explored to date. We expressed avocado PaDGAT1, PaDGAT2, and PaPDAT1 in bakers' yeast and leaves of Nicotiana benthamiana. PaDGAT1 complemented the TAG biosynthesis deficiency in the quadruple mutant yeast strain H1246, and substantially elevated total cellular lipid content. In vitro enzyme assays showed that PaDGAT1 prefers oleic acid compared to palmitic acid (16:0). Both PaDGAT1 and PaPDAT1 increased the lipid content and elevated oleic acid levels when expressed independently or together, transiently in N. benthamiana leaves. These results indicate that PaDGAT1 and PaPDAT1 prefer oleate-containing substrates, and their coordinated expression likely contributes to sustained TAG synthesis that is enriched in oleic acid. This study establishes a knowledge base for future metabolic engineering studies focused on exploitation of the biochemical properties of PaDGAT1 and PaPDAT1.

A role of salt bridges in mediating drug potency: A lesson from the N-myristoyltransferase inhibitors

The salt bridge is the strongest non-covalent interaction in nature and is known to participate in protein folding, protein-protein interactions, and molecular recognition. However, the role of salt bridges in the context of drug design has remained not well understood. Here, we report that a common feature in the mechanism of inhibition of the N-myristoyltransferases (NMT), promising targets for the treatment of protozoan infections and cancer, is the formation of a salt bridge between a positively charged chemical group of the small molecule and the negatively charged C-terminus of the enzyme. Substituting the inhibitor positively charged amine group with a neutral methylene group prevents the formation of the salt bridge and leads to a dramatic activity loss. Molecular dynamics simulations have revealed that salt bridges stabilize the NMT-ligand complexes by functioning as molecular clips that stabilize the conformation of the protein structure. As such, the creation of salt bridges between the ligands and their protein targets may find an application as a valuable tool in rational drug design.

Dynamic Evolution of Avian RNA Virus Sensors: Repeated Loss of RIG-I and RIPLET

Retinoic acid-inducible gene I (RIG-I) and melanoma differentiation-associated protein 5 (MDA5) are key RNA virus sensors belonging to the RIG-I-like receptor (RLR) family. The activation of the RLR inflammasome leads to the establishment of antiviral state, mainly through interferon-mediated signaling. The evolutionary dynamics of RLRs has been studied mainly in mammals, where rare cases of RLR gene losses were described. By in silico screening of avian genomes, we previously described two independent disruptions of MDA5 in two bird orders. Here, we extend this analysis to approximately 150 avian genomes and report 16 independent evolutionary events of RIG-I inactivation. Interestingly, in almost all cases, these inactivations are coupled with genetic disruptions of RIPLET/RNF135, an ubiquitin ligase RIG-I regulator. Complete absence of any detectable RIG-I sequences is unique to several galliform species, including the domestic chicken (Gallus gallus). We further aimed to determine compensatory evolution of MDA5 in RIG-I-deficient species. While we were unable to show any specific global pattern of adaptive evolution in RIG-I-deficient species, in galliforms, the analyses of positive selection and surface charge distribution support the hypothesis of some compensatory evolution in MDA5 after RIG-I loss. This work highlights the dynamic nature of evolution in bird RNA virus sensors.

Ion-specific binding of cations to the carboxylate and of anions to the amide of alanylalanine

Studies of ion-specific effects on oligopeptides have aided our understanding of Hofmeister effects on proteins, yet the use of different model peptides and different experimental sensitivities have led to conflicting conclusions. To resolve these controversies, we study a small model peptide, L-Alanyl-L-alanine (2Ala), carrying all fundamental chemical protein motifs: C-terminus, amide bond, and N-terminus. We elucidate the effect of GdmCl, LiCl, KCl, KI, and KSCN by combining dielectric relaxation, nuclear magnetic resonance (¹H-NMR), and (two-dimensional) infrared spectroscopy. Our dielectric results show that all ions reduce the rotational mobility of 2Ala, yet the magnitude of the reduction is larger for denaturing cations than for anions. The NMR chemical shifts of the amide group are particularly sensitive to denaturing anions, indicative of anion-amide interactions. Infrared experiments reveal that LiCl alters the spectral homogeneity and dynamics of the carboxylate, but not the amide group. Interaction of LiCl with the negatively charged pole of 2Ala, the COO^- group, can explain the marked cationic effect on dipolar rotation, while interaction of anions between the poles, at the amide, only weakly perturbs dipolar dynamics. As such, our results provide a unifying view on ions' preferential interaction sites at 2Ala and help rationalize Hofmeister effects on proteins.

Enhancement of the polyethylene terephthalate and mono-(2-hydroxyethyl) terephthalate degradation activity of Ideonella sakaiensis PETase by an electrostatic interaction-based strategy

The serious environmental pollution that came up with the continuously growing demand for polyethylene terephthalate (PET) has attracted global concern. The IsPETase which has shown the highest PET degradation activity under ambient temperature is a promising enzyme for PET biodegradation, while poor thermostability limited its practical application. Herein, an electrostatic interaction-based strategy was applied for rational design of IsPETase towards enhanced thermostability. The IsPETaseI139R variant displayed the highest Tm value of 56.4 °C and 3.6-times higher PET degradation activity. Molecular simulations demonstrated that the introduction of salt bridges stabilized the local structures, resulting in robust thermostability. Meanwhile, the IsPETaseS92K/D157E/R251A not only exhibited higher thermostability but also showed a 1.74-fold kcat increase towards mono-(2-hydroxyethyl) terephthalate, which ultimately achieved PET depolymerization to complete monomer TPA. Collectively, the electrostatic interaction-based strategy, together with the derived IsPETase variants, could help promote the bio-recycle of PET, reducing the severe global burden of PET waste.

Analysis of Protein-Protein Interactions for Intermolecular Bond Prediction

Protein-protein interactions often involve a complex system of intermolecular interactions between residues and atoms at the binding site. A comprehensive exploration of these interactions can help reveal key residues involved in protein-protein recognition that are not obvious using other protein analysis techniques. This paper presents and extends DiffBond, a novel method for identifying and classifying intermolecular bonds while applying standard definitions of bonds in chemical literature to explain protein interactions. DiffBond predicted intermolecular bonds from four protein complexes: Barnase-Barstar, Rap1a-raf, SMAD2-SMAD4, and a subset of complexes formed from three-finger toxins and nAChRs. Based on validation through manual literature search and through comparison of two protein complexes from the SKEMPI dataset, DiffBond was able to identify intermolecular ionic bonds and hydrogen bonds with high precision and recall, and identify salt bridges with high precision. DiffBond predictions on bond existence were also strongly correlated with observations of Gibbs free energy change and electrostatic complementarity in mutational experiments. DiffBond can be a powerful tool for predicting and characterizing influential residues in protein-protein interactions, and its predictions can support research in mutational experiments and drug design.

Machine Learning Enabled Structure-Based Drug Repurposing Approach to Identify Potential CYP1B1 Inhibitors

Drug-metabolizing enzyme (DME)-mediated pharmacokinetic resistance of some clinically approved anticancer agents is one of the main reasons for cancer treatment failure. In particular, some commonly used anticancer medicines, including docetaxel, tamoxifen, imatinib, cisplatin, and paclitaxel, are inactivated by CYP1B1. Currently, no approved drugs are available to treat this CYP1B1-mediated inactivation, making the pharmaceutical industries strive to discover new anticancer agents. Because of the extreme complexity and high risk in drug discovery and development, it is worthwhile to come up with a drug repurposing strategy that may solve the resistance problem of existing chemotherapeutics. Therefore, in the current study, a drug repurposing strategy was implemented to find the possible CYP1B1 inhibitors using machine learning (ML) and structure-based virtual screening (SB-VS) approaches. Initially, three different ML models were developed such as support vector machines (SVMs), random forest (RF), and artificial neural network (ANN); subsequently, the best-selected ML model was employed for virtual screening of the selleckchem database to identify potential CYP1B1 inhibitors. The inhibition potency of the obtained hits was judged by analyzing the crucial active site amino acid interactions against CYP1B1. After a thorough assessment of docking scores, binding affinities, as well as binding modes, four compounds were selected and further subjected to in vitro analysis. From the in vitro analysis, it was observed that chlorprothixene, nadifloxacin, and ticagrelor showed promising inhibitory activity toward CYP1B1 in the IC₅₀ range of 0.07-3.00 μM. These new chemical scaffolds can be explored as adjuvant therapies to address CYP1B1-mediated drug-resistance problems.

Not sorcery after all: Roles of multiple charged residues in membrane insertion of gasdermin-A3

Gasdermins execute programmatory cell death, known as pyroptosis, by forming medium-sized membrane pores. Recently, the molecular structure of those pores as well as the diversity in their shape and size have been revealed by cryoTEM and atomic force microscopy, respectively. Even though a growth of smaller to larger oligomers and reshaping from slits to rings could be documented, the initiation of the gasdermin pore formation remains a mystery. In one hypothesis, gasdermin monomers insert into membranes before associating into oligomeric pores. In the other hypothesis, gasdermin oligomers preassemble on the membrane surface prior to membrane insertion. Here, by studying the behavior of monomeric membrane-inserted gasdermin-A3 (GSDMA3), we unveil that a monomeric gasdermin prefers the membrane-adsorbed over the membrane-inserted state. Our results thus support the hypothesis of oligomers preassembling on the membrane surface before membrane penetration. At the same time, our simulations of small membrane-inserted arcs of GSDMA3 suggest that the inserting oligomer can be small and does not have to comprise a full ring of approximately 26–30 subunits. Moreover, our simulations have revealed an astonishingly large impact of salt-bridge formation and protein surroundings on the transmembrane passage of charged residues, reducing the energetic cost by up to 53% as compared to their free forms. The here observed free energy barrier of mere 5.6 kcal/mol for the membrane insertion of monomeric GSDMA3 explains the surprising ability of gasdermins to spontaneously self-insert into cellular membranes.

Evaluation of immune evasion in SARS-CoV-2 Delta and Omicron variants

Emerging SARS-CoV-2 variants with higher transmissibility and immune escape remain a persistent threat across the globe. This is evident from the recent outbreaks of the Delta (B.1.617.2) and Omicron variants. These variants have originated from different continents and spread across the globe. In this study, we explored the genomic and structural basis of these variants for their lineage defining mutations of the spike protein through computational analysis, protein modeling, and molecular dynamic (MD) simulations. We further experimentally validated the importance of these deletion mutants for their immune escape using a pseudovirus-based neutralization assay, and an antibody (4A8) that binds directly to the spike protein's NTD. Delta variant with the deletion and mutations in the NTD revealed a better rigidity and reduced flexibility as compared to the wild-type spike protein (Wuhan isolate). Furthermore, computational studies of 4A8 monoclonal antibody (mAb) revealed a reduced binding of Delta variant compared to the wild-type strain. Similarly, the MD simulation data and virus neutralization assays revealed that the Omicron also exhibits immune escape, as antigenic beta-sheets appear to be disrupted. The results of the present study demonstrate the higher possibility of immune escape and thereby achieved better fitness advantages by the Delta and Omicron variants, which warrants further demonstrations through experimental evidences. Our study, based on in-silico computational modelling, simulations, and pseudovirus-based neutralization assay, highlighted and identified the probable mechanism through which the Delta and Omicron variants are more pathogenically evolved with higher transmissibility as compared to the wild-type strain.

Accurate Description of Solvent-Exposed Salt Bridges with a Non-polarizable Force Field Incorporating Solvent Effects

The strength of salt bridges resulting from the interaction of cations and anions is modulated by their environment. However, polarization of the solvent molecules by the charged moieties makes the accurate description of cation-anion interactions in an aqueous solution by means of a pairwise additive potential energy function and classical combination rules particularly challenging. In this contribution, aiming at improving the representation of solvent-exposed salt-bridge interactions with an all-atom non-polarizable force field, we put forth here a parametrization strategy. First, the interaction of a cation and an anion is characterized by hybrid quantum mechanical/molecular mechanics (QM/MM) potential of mean force (PMF) calculations, whereby constantly exchanging solvent molecules around the ions are treated at the quantum mechanical level. The Lennard-Jones (LJ) parameters describing the salt-bridge ion pairs are then optimized to match the reference QM/MM PMFs through the so-called nonbonded FIX, or NBFIX, feature of the CHARMM force field. We apply the new set of parameters, coined CHARMM36m-SBFIX, to the calculation of association constants for the ammonium-acetate and guanidinium-acetate complexes, the osmotic pressures for glycine zwitterions, guanidinium, and acetate ions, and to the simulation of both folded and intrinsically disordered proteins. Our findings indicate that CHARMM36m-SBFIX improves the description of solvent-exposed salt-bridge interactions, both structurally and thermodynamically. However, application of this force field to the standard binding free-energy calculation of a protein-ligand complex featuring solvent-excluded salt-bridge interactions leads to a poor reproduction of the experimental value, suggesting that the parameters optimized in an aqueous solution cannot be readily transferred to describe solvent-excluded salt-bridge interactions. Put together, owing to their sensitivity to the environment, modeling salt-bridge interactions by means of a single, universal set of LJ parameters remains a daunting theoretical challenge.

Improvement on Thermostability of Pectate Lyase and Its Potential Application to Ramie Degumming

In order to obtain a thermostable pectate lyase for ramie degumming, a rational design based on structural analysis was carried out on a novel pectate lyase (Pel419) derived from the Dickeya Dadantii DCE-01 for high-efficiency ramie degumming. A total of five potential amino acid sites were chosen to replace residues. Then, the mutant enzymes were subjected to the heterologous expressions in Escherichia coli and their enzymatic characteristics were determined. The optimal reaction temperature for the five mutants kept consistent with that for the wild type. The enzyme activity and thermal stability of mutant V52A were significantly improved. Meanwhile, the weight loss rate obtained by V52A with the best enzymatic characteristics in the ramie degumming process at 50 °C is comparable with that obtained by commercial cotton-ramie processing pectinases, indicating that V52A was a potential industrial enzyme that could be applied to large-scale ramie degumming. In this study, the biological functions of conservative residues of Pel419 were preliminarily explored. The mutant V52A with both enzymatic activity and improved heat resistance was acquired, providing a superior material for developing enzyme preparations of ramie degumming, and rendering an effective method for the rational design aiming to improve the thermostability of pectate lyase.

SpySwitch enables pH- or heat-responsive capture and release for plug-and-display nanoassembly

Proteins can be empowered via SpyTag for anchoring and nanoassembly, through covalent bonding to SpyCatcher partners. Here we generate a switchable version of SpyCatcher, allowing gentle purification of SpyTagged proteins. We introduce numerous histidines adjacent to SpyTag’s binding site, giving moderate pH-dependent release. After phage-based selection, our final SpySwitch allows purification of SpyTag- and SpyTag003-fusions from bacterial or mammalian culture by capture at neutral pH and release at pH 5, with purity far beyond His-tag methods. SpySwitch is also thermosensitive, capturing at 4 °C and releasing at 37 °C. With flexible choice of eluent, SpySwitch-purified proteins can directly assemble onto multimeric scaffolds. 60-mer multimerization enhances immunogenicity and we use SpySwitch to purify receptor-binding domains from SARS-CoV-2 and 11 other sarbecoviruses. For these receptor-binding domains we determine thermal resilience (for mosaic vaccine development) and cross-recognition by antibodies. Antibody EY6A reacts across all tested sarbecoviruses, towards potential application against new coronavirus pandemic threats.

The SpyCatcher-SpyTag system allows protein anchoring and nanoassembly. Here, the authors engineer SpySwitch, a dually switchable Catcher which allows gentle purification of SpyTagged proteins prior to downstream applications such as the assembly of virus-like particles.

Accessibility of ENaC extracellular domain central core residues

The epithelial Na+ channel (ENaC)/degenerin family has a similar extracellular architecture, where specific regulatory factors interact and alter channel gating behavior. The extracellular palm domain serves as a key link to the channel pore. In this study, we used cysteine-scanning mutagenesis to assess the functional effects of Cys-modifying reagents on palm domain β10 strand residues in mouse ENaC. Of the 13 ENaC α subunit mutants with Cys substitutions examined, only mutants at sites in the proximal region of β10 exhibited changes in channel activity in response to methanethiosulfonate reagents. Additionally, Cys substitutions at three proximal sites of β and γ subunit β10 strands also rendered mutant channels methanethiosulfonate-responsive. Moreover, multiple Cys mutants were activated by low concentrations of thiophilic Cd2+. Using the Na+ self-inhibition response to assess ENaC gating behavior, we identified four α, two β, and two γ subunit β10 strand mutations that changed the Na+ self-inhibition response. Our results suggest that the proximal regions of β10 strands in all three subunits are accessible to small aqueous compounds and Cd2+ and have a role in modulating ENaC gating. These results are consistent with a structural model of mouse ENaC that predicts the presence of aqueous tunnels adjacent to the proximal part of β10 and with previously resolved structures of a related family member where palm domain structural transitions were observed with channels in an open or closed state.

Triptorelin-functionalized PEG-coated biosynthesized gold nanoparticles: Effects of receptor-ligand interactions on adhesion to triple negative breast cancer cells

This paper presents the results of an experimental and computational study of the adhesion of triptorelin-conjugated PEG-coated biosynthesized gold nanoparticles (GNP-PEG-TRP) to triple-negative breast cancer (TNBC) cells. The adhesion is studied at the nanoscale using a combination of atomic force microscopy (AFM) experiments and molecular dynamics (MD) simulations. The AFM measurements showed that the triptorelin-functionalized gold nanoparticles (GNP-TRP and GNP-PEG-TRP) have higher adhesion to triple-negative breast cancer cells (TNBC) than non-tumorigenic breast cells. The increased adhesion of GNP-TRP and GNP-PEG-TRP to TNBC is also attributed to the overexpression of LHRH receptors on the surfaces of both TNBC. Finally, the molecular dynamics model reveals insights into the effects of receptor density, molecular configuration, and receptor-ligand docking characteristics on the interactions of triptorelin-functionalized PEG-coated gold nanoparticles with TNBC. A three to nine-fold increase in the adhesion is predicted between triptorelin-functionalized PEG-coated gold nanoparticles and TNBC cells. The implications of the results are then discussed for the specific targeting of TNBC.

Structural basis of TRPV5 regulation by physiological and pathophysiological modulators

Transient receptor potential vanilloid 5 (TRPV5) is a kidney-specific Ca2+-selective ion channel that plays a key role in Ca2+ homeostasis. The basal activity of TRPV5 is balanced through activation by phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) and inhibition by Ca2+-bound calmodulin (CaM). Parathyroid hormone (PTH), the key extrinsic regulator of Ca2+ homeostasis, increases the activity of TRPV5 via protein kinase A (PKA)-mediated phosphorylation. Metabolic acidosis leads to reduced TRPV5 activity independent of PTH, causing hypercalciuria. Using cryoelectron microscopy (cryo-EM), we show that low pH inhibits TRPV5 by precluding PI(4,5)P2 activation. We capture intermediate conformations at low pH, revealing a transition from open to closed state. In addition, we demonstrate that PI(4,5)P2 is the primary modulator of channel gating, yet PKA controls TRPV5 activity by preventing CaM binding and channel inactivation. Our data provide detailed molecular mechanisms for regulation of TRPV5 by two key extrinsic modulators, low pH and PKA.

Molecular Dynamics and MM-PBSA Analysis of the SARS-CoV-2 Gamma Variant in Complex with the hACE-2 Receptor

The SARS-CoV-2 virus, since its appearance in 2019, has caused millions of cases and deaths. To date, there is no effective treatment or a vaccine that is fully protective. Despite the efforts made by governments and health institutions around the globe to control its propagation, the evolution of the virus has accelerated, diverging into hundreds of variants. However, not all of them are variants of concern (VoC’s). VoC’s have appeared in different regions and throughout the two years of the pandemic they have spread around the world. Specifically, in South America, the gamma variant (previously known as P.1) appeared in early 2021, bringing with it a second wave of infections. This variant contains the N501Y, E484K and K417T mutations in the receptor binding domain (RBD) of the spike protein. Although these mutations have been described experimentally, there is still no clarity regarding their role in the stabilization of the complex with the human angiotensin converting enzyme 2 (hACE-2) receptor. In this article we dissect the influence of mutations on the interaction with the hACE-2 receptor using molecular dynamics and estimations of binding affinity through a screened version of the molecular mechanics Poisson Boltzmann surface area (MM-PBSA) and interaction entropy. Our results indicate that mutations E484K and K417T compensate each other in terms of binding affinity, while the mutation N501Y promotes a more convoluted effect. This effect consists in the adoption of a cis configuration in the backbone of residue Y495 within the RBD, which in turn promotes polar interactions with the hACE-2 receptor. These results not only correlate with experimental observations and complement previous knowledge, but also expose new features associated with the specific contribution of concerned mutations. Additionally, we propose a recipe to assess the residue-specific contribution to the interaction entropy.

Observation of Arginine Side-Chain Motions Coupled to the Global Conformational Exchange Process in Deubiquitinase A

Coupled motions have been demonstrated to be functionally important in a number of enzymes. Noncovalent side-chain interactions play essential roles in coordinating the motions across different structural elements in a protein. However, most of the dynamic studies of proteins are focused on backbone amides or methyl groups in the side chains and little is known about the polar and charged side chains. We have previously characterized the conformational dynamics of deubiquitinase A (DUBA), an isopeptidase, on the microsecond-to-millisecond (μs-ms) time scales with the amide ¹H Carr-Purcell-Meiboom-Gill (CPMG) experiment. We detected a global conformational exchange process on a time scale of approximately 200 μs, which involves most of the structural elements in DUBA, including the active site and the substrate binding interface. Here, we extend our previous study on backbone amides to the arginine side-chain N^ε-H^ε groups using a modified ¹H CPMG pulse sequence that can efficiently detect both backbone amide and arginine side-chain N^ε-H^ε signals in a single experiment. We found that the side chains of three arginines display motions on the same time scale as the backbone amides. Mutations of two of the three arginines to alanines result in a decrease in enzyme activity. One of these two arginines is located in a loop involved in substrate binding. This loop is not visible in the backbone amide-detected experiments due to excess line broadening induced by motions on the μs-ms time scales. These results clearly demonstrate that the motions of some arginine side chains are coupled to the global conformational exchange process and provide an additional probe for motions in a functionally important loop that did not yield visible backbone amide signals, suggesting the value of side-chain experiments on DUBA. The modified ¹H CPMG pulse sequence allows the simultaneous characterization of backbone and arginine side-chain dynamics without any increase in data acquisition time and can be applied to the dynamic studies of any protein that displays measurable amide ¹H relaxation dispersion.

A conserved sequence in the small intracellular loop of tetraspanins forms an M-shaped inter-helix turn

Tetraspanins are a family of small proteins with four transmembrane segments (TMSs) playing multiple roles in human physiology. Nevertheless, we know little about the factors determining their structure. In the study at hand, we focus on the small intracellular loop (SIL) between TMS2 and TMS3. There we have identified a conserved five amino acid core region with three charged residues forming an M-shaped backbone, which we call M-motif. The M´s plane runs parallel to the membrane surface and the central amino acid constitutes the inter-helix turning point. At the second position of the M-motif, in tetraspanin crystal structures we identified a glutamate oriented towards a lysine in the juxtamembrane region of TMS1. Using Tspan17 as example, we find that by mutating either the glutamate or juxtamembrane-lysine, but not upon glutamate/lysine swapping, expression level, maturation and ER-exit are reduced. We conclude that the SIL is more than a short linking segment but propose it is involved in shaping the tertiary structure of tetraspanins.

Potential effects and relevant lead compounds of Vigna mungo (L.) Hepper seeds against bacterial infection, helminthiasis, thrombosis and neuropharmacological disorders

Multidrug-resistant bacterial infections, helminthiasis, thrombosis, anxiety and insomnia are some of the major global health concerns. Vigna mungo (L.) Hepper (VM) has been used traditionally to treat microbial infection, helminthic disorder, schizophrenia, memory loss, and blood circulatory problem. This research aims to discover antibacterial, anthelmintic, thrombolytic and neuropharmacological effects of the methanol extract of Vigna mungo seeds (MESVM), and also in-silico prediction of relevant lead compounds by molecular docking and ADME/T analysis. The crude extracts and subsequent fractions of MESVM were investigated for antibacterial activity by disc diffusion method, anthelmintic activity by paralysis and death test on earthworms, and thrombolytic activity by in vitro blood clot dissolution test. Open-field test and elevated plus maze test were performed for evaluating anxiolytic activity of the extracts. Using molecular docking, ligand poses of selected VM seeds’ phytoconstituents were predicted targeting tubulin, GlcN-6-P synthase, and human tissue plasminogen activator proteins for anthelmintic, antibacterial, and thrombolytic activity, respectively. In the antibacterial activity test, the MESVM at 10000 μg/mL concentration created highest and significant (P < 0.001) zone of inhibition against Staphylococcus aureus (15.42 mm) and Escherichia coli (12 mm) compared with tetracycline. The MESVM exhibited remarkable anthelmintic activity at 50 mg/mL concentration with 35.4 min paralysis time, 75.2 min death time and were closer to the durations of standard drug albendazole. No test extract showed anxiolytic activity. In thrombolytic activity test, all concentrations of MESVM produced clot lytic activity with high significance (P < 0.001) in comparison with the blank. In docking, 2′-hydroxygenistein, cyclokievitone hydrate, and aureol displayed maximum affinity to the target proteins for anthelmintic, antibacterial, and thrombolytic activity, respectively. This research revealed that the MESVM demonstrated potential anthelmintic, antibacterial and thrombolytic effects that confirmed the folkloric uses of VM and the found relevant lead compounds might be further optimized in future drug development.

Importance of Binding Site Hydration and Flexibility Revealed When Optimizing a Macrocyclic Inhibitor of the Keap1-Nrf2 Protein-Protein Interaction

Upregulation of the transcription factor Nrf2 by inhibition of the interaction with its negative regulator Keap1 constitutes an opportunity for the treatment of disease caused by oxidative stress. We report a structurally unique series of nanomolar Keap1 inhibitors obtained from a natural product-derived macrocyclic lead. Initial exploration of the structure-activity relationship of the lead, followed by structure-guided optimization, resulted in a 100-fold improvement in inhibitory potency. The macrocyclic core of the nanomolar inhibitors positions three pharmacophore units for productive interactions with key residues of Keap1, including R415, R483, and Y572. Ligand optimization resulted in the displacement of a coordinated water molecule from the Keap1 binding site and a significantly altered thermodynamic profile. In addition, minor reorganizations of R415 and R483 were accompanied by major differences in affinity between ligands. This study therefore indicates the importance of accounting both for the hydration and flexibility of the Keap1 binding site when designing high-affinity ligands.

Highly Similar Sequence and Structure Yet Different Biophysical Behavior: A Computational Study of Two Triosephosphate Isomerases

Homodimeric triosephosphate isomerases (TIMs) from Trypanosoma cruzi (TcTIM) and Trypanosoma brucei (TbTIM) have markedly similar amino-acid sequences and three-dimensional structures. However, several of their biophysical parameters, such as their susceptibility to sulfhydryl agents and their reactivation speed after being denatured, have significant differences. The causes of these differences were explored with microsecond-scale molecular dynamics (MD) simulations of three different TIM proteins: TcTIM, TbTIM, and a chimeric protein, Mut1. We examined their electrostatic interactions and explored the impact of simulation length on them. The same salt bridge between catalytic residues Lys 14 and Glu 98 was observed in all three proteins, but key differences were found in other interactions that the catalytic amino acids form. In particular, a cation-π interaction between catalytic amino acids Lys 14 and His 96 and both a salt bridge and a hydrogen bond between catalytic Glu 168 and residue Arg 100 were only observed in TcTIM. Furthermore, although TcTIM forms less hydrogen bonds than TbTIM and Mut1, its hydrogen bond network spans almost the entire protein, connecting the residues in both monomers. This work provides new insight into the mechanisms that give rise to the different behavior of these proteins. The results also show the importance of long simulations.

Enhanced control of bioactivity of tissue plasminogen activator (tPA) through domain-directed enzymatic oxidation of terminal galactose

Introduction: In targeted enzyme prodrug constructs, it is critical to control the bioactivity of the drug in its prodrug form. The preparation of such constructs often involves conjugation reactions directed to functional groups on amino acid side chains of the protein, which result in random conjugation and incomplete control of bioactivity of a prodrug, which may result in significant nontarget effect. Thus, more specific method of modification is desired. If the drug is a glycoprotein, enzymatic oxidation may offer an alternative approach for therapeutic glycoproteins. Methods: Tissue plasminogen activator (tPA), a model glycoprotein enzyme, was treated with galactose oxidase (GO) and horseradish peroxidase, followed by thiolation reaction and conjugation with low molecular weight heparin (LMWH). The LMWH-tPA conjugate was isolated by ion-exchange chromatography followed by centrifugal filtration. The conjugate was characterized for its fibrinolytic activity and for its plasminogen activation through an indirect amidolytic assay with a plasmin-specific substrate S-2251 when LMWH-tPA conjugate is complexed with protamine-albumin conjugate, followed by triggered activation in the presence of heparin. Results: LMWH-tPA conjugate prepared via enzymatic oxidation retained ~95% of its fibrinolytic activity with respect to native tPA. Upon complexation with protamine-albumin conjugate, the activity of LMWH-tPA was effectively inhibited (~90%) whereas the LMWH-tPA prepared by random thiolation exhibited ~55% inhibition. Addition of heparin fully generated the activities of both conjugates. Conclusion: The tPA was successfully modified via enzymatic oxidation by GO, resulting in enhanced control of its activity in the prodrug construct. This approach can be applied to other therapeutic glycoproteins.