Modeling protein quaternary structure of homo- and hetero-oligomers beyond binary interactions by homology

Recent advances and challenges in protein complex model accuracy estimation

Estimation of model accuracy plays a crucial role in protein structure prediction, aiming to evaluate the quality of predicted protein structure models accurately and objectively. This process is not only key to screening candidate models that are close to the real structure, but also provides guidance for further optimization of protein structures. With the significant advancements made by AlphaFold2 in monomer structure, the problem of single-domain protein structure prediction has been widely solved. Correspondingly, the importance of assessing the quality of single-domain protein models decreased, and the research focus has shifted to estimation of model accuracy of protein complexes. In this review, our goal is to provide a comprehensive overview of the reference and statistical metrics, as well as representative methods, and the current challenges within four distinct facets (Topology Global Score, Interface Total Score, Interface Residue-Wise Score, and Tertiary Residue-Wise Score) in the field of complex EMA.

New bacteriophage-derived lysins, LysJ and LysF, with the potential to control Bacillus anthracis

Bacillus anthracis is an etiological agent of anthrax, a severe zoonotic disease that can be transmitted to people and cause high mortalities. Bacteriophages and their lytic enzymes, endolysins, have potential therapeutic value in treating infections caused by this bacterium as alternatives or complements to antibiotic therapy. They can also be used to identify and detect B. anthracis. Endolysins of two B. anthracis Wbetavirus phages, J5a and F16Ba which were described by us recently, differ significantly from the best-known B. anthracis phage endolysin PlyG from Wbetavirus genus bacteriophage Gamma and a few other Wbetavirus genus phages. They are larger than PlyG (351 vs. 233 amino acid residues), contain a signal peptide at their N-termini, and, by prediction, have a different fold of cell binding domain suggesting different structural basis of cell epitope recognition. We purified in a soluble form the modified versions of these endolysins, designated by us LysJ and LysF, respectively, and depleted of signal peptides. Both modified endolysins could lyse the B. anthracis cell wall in zymogram assays. Their activity against the living cells of B. anthracis and other species of Bacillus genus was tested by spotting on the layers of bacteria in soft agar and by assessing the reduction of optical density of bacterial suspensions. Both methods proved the effectiveness of LysJ and LysF in killing the anthrax bacilli, although the results obtained by each method differed. Additionally, the lytic efficiency of both proteins was different, which apparently correlates with differences in their amino acid sequence. KEY POINTS: • LysJ and LysF are B. anthracis-targeting lysins differing from lysins studied so far • LysJ and LysF could be overproduced in E. coli in soluble and active forms • LysJ and LysF are active in killing cells of B. anthracis virulent strains.

Temperature effect in the inhibition of PLA2 activity of Bothrops brazili venom by Rosmarinic and Chlorogenic acids, experimental and computational approaches

Myotoxicity caused by snakebite envenoming emerges as one of the main problems of ophidic accidents as it is not well neutralized by the current serum therapy. A promising alternative is to search for efficient small molecule inhibitors that can act against multiple venom components. Phospholipase A2 (PLA2) is frequently found in snake venom and is usually associated with myotoxicity. Thus it represents an excellent target for the search of new treatments. This work reports the effect of temperature in the inhibition of catalytic properties of PLA2 from Bothrops brazili venom by Rosmarinic (RSM) and Chlorogenic (CHL) acids through experimental and computational approaches. Three temperatures were evaluated (25, 37 and 50 °C). In the experimental section, enzymatic assays showed that RSM is a better inhibitor in all three temperatures. At 50 °C, the inhibition efficiency decayed significantly for both acids. Docking studies revealed that both ligands bind to the hydrophobic channel of the protein dimer where the phospholipid binds in the catalytic process, interacting with several functional residues. In this context, RSM presents better interaction energies due to stronger interactions with chain B of the dimer. Molecular dynamics simulations showed that RSM can establish selective interactions with ARG112B of PLA2, which is located next to residues of the putative Membrane Disruption Site in PLA2-like structures. The affinity of RSM and CHL acids towards PLA2 is mainly driven by electrostatic interactions, especially salt bridge interactions established with residues ARG33B (for CHL) and ARG112B (RSM) and hydrogen bonds with residue ASP89A. The inability of CHL to establish a stable interaction with ARG112B was identified as the reason for its lower inhibition efficiency compared to RSM at the three temperatures. Furthermore, extensive structural analysis was performed to explain the lower inhibition efficiency at 50 °C for both ligands. The analysis performed in this work provides important information for the future design of new inhibitors.Communicated by Ramaswamy H. Sarma.

EGG: Accuracy Estimation of Individual Multimeric Protein Models Using Deep Energy-Based Models and Graph Neural Networks

Reliable and accurate methods of estimating the accuracy of predicted protein models are vital to understanding their respective utility. Discerning how the quaternary structure conforms can significantly improve our collective understanding of cell biology, systems biology, disease formation, and disease treatment. Accurately determining the quality of multimeric protein models is still computationally challenging, as the space of possible conformations is significantly larger when proteins form in complex with one another. Here, we present EGG (energy and graph-based architectures) to assess the accuracy of predicted multimeric protein models. We implemented message-passing and transformer layers to infer the overall fold and interface accuracy scores of predicted multimeric protein models. When evaluated with CASP15 targets, our methods achieved promising results against single model predictors: fourth and third place for determining the highest-quality model when estimating overall fold accuracy and overall interface accuracy, respectively, and first place for determining the top three highest quality models when estimating both overall fold accuracy and overall interface accuracy.

Characterization of Sf9 cell clones with differential susceptibilities to Sf-rhabdovirus X+3.7 and Sf-rhabdovirus X- replication

Our laboratory previously discovered a novel rhabdovirus in the Spodoptera frugiperda Sf9 insect cell line that was designated as Sf-rhabdovirus. Using limiting dilution, this cell line was found to be a mixed population of cells infected by Sf-rhabdovirus variants containing either the full length X accessory gene with a 3.7 kb internal duplication (designated as Sf-rhabdovirus X+3.7) or lacking the duplication and part of the X gene (designated as Sf-rhabdovirus X-), and cells that were negative for Sf-rhabdovirus. In this paper, we found that the Sf-rhabdovirus negative cell clones had sub-populations with different susceptibilities to the replication of Sf-rhabdovirus X+3.7 and X- variants: cell clone Sf9-13F12 was more sensitive to replication by both virus variants compared to Sf9-3003; moreover, Sf9-3003 showed more resistance to X+3.7 replication than to X- replication. RNA-Seq analysis indicated significant differentially expressed genes in the Sf9-13F12 and Sf9-3003 cell clones further supporting that distinct sub-populations of virus-negative cells co-exist in the parent Sf9 cell line.

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.

A Survey of Deep Learning Methods for Estimating the Accuracy of Protein Quaternary Structure Models

The quality prediction of quaternary structure models of a protein complex, in the absence of its true structure, is known as the Estimation of Model Accuracy (EMA). EMA is useful for ranking predicted protein complex structures and using them appropriately in biomedical research, such as protein-protein interaction studies, protein design, and drug discovery. With the advent of more accurate protein complex (multimer) prediction tools, such as AlphaFold2-Multimer and ESMFold, the estimation of the accuracy of protein complex structures has attracted increasing attention. Many deep learning methods have been developed to tackle this problem; however, there is a noticeable absence of a comprehensive overview of these methods to facilitate future development. Addressing this gap, we present a review of deep learning EMA methods for protein complex structures developed in the past several years, analyzing their methodologies, data and feature construction. We also provide a prospective summary of some potential new developments for further improving the accuracy of the EMA methods.

Different receptor models show differences in ligand binding strength and location: a computational drug screening for the tick-borne encephalitis virus

The tick-borne encephalitis virus (TBE) is a neurotrophic disease that has spread more rapidly throughout Europe and Asia in the past few years. At the same time, no cure or specific therapy is known to battle the illness apart from vaccination. To find a pharmacologically relevant drug, a computer-aided drug screening was initiated. Such a procedure probes a possible binding of a drug to the RNA Polymerase of TBE. The crystal structure of the receptor, however, includes missing and partially modeled regions, which rendered the structure incomplete and of questionable use for a thorough drug screening procedure. The quality of the receptor model was addressed by studying three putative structures created. We show that the choice of receptor models greatly influences the binding affinity of potential drug molecules and that the binding location could also be significantly impacted. We demonstrate that some drug candidates are unsuitable for one model but show decent results for another. Without any prejudice on the three employed receptor models, the study reveals the imperative need to investigate the receptor structure before drug binding is probed whether experimentally or computationally.

Copper Resistance Mechanism and Copper Response Genes in Corynebacterium crenatum

Heavy metal resistance mechanisms and heavy metal response genes are crucial for microbial utilization in heavy metal remediation. Here, Corynebacterium crenatum was proven to possess good tolerance in resistance to copper. Then, the transcriptomic responses to copper stress were investigated, and the vital pathways and genes involved in copper resistance of C. crenatum were determined. Based on transcriptome analysis results, a total of nine significantly upregulated DEGs related to metal ion transport were selected for further study. Among them, GY20_RS0100790 and GY20_RS0110535 belong to transcription factors, and GY20_RS0110270, GY20_RS0100790, and GY20_RS0110545 belong to copper-binding peptides. The two transcription factors were studied for the function of regulatory gene expression. The three copper-binding peptides were displayed on the C. crenatum surface for a copper adsorption test. Furthermore, the nine related metal ion transport genes were deleted to investigate the effect on growth in copper stress. This investigation provided the basis for utilizing C. crenatum in copper bioremediation.

Novel benzimidazole angiotensin receptor blockers with anti-SARS-CoV-2 activity equipotent to that of nirmatrelvir: computational and enzymatic studies

BACKGROUND

Hypertension worsens outcomes in SARS-CoV-2 patients. Sartans, a type of antihypertensive angiotensin receptor blocker-(ARB), reduce COVID-19 morbidity and mortality by targeting angiotensin-converting enzyme-2 (ACE2). This study aimed to evaluate the antiviral and antihypertensive effects of nirmatrelvir, commercial sartans (candesartan, losartan, and losartan carboxylic (Exp3174)), and newly synthesized sartans (benzimidazole-N-biphenyl carboxyl (ACC519C) and benzimidazole-N-biphenyl tetrazole (ACC519T)), compared to nirmatrelvir, the antiviral component of Paxlovid.

RESEARCH DESIGN AND METHODS

Surface plasmon resonance (SPR) and enzymatic studies assessed drug effects on ACE2. Antiviral abilities were tested with SARS-CoV-2-infected Vero E6 cells, and antihypertensive effects were evaluated using angiotensin II-contracted rabbit iliac arteries.

RESULTS

Benzimidazole-based candesartan and ACC519C showed antiviral activity comparable to nirmatrelvir (95% inhibition). Imidazole-based losartan, Exp3174, and ACC519T were less potent (75%-80% and 50%, respectively), with Exp3174 being the least effective. SPR analysis indicated high sartans-ACE2 binding affinity. Candesartan and nirmatrelvir combined had greater inhibitory and cytopathic effects (3.96%) than individually (6.10% and 5.08%). ACE2 enzymatic assays showed varying effects of novel sartans on ACE2. ACC519T significantly reduced angiotensin II-mediated contraction, unlike nirmatrelvir and ACC519T(2).

CONCLUSION

This study reports the discovery of a new class of benzimidazole-based sartans that significantly inhibit SARS-CoV-2, likely due to their interaction with ACE2.

Secretory Production of the Hericium erinaceus Laccase from Saccharomyces cerevisiae

Mushroom laccases play a crucial role in lignin depolymerization, one of the most critical challenges in lignin utilization. Importantly, laccases can utilize a wide range of substrates, such as toxicants and antibiotics. This study isolated a novel laccase, named HeLac4c, from endophytic white-rot fungi Hericium erinaceus mushrooms. The cDNAs for this enzyme were 1569 bp in length and encoded a protein of 523 amino acids, including a 20 amino-acid signal peptide. Active extracellular production of glycosylated laccases from Saccharomyces cerevisiae was successfully achieved by selecting an optimal translational fusion partner. We observed that 5 and 10 mM Ca2+, Zn2+, and K+ increased laccase activity, whereas 5 mM Fe2+ and Al3+ inhibited laccase activity. The laccase activity was inhibited by the addition of low concentrations of sodium azide and L-cysteine. The optimal pH for the 2,2'-Azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) diammonium salt was 4.4. Guaiacylglycerol-β-guaiacyl ether, a lignin model compound, was polymerized by the HeLac4c enzyme. These results indicated that HeLac4c is a novel oxidase biocatalyst for the bioconversion of lignin into value-added products for environmental biotechnological applications.

Establishing comprehensive quaternary structural proteomes from genome sequence

A critical body of knowledge has developed through advances in protein microscopy, protein-fold modeling, structural biology software, availability of sequenced bacterial genomes, large-scale mutation databases, and genome-scale models. Based on these recent advances, we develop a computational framework that; i) identifies the oligomeric structural proteome encoded by an organism's genome from available structural resources; ii) maps multi-strain alleleomic variation, resulting in the structural proteome for a species; and iii) calculates the 3D orientation of proteins across subcellular compartments with residue-level precision. Using the platform, we; iv) compute the quaternary E. coli K-12 MG1655 structural proteome; v) use a dataset of 12,000 mutations to build Random Forest classifiers that can predict the severity of mutations; and, in combination with a genome-scale model that computes proteome allocation, vi) obtain the spatial allocation of the E. coli proteome. Thus, in conjunction with relevant datasets and increasingly accurate computational models, we can now annotate quaternary structural proteomes, at genome-scale, to obtain a molecular-level understanding of whole-cell functions.

Significance

Advancements in experimental and computational methods have revealed the shapes of multi-subunit proteins. The absence of a unified platform that maps actionable datatypes onto these increasingly accurate structures creates a barrier to structural analyses, especially at the genome-scale. Here, we describe QSPACE, a computational annotation platform that evaluates existing resources to identify the best-available structure for each protein in a user's query, maps the 3D location of actionable datatypes ( e.g. , active sites, published mutations) onto the selected structures, and uses third-party APIs to determine the subcellular compartment of all amino acids of a protein. As proof-of-concept, we deployed QSPACE to generate the quaternary structural proteome of E. coli MG1655 and demonstrate two use-cases involving large-scale mutant analysis and genome-scale modelling.

A Structural In Silico Analysis of the Immunogenicity of L-Asparaginase from Penicillium cerradense

L-asparaginase is an essential drug used to treat acute lymphoid leukemia (ALL), a cancer of high prevalence in children. Several adverse reactions associated with L-asparaginase have been observed, mainly caused by immunogenicity and allergenicity. Some strategies have been adopted, such as searching for new microorganisms that produce the enzyme and applying protein engineering. Therefore, this work aimed to elucidate the molecular structure and predict the immunogenic profile of L-asparaginase from Penicillium cerradense, recently revealed as a new fungus of the genus Penicillium and producer of the enzyme, as a motivation to search for alternatives to bacterial L-asparaginase. In the evolutionary relationship, L-asparaginase from P. cerradense closely matches Aspergillus species. Using in silico tools, we characterized the enzyme as a protein fragment of 378 amino acids (39 kDa), including a signal peptide containing 17 amino acids, and the isoelectric point at 5.13. The oligomeric state was predicted to be a homotetramer. Also, this L-asparaginase presented a similar immunogenicity response (T- and B-cell epitopes) compared to Escherichia coli and Dickeya chrysanthemi enzymes. These results suggest a potentially useful L-asparaginase, with insights that can drive strategies to improve enzyme production.

Origin of Enantioselectivity in Engineered Cytochrome c-Catalyzed Carbon-Radical FePP Hydrolysis Revealed Using QM/MM (ABEEM Polarizable Force Field) and MD Simulations

The origin of highly efficient asymmetric aminohydroxylation of styrene catalyzed by engineered cytochrome c is investigated by the developed Atom-Bond Electronegativity Equalization Method polarizable force field (ABEEM PFF), which is a combined outcome of electronic and steric effects. Model molecules were used to establish the charge parameters of the ABEEM PFF, for which the bond-stretching and angle-bending parameters were obtained by using a combination of modified Seminario and scan methods. The interactions between carbon-radical Fe-porphyrin (FePP) and waters are simulated by molecular dynamics, which shows a clear preference for the pre-R over the pre-S. This preference is attributed to the hydrogen-bond between the mutated 100S and 101P residues as well as van der Waals interactions, enforcing a specific conformation of the carbon-radical FePP complex within the binding pocket. Meanwhile, the hydrogen-bond between water and the nitrogen atom in the active intermediate dictates the stereochemical outcome. Quantum mechanics/molecular mechanics (QM/MM (ABEEM PFF)) and free-energy perturbation calculations elucidate that the 3RTS is characterized by sandwich-like structure among adjacent amino acid residues, which exhibits greater stability than crowed arrangement in 3STS and enables the R enantiomer to form more favorably. Thus, this study provides mechanistic insight into the catalytic reaction of hemoproteins.

Generalized biomolecular modeling and design with RoseTTAFold All-Atom

Deep-learning methods have revolutionized protein structure prediction and design but are presently limited to protein-only systems. We describe RoseTTAFold All-Atom (RFAA), which combines a residue-based representation of amino acids and DNA bases with an atomic representation of all other groups to model assemblies that contain proteins, nucleic acids, small molecules, metals, and covalent modifications, given their sequences and chemical structures. By fine-tuning on denoising tasks, we developed RFdiffusion All-Atom (RFdiffusionAA), which builds protein structures around small molecules. Starting from random distributions of amino acid residues surrounding target small molecules, we designed and experimentally validated, through crystallography and binding measurements, proteins that bind the cardiac disease therapeutic digoxigenin, the enzymatic cofactor heme, and the light-harvesting molecule bilin.

The structure assessment web server: for proteins, complexes and more

The 'structure assessment' web server is a one-stop shop for interactive evaluation and benchmarking of structural models of macromolecular complexes including proteins and nucleic acids. A user-friendly web dashboard links sequence with structure information and results from a variety of state-of-the-art tools, which facilitates the visual exploration and evaluation of structure models. The dashboard integrates stereochemistry information, secondary structure information, global and local model quality assessment of the tertiary structure of comparative protein models, as well as prediction of membrane location. In addition, a benchmarking mode is available where a model can be compared to a reference structure, providing easy access to scores that have been used in recent CASP experiments and CAMEO. The structure assessment web server is available at https://swissmodel.expasy.org/assess.

Engineering an artificial catch bond using mechanical anisotropy

Catch bonds are a rare class of protein-protein interactions where the bond lifetime increases under an external pulling force. Here, we report how modification of anchor geometry generates catch bonding behavior for the mechanostable Dockerin G:Cohesin E (DocG:CohE) adhesion complex found on human gut bacteria. Using AFM single-molecule force spectroscopy in combination with bioorthogonal click chemistry, we mechanically dissociate the complex using five precisely controlled anchor geometries. When tension is applied between residue #13 on CohE and the N-terminus of DocG, the complex behaves as a two-state catch bond, while in all other tested pulling geometries, including the native configuration, it behaves as a slip bond. We use a kinetic Monte Carlo model with experimentally derived parameters to simulate rupture force and lifetime distributions, achieving strong agreement with experiments. Single-molecule FRET measurements further demonstrate that the complex does not exhibit dual binding mode behavior at equilibrium but unbinds along multiple pathways under force. Together, these results show how mechanical anisotropy and anchor point selection can be used to engineer artificial catch bonds.

Functional characterization and structural insights of three cutinases from the ascomycete Fusarium verticillioides

Cutinases are serine esterases that belong to the α/β hydrolases superfamily. The natural substrates for these enzymes are cutin and suberin, components of the plant cuticle, the first barrier in the defense system against pathogen invasion. It is well-reported that plant pathogens produce cutinases to facilitate infection. Fusarium verticillioides, one important corn pathogens, is an ascomycete upon which its cutinases are poorly explored. Consequently, the objective of this study was to perform the biochemical characterization of three precursor cutinases (FvCut1, FvCut2, and FvCut3) from F. verticillioides and to obtain structural insights about them. The cutinases were produced in Escherichia coli and purified. FvCut1, FvCut2, and FvCut3 presented optimal temperatures of 20, 40, and 35 °C, and optimal pH of 9, 7, and 8, respectively. Some chemicals stimulated the enzymatic activity. The kinetic parameters revealed that FvCut1 has higher catalytic efficiency (Kcat/Km) in the p-nitrophenyl-butyrate (p-NPB) substrate. Nevertheless, the enzymes were not able to hydrolyze polyethylene terephthalate (PET). Furthermore, the three-dimensional models of these enzymes showed structural differences among them, mainly FvCut1, which presented a narrower opening cleft to access the catalytic site. Therefore, our study contributes to exploring the diversity of fungal cutinases and their potential biotechnological applications.

2,4-disubstituted 6-fluoroquinolines as potent antiplasmodial agents: QSAR, homology modeling, molecular docking and ADMET studies

Objective

This work was designed to study 2,4-disubstituted 6-fluoroquinolines as antiplasmodial agents by using in silico techniques, to aid in the design of novel analogs with high potency against malaria and high inhibition of Plasmodium falciparum translation elongation factor 2 (PfeEF2), a novel drug target.

Methods

Quantitative structure-activity relationships (QSAR) of 2,4-disubstituted 6-fluoroquinolines were studied with the genetic function approximation technique in Material Studio software. The 3D structure of PfeEF2 was modeled in the SWISS-MODEL workspace through homology modeling. A molecular docking study of the modeled PfeEF2 and 2,4-disubstituted 6-fluoroquinolines was conducted with Autodock Vina in Pyrx software. Furthermore, the in silico pharmacokinetic properties of selected compounds were investigated.

Results

A robust, reliable and predictive QSAR model was developed that related the chemical structures of 2,4-disubstituted 6-fluoroquinolines to their antiplasmodium activities. The model had an internal squared correlation coefficient R2 of 0.921, adjusted squared correlation coefficient R2adj of 0.878, leave-one-out cross-validation coefficient Q2cv of 0.801 and predictive squared correlation coefficient R2pred of 0.901. The antiplasmodium activity of 6-fluoroquinolines was found to depend on the n5Ring, GGI9, TDB7u, TDB8u and RDF75i physicochemical properties: n5Ring, TDB8u and RDF75i were positively associated, whereas GGI9 and TDB7u were negatively associated, with the antiplasmodium activity of the compounds. Stable complexes formed between the compounds and modeled PfeEF2, with binding affinity ranging from -8.200 to -10.700 kcal/mol. Compounds 5, 11, 16, 22 and 24 had better binding affinities than quinoline-4-carboxamide (DDD107498), as well as good pharmacokinetic properties, and therefore may be better inhibitors of this novel target.

Conclusion

QSAR and docking studies provided insight into designing novel 2,4-disubstituted 6-fluoroquinolines with high antiplasmodial activity and good structural properties for inhibiting a novel antimalarial drug target.

In vitro and in silico studies of enterobactin-inspired Ciprofloxacin and Fosfomycin first generation conjugates on the antibiotic resistant E. coli OQ866153

BACKGROUND

The emergence of antimicrobial resistance in bacterial pathogens is a growing concern worldwide due to its impact on the treatment of bacterial infections. The "Trojan Horse" strategy has been proposed as a potential solution to overcome drug resistance caused by permeability issues.

OBJECTIVE

The objective of our research was to investigate the bactericidal activity and mechanism of action of the "Trojan Horse" strategy using enterobactin conjugated with Ciprofloxacin and Fosfomycin against the antibiotic-resistant Escherichia coli strain OQ866153.

METHODOLOGY

Enterobactin, a mixed ligand of E. coli OQ866153, was conjugated with Ciprofloxacin and Fosfomycin individually to aid active absorption via specific enterobactin binding proteins (FepABCDG). The effectiveness of the conjugates was assessed by measuring their bactericidal activity against E. coli OQ866153, as well as their ability to inhibit DNA gyrase enzyme and biofilm formation.

RESULTS

The Fe+3-enterobactin-Ciprofloxacin conjugate effectively inhibited the DNA gyrase enzyme (Docking score = -8.597 kcal/mol) and resulted in a lower concentration (25 μg/ml) required to eliminate supercoiled DNA plasmids compared to the parent drug (35 μg/ml; Docking score = -6.264 kcal/mol). The Fe+3-Enterobactin-Fosfomycin conjugate showed a higher inhibition percentage (100%) of biofilm formation compared to Fosfomycin (21.58%) at a concentration of 2 mg/ml, with docking scores of -5.481 and -3.756 kcal/mol against UDP-N acetylglucosamine 1-carboxyvinyltransferase MurA.

CONCLUSION

The findings of this study suggest that the "Trojan Horse" strategy using enterobactin conjugated with Ciprofloxacin and Fosfomycin can effectively overcome permeability issues caused by efflux proteins and enhance the bactericidal activity of these drugs against antibiotic-resistant strains of E. coli.

An artificial metalloenzyme that can oxidize water photocatalytically: design, synthesis, and characterization

In nature, light-driven water oxidation (WO) catalysis is performed by photosystem II via the delicate interplay of different cofactors positioned in its protein scaffold. Artificial systems for homogeneous photocatalytic WO are based on small molecules that often have limited solubility in aqueous solutions. In this work, we alleviated this issue and present a cobalt-based WO-catalyst containing artificial metalloenzyme (ArM) that is active in light-driven, homogeneous WO catalysis in neutral-pH aqueous solutions. A haem-containing electron transfer protein, cytochrome B5 (CB5), served to host a first-row transition-metal-based WO catalyst, CoSalen (CoIISalen, where H2Salen = N,N'-bis(salicylidene)ethylenediamine), thus producing an ArM capable of driving photocatalytic WO. The CoSalen ArM formed a water-soluble pre-catalyst in the presence of [Ru(bpy)3](ClO4)2 as photosensitizer and Na2S2O8 as the sacrificial electron acceptor, with photocatalytic activity similar to that of free CoSalen. During photocatalysis, the CoSalen-protein interactions were destabilized, and the protein partially unfolded. Rather than forming tens of nanometer sized CoOx nanoparticles as free CoSalen does under photocatalytic WO conditions, the CB5 : CoSalen ArM showed limited protein cross-linking and remained soluble. We conclude that a weak, dynamic interaction between a soluble cobalt species and apoCB5 was formed, which generated a catalytically active adduct during photocatalysis. A detailed analysis was performed on protein stability and decomposition processes during the harsh oxidizing reaction conditions of WO, which will serve for the future design of WO ArMs with improved activity and stability.

Exploring the functional profiles of odorant binding proteins crucial for sensing key odorants in the new leaves of coconut palms in Rhynchophorus ferrugineus

The red palm weevil (RPW), Rhynchophorus ferrugineus (Curculionidae: Coleoptera) is a highly destructive global pest of coconut trees, with a preference for laying its eggs on new leaves. Females can identify where to lay eggs by using their sense of smell to detect specific odorants found in new leaves. In this study, we focused on the two odorants commonly found in new leaves by GC-MS: trans, trans-2,4-nonadienal and trans-2-nonenal. Our behavioral assays demonstrated a significant attraction of females to both of these odorants, with their electrophysiological responses being dose-dependent. Furthermore, we examined the expression patterns induced by these odorants in eleven RferOBP genes. Among them, RferOBP3 and RferOBP1768 exhibited the most significant and simultaneous upregulation. To further understand the role of these two genes, we conducted experiments with females injected with OBP-dsRNA. This resulted in a significant decrease in the expression of RferOBP3 and RferOBP1768, as well as impaired the perception of the two odorants. A fluorescence competitive binding assay also showed that both RferOBPs strongly bound to the odorants. Additionally, sequence analysis revealed that these two RferOBPs belong to the Minus-C family and possess four conserved cysteines. Molecular docking simulations showed strong interactions between these two RferOBPs and the odorant molecules. Overall, our findings highlight the crucial role of RferOBP3 and RferOBP1768 in the olfactory perception of the key odorants in coconut palm new leaves. This knowledge significantly improves our understanding of how RPW females locate sites for oviposition and lays the foundation for future research on the development of environmentally friendly pest attractants.

Identification and mechanism determination of the efflux pump subunit amrB gene mutations linked to gentamicin susceptibility in clinical Burkholderia pseudomallei from Malaysian Borneo

The bacterium Burkholderia pseudomallei is typically resistant to gentamicin but rare susceptible strains have been isolated in certain regions, such as Thailand and Sarawak, Malaysia. Recently, several amino acid substitutions have been reported in the amrB gene (a subunit of the amrAB-oprA efflux pump gene) that confer gentamicin susceptibility. However, information regarding the mechanism of the substitutions conferring the susceptibility is lacking. To understand the mechanism of amino acid substitution that confers susceptibility, this study identifies the corresponding mutations in clinical gentamicin-susceptible B. pseudomallei isolates from the Malaysian Borneo (n = 46; Sarawak: 5; Sabah: 41). Three phenotypically confirmed gentamicin-susceptible (GENs) strains from Sarawak, Malaysia, were screened for mutations in the amrB gene using gene sequences of gentamicin-resistant (GENr) strains (QEH 56, QEH 57, QEH20, and QEH26) and publicly available sequences (AF072887.1 and BX571965.1) as the comparator. The effect of missense mutations on the stability of the AmrB protein was determined by calculating the average energy change value (ΔΔG). Mutagenesis analysis identified a polymorphism-associated mutation, g.1056 T > G, a possible susceptible-associated in-frame deletion, Delta V412, and a previously confirmed susceptible-associated amino acid substitution, T368R, in each of the three GENs isolates. The contribution of Delta V412 needs further confirmation by experimental mutagenesis analysis. The mechanism by which T368R confers susceptibility, as elucidated by in silico mutagenesis analysis using AmrB-modeled protein structures, is proposed to be due to the location of T368R in a highly conserved region, rather than destabilization of the AmrB protein structure.

The Chemokine CXCL14 as a Potential Immunotherapeutic Agent for Cancer Therapy

There is great enthusiasm toward the development of novel immunotherapies for the treatment of cancer, and given their roles in immune system regulation, chemokines stand out as promising candidates for use in new cancer therapies. Many previous studies have shown how chemokine signaling pathways could be targeted to halt cancer progression. We and others have revealed that the chemokine CXCL14 promotes antitumor immune responses, suggesting that CXCL14 may be effective for cancer immunotherapy. However, it is still unknown what mechanism governs CXCL14-mediated antitumor activity, how to deliver CXCL14, what dose to apply, and what combinations with existing therapy may boost antitumor immune responses in cancer patients. Here, we provide updates on the role of CXCL14 in cancer progression and discuss the potential development and application of CXCL14 as an immunotherapeutic agent.

3,4-Dichlorophenylacetic acid acts as an auxin analog and induces beneficial effects in various crops

Auxins and their analogs are widely used to promote root growth, flower and fruit development, and yield in crops. The action characteristics and application scope of various auxins are different. To overcome the limitations of existing auxins, expand the scope of applications, and reduce side effects, it is necessary to screen new auxin analogs. Here, we identified 3,4-dichlorophenylacetic acid (Dcaa) as having auxin-like activity and acting through the auxin signaling pathway in plants. At the physiological level, Dcaa promotes the elongation of oat coleoptile segments, the generation of adventitious roots, and the growth of crop roots. At the molecular level, Dcaa induces the expression of auxin-responsive genes and acts through auxin receptors. Molecular docking results showed that Dcaa can bind to auxin receptors, among which TIR1 has the highest binding activity. Application of Dcaa at the root tip of the DR5:GUS auxin-responsive reporter induces GUS expression in the root hair zone, which requires the PIN2 auxin efflux carrier. Dcaa also inhibits the endocytosis of PIN proteins like other auxins. These results provide a basis for the application of Dcaa in agricultural practices.

Antiviral fibrils of self-assembled peptides with tunable compositions

The lasting threat of viral pandemics necessitates the development of tailorable first-response antivirals with specific but adaptive architectures for treatment of novel viral infections. Here, such an antiviral platform has been developed based on a mixture of hetero-peptides self-assembled into functionalized β-sheets capable of specific multivalent binding to viral protein complexes. One domain of each hetero-peptide is designed to specifically bind to certain viral proteins, while another domain self-assembles into fibrils with epitope binding characteristics determined by the types of peptides and their molar fractions. The self-assembled fibrils maintain enhanced binding to viral protein complexes and retain high resilience to viral mutations. This method is experimentally and computationally tested using short peptides that specifically bind to Spike proteins of SARS-CoV-2. This platform is efficacious, inexpensive, and stable with excellent tolerability.

Human Hsp70 Substrate-Binding Domains Recognize Distinct Client Proteins

The 13 Hsp70 proteins in humans act on unique sets of substrates with diversity often being attributed to J-domain-containing protein (Hsp40 or JDP) cofactors. We were therefore surprised to find drastically different binding affinities for Hsp70-peptide substrates, leading us to probe substrate specificity among the 8 canonical Hsp70s from humans. We used peptide arrays to characterize Hsp70 binding and then mined these data using machine learning to develop an algorithm for isoform-specific prediction of Hsp70 binding sequences. The results of this algorithm revealed recognition patterns not predicted based on local sequence alignments. We then showed that none of the human isoforms can complement heat-shocked DnaK knockout Escherichia coli cells. However, chimeric Hsp70s consisting of the human nucleotide-binding domain and the substrate-binding domain of DnaK complement during heat shock, providing further evidence in vivo of the divergent function of the Hsp70 substrate-binding domains. We also demonstrated that the differences in heat shock complementation among the chimeras are not due to loss of DnaJ binding. Although we do not exclude JDPs as additional specificity factors, our data demonstrate substrate specificity among the Hsp70s, which has important implications for inhibitor development in cancer and neurodegeneration.

Optimizing ErCas12a for efficient gene editing in Arabidopsis thaliana

The ErCas12a nuclease, also known as MAD7, is part of a CRISPR/Cas system from Eubacterium rectale and distantly related to Cas12a nucleases. As it shares only 31% sequence homology with the commonly used AsCas12a, its intellectual property may not be covered by the granted patent rights for Cas12a nucleases. Thus, ErCas12a became an attractive alternative for practical applications. However, the editing efficiency of ErCas12a is strongly target sequence- and temperature-dependent. Therefore, optimization of the enzyme activity through protein engineering is especially attractive for its application in plants, as they are cultivated at lower temperatures. Based on the knowledge obtained from the optimization of Cas12a nucleases, we opted to improve the gene editing efficiency of ErCas12a by introducing analogous amino acid exchanges. Interestingly, neither of these mutations analogous to those in the enhanced or Ultra versions of AsCas12a resulted in significant editing enhancement of ErCas12a in Arabidopsis thaliana. However, two different mutations, V156R and K172R, in putative alpha helical structures of the enzyme showed a detectable improvement in editing. By combining these two mutations, we obtained an improved ErCas12a (imErCas12a) variant, showing several-fold increase in activity in comparison to the wild-type enzyme in Arabidopsis. This variant yields strong editing efficiencies at 22 °C which could be further increased by raising the cultivation temperature to 28 °C and even enabled editing of formerly inaccessible targets. Additionally, no enhanced off-site activity was detected. Thus, imErCas12a is an economically attractive and efficient alternative to other CRISPR/Cas systems for plant genome engineering.

Bioactive metabolites of licorice and thyme as potential inhibitors of Cox1 enzyme of phytopathogens of Capsicum annuum L.: In-silico approaches

Cytochrome c oxidase subunit 1 (Cox1), a key enzyme, has a crucial role in cellular respiration in eukaryotes and prokaryotes. Generally, respiratory inhibitors are considered one of the types of chemical pesticides. Thyme oil and licorice aqueous extract have been reported to have antifungal activities against fungal phytopathogens of Capsicum annuum L., i.e., Colletotrichum capsici, Fusarium oxysporum, and Pythium aphanidermatum. The present study focuses on identifying the key bioactive molecules of thyme and licorice botanicals inhibiting the activity of the Cox1 enzymes of the above mentioned phytopathogens, employing the in-silico approach. From a wide range of bioactive molecules screened, the molecular docking indicated trans-carveol, carvacrol, kaempferol 3-rhamnoside 7-xyloside, kaempferitrin, and astragalin 7-rhamnoside as the potential inhibitors for Cox1 of C. capsici, β-Caryophyllene, Caryophyllene acetate, hispaglabridin A, kaempferol 3-rhamnoside 7-xyloside and licorice glycoside A for Cox1 of F. oxysporum and (+)-Longifolen, Caryophyllene acetate, Hispaglabridin A, Neoliquiritin 2''-apioside and Licorice-saponin A3 for Cox1 of P. aphanidermatum. Most of the top-scoring bioactive molecules exhibited higher binding affinity with the targets than the chemical compound, i.e., carbendazim. Density functional theory (DFT) analysis confirmed the reactivity of the top-docked compounds. Molecular dynamic simulations confirmed the stability of docked complexes when evaluated through multiple descriptors. Additionally, MM/PBSA analysis supported the findings, indicating the spontaneous binding of the enzymes to the screened ligands. ADMET analysis revealed the safety of the selected bioactive compounds. The present findings could be useful in developing biopesticidal formulations as efficient and sustainable alternatives to chemical pesticides.Communicated by Ramaswamy H. Sarma.

Synergistic and antagonistic interactions of oxybenzone and ocean acidification: new insight into vulnerable cellular processes in non-calcifying anthozoans

Cnidarians face significant threats from ocean acidification (OA) and anthropogenic pollutants such as oxybenzone (BP-3). The convergence of threats from multiple stressors is an important area to investigate because of potential significant synergistic or antagonistic interactions. Real-time quantitative PCR was performed to characterize the expression profiles of twenty-two genes of interest (GOI) in sea anemones (Exaiptasia diaphana) exposed to one of four treatments: 1) 96 h of OA conditions followed by a 4 h exposure to 20 ppb BP-3; 2) Exposure to 4 h 20 ppb BP-3 without 96 h of OA; 3) Exposure to 96 h of OA alone; or 4) laboratory conditions with no exposure to BP-3 and/or OA. These 22 GOIs represent cellular processes associated with proton-dependent transport, sodium-dependent transport, metal cation binding/transport, extracellular matrix, amino acid metabolism/transport, immunity, and/or steroidogenesis. These 22 GOIs provide new insight into vulnerable cellular processes in non-calcifying anthozoans exposed to OA and BP-3. Expression profiles were categorized as synergistic, antagonistic, or additive of BP-3 in the presence of OA. Two GOIs were synergistic. Fifteen GOIs were antagonistic and the remaining five GOIs were additive in response to BP-3 in acidified seawater. A subset of these GOIs appear to be candidate biomarkers for future in situ investigations. In human health, proton-dependent monocarboxylate transporters (MCTs) are promising pharmacological targets and recognized as potential biomarkers. By comparison, these same MCTs appear to be targets of xenobiotic chemical pollutants in cnidarian physiology. In the presence of BP-3, a network of collagen synthesis genes are upregulated and antagonistic in their expression profiles. Cytochrome b561 is a critical protein required for collagen synthesis and in silico modeling demonstrates BP-3 binds in the pocket of cytochrome b561. Understanding the underlying molecular mechanisms of "drug-like" compounds such as BP-3 may lead to a more comprehensive interpretation of transcriptional expression profiles. The collective antagonistic responses of GOIs associated with collagen synthesis strongly suggests these GOIs should be considered candidate biomarkers of effect. GOIs with synergistic and additive responses represent candidate biomarkers of exposure. Results show the effects of OA and BP-3 are interactive with respect to their impact on cnidarians. This investigation offers mechanistic data that supports the expression profiles and underpins higher order physiological responses.

Matrix Metalloproteinase-Responsive Delivery of PEGylated Fibroblast Growth Factor 2

Attachment of polyethylene glycol (PEG) chains is a common, well-studied, and Food and Drug Administration-approved method to address the pharmacokinetic challenges of therapeutic proteins. Occasionally, PEGylation impairs the activity of pharmacodynamics (PD). To overcome this problem, disease-relevant cleavable linkers between the polymer and the therapeutic protein can unleash full PD by de-PEGylating the protein at its target site. In this study, we engineered a matrix metalloproteinase (MMP)-responsive fibroblast growth factor 2 (FGF-2) mutant that was site-specifically extended with a PEG polymer chain. Using bioinspired strategies, the bioconjugate was designed to release the native protein at the desired structure/environment with preservation of the proliferative capacity in vitro on NIH3T3 cells. In vivo, hepatic exposure was diminished but not its renal distribution over time compared to unconjugated FGF-2. By releasing the growth factor from the PEG polymer in response to MMP cleavage, restored FGF-2 may enter hard-to-reach tissues and activate cell surface receptors or nuclear targets.

Whole-exome sequencing reveals a likely pathogenic LMNA variant causing hypertrophic cardiomyopathy

OBJECTIVE

We studied the clinical and molecular features of a family with hypertrophic cardiomyopathy (HCM).

BACKGROUND

A very heterogeneous disease affecting the heart muscle, HCM is mostly caused by variants in the proteins of sarcomeres. The detection of HCM pathogenic variants can affect the handling of patients and their families.

METHODS

Whole-exome sequencing (WES) was performed to assess the genetic cause(s) of HCM in a consanguineous Iranian family.

RESULTS

Missense likely pathogenic variant c.1279C>T (p.Arg427Cys) within exon 7 of the LMNA gene (NM_170707) was found. The segregations were confirmed by polymerase chain reaction-based Sanger sequencing.

CONCLUSIONS

Variant c.1279C>T (p.Arg427Cys) in the LMNA gene seemed to have been the cause of HCM in the family. A few LMNA gene variants related to HCM phenotypes have been recognized so far. Identifying HCM genetic basis confers significant opportunities to understand how the disease can develop and, by extension, how this progression can be arrested. Our study supports WES effectiveness for first-tier variant screening of HCM in a clinical setting.

Phylogenetic conservation of Trop-2 across species-rodent and primate genomics model anti-Trop-2 therapy for pre-clinical benchmarks

A phylogenetic conservation analysis of Trop-2 across vertebrate species showed a high degree of sequence conservation, permitting to explore multiple models as pre-clinical benchmarks. Sequence divergence and incomplete conservation of expression patterns were observed in mouse and rat. Primate Trop-2 sequences were found to be 95%-100% identical to the human sequence. Comparative three-dimension primate Trop-2 structures were obtained with AlphaFold and homology modeling. This revealed high structure conservation of Trop-2 (0.66 ProMod3 GMQE, 0.80-0.86 ± 0.05 QMEANDisCo scores), with conservative amino acid changes at variant sites. Primate TACSTD2/TROP2 cDNAs were cloned and transfectants for individual ORF were shown to be efficiently recognized by humanized anti-Trop-2 monoclonal antibodies (Hu2G10, Hu2EF). Immunohistochemistry analysis of Macaca mulatta (rhesus monkey) tissues showed Trop-2 expression patterns that closely followed those in human tissues. This led us to test Trop-2 targeting in vivo in Macaca fascicularis (cynomolgus monkey). Intravenously injected Hu2G10 and Hu2EF were well tolerated from 5 to 10 mg/kg. Neither neurological, respiratory, digestive, urinary symptoms, nor biochemical or hematological toxicities were detected during 28-day observation. Blood serum pharmacokinetic (PK) studies were conducted utilizing anti-idiotypic antibodies in capture-ELISA assays. Hu2G10 (t1/2 = 6.5 days) and Hu2EF (t1/2 = 5.5 days) were stable in plasma, and were detectable in the circulation up to 3 weeks after the infusion. These findings validate primates as reliable models for Hu2G10 and Hu2EF toxicity and PK, and support the use of these antibodies as next-generation anti-Trop-2 immunotherapy tools.

Schistosoma mansoni vaccine candidates identified by unbiased phage display screening in self-cured rhesus macaques

Schistosomiasis, a challenging neglected tropical disease, affects millions of people worldwide. Developing a prophylactic vaccine against Schistosoma mansoni has been hindered by the parasite's biological complexity. In this study, we utilized the innovative phage-display immunoprecipitation followed by a sequencing approach (PhIP-Seq) to screen the immune response of 10 infected rhesus macaques during self-cure and challenge-resistant phases, identifying vaccine candidates. Our high-throughput S. mansoni synthetic DNA phage-display library encoded 99.6% of 119,747 58-mer peptides, providing comprehensive coverage of the parasite's proteome. Library screening with rhesus macaques' antibodies, from the early phase of establishment of parasite infection, identified significantly enriched epitopes of parasite extracellular proteins known to be expressed in the digestive tract, shifting towards intracellular proteins during the late phase of parasite clearance. Immunization of mice with a selected pool of PhIP-Seq-enriched phage-displayed peptides from MEG proteins, cathepsins B, and asparaginyl endopeptidase significantly reduced worm burden in a vaccination assay. These findings enhance our understanding of parasite-host immune responses and provide promising prospects for developing an effective schistosomiasis vaccine.

PLAT domain protein 1 (PLAT1/PLAFP) binds to the Arabidopsis thaliana plasma membrane and inserts a lipid

Robust agricultural yields depend on the plant's ability to fix carbon amid variable environmental conditions. Over seasonal and diurnal cycles, the plant must constantly adjust its metabolism according to available resources or external stressors. The metabolic changes that a plant undergoes in response to stress are well understood, but the long-distance signaling mechanisms that facilitate communication throughout the plant are less studied. The phloem is considered the predominant conduit for the bidirectional transport of these signals in the form of metabolites, nucleic acids, proteins, and lipids. Lipid trafficking through the phloem in particular attracted our attention due to its reliance on soluble lipid-binding proteins (LBP) that generate and solubilize otherwise membrane-associated lipids. The Phloem Lipid-Associated Family Protein (PLAFP) from Arabidopsis thaliana is generated in response to abiotic stress as is its lipid-ligand phosphatidic acid (PA). PLAFP is proposed to transport PA through the phloem in response to drought stress. To understand the interactions between PLAFP and PA, nearly 100 independent systems comprised of the protein and one PA, or a plasma membrane containing varying amounts of PA, were simulated using atomistic classical molecular dynamics methods. In these simulations, PLAFP is found to bind to plant plasma membrane models independent of the PA concentration. When bound to the membrane, PLAFP adopts a binding pose where W41 and R82 penetrate the membrane surface and anchor PLAFP. This triggers a separation of the two loop regions containing W41 and R82. Subsequent simulations indicate that PA insert into the β-sandwich of PLAFP, driven by interactions with multiple amino acids besides the W41 and R82 identified during the insertion process. Fine-tuning the protein-membrane and protein-PA interface by mutating a selection of these amino acids may facilitate engineering plant signaling processes by modulating the binding response.

Intrinsic disorder in the dynamic evolution of structure, stability, and flexibility of potyviral VLP assemblies: A computational study

An all-atom Molecular Dynamics (MD) study was applied to three viral nanoparticles (VLPs) of Turnip mosaic virus (TuMV), a potyvirus: the particles genetically functionalized with two peptides, VIP (human vasoactive intestinal peptide) and VEGFR (peptide derived from the human receptor 3 of the vascular endothelial growth factor), and the non-functionalized VLP. Previous experimental results showed that VIP-VLP was the only construct of the three that was not viable. VLPs subjected to our MD study were modeled by four complete turns of the particle involving 35 subunits of the coat protein (CP). The MD simulations showed differences in structures and interaction energies associated to the crucial contribution of the disordered N-terminal arms of CP to the global stability of the particle. These differences suggested an overall stability greater in VEGFR-VLP and smaller in VIP-VLP as compared to the unfunctionalized VLP. Our novel MD study of potyviral VLPs revealed essential clues about structure and interactions of these assembled protein particles and suggests that the computational prediction of the viability of VLPs can be a valuable contribution in the field of viral nanobiotechnology.

A single-domain intrabody targeting the follicle-stimulating hormone receptor impacts FSH-induced G protein-dependent signalling

Intracellular variable fragments of heavy-chain antibody from camelids (intra-VHH) have been successfully used as chaperones to solve the 3D structure of active G protein-coupled receptors bound to their transducers. However, their effect on signalling has been poorly explored, although they may provide a better understanding of the relationships between receptor conformation and activity. Here, we isolated and characterized iPRC1, the first intra-VHH recognizing a member of the large glycoprotein hormone receptor family, the follicle-stimulating hormone receptor (FSHR). This intra-VHH recognizes the FSHR third intracellular loop and decreases cAMP production in response to FSH, without altering Gαs recruitment. Hence, iPRC1 behaves as an allosteric modulator and provides a new tool to complete structure/activity studies performed thus far on this receptor.

Molecular cloning, overexpression, characterization, and In silico modelling analysis of a novel GDSL autotransporter-dependent outer membrane lipase (OML) of Pseudomonas guariconensis

The outer membrane lipase (oml) gene, encoding a novel autotransporter-dependent lipase from Pseudomonas guariconensis, was cloned and sequenced. The oml gene has an open reading frame of 1866 bp. It encodes the 621 amino acid autotransporter-dependent GDSL lipase (OML), which has the highest sequence similarity (64.08 %) with the EstA of Pseudomonas aeruginosa (PDB:3kvn.1. A). OML was expressed and purified, which showed a purified band of approximately 70 kDa. The purified enzyme showed maximum activity at pH 9 and 40 °C. Substrate specificity studies and kinetic study by Lineweaver-Burk plot of purified OML showed Km of 1.27 mM and Vmax of 333.33 U/mL with p-nitrophenyl palmitate. The purified enzyme showed good stability in the presence of hexane, methanol, and ethanol, while the presence of the metal ion Mg2+ showed maximum lipase activity. Bioinformatics analysis supported the in vitro findings by predicting enzyme substrate specificity towards long-chain fatty acids and fatty acids with shorter chain lengths. The stability of the interaction of the protein-ligand complex (OML-ricinoleic acid) was confirmed using MDS and castor oil bioconversion using purified OML was confirmed using High-Performance Liquid Chromatography (HPLC).

Outer membrane β-barrel structure prediction through the lens of AlphaFold2

Most proteins found in the outer membrane of gram-negative bacteria share a common domain: the transmembrane β-barrel. These outer membrane β-barrels (OMBBs) occur in multiple sizes and different families with a wide range of functions evolved independently by amplification from a pool of homologous ancestral ββ-hairpins. This is part of the reason why predicting their three-dimensional (3D) structure, especially by homology modeling, is a major challenge. Recently, DeepMind's AlphaFold v2 (AF2) became the first structure prediction method to reach close-to-experimental atomic accuracy in CASP even for difficult targets. However, membrane proteins, especially OMBBs, were not abundant during their training, raising the question of how accurate the predictions are for these families. In this study, we assessed the performance of AF2 in the prediction of OMBBs and OMBB-like folds of various topologies using an in-house-developed tool for the analysis of OMBB 3D structures, and barrOs. In agreement with previous studies on other membrane protein classes, our results indicate that AF2 predicts transmembrane β-barrel structures at high accuracy independently of the use of templates, even for novel topologies absent from the training set. These results provide confidence on the models generated by AF2 and open the door to the structural elucidation of novel transmembrane β-barrel topologies identified in high-throughput OMBB annotation studies or designed de novo.

Identification of two novel and one rare mutation in DYRK1A and prenatal diagnoses in three Chinese families with intellectual Disability-7

Background and purpose: Intellectual disability-7 (MRD7) is a subtype disorder of intellectual disability (MRD) involving feeding difficulties, hypoactivity, and febrile seizures at an age of early onset, then progressive intellectual and physical development deterioration. We purposed to identify the underlying causative genetic factors of three individuals in each Chinese family who presented with symptoms of intellectual disability and facial dysmorphic features. We provided prenatal diagnosis for the three families and genetic counseling for the prevention of this disease. Methods: We collected retrospective clinical diagnostic evidence for the three probands in our study, which included magnetic resonance imaging (MRI), computerized tomography (CT), electroencephalogram (EEG), and intelligence tests for the three probands in our study. Genetic investigation of the probands and their next of kin was performed by Trio-whole exome sequencing (WES). Sanger sequencing or quantitative PCR technologies were then used as the next step to verify the variants confirmed with Trio-WES for the three families. Moreover, we performed amniocentesis to explore the state of the three pathogenic variants in the fetuses by prenatal molecular genetic diagnosis at an appropriate gestational period for the three families. Results: The three probands and one fetus were clinically diagnosed with microcephaly and exhibited intellectual developmental disability, postnatal feeding difficulties, and facial dysmorphic features. Combining probands' clinical manifestations, Trio-WES uncovered the three heterozygous variants in DYRK1A: a novel variant exon3_exon4del p.(Gly4_Asn109del), a novel variant c.1159C>T p.(Gln387*), and a previously presented but rare pathogenic variant c.1309C>T p.(Arg437*) (NM_001396.5) in three families, respectively. In light of the updated American College of Medical Genetic and Genomics (ACMG) criterion, the variant of exon3_exon4del and c.1159C>T were both classified as likely pathogenic (PSV1+PM6), while c1309C>T was identified as pathogenic (PVS1+PS2_Moderate+PM2). Considering clinical features and molecular testimony, the three probands were confirmed diagnosed with MRD7. These three discovered variants were considered as the three causal mutations for MRD7. Prenatal diagnosis detected the heterozygous dominant variant of c.1159C>T p.(Gln387*) in one of the fetuses, indicating a significant probability of MRD7, subsequently the gestation was intervened by the parents' determination and professional obstetrical operation. On the other side, prenatal molecular genetic testing revealed wild-type alleles in the other two fetuses, and their parents both decided to sustain the gestation. Conclusion: We identified two novel and one rare mutation in DYRK1A which has broadened the spectrum of DYRK1A and provided evidence for the diagnosis of MRD7 at the molecular level. Besides, this study has supported the three families with MRD7 to determine the causative genetic factors efficiently and provide concise genetic counseling for the three families by using Trio-WES technology.

The Antiviral Activity of the Lectin Griffithsin against SARS-CoV-2 Is Enhanced by the Presence of Structural Proteins

Although COVID-19 transmission has been reduced by the advent of vaccinations and a variety of rapid monitoring techniques, the SARS-CoV-2 virus itself has shown a remarkable ability to mutate and persist. With this long track record of immune escape, researchers are still exploring prophylactic treatments to curtail future SARS-CoV-2 variants. Specifically, much focus has been placed on the antiviral lectin Griffithsin in preventing spike protein-mediated infection via the hACE2 receptor (direct infection). However, an oft-overlooked aspect of SARS-CoV-2 infection is viral capture by attachment receptors such as DC-SIGN, which is thought to facilitate the initial stages of COVID-19 infection in the lung tissue (called trans-infection). In addition, while immune escape is dictated by mutations in the spike protein, coronaviral virions also incorporate M, N, and E structural proteins within the particle. In this paper, we explored how several structural facets of both the SARS-CoV-2 virion and the antiviral lectin Griffithsin can affect and attenuate the infectivity of SARS-CoV-2 pseudovirus. We found that Griffithsin was a better inhibitor of hACE2-mediated direct infection when the coronaviral M protein is present compared to when it is absent (possibly providing an explanation regarding why Griffithsin shows better inhibition against authentic SARS-CoV-2 as opposed to pseudotyped viruses, which generally do not contain M) and that Griffithsin was not an effective inhibitor of DC-SIGN-mediated trans-infection. Furthermore, we found that DC-SIGN appeared to mediate trans-infection exclusively via binding to the SARS-CoV-2 spike protein, with no significant effect observed when other viral proteins (M, N, and/or E) were present. These results provide etiological data that may help to direct the development of novel antiviral treatments, either by leveraging Griffithsin binding to the M protein as a novel strategy to prevent SARS-CoV-2 infection or by narrowing efforts to inhibit trans-infection to focus on DC-SIGN binding to SARS-CoV-2 spike protein.

Evolution of high-molecular-mass hyaluronic acid is associated with subterranean lifestyle

Hyaluronic acid is a major component of extracellular matrix which plays an important role in development, cellular response to injury and inflammation, cell migration, and cancer. The naked mole-rat (Heterocephalus glaber) contains abundant high-molecular-mass hyaluronic acid in its tissues, which contributes to this species' cancer resistance and possibly to its longevity. Here we report that abundant high-molecular-mass hyaluronic acid is found in a wide range of subterranean mammalian species, but not in phylogenetically related aboveground species. These subterranean mammalian species accumulate abundant high-molecular-mass hyaluronic acid by regulating the expression of genes involved in hyaluronic acid degradation and synthesis and contain unique mutations in these genes. The abundant high-molecular-mass hyaluronic acid may benefit the adaptation to subterranean environment by increasing skin elasticity and protecting from oxidative stress due to hypoxic conditions. Our work suggests that high-molecular-mass hyaluronic acid has evolved with subterranean lifestyle.

Morphological and molecular phylogenetic characterization of Sarcocystis kani sp. nov. and other novel, closely related Sarcocystis spp. infecting small mammals and colubrid snakes in Asia

We investigated the morphology and phylogenetic relationships of novel and previously recognized Sarcocystis spp. infecting small mammals and colubrid snakes in Asia. The nuclear 18S rRNA and mitochondrial cox1 of Sarcocystis sp.1 from mangrove snakes (Boiga dendrophila) in Thailand and Sarcocystis sp.2 from a ricefield rat (Rattus argentiventer) in Sumatra were partially sequenced. Sporocysts of Sarcocystis sp.1 induced development of sarcocysts in experimentally infected rats, which showed a unique ultrastructure that was observed previously by S.P. Kan in rats from Malaysia; therefore, we describe this species as Sarcocystis kani sp. nov. Its integration into the 18S rRNA phylogeny of Sarcocystis spp. cycling between small mammals and colubrid snakes helped clarify relationships among the so-called S. zuoi-complex of molecularly cryptic species: Sarcocystis kani sp. nov., S. sp.2, S. attenuati, S. scandentiborneensis, and S. zuoi were all included in this clade. Tree topology was resolved into dichotomies congruent with the morphological disparities between the taxa. However, cox1 gene sequencing (including newly sequenced S. singaporensis and S. zamani) revealed that Sarcocystis kani, S. attenuati, and S. scandentiborneensis were identical suggesting a recent, common ancestry. To identify other distinctive features, lineage-specific molecular patterns within both genes were examined revealing that all 18S rRNA sequences of the S. zuoi - complex possess a unique, 7-nt long motif in helix 38 of domain V7 that was different in S. clethrionomyelaphis which branched off basally from the complex. Three-dimensional homology modelling of COX1 protein structure identified amino acid substitutions within the barcode area specific for the S. zuoi-complex and substantial divergence in structurally important amino acids between Sarcocystis species of snakes as definitive hosts and other lineages of the Sarcocystidae. We discuss the utility of selected genes for species delimitation of the Sarcocystis spp. under investigation, which probably evolved during recent radiations of their intermediate and definitive hosts.

Improved multimer prediction using massive sampling with AlphaFold in CASP15

AlphaFold2 has revolutionized structure prediction by achieving high accuracy comparable to experimentally determined structures. However, there is still room for improvement, especially for challenging cases like multimers. A key to the success of AlphaFold is its ability to assess and rank its own predictions. Our basic idea for the Wallner group in CASP15 was to exploit this excellent scoring function in AlphaFold by massive sampling. To achieve this goal, we conducted AlphaFold runs using six different settings, using templates, without templates, and with an increased number of recycles for both multimer v1 and v2 weights. In all instances, we enabled dropout layers during inference, allowing for sampling of uncertainty and enhancing the diversity of the generated models. In total, 274 289 models were generated for the 38 targets in CASP15, with a median of 4810 models per target. Of these 38 targets, 10 were high quality, 11 were medium quality, 11 were acceptable, and only 6 were incorrect. The improvement over the baseline method, NBIS-AF2-multimer, is substantial, with the mean DockQ increasing from 0.43 to 0.56, with several targets showing a DockQ score increase of +0.6 units. Remarkable, considering Wallner and NBIS-AF2-multimer were using identical input data. The success can be attributed to the diversified sampling using dropout with different settings and, in particular, the use of multimer v1, which is much more susceptible to sampling compared with v2. The method is available here: http://wallnerlab.org/AFsample/.

Estimating protein complex model accuracy based on ultrafast shape recognition and deep learning in CASP15

This article reports and analyzes the results of protein complex model accuracy estimation by our methods (DeepUMQA3 and GraphGPSM) in the 15th Critical Assessment of techniques for protein Structure Prediction (CASP15). The new deep learning-based multimeric complex model accuracy estimation methods are proposed based on the ensemble of three-level features coupling with deep residual/graph neural networks. For the input multimeric complex model, we describe it from three levels: overall complex features, intra-monomer features, and inter-monomer features. We designed an overall ultrafast shape recognition (USR) to characterize the relationship between local residues and the overall complex topology, and an inter-monomer USR to characterize the relationship between the residues of one monomer and the topology of other monomers. DeepUMQA3 (Group name: GuijunLab-RocketX) ranked first in the interface residue accuracy estimation of CASP15. The Pearson correlation between the interface residue Local Distance Difference Test (lDDT) predicted by DeepUMQA3 and the real lDDT is 0.570, the only method that exceeds 0.5. Among the top 5 methods, DeepUMQA3 achieved the highest Pearson correlation of lDDT on 25 out of 39 targets. GraphGPSM (Group name: GuijunLab-PAthreader) has TM-score Pearson correlations greater than 0.9 on 14 targets, showing a good ability to estimate the overall fold accuracy. The DeepUMQA3 server is available at http://zhanglab-bioinf.com/DeepUMQA/ and the GraphGPSM server is available at http://zhanglab-bioinf.com/GraphGPSM/.

Assessment of the assessment-All about complexes

Predicting model quality is a fundamental component of any modeling procedure, and blind assessment of these methods constitutes a crucial aspect of the Critical Assessment of Protein Structure Prediction (CASP) experiment. Historically, the main focus was on assessing methods that predict global and per-residue accuracies in tertiary structure models. This focus shifted with the community's increased efforts in modeling complexes and assemblies. We asked the community to process the models from the CASP15 assembly category and provide estimates of the accuracy of the predicted quaternary structure, both globally and at the local interface level. Besides identifying remarkable accuracy of modeling groups in assessing their own predictions, we set up a benchmarking pipeline to highlight different aspects of quaternary structure models and introduced a simple consensus EMA method as baseline. While participating methods showed commendable performance, the baseline was difficult to surpass. It is important to point out that prediction performance varies for the individual CASP targets, highlighting potential areas of improvement and challenges ahead.

Glycoside hydrolases in the biodegradation of lignocellulosic biomass

Lignocellulose is a plentiful and intricate biomass substance made up of cellulose, hemicellulose, and lignin. Cellulose and hemicellulose are polysaccharides characterized by different compositions and degrees of polymerization. As renewable resources, their applications are eco-friendly and can help reduce reliance on petrochemical resources. This review aims to illustrate cellulose, hemicellulose, and their structures and hydrolytic enzymes. To obtain desirable enzyme sources for the high hydrolysis of lignocellulose, highly stable, efficient and thermophilic enzyme sources, and new technologies, such as rational design and machine learning, have been introduced in detail. Generally, the efficient biodegradation of abundant natural biomass into fermentable sugars or other intermediates has great potential in practical applications.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03819-1.

Cooperative and structural relationships of the trimeric Spike with infectivity and antibody escape of the strains Delta (B.1.617.2) and Omicron (BA.2, BA.5, and BQ.1)

Herein, we conducted simulations of trimeric Spike from several SARS-CoV-2 variants of concern (Delta and Omicron sub-variants BA.2, BA.5, and BQ.1) and investigated the mechanisms by which specific mutations confer resistance to neutralizing antibodies. We observed that the mutations primarily affect the cooperation between protein domains within and between protomers. The substitutions K417N and L452R expand hydrogen bonding interactions, reducing their interaction with neutralizing antibodies. By interacting with nearby residues, the K444T and N460K mutations in the SpikeBQ.1 variant potentially reduces solvent exposure, thereby promoting resistance to antibodies. We also examined the impact of D614G, P681R, and P681H substitutions on Spike protein structure that may be related to infectivity. The D614G substitution influences communication between a glycine residue and neighboring domains, affecting the transition between up- and -down RBD states. The P681R mutation, found in the Delta variant, enhances correlations between protein subunits, while the P681H mutation in Omicron sub-variants weakens long-range interactions that may be associated with reduced fusogenicity. Using a multiple linear regression model, we established a connection between inter-protomer communication and loss of sensitivity to neutralizing antibodies. Our findings underscore the importance of structural communication between protein domains and provide insights into potential mechanisms of immune evasion by SARS-CoV-2. Overall, this study deepens our understanding of how specific mutations impact SARS-CoV-2 infectivity and shed light on how the virus evades the immune system.

New prediction categories in CASP15

Prediction categories in the Critical Assessment of Structure Prediction (CASP) experiments change with the need to address specific problems in structure modeling. In CASP15, four new prediction categories were introduced: RNA structure, ligand-protein complexes, accuracy of oligomeric structures and their interfaces, and ensembles of alternative conformations. This paper lists technical specifications for these categories and describes their integration in the CASP data management system.

A computational simulation appraisal of banana lectin as a potential anti-SARS-CoV-2 candidate by targeting the receptor-binding domain

BACKGROUND

The ongoing concern surrounding coronavirus disease 2019 (COVID-19) primarily stems from continuous mutations in the genome of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), leading to the emergence of numerous variants. The receptor-binding domain (RBD) in the S1 subunit of the S protein of the virus plays a crucial role in recognizing the host's angiotensin-converting enzyme 2 (hACE2) receptor and facilitating cell membrane fusion processes, making it a potential target for preventing viral entrance into cells. This research aimed to determine the potential of banana lectin (BanLec) proteins to inhibit SARS-CoV-2 attachment to host cells by interacting with RBD through computational modeling.

MATERIALS AND METHODS

The BanLecs were selected through a sequence analysis process. Subsequently, the genes encoding BanLec proteins were retrieved from the Banana Genome Hub database. The FGENESH online tool was then employed to predict protein sequences, while web-based tools were utilized to assess the physicochemical properties, allergenicity, and toxicity of BanLecs. The RBDs of SARS-CoV-2 were modeled using the SWISS-MODEL in the following step. Molecular docking procedures were conducted with the aid of ClusPro 2.0 and HDOCK web servers. The three-dimensional structures of the docked complexes were visualized using PyMOL. Finally, molecular dynamics simulations were performed to investigate and validate the interactions of the complexes exhibiting the highest interactions, facilitating the simulation of their dynamic properties.

RESULTS

The BanLec proteins were successfully modeled based on the RNA sequences from two species of banana (Musa sp.). Moreover, an amino acid modification in the BanLec protein was made to reduce its mitogenicity. Theoretical allergenicity and toxicity predictions were conducted on the BanLecs, which suggested they were likely non-allergenic and contained no discernible toxic domains. Molecular docking analysis demonstrated that both altered and wild-type BanLecs exhibited strong affinity with the RBD of different SARS-CoV-2 variants. Further analysis of the molecular docking results showed that the BanLec proteins interacted with the active site of RBD, particularly the key amino acids residues responsible for RBD's binding to hACE2. Molecular dynamics simulation indicated a stable interaction between the Omicron RBD and BanLec, maintaining a root-mean-square deviation (RMSD) of approximately 0.2 nm for a duration of up to 100 ns. The individual proteins also had stable structural conformations, and the complex demonstrated a favorable binding-free energy (BFE) value.

CONCLUSIONS

These results confirm that the BanLec protein is a promising candidate for developing a potential therapeutic agent for combating COVID-19. Furthermore, the results suggest the possibility of BanLec as a broad-spectrum antiviral agent and highlight the need for further studies to examine the protein's safety and effectiveness as a potent antiviral agent.