A Global Ramachandran Score Identifies Protein Structures with Unlikely Stereochemistry

Computer Integrated Dominant Epitopes Evoke Protective Immune Response Against Streptococcus pneumoniae

Streptococcus pneumoniae is a gram-positive bacterium responsible for various diseases like pneumonia, acute otitis media, sinusitis, meningitis and bacteraemia. These diseases cause a significant amount of morbidity and mortality. Although polysaccharide vaccines are available, the protection provided by these vaccines is serotype-dependent and not enough in sensitive populations like children and older people. Designing a subunit vaccine by using proteins that are responsible for the pathogenesis of diseases can provide better protection against bacterial infections. In this study, we present the design of a novel multi-epitope vaccine against Streptococcus pneumoniae using an immunoinformatic approach. More than 1170 epitopes were identified against B cells, cytotoxic T lymphocytes and helper T lymphocytes from more than 60 pneumococcal proteins. Epitopes were further screened, and potential epitopes were selected for vaccine development. Seven different vaccine combinations that harbour the 15 dominant B-cell, cytotoxic T cell and helper T cell epitopes were evaluated with linker and β-defensin adjuvant to finalise the best vaccine construct. Bioinformatics tools were used to analyse the construct's physicochemical properties, secondary and tertiary structures, allergenicity, antigenicity and immunogenicity. Docking studies with the TLR-4 receptor and molecular dynamics simulations indicated strong binding affinity and stability. In silico immune response simulations predicted robust IgG immune response generation and observed more than 200 000 IgG1 + IgG2 counts per mL. Similarly, cell-mediated immunity was also enhanced by the designed vaccine construct. The construct was codon-optimised and cloned in silico for expression in Escherichia coli. These findings suggest that the construct is a promising candidate for further experimental validation.

Molecular Dynamic Study on the Structure and Thermal Stability of Mutant Pediocin PA-1 Peptides Engineered with Cysteine Substitutions

Pediocin and analogous bacteriocins, valued for thermal stability, serve as versatile antimicrobials in the food sector. Improving their resilience at high temperatures and deriving derivatives not only benefit food production but also offer broad-spectrum antimicrobial potential in pharmaceuticals, spanning treatments for peptic ulcers, women's health, and novel anticancer agents. The study aims to create mutant peptides capable of establishing a third disulfide bond or enhanced through cysteine substitutions. This involves introducing additional Cys residues into the inherent structure of pediocin PA-1 to facilitate disulfide bond formation. Five mutants (Mut 1-5) were systematically generated with double Cys substitutions and assessed for thermal stability through MD simulations across temperatures (298-394 K). The most robust mutants (Mut 1, Mut 4-5) underwent extended analysis via MD simulations, comparing their structural stability, secondary structure, and surface accessibility to the reference Pediocin PA-1 molecule. This comprehensive assessment aims to understand how Cys substitutions influence disulfide bonds and the overall thermal stability of the mutant peptides. In silico analysis indicated that Mut 1 and Mut 5, along with the reference structure, lose their helical structure and one natural disulfide bond at high temperatures, and may impacting antimicrobial activity. Conversely, Mut 4 retained its helical structure and exhibited thermal stability similar to Pediocin PA-1. Pending further experimental validation, this study implies Mut 4 may have high stability and exceptional resistance to high temperatures, potentially serving as an effective antimicrobial alternative.

Transport of phenoxyacetic acid herbicides by PIN-FORMED auxin transporters

Auxins are a group of phytohormones that control plant growth and development. Their crucial role in plant physiology has inspired development of potent synthetic auxins that can be used as herbicides. Phenoxyacetic acid derivatives are a widely used group of auxin herbicides in agriculture and research. Despite their prevalence, the identity of the transporters required for distribution of these herbicides in plants is both poorly understood and the subject of controversial debate. Here we show that PIN-FORMED auxin transporters transport a range of phenoxyacetic acid herbicides across the membrane. We go on to characterize the molecular determinants of substrate specificity using a variety of different substrates as well as protein mutagenesis to probe the binding site. Finally, we present cryogenic electron microscopy structures of Arabidopsis thaliana PIN8 bound to either 2,4-dichlorophenoxyacetic acid or 4-chlorophenoxyacetic acid. These structures represent five key states from the transport cycle, allowing us to describe conformational changes associated with the transport cycle. Overall, our results reveal that phenoxyacetic acid herbicides use the same export machinery as endogenous auxins and exemplify how transporter binding sites undergo transformations that dictate substrate specificity. These results provide a foundation for future development of novel synthetic auxins and for precision breeding of herbicide-resistant crop plants.

Functional characterization and structural prediction of hypothetical proteins in monkeypox virus and identification of potential inhibitors

The excessive activation of the monkeypox virus (MPXV-Congo_8-156) is linked to various skin and respiratory disorders such as rashes, fluid-filled blisters, swollen lymph nodes and encephalitis (inflammation of the brain), highlighting MPXV-Congo_8-156 as a promising target for drug intervention. Despite the effectiveness of Cidofovir, in inhibiting MPXV activity, its limited ability to penetrate the skin and its strong side effects restrict its application. To address this challenge, we screened 500 compounds capable of penetrating the skin and gastrointestinal tract to identify potent MPXV inhibitors. Various characterization schemes and structural models of MPXV-Congo_8-156 were explored with bioinformatics tools like PROTPARAM, SOPMA, SWISS-MODEL and PROCHECK. Using molecular docking in PyRx, we evaluated the binding affinities of these compounds with MPXV-Congo_8-156 and identified the top five candidates ranging from - 9.2 to - 8.8 kcal/mol. ADMET analysis indicated that all five compounds were safer alternatives, showing no AMES toxicity or carcinogenicity in toxicological assessments. Molecular dynamics (MD) simulations, conducted for 100 ns each, confirmed the docking interactions of the top five compounds alongside the control (Cidofovir), validating their potential as MPXV inhibitors. The compounds with PubChem CID numbers 4061636, 4422538, 3583576, 4856107 and 4800629 demonstrated strong support in terms of root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), radius of gyration (Rg), solvent-accessible surface area (SASA) value, hydrogen bond analysis, and Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) analysis. Thus, our investigation identified these five compounds as promising inhibitors of MPXV, offering potential therapeutic avenues. However, further in vivo studies are necessary to validate our findings.

Molecular dynamic simulations to assess the structural variability of ClpV from Enterobacter cloacae

The Enterobacter cloacae complex (ECC) consists of six Enterobacter species (E. cloacae, hormaechei, kobei, ludwigii, nimipressuralis and asburiae) that have emerged as nosocomial pathogens of interest, with E. cloacae and Enterobacter hormachei being the most frequently isolated ECC species in human clinical specimens and intensive care unit (ICU) patients. Many nosocomial outbreaks of E. cloacae have been related to transmission through contaminated surgical equipment and operative cleaning solutions. As this pathogen evades the action of antibiotics, it is important to find alternative targets to limit the devastating effects of these pathogens. ClpV is a Clp ATPase which dissociates and recycles the contracted sheath of the bacterial type VI secretion system (T6SS), thereby regulating bacterial populations and facilitating environmental colonization. Seventy-one Enterobacter strains were mined for Clp ATPase proteins. All the investigated strains contained ClpA, ClpB, ClpX and ClpV while only 20% contained ClpK. All the investigated strains contained more than one ClpV protein, and the ClpV proteins showed significant variations. Three ClpV proteins from E. cloacae strain E3442 were then investigated to determine the structural difference between each protein. Homology modelling showed the proteins to be structurally similar to each other, however the physicochemical characteristics of the proteins vary. Additionally, physicochemical analysis and molecular dynamic simulations showed that the proteins were highly dynamic and not significantly different from each other. Further investigation of the proteins in silico and in vitro in the presence and absence of various ligands and proteins could be performed to determine whether the proteins all interact with their surroundings in the same manner. This would allow one to determine why multiple homologs of the same protein are expressed by pathogens.

Anxiogenic-like effects of coumarin, possibly through the GABAkine interaction pathway: Animal studies with in silico approaches

BACKGROUND

Anxiety disorder is the most common mental illness and a major contributor to impairment. Thus, there is an urgent need to find novel lead compounds to mitigate anxiety. It is widely recognized that the neurobiology of anxiety-related behavior involves GABAergic systems.

OBJECTIVES

This research aimed to examine the anxiogenic action of coumarin (CMN), a natural benzopyrone derived from plants, and determine its underlying mechanism through in vivo and in silico investigations.

METHODS

This was accomplished by using a variety of behavioral procedures, including open field, swing, hole cross, and light-dark tests, on male and female Swiss albino mice that had been orally administered three experimental doses of CMN (1, 2, and 4 mg/kg). The CMN group was also examined with the GABAA receptor agonist diazepam (DZP, 2 mg/kg) and flumazenil antagonist (FLU, 0.1 mg/kg). Furthermore, CMN and standards were subjected to a molecular docking analysis to determine their binding affinities for the GABAA receptor subunits (α1, α4, β2, γ2, and δ). Several software programs were used to visualize the ligand-receptor interaction and analyze the pharmacokinetic profile.

RESULTS

Compared to typical treatments, our results show that CMN (1 mg/kg) significantly (p < 0.05) increases the locomotor activity of animals. Furthermore, CMN exerted the highest binding affinity (-6.5 kcal/mol) with the GABA-α1 receptor compared to conventional DZP. Along with FLU, CMN displayed several hydrophobic and hydrogen bonds with GABAA receptor subunits. The pharmacokinetic and drug-like properties of CMN are also remarkable. In animal studies, CMN worked synergistically with FLU to provide anxiogenic-like effects.

CONCLUSION

We conclude that, based on in vivo and in silico data, CMN, alone or in combination with FLU, may be employed in future neurological clinical studies. However, further research is needed to confirm this behavioral activity and elucidate the possible mechanism of action.

Phytonutrients as a Defensive Barrier Against G Ectodomain Fusion in Chandipura Virus Infection

Viruses, microscopic menace that transcends time leaving its mark on every era have been silent predators since the dawn of civilization, evolving with us and shaping our history. Chandipura virus (CHPV), a potent member of the Rhabdoviridae family poses a significant threat in India with rapid neuroinvasive potential leading to fatal encephalitis, particularly in children. Given the scarcity of research, our study consolidates critical information regarding its lifecycle, fusion process, and reviewed the LRP1 and GRP78 as CHPV target receptors. With no FDA-approved drugs currently available for CHPV prevention, our research focuses on identifying potential molecules that can disrupt the virus at its most critical juncture, the fusion stage. The results derived from compounds screening indicated Silibinin, 3-(2,3-Dihydroxy-3-Methylbutyl)-6-Hydroxy-2-[(1E,5E)-3,4,10-Trihydroxyundeca-1,5-Dienyl] Benzaldehyde, Budmunchiamine L5, and L4 as a leading molecule may efficaciously inhibit G ectodomain fusion. By analyzing pharmacokinetic properties through radar graph, outcomes support the nomination of four compounds as potential inhibitory molecules and ensure they possess the optimal balance of drug-like characteristics. Working with the CHPV presents significant challenges, making the in silico parameters crucial in guiding future research. Our study sought to pioneer the discovery of therapeutic molecules against the CHPV, providing a foundational framework for developing effective antiviral strategies.

Anxiolytic Efficacy of Indirubin: In Vivo Approach Along with Receptor Binding Profiling and Molecular Interaction with GABAergic Pathways

Anxiety is a natural response to stress, characterized by feelings of worry, fear, or unease. The current research was conducted to investigate the anxiolytic effect of indirubin (IND) in different behavioral paradigms in Swiss albino mice. To observe the animal's behavioural response to assess anxiolytic activity, different tests were performed, such as the open-field (square cross, grooming, and rearing), swing, dark-light, and hole cross tests. The experimental mice were administered IND (5 and 10 mg/kg, p.o.), where diazepam (DZP) and vehicle were used as positive and negative controls, respectively. In addition, a combination treatment (DZP+IND-10) was provided to the animals to determine the modulatory effect of IND on DZP. Molecular docking approach was also conducted to determine the binding energy of IND with the GABAA receptor (α2 and α3 subunits) and pharmacokinetics were also estimated. The findings revealed that IND dose-dependently significantly (p<0.05) reduced the animal's movement exerting calming behavior like DZP. IND also demonstrated the highest docking score (-7.7 kcal/mol) against the α3 subunit, while DZP showed a lower docking value (-6.4 kcal/mol) than IND. The ADMET analysis revealed that IND has proper drug-likeness and pharmacokinetic characteristics. In conclusion, IND exerted anxiolytic effects through GABAergic Pathways.

Molecular dynamics simulation based prediction of T-cell epitopes for the production of effector molecules for liver cancer immunotherapy

Liver cancer is the sixth most frequent malignancy and the fourth major cause of deaths worldwide. The current treatments are only effective in early stages of cancer. To overcome the therapeutic challenges and exploration of immunotherapeutic options, broad spectral therapeutic vaccines could have significant impact. Based on immunoinformatic and integrated machine learning tools, we predicted the potential therapeutic vaccine candidates of liver cancer. In this study, machine learning and MD simulation-based approach are effectively used to design T-cell epitopes that aid the immune system against liver cancer. Antigenicity, molecular weight, subcellular localization and expression site predictions were used to shortlist liver cancer associated proteins including AMBP, CFB, CDHR5, VTN, APOBR, AFP, SERPINA1 and APOE. We predicted CD8+ T-cell epitopes of these proteins containing LGEGATEAE, LLYIGKDRK, EDIGTEADV, QVDAAMAGR, HLEARKKSK, HLCIRHEMT, LKLSKAVHK, EQGRVRAAT and CD4+ T-cell epitopes of VLGEGATEA, WVTKQLNEI, VEEDTKVNS, FTRINCQGK, WGILGREEA, LQDGEKIMS, VKFNKPFVF, VRAATVGSL. We observed the substantial physicochemical properties of these epitopes with a significant binding affinity with MHC molecules. A polyvalent construct of these epitopes was designed using suitable linkers and adjuvant indicated significant binding energy (>-10.5 kcal/mol) with MHC class-I and II molecule. Based on in silico cloning, we found the considerable compatibility of this polyvalent construct with the E. coli expression system and the efficiency of its translation in host. The system-level and machine learning based cross validations showed the possible effect of these T-cell epitopes as potential vaccine candidates for the treatment of liver cancer.

Omics-centric evidences of fipronil biodegradation by Rhodococcus sp. FIP_B3

The widespread use of the pesticide fipronil in domestic and agriculture sectors has resulted in its accumulation across the environment. Its use to assure food security has inadvertently affected soil microbiome composition, fertility and, ultimately, human health. Degradation of residual fipronil present in the environment using specific microbial species is a promising strategy for its removal. The present study delves into the omics approach for fipronil biodegradation using the native bacterium Rhodococcus sp. FIP_B3. It has been observed that within 40 days, nearly 84% of the insecticide gets degraded. The biodegradation follows a pseudo-first-order kinetics (k = 0.0197/d with a half-life of ∼11 days). Whole genome analysis revealed Cytochrome P450 monooxygenase, peroxidase-related enzyme, haloalkane dehalogenase, 2-nitropropane dioxygenase, and aconitate hydratase are involved in the degradation process. Fipronil-sulfone, 5-amino-1-(2-chloro-4-(trifluoromethyl)phenyl)-4- ((trifluoromethyl)sulfonyl)-1H-pyrazole-3-carbonitrile, (E)-5-chloro-2-oxo-3- (trifluoromethyl)pent-4-enoic acid, 4,4,4-trifluoro-2-oxobutanoic acid, and 3,3,3- trifluoropropanoic acid were identified as the major metabolites that support the bacterial degradation of fipronil. In-silico molecular docking and molecular dynamic simulation-based analyses of degradation pathway intermediates with their respective enzymes have indicated stable interactions with significant binding energies (-5.9 to -9.7 kcal/mol). These results have provided the mechanistic cause of the elevated potential of Rhodococcus sp. FIP_B3 for fipronil degradation and will be advantageous in framing appropriate strategies for the bioremediation of fipronil-contaminated environment.

Binding properties of olfactory proteins to host volatiles, free fatty acids and cuticular hydrocarbons in the termite Reticulitermes aculabialis

As eusocial insects prevalent in tropical and subtropical regions, termites are characterized by highly organized behaviors and exceptional adaptability, rooted in caste differentiation and chemical communication. These traits make them excellent models for studying insect social structures and ecological interactions. Investigating how termites use chemical signals to perceive and respond to their environment provides insights into their coordination and adaptation within complex ecosystems. This study delved into the chemosensory mechanisms of Reticulitermes aculabialis, examining the interactions of four olfactory proteins with 70 ligands, including host volatiles, cuticular hydrocarbons (CHCs), and free fatty acids (FFAs). Molecular docking simulations revealed varied affinities of the olfactory proteins for long-chain hydrocarbons (n-C23 to n-C28), suggesting a nuanced chemical communication system through specific hydrocarbon detection. RacuCSP1 and RacuCSP2 exhibited specific binding to linoleic acid and undecanoic acid, respectively, highlighting the significance of FFAs in the physiological and behavioral processes of termites. The four olfactory proteins showed a strong affinity for longifolene in fluorescence competitive binding experiments. Notably, RacuOBPs exhibited unique affinities for terpenoid volatiles such as β-lonone and neocembrene, while RacuCSPs specifically bound with terpenoids like 3-carene, myrtenol, α-pinene oxide and β-pinene indicating their critical roles in host detection. Behavioral observations following gene silencing revealed that RacuOBP5 was essential for recognizing longifolene and α-lonone recognition, while RacuCSP1 was key for detecting α-pinene in termites. These findings enhance our understanding of the termite chemosensory system and offer insights for developing precise pest management strategies.

In silico design and assessment of a multi-epitope peptide vaccine against multidrug-resistant Acinetobacter baumannii

Acinetobacter baumannii, an opportunistic and notorious nosocomial pathogen, is responsible for many infections affecting soft tissues, skin, lungs, bloodstream, and urinary tract, accounting for more than 722,000 cases annually. Despite the numerous advancements in therapeutic options, no approved vaccine is currently available for this particular bacterium. Consequently, this study focused on creating a rational vaccine design using bioinformatics tools. Three outer membrane proteins with immunogenic potential and properties of good vaccine candidates were used to select epitopes based on good antigenic properties, non-allergenicity, high binding scores, and a low IC50 value. A multi-epitope peptide (MEP) construct was created by sequentially linking the epitopes using suitable linkers. ClusPro 2.0 and C-ImmSim web servers were used for docking analysis with TLR2/TLR4 and immune response respectively. The Ramachandran plot showed an accurate model of the MEP with 100% residue in the most favored and allowed regions. The construct was highly antigenic, stable, non-allergenic, non-toxic, and soluble, and showed maximum population coverage. Additionally, molecular docking demonstrated strong binding between the designed MEP vaccine and TLR2/TLR4. In silico immunological simulations showed significant increases in T-cell and B-cell populations. Finally, codon optimization and in silico cloning were conducted using the pET-28a (+) plasmid vector to evaluate the efficiency of the expression of vaccine peptide in the host organism (Escherichia coli). This designed MEP vaccine would support and accelerate the laboratory work to develop a potent vaccine targeting MDR Acinetobacter baumannii.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-024-00292-3.

Investigating the effects of the ARG258HIS mutation on RAD51C in inherited Fanconi Anemia and cancer disease

Fanconi anemia is a rare chromosomal instability disorder associated with developmental abnormalities, bone marrow failure, and a heightened susceptibility to leukemia and other cancers. It is an autosomal recessive genetic disorder, necessitating both parents to carry the faulty gene. Diagnostic methods include blood tests, chromosome breakage assessments, and genetic testing. While there is no cure, treatments encompass blood transfusions, bone marrow transplants, and gene therapy, with patients requiring regular check-ups, supportive care, and cancer screening to enhance their quality of life. In this study, we identify a specific substitution (R258H) targeting the crucial binding site of the alpha-helix region in RAD51C. This substitution induces structural disorder in distinct regions, as indicated by the near absence of electron density for multiple amino acids. Intriguingly, these disordered regions do not follow a continuous sequence from the mutation site and extend across domain boundaries. We utilized computational prediction algorithms and Molecular Dynamics (MD) simulations to model RAD51C and its mutation (R258H) structurally. These simulations highlighted alterations in conformational dynamics, the Free Energy Landscape (FEL), and intrinsic molecular motions induced by the mutation, suggesting structural destabilization that could disrupt its function. This observed destabilization in RAD51C due to mutations offers valuable insights that may serve as diagnostic markers for individuals carrying these mutations, particularly in Fanconi anemia.

Antifungal activity of citronellal against Trichophyton rubrum and its predictive mechanism of action by CYP51 inhibition through molecular docking

The present study aimed to investigate the antifungal activity of citronellal (CIT) against clinical isolates of T. rubrum and to show the possible mechanism of action involved. The antifungal potential of CIT was evaluated from the Minimum Inhibitory Concentration (MIC), Minimum Fungicide Concentration (MFC) and assays with ergosterol and sorbitol, to elucidate the possible mechanisms of action, and molecular docking. MIC and MFC values ranged from 4 to 512 µg/mL. Regarding the mechanism of action, the monoterpene demonstrated interaction with fungal ergosterol. In addition, it is possible to observe that CIT acts on crucial enzymes for the biosynthesis and maintenance of the fungal cell membrane, due to the ability of the monoterpene to bind to CYP51. The results obtained in this research demonstrate that CIT has the potential to become, in the future, a product for the treatment of dermatophytosis.

In Silico Investigation of the Interactions Between Cotton Leaf Curl Multan Virus Proteins and the Transcriptional Gene Silencing Factors of Gossypium hirsutum L

The highly dynamic nature of the Cotton leaf curl virus (CLCuV) complex (causing Cotton leaf curl disease, a significant global threat to cotton) presents a formidable challenge in unraveling precise molecular mechanisms governing viral-host interactions. To address this challenge, the present study investigated the molecular interactions of 6 viral proteins (Rep, TrAP, C4, C5, V2, and βC1) with 18 cotton Transcriptional Gene Silencing (TGS) proteins. Protein-protein dockings conducted for different viral-host protein pairs using Clustered Protein Docking (ClusPro) and Global RAnge Molecular Matching (GRAMM) (216 docking runs), revealed variable binding energies. The interacting pairs with the highest binding affinities were further scrutinized using bioCOmplexes COntact MAPS (COCOMAPS) server, which revealed robust binding of three viral proteins- TrAP, C4, and C5 with 14 TGS proteins, identifying several novel interactions (not reported yet by earlier studies), such as TrAP targeting DCL3, HDA6, and SUVH6; C4 targeting RAV2, CMT2, and DMT1; and C5 targeting CLSY1, RDR1, RDR2, AGO4, SAMS, and SAHH. Visualizing these interactions in PyMol provided a detailed insight into interacting regions. Further assessment of the impact of 18 variants of the C4 protein on interaction with CMT2 revealed no correlation between sequence variation and docking energies. However, conserved residues in the C4 binding regions emerged as potential targets for disrupting viral integrity. Hence, this study provides valuable insights into the viral-host interplay, advancing our understanding of Cotton leaf curl Multan virus pathogenicity and opening novel avenues for devising various antiviral strategies by targeting the host-viral interacting regions after experimental validation.

Structural basis for α-tubulin-specific and modification state-dependent glutamylation

Microtubules have spatiotemporally complex posttranslational modification patterns. Tubulin tyrosine ligase-like (TTLL) enzymes introduce the most prevalent modifications on α-tubulin and β-tubulin. How TTLLs specialize for specific substrate recognition and ultimately modification-pattern generation is largely unknown. TTLL6, a glutamylase implicated in ciliopathies, preferentially modifies tubulin α-tails in microtubules. Cryo-electron microscopy, kinetic analysis and single-molecule biochemistry reveal an unprecedented quadrivalent recognition that ensures simultaneous readout of microtubule geometry and posttranslational modification status. By binding to a β-tubulin subunit, TTLL6 modifies the α-tail of the longitudinally adjacent tubulin dimer. Spanning two tubulin dimers along and across protofilaments (PFs) ensures fidelity of recognition of both the α-tail and the microtubule. Moreover, TTLL6 reads out and is stimulated by glutamylation of the β-tail of the laterally adjacent tubulin dimer, mediating crosstalk between α-tail and β-tail. This positive feedback loop can generate localized microtubule glutamylation patterns. Our work uncovers general principles that generate tubulin chemical and topographic complexity.

Structural investigation, computational analysis, and theoretical cryoprotectant approach of antifreeze protein type IV mutants

Antifreeze proteins (AFPs) have unique features to sustain life in sub-zero environments due to ice recrystallization inhibition (IRI) and thermal hysteresis (TH). AFPs are in demand as agents in cryopreservation, but some antifreeze proteins have low levels of activity. This research aims to improve the cryopreservation activity of an AFPIV. In this in silico study, the helical peptide afp1m from an Antarctic yeast AFP was modeled into a sculpin AFPIV, to replace each of its four α-helices in turn, using various computational tools. Additionally, a new linker between the first two helices of AFPIV was designed, based on a flounder AFPI, to boost the ice interaction activity of the mutants. Bioinformatics tools such as ExPASy Prot-Param, Pep-Wheel, SOPMA, GOR IV, Swiss-Model, Phyre2, MODFOLD, MolPropity, and ProQ were used to validate and analyze the structural and functional properties of the model proteins. Furthermore, to evaluate the AFP/ice interaction, molecular dynamics (MD) simulations were executed for 20, 100, and 500 ns at various temperatures using GROMACS software. The primary, secondary, and 3D modeling analysis showed the best model for a redesigned antifreeze protein (AFP1mb, with afp1m in place of the fourth AFPIV helix) with a QMEAN (Swiss-Model) Z score value of 0.36, a confidence of 99.5%, a coverage score of 22%, and a p value of 0.01. The results of the MD simulations illustrated that AFP1mb had more rigidity and better ice interactions as a potential cryoprotectant than the other models; it also displayed enhanced activity in limiting ice growth at different temperatures.

A pan-genomic analysis based multi-epitope vaccine development by targeting Stenotrophomonas maltophilia using reverse vaccinology method: an in-silico approach

Antibiotic resistance in bacteria leads to high mortality rates and healthcare costs, a significant concern for public health. A colonizer of the human respiratory system, Stenotrophomonas maltophilia is frequently associated with hospital-acquired infections in individuals with cystic fibrosis, cancer, and other chronic illnesses. The importance of this study is underscored by its capacity to meet the critical demand for effective preventive strategies against this pathogen, particularly among susceptible groups of cystic fibrosis and those undergoing cancer treatment. In this study, we engineered a multi-epitope vaccine targeting S. maltophilia through genomic analysis, reverse vaccination strategies, and immunoinformatic techniques by examining a total of 81 complete genomes of S. maltophilia strains. Our investigation revealed 1945 core protein-coding genes alongside their corresponding proteomic sequences, with 191 of these genes predicted to exhibit virulence characteristics. Out of the filtered proteins, three best antigenic proteins were selected for epitope prediction while seven epitopes each from CTL, HTL, and B cell were chosen for vaccine development. The vaccine was refined and validated, showing highly antigenic and desirable physicochemical features. Molecular docking assessments revealed stable binding with TLR-4. Molecular dynamic simulation demonstrated stable dynamics with minor alterations. The originality of this investigation is rooted in the thorough techniques aimed at designing a vaccine that directly targets S. maltophilia, a microorganism of considerable clinical relevance that currently lacks an available vaccine. This study not only responds to a pressing public health crisis but also lays the groundwork for subsequent research endeavors focused on the prevention of S. maltophilia outbreaks. Further evidence from studies in mice models is needed to confirm immune protection against S. maltophilia.

Antiviral Potential of Chiococca alba (L.) Hitchc. Plant Extracts Against Chikungunya and Mayaro Viruses

Chikungunya and Mayaro fevers are viral infectious diseases characterized by fever and arthralgia, for which there are currently no effective vaccines or treatments. The urgent need for novel antiviral agents against Chikungunya virus (CHIKV) and Mayaro virus (MAYV) has led to interest in plant-based compounds that can disrupt the viral replication cycle. Chiococca alba (L.) Hitchc., a Neotropical plant traditionally used by Yucatec Maya healers as an antipyretic and antirheumatic, may hold potential as a source of antiviral agents. This study aimed to evaluate the antiviral potential of C. alba methanolic extracts (CAH21 and CAH24) against CHIKV and MAYV through preliminary in vitro and in silico analyses. The cytotoxicity of two methanolic extracts from C. alba roots was assessed in Vero cells using the neutral red assay, and their viral activity was determined via plaque assay post-treatment. Given the observed antiviral effects, we used computational predictions to explore interactions between the multifunctional nsP2 proteases and secondary metabolites identified in C. alba extracts. The metabolites were identified using high-performance liquid chromatography (HPLC) and gas chromatography-mass spectrometry (GC-MS). Phytochemical analysis revealed the presence of flavonoids, coumarins, and phenolic acids in the C. alba extracts. In vitro assays demonstrated that both extracts inhibited over 70% of activity against CHIKV and MAYV at a concentration of 60 µg/mL. In silico predictions suggested that the flavonoids naringin and vitexin had the highest affinity for the nsP2 proteases of CHIKV and MAYV, indicating their potential as viral inhibitors. Our findings revealed that C. alba extract represents a promising source of novel antiviral compounds.

Effect of cysteine mutation at Ca2+ coordinating residues to the autolysis, folding and hydrophobicity of full length and mature Rand protease: molecular dynamics simulation and essential dynamics

Rand protease is a serine protease that shared common characteristics with members of the MEROPS S8 subtilisin family. It is thermostable, highly stable in organic solvent and broad in specificity. Many structures of homologous protein solved by X-ray crystallography and NMR have been deposited to Protein Data Bank (PDB) which allowed this study to rely on structure prediction by deep learning to build three-dimensional (3D) structure of full length and mature Rand protease (flRP and mRP). In silico cysteine mutation to 7 predicted high affinity Ca2+ coordinating residues were introduced, and the mutants were subjected to molecular dynamics simulation to study its effect on flRP and mRP. MD simulation showed a marked increase in flexibility of the pro-peptide segment indicating the impact of single cysteine substitution at high affinity Ca2+ coordinating residues to autolysis of flRP. MD simulation for mRP reaffirmed the role of Ca2+ coordinating sites in providing stability to Rand protease. In addition, these residues also affect the autolysis, folding and hydrophobicity of RP. Essential dynamics observed large contribution of the first few eigenvectors of flRP, mRP and their high affinity Ca2+ coordinating residues mutants to the TMSF values which indicates that these values account for a large portion of the overall atomic fluctuations. These results have given a more comprehensive understanding on the role of cysteine substituted Ca2+ coordinating surface loop to the structure of flRP and mRP which are important in contributing to the structural stability of subtilisin.Communicated by Ramaswamy H. Sarma.

Bimodal substrate binding in the active site of the glycosidase BcX

Bacillus circulans xylanase (BcX) from the glycoside hydrolase family 11 degrades xylan through a retaining, double-displacement mechanism. The enzyme is thought to hydrolyze glycosidic bonds in a processive manner and has a large, active site cleft, with six subsites allowing the binding of six xylose units. Such an active site architecture suggests that oligomeric xylose substrates can bind in multiple ways. In the crystal structure of the catalytically inactive variant BcX E78Q, the substrate xylotriose is observed in the active site, as well as bound to the known secondary binding site and a third site on the protein surface. Nuclear magnetic resonance (NMR) titrations with xylose oligomers of different lengths yield nonlinear chemical shift trajectories for active site nuclei resonances, indicative of multiple binding orientations for these substrates for which binding and dissociation are in fast exchange on the NMR timescale, exchanging on the micro- to millisecond timescale. Active site binding can be modeled with a 2 : 1 model with dissociation constants in the low and high millimolar range. Extensive mutagenesis of active site residues indicates that tight binding occurs in the glycon binding site and is stabilized by Trp9 and the thumb region. Mutations F125A and W71A lead to large structural rearrangements. Binding at the glycon site is sensed throughout the active site, whereas the weak binding mostly affects the aglycon site. The interactions with the two active site locations are largely independent of each other and of binding at the secondary binding site.

Physiological modeling of the metaverse of the Mycobacterium tuberculosis β-CA inhibition mechanism

Tuberculosis (TB) is an infectious disease that primarily affects the lungs of humans and accounts for Mycobacterium tuberculosis (Mtb) bacteria as the etiologic agent. In this study, we introduce a computational framework designed to identify the important chemical features crucial for the effective inhibition of Mtb β-CAs. Through applying a mechanistic model, we elucidated the essential features pivotal for robust inhibition. Using this model, we engineered molecules that exhibit potent inhibitory activity and introduce relevant novel chemistry. The designed molecules were prioritized for synthesis based on their predicted pKi values via the QSAR (Quantitative Structure-Activity Relationship) model. All the rationally designed and synthesized compounds were evaluated in vitro against different carbonic anhydrase isoforms expressed from the pathogen Mtb; moreover, the off-target and widely human-expressed CA I and II were also evaluated. Among the reported derivatives, 2, 4, and 5 demonstrated the most valuable in vitro activity, resulting in promising candidates for the treatment of TB infection. All the synthesized molecules exhibited favorable pharmacokinetic and toxicological profiles based on in silico predictions. Docking analysis confirmed that the zinc-binding groups bind effectively into the catalytic triad of the Mtb β-Cas, supporting the in vitro outcomes with these binding interactions. Furthermore, molecules with good prediction accuracies according to previously established mechanistic and QSAR models were utilized to delve deeper into the realm of systems biology to understand their mechanism in combating tuberculotic pathogenesis. The results pointed to the key involvement of the compounds in modulating immune responses via NF-κβ1, SRC kinase, and TNF-α to modulate granuloma formation and clearance via T cells. This dual action, in which the pathogen's enzyme is inhibited while modulating the human immune machinery, represents a paradigm shift toward more effective and comprehensive treatment approaches for combating tuberculosis.

Exploration of dynamic interaction between β-lactoglobulin and casein micelles during UHT milk process

The protein aggregation induced by UHT treatment shortens the shelf life of UHT milk. However, the mechanism of β-Lg induced casein micelle aggregation remains unclear. Herein, the dynamic interaction between β-Lg and casein micelles during UHT processing was investigated by experimental techniques and molecular dynamics simulations. Results showed that β-Lg decreased the stability of casein micelles, increased their size and zeta potential. Raman and FTIR spectra analysis suggested that hydrogen and disulfide bonds facilitated their interaction. Cryo-TEM showed that the formation of the casein micelle/β-Lg complex involved rigid binding, flexible linking, and severe cross-linking aggregation during UHT processing. SAXS and MST demonstrated β-Lg bound to κ-casein on micelle surfaces with a dissociation constant (Kd) of 3.84 ± 1.14 μm. Molecular docking and dynamic simulations identified the interacting amino acid residues and clarified that electrostatic and van der Waals forces drove the interaction. UHT treatment increased hydrogen bonds and decreased total binding energy. The non-covalent binding promoted the formation of disulfide bonds between β-Lg and casein micelles under heat treatment. Ultimately, it was concluded that non-covalent interaction and disulfide bonding resulted in casein micelle/β-Lg aggregates. These findings provided scientific insights into protein aggregation in UHT milk.

Co-localization and co-expression of Olfml3 with Iba1 in brain of mice

Olfml3 is a microglia-specific protein whose role in neuroinflammation is elusive. In silico analysis was conducted to characterize the Olfml3 protein, followed by molecular docking and MD simulation to check possible interaction with Iba1. Further, expression and co-localization analysis was performed in the LPS-induced neuroinflammatory mice brains. Results suggest that Olfml3 physically interacts with Iba1. Olfml3 and Iba1 expression increases during neuroinflammation in mice brains. Olfml3 was observed to co-localize with Iba1, and the number of Olfml3 and Iba1 dual-positive cells increased in the brain of the neuroinflammatory mice model. Thus, Olfml3 could potentially participate in microglia functions by interacting with Iba1.

In silico evidence of bitopertin's broad interactions within the SLC6 transporter family

The Glycine Transporter Type 1 (GlyT1) significantly impacts central nervous system functions, influencing glycinergic and glutamatergic neurotransmission. Bitopertin, the first GlyT1 inhibitor in clinical trials, was developed for schizophrenia treatment but showed limited efficacy. Despite this, bitopertin's repositioning could advance treating various pathologies. This study aims to understand bitopertin's mechanism of action using computational methods, exploring off-target effects, and providing a comprehensive pharmacological profile. Similarity Ensemble Approach (SEA) and SwissTargetPrediction initially predicted targets, followed by molecular modeling on SWISS-MODEL and GalaxyWeb servers. Binding sites were identified using PrankWeb, and molecular docking was performed with DockThor and GOLD software. Molecular dynamics analyses were conducted on the Visual Dynamics platform. Reverse screening on SEA and SwissTargetPrediction identified GlyT1 (SLC6A9), GlyT2 (SLC6A5), PROT (SLC6A7), and DAT (SLC6A3) as potential bitopertin targets. Homology modeling on SwissModel generated high-resolution models, optimized further on GalaxyWeb. PrankWeb identified similar binding sites in GlyT1, GlyT2, PROT, and DAT, indicating potential interaction. Docking studies suggested bitopertin's interaction with GlyT1 and proximity to GlyT2 and PROT. Molecular dynamics confirmed docking results, highlighting bitopertin's target stability beyond GlyT1. The study concludes that bitopertin potentially interacts with multiple SLC6 family targets, indicating a broader pharmacological property.

In silico and experimental characterization of a new polyextremophilic subtilisin-like protease from Microbacterium metallidurans and its application as a laundry detergent additive

Considering the current growing interest in new and improved enzymes for use in a variety of applications, the present study aimed to characterize a novel detergent-stable serine alkaline protease from the extremophilic actinobacterium Microbacterium metallidurans TL13 (MmSP) using a combined in silico and experimental approach. The MmSP showed a close phylogenetic relationship with high molecular weight S8 peptidases of Microbacterium species. Moreover, its physical and chemical parameters computed using Expasy's ProtParam tool revealed that MmSP is hydrophilic, halophilic and thermo-alkali stable. 3D structure modelling and functional prediction of TL13 serine protease resulted in the detection of five characteristic domains: [catalytic subtilase domain, fibronectin (Fn) type-III domain, peptidase inhibitor I9, protease-associated (PA) domain and bacterial Ig-like domain (group 3)], as well as the three amino acid residues [aspartate (D182), histidine (H272) and serine (S604)] in the catalytic subtilase domain. The extremophilic strain TL13 was tested for protease production using agricultural wastes/by-products as carbon substrates. Maximum enzyme activity (390 U/gds) was obtained at 8th day fermentation on potato peel medium. Extracellular extract was concentrated and partially purified using ammonium sulfate precipitation methodology (1.58 folds purification fold). The optimal pH, temperature and salinity of MmSP were 9, 60 °C and 1 M NaCl, respectively. The MmSP protease showed broad pH stability, thermal stability, salt tolerance and detergent compatibility. In order to achieve the maximum stain removal efficacy by the TL 13 serine protease, the operation conditions were optimized using a Box-Behnken Design (BBD) with four variables, namely, time (15-75 min), temperature (30-60 °C), MmSP enzyme concentration (5-10 U/mL) and pH (7-11). The maximum stain removal yield (95 ± 4%) obtained under the optimal enzymatic operation conditions (treatment with 7.5 U/mL of MmSP during 30 min at 32 °C and pH9) was in good agreement with the value predicted by the regression model (98 ± %), which prove the validity of the fitted model. In conclusion, MmSP appears to be a good candidate for industrial applications, particularly in laundry detergent formulations, due to its high hydrophilicity, alkali-halo-stability, detergent compatibility and stain removal efficiency.

Transport of herbicides by PIN-FORMED auxin transporters

Auxins are a group of phytohormones that control plant growth and development 1. Their crucial role in plant physiology has inspired development of potent synthetic auxins that can be used as herbicides 2. Phenoxyacetic acid derivatives are a widely used group of auxin herbicides in agriculture and research. Despite their prevalence, the identity of the transporters required for distribution of these herbicides in plants is both poorly understood and the subject of controversial debate 3,4. Here we show that PIN-FORMED auxin transporters transport a range of phenoxyacetic acid herbicides across the membrane and we characterize the molecular determinants of this process using a variety of different substrates as well as protein mutagenesis to control substrate specificity. Finally, we present Cryo-EM structures of Arabidopsis thaliana PIN8 with 2,4-dichlorophenoxyacetic acid (2,4-D) or 4-chlorophenoxyacetic acid (4-CPA) bound. These structures represent five key states from the transport cycle, allowing us to describe conformational changes associated with substrate binding and transport across the membrane. Overall, our results reveal that phenoxyacetic acid herbicides use the same export machinery as endogenous auxins and exemplify how transporter binding sites undergo transformations that dictate substrate specificity. These results enable development of novel synthetic auxins and for guiding precision breeding of herbicide resistant crop plants.

Geraniol Potentiates the Effect of Fluconazole against Planktonic and Sessile Cells of Azole-Resistant Candida tropicalis: In Vitro and In Vivo Analyses

Candida tropicalis is regarded as an opportunistic pathogen, causing diseases ranging from superficial infections to life-threatening disseminated infections. The ability of this yeast to form biofilms and develop resistance to antifungals represents a significant therapeutic challenge. Herein, the effect of geraniol (GER), alone and combined with fluconazole (FLZ), was evaluated in the planktonic and sessile cells of azole-resistant C. tropicalis. GER showed a time-dependent fungicidal effect on the planktonic cells, impairing the cell membrane integrity. Additionally, GER inhibited the rhodamine 6G efflux, and the molecular docking analyzes supported the binding affinity of GER to the C. tropicalis Cdr1 protein. GER exhibited a synergism with FLZ against the planktonic and sessile cells, inhibiting the adhesion of the yeast cells and the viability of the 48-h biofilms formed on abiotic surfaces. C. tropicalis biofilms treated with GER, alone or combined with FLZ, displayed morphological and ultrastructural alterations, including a decrease in the stacking layers and the presence of wilted cells. Moreover, neither GER alone nor combined with FLZ caused toxicity, and both treatments prolonged the survival of the Galleria mellonella larvae infected with azole-resistant C. tropicalis. These findings indicate that the combination of GER and FLZ may be a promising strategy to control azole-resistant C. tropicalis infections.

The Mechanism of Action of L-Tyrosine Derivatives against Chikungunya Virus Infection In Vitro Depends on Structural Changes

Although the disease caused by chikungunya virus (CHIKV) is of great interest to public health organizations around the world, there are still no authorized antivirals for its treatment. Previously, dihalogenated anti-CHIKV compounds derived from L-tyrosine (dH-Y) were identified as being effective against in vitro infection by this virus, so the objective of this study was to determine the mechanisms of its antiviral action. Six dH-Y compounds (C1 to C6) dihalogenated with bromine or chlorine and modified in their amino groups were evaluated by different in vitro antiviral strategies and in silico tools. When the cells were exposed before infection, all compounds decreased the expression of viral proteins; only C4, C5 and C6 inhibited the genome; and C1, C2 and C3 inhibited infectious viral particles (IVPs). Furthermore, C1 and C3 reduce adhesion, while C2 and C3 reduce internalization, which could be related to the in silico interaction with the fusion peptide of the E1 viral protein. Only C3, C4, C5 and C6 inhibited IVPs when the cells were exposed after infection, and their effect occurred in late stages after viral translation and replication, such as assembly, and not during budding. In summary, the structural changes of these compounds determine their mechanism of action. Additionally, C3 was the only compound that inhibited CHIKV infection at different stages of the replicative cycle, making it a compound of interest for conversion as a potential drug.

AQuaRef: Machine learning accelerated quantum refinement of protein structures

Cryo-EM and X-ray crystallography provide crucial experimental data for obtaining atomic-detail models of biomacromolecules. Refining these models relies on library-based stereochemical restraints, which, in addition to being limited to known chemical entities, do not include meaningful noncovalent interactions relying solely on nonbonded repulsions. Quantum mechanical (QM) calculations could alleviate these issues but are too expensive for large molecules. We present a novel AI-enabled Quantum Refinement (AQuaRef) based on AIMNet2 neural network potential mimicking QM at substantially lower computational costs. By refining 41 cryo-EM and 30 X-ray structures, we show that this approach yields atomic models with superior geometric quality compared to standard techniques, while maintaining an equal or better fit to experimental data.

Tubulin code eraser CCP5 binds branch glutamates by substrate deformation

Microtubule function is modulated by the tubulin code, diverse posttranslational modifications that are altered dynamically by writer and eraser enzymes1. Glutamylation-the addition of branched (isopeptide-linked) glutamate chains-is the most evolutionarily widespread tubulin modification2. It is introduced by tubulin tyrosine ligase-like enzymes and erased by carboxypeptidases of the cytosolic carboxypeptidase (CCP) family1. Glutamylation homeostasis, achieved through the balance of writers and erasers, is critical for normal cell function3-9, and mutations in CCPs lead to human disease10-13. Here we report cryo-electron microscopy structures of the glutamylation eraser CCP5 in complex with the microtubule, and X-ray structures in complex with transition-state analogues. Combined with NMR analysis, these analyses show that CCP5 deforms the tubulin main chain into a unique turn that enables lock-and-key recognition of the branch glutamate in a cationic pocket that is unique to CCP family proteins. CCP5 binding of the sequences flanking the branch point primarily through peptide backbone atoms enables processing of diverse tubulin isotypes and non-tubulin substrates. Unexpectedly, CCP5 exhibits inefficient processing of an abundant β-tubulin isotype in the brain. This work provides an atomistic view into glutamate branch recognition and resolution, and sheds light on homeostasis of the tubulin glutamylation syntax.

A computational approach to identify phytochemicals as potential inhibitor of acetylcholinesterase: Molecular docking, ADME profiling and molecular dynamics simulations

Inhibition of acetylcholinesterase (AChE) is a crucial target in the treatment of Alzheimer's disease (AD). Common anti-acetylcholinesterase drugs such as Galantamine, Rivastigmine, Donepezil, and Tacrine have significant inhibition potential. Due to side effects and safety concerns, we aimed to investigate a wide range of phytochemicals and structural analogues of these compounds. Compounds similar to the established drugs, and phytochemicals were investigated as potential inhibitors for AChE in treating AD. A total of 2,270 compound libraries were generated for further analysis. Initial virtual screening was performed using Pyrx software, resulting in 638 molecules showing higher binding affinities compared to positive controls Tacrine (-9.0 kcal/mol), Donepezil (-7.3 kcal/mol), Galantamine (-8.3 kcal/mol), and Rivastigmine (-6.4 kcal/mol). Subsequently, ADME properties were assessed, including blood-brain barrier permeability and Lipinski's rule of five violations, leading to 88 compounds passing the ADME analysis. Among the rivastigmine analogous, [3-(1-methylpiperidin-2-yl)phenyl] N,N-diethylcarbamate showed interaction with Tyr123, Tyr336, Tyr340, Phe337, Trp285 residues of AChE. Tacrine similar compounds, such as 4-amino-2-styrylquinoline, exhibited bindings with Tyr123, Phe337, Tyr336, Trp285, Trp85, Gly119, and Gly120 residues. A phytocompound (bisdemethoxycurcumin) showed interaction with Trp285, Tyr340, Trp85, Tyr71, and His446 residues of AChE with favourable binding. These findings underscore the potential of these compounds as novel inhibitors of AChE, offering insights into alternative therapeutic avenues for AD. A 100ns simulation analysis confirmed the stability of protein-ligand complex based on the RMSD, RMSF, ligand properties, PCA, DCCM and MMGBS parameters. The investigation suggested 3 ligands as a potent inhibitor of AChE which are [3-(1-methylpiperidin-2-yl)phenyl] N,N-diethylcarbamate, 4-Amino-2-styrylquinoline and bisdemethoxycurcumin. Furthermore, investigation, including in-vitro and in-vivo studies, is needed to validate the efficacy, safety profiles, and therapeutic potential of these compounds for AD treatment.

Unlocking the potential of feruloyl esterase from Myceliophthora verrucosa: a key player in efficient conversion of biorefinery-relevant pretreated rice straw

The lignocellulolytic accessory enzyme, Feruloyl esterase C (FE_5DR), encoded in the genome of thermotolerant Myceliophthora verrucosa was successfully cloned and heterologously expressed in Pichia pastoris. The expressed FE_5DR was purified using UNOsphere™ Q anion exchange chromatography column, exhibiting a homogeneous band of ~ 39 kDa. Its optimum temperature was determined to be 60 °C, with an optimal pH of 6.0. Additionally, the enzyme activity of FE_5DR was significantly enhanced by preincubation in a buffer containing Mg2+, Cu2+ and Ca2 metal ions. Enzyme kinetic parameters, computed from double reciprocal Lineweaver-Burk plots, yielded observed Vmax and Km values of 0.758 U/mg and 0.439 mM, respectively. Furthermore, the potential of custom-made cocktails comprising FE_5DR and benchmark cellulase derived from the developed mutant strain of Aspergillus allahabadii MAN 40, as well as the biorefinery-relevant lignocellulolytic enzyme Cellic CTec 3, resulted in improved saccharification of unwashed acid pretreated (UWAP) rice straw slurry and mild alkali deacetylated (MAD) rice straw when compared to benchmark MAN 40 and Cellic CTec 3.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04013-7.

Cordycepin Triphosphate as a Potential Modulator of Cellular Plasticity in Cancer via cAMP-Dependent Pathways: An In Silico Approach

Cordycepin, or 3'-deoxyadenosine, is an adenosine analog with a broad spectrum of biological activity. The key structural difference between cordycepin and adenosine lies in the absence of a hydroxyl group at the 3' position of the ribose ring. Upon administration, cordycepin can undergo an enzymatic transformation in specific tissues, forming cordycepin triphosphate. In this study, we conducted a comprehensive analysis of the structural features of cordycepin and its derivatives, contrasting them with endogenous purine-based metabolites using chemoinformatics and bioinformatics tools in addition to molecular dynamics simulations. We tested the hypothesis that cordycepin triphosphate could bind to the active site of the adenylate cyclase enzyme. The outcomes of our molecular dynamics simulations revealed scores that are comparable to, and superior to, those of adenosine triphosphate (ATP), the endogenous ligand. This interaction could reduce the production of cyclic adenosine monophosphate (cAMP) by acting as a pseudo-ATP that lacks a hydroxyl group at the 3' position, essential to carry out nucleotide cyclization. We discuss the implications in the context of the plasticity of cancer and other cells within the tumor microenvironment, such as cancer-associated fibroblast, endothelial, and immune cells. This interaction could awaken antitumor immunity by preventing phenotypic changes in the immune cells driven by sustained cAMP signaling. The last could be an unreported molecular mechanism that helps to explain more details about cordycepin's mechanism of action.

Conformational changes in the Niemann-Pick type C1 protein NCR1 drive sterol translocation

The membrane protein Niemann-Pick type C1 (NPC1, named NCR1 in yeast) is central to sterol homeostasis in eukaryotes. Saccharomyces cerevisiae NCR1 is localized to the vacuolar membrane, where it is suggested to carry sterols across the protective glycocalyx and deposit them into the vacuolar membrane. However, documentation of a vacuolar glycocalyx in fungi is lacking, and the mechanism for sterol translocation has remained unclear. Here, we provide evidence supporting the presence of a glycocalyx in isolated S. cerevisiae vacuoles and report four cryo-EM structures of NCR1 in two distinct conformations, named tense and relaxed. These two conformations illustrate the movement of sterols through a tunnel formed by the luminal domains, thus bypassing the barrier presented by the glycocalyx. Based on these structures and on comparison with other members of the Resistance-Nodulation-Division (RND) superfamily, we propose a transport model that links changes in the luminal domains with a cycle of protonation and deprotonation within the transmembrane region of the protein. Our model suggests that NPC proteins work by a generalized RND mechanism where the proton motive force drives conformational changes in the transmembrane domains that are allosterically coupled to luminal/extracellular domains to promote sterol transport.

Conserved proline residues prevent dimerization and aggregation in the β-lactamase BlaC

Evolution leads to conservation of amino acid residues in protein families. Conserved proline residues are usually considered to ensure the correct folding and to stabilize the three-dimensional structure. Surprisingly, proline residues that are highly conserved in class A β-lactamases were found to tolerate various substitutions without large losses in enzyme activity. We investigated the roles of three conserved prolines at positions 107, 226, and 258 in the β-lactamase BlaC from Mycobacterium tuberculosis and found that mutations can lead to dimerization of the enzyme and an overall less stable protein that is prone to aggregate over time. For the variant Pro107Thr, the crystal structure shows dimer formation resembling domain swapping. It is concluded that the proline substitutions loosen the structure, enhancing multimerization. Even though the enzyme does not lose its properties without the conserved proline residues, the prolines ensure the long-term structural integrity of the enzyme.

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

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

AutoEpiCollect, a Novel Machine Learning-Based GUI Software for Vaccine Design: Application to Pan-Cancer Vaccine Design Targeting PIK3CA Neoantigens

Previous epitope-based cancer vaccines have focused on analyzing a limited number of mutated epitopes and clinical variables preliminarily to experimental trials. As a result, relatively few positive clinical outcomes have been observed in epitope-based cancer vaccines. Further efforts are required to diversify the selection of mutated epitopes tailored to cancers with different genetic signatures. To address this, we developed the first version of AutoEpiCollect, a user-friendly GUI software, capable of generating safe and immunogenic epitopes from missense mutations in any oncogene of interest. This software incorporates a novel, machine learning-driven epitope ranking method, leveraging a probabilistic logistic regression model that is trained on experimental T-cell assay data. Users can freely download AutoEpiCollectGUI with its user guide for installing and running the software on GitHub. We used AutoEpiCollect to design a pan-cancer vaccine targeting missense mutations found in the proto-oncogene PIK3CA, which encodes the p110ɑ catalytic subunit of the PI3K kinase protein. We selected PIK3CA as our gene target due to its widespread prevalence as an oncokinase across various cancer types and its lack of presence as a gene target in clinical trials. After entering 49 distinct point mutations into AutoEpiCollect, we acquired 361 MHC Class I epitope/HLA pairs and 219 MHC Class II epitope/HLA pairs. From the 49 input point mutations, we identified MHC Class I epitopes targeting 34 of these mutations and MHC Class II epitopes targeting 11 mutations. Furthermore, to assess the potential impact of our pan-cancer vaccine, we employed PCOptim and PCOptim-CD to streamline our epitope list and attain optimized vaccine population coverage. We achieved a world population coverage of 98.09% for MHC Class I data and 81.81% for MHC Class II data. We used three of our predicted immunogenic epitopes to further construct 3D models of peptide-HLA and peptide-HLA-TCR complexes to analyze the epitope binding potential and TCR interactions. Future studies could aim to validate AutoEpiCollect's vaccine design in murine models affected by PIK3CA-mutated or other mutated tumor cells located in various tissue types. AutoEpiCollect streamlines the preclinical vaccine development process, saving time for thorough testing of vaccinations in experimental trials.

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.

Analysis of SIKVAV's receptor affinity, pharmacokinetics, and pharmacological characteristics: a matrikine with potent biological function

Matrikines are biologically active peptides generated from fragments fragmentation of extracellular matrix components (ECM) that are functionally distinct from the original full-length molecule. The active matricryptic sites can be unmasked by ECM components enzymatic degradation or multimerization, heterotypic binding, adsorption to other molecules, cell-mediated mechanical forces, exposure to reactive oxygen species, ECM denaturation, and others. Laminin α1-derived peptide (SIKVAV) is a bioactive peptide derived from laminin-111 that participates in tumor development, cell proliferation, angiogenesis in various cell types. SIKVAV has also a potential pharmaceutical activity that may be used for tissue regeneration and bioengineering in Alzheimer's disease and muscular dystrophies. In this work, we made computational analyzes of SIKVAV regarding the ADMET panel, that stands for Administration, Distribution, Metabolism, Excretion, and Toxicity. Docking analyzes using the α3β1 and α6β1 integrin receptors were performed to fill in the gaps in the SIKVAV's signaling pathway and coupling tests showed that SIKVAV can interact with both receptors. Moreover, there is no indication of cytotoxicity, mutagenic or carcinogenic activity, skin or oral sensitivity. Our analysis suggests that SIKVAV has a high probability of interacting with peroxisome proliferator-activated receptor-gamma (NR-PPAR-γ), which has anti-inflammatory activity. The results of bioinformatics can help understand the participation of SIKVAV in homeostasis and influence the understanding of how this peptide can act as a biological asset in the control of dystrophies, neurodegenerative diseases, and tissue engineering.

Combined in vitro/in silico approaches, molecular dynamics simulations and safety assessment of the multifunctional properties of thymol and carvacrol: A comparative insight

Bioactive compounds derived from medicinal plants have acquired immense attentiveness in drug discovery and development. The present study investigated in vitro and predicted in silico the antibacterial, antifungal, and antiviral properties of thymol and carvacrol, and assessed their safety. The performed microbiological assays against Pseudomonas aeruginosa, Escherichia coli, Salmonella enterica Typhimurium revealed that the minimal inhibitory concentration values ranged from (0.078 to 0.312 mg/mL) and the minimal fungicidal concentration against Candida albicans was 0.625 mg/mL. Molecular docking simulations, stipulated that these compounds could inhibit bacterial replication and transcription functions by targeting DNA and RNA polymerases receptors with docking scores varying between (-5.1 to -6.9 kcal/mol). Studied hydroxylated monoterpenes could hinder C. albicans growth by impeding lanosterol 14α-demethylase enzyme and showed a (ΔG=-6.2 and -6.3 kcal/mol). Computational studies revealed that thymol and carvacrol could target the SARS-Cov-2 spike protein of the Omicron variant RBD domain. Molecular dynamics simulations disclosed that these compounds have a stable dynamic behavior over 100 ns as compared to remdesivir. Chemo-computational toxicity prediction using Protox II webserver indicated that thymol and carvacrol could be safely and effectively used as drug candidates to tackle bacterial, fungal, and viral infections as compared to chemical medication.

Immunoinformatics Prediction and Protective Efficacy of Vaccine Candidate PiuA-PlyD4 Against Streptococcus Pneumoniae

Purpose

This study was designed to evaluate the immune protective efficacy of the novel Streptococcus pneumoniae (S. pneumoniae) protein vaccine PiuA-PlyD4 through immunoinformatics prediction and in vitro and in vivo experiments.

Methods

In this study, we conducted immunoinformatics prediction and protection analysis on the fusion protein PiuA-PlyD4. The epitope composition of the vaccine was analyzed based on the prediction of B-cell and helper T-cell epitopes. Meanwhile, the molecular docking of PiuA and TLR2/4 was simulated. After immunizing C57BL/6 mice with the prepared vaccine, the biological safety, immunogenicity and conservation were evaluated. By constructing different infection models and from the aspects of adhesion inhibition and cytokines, the protective effect of the fusion protein vaccine PiuA-PlyD4 on S. pneumoniae infection was explored.

Results

PiuA-PlyD4 has abundant B-cell and helper T-cell epitopes and shows a high antigenicity score and structural stability. Molecular docking analysis suggested the potential interaction between PiuA and TLR2/4. The specific antibody titer of fusion protein antiserum was as high as (7.81±2.32) ×105. The protective effect of the immunized mice on nasal and lung colonization was significantly better than that of the control group, and the survival rate against S. pneumoniae infection of serotype 3 reached 50%. Cytokine detection showed that the humoral immune response, Th1, Th2 and Th17 cellular immune pathways were all involved in the process.

Conclusion

The study indicates that PiuA-PlyD4, whether the results are predicted by immunoinformatics or experimentally validated in vivo and in vitro, has good immunogenicity and immunoreactivity and can provide effective protection against S. pneumoniae infection. Therefore, it can be considered a promising prophylactic vaccine candidate for S. pneumoniae.

Quality Assessment of Selected Protein Structures Derived from Homology Modeling and AlphaFold

With technology advancing, many prediction algorithms have been developed to facilitate the modeling of inherently dynamic and flexible macromolecules such as proteins. Improvements in the prediction of protein structures have attracted a great deal of attention due to the advantages they offer, e.g., in drug design. While trusted experimental methods, such as X-ray crystallography, NMR spectroscopy, and electron microscopy, are preferred structure analysis techniques, in silico approaches are also being widely used. Two computational methods, which are on opposite ends of the spectrum with respect to their modus operandi, i.e., homology modeling and AlphaFold, have been established to provide high-quality structures. Here, a comparative study of the quality of structures either predicted by homology modeling or by AlphaFold is presented based on the characteristics determined by experimental studies using structure validation servers to fulfill the purpose. Although AlphaFold is able to predict high-quality structures, high-confidence parts are sometimes observed to be in disagreement with experimental data. On the other hand, while the structures obtained from homology modeling are successful in incorporating all aspects of the experimental structure used as a template, this method may struggle to accurately model a structure in the absence of a suitable template. In general, although both methods produce high-quality models, the criteria by which they are superior to each other are different and thus discussed in detail.

Unraveling Extremely Damaging IRAK4 Variants and Their Potential Implications for IRAK4 Inhibitor Efficacy

Interleukin-1-receptor-associated kinase 4 (IRAK4) possesses a crucial function in the toll-like receptor (TLR) signaling pathway, and the dysfunction of this molecule could lead to various infectious and immune-related diseases in addition to cancers. IRAK4 genetic variants have been linked to various types of diseases. Therefore, we conducted a comprehensive analysis to recognize the missense variants with the most damaging impacts on IRAK4 with the employment of diverse bioinformatics tools to study single-nucleotide polymorphisms' effects on function, stability, secondary structures, and 3D structure. The residues' location on the protein domain and their conservation status were investigated as well. Moreover, docking tools along with structural biology were engaged in analyzing the SNPs' effects on one of the developed IRAK4 inhibitors. By analyzing IRAK4 gene SNPs, the analysis distinguished ten variants as the most detrimental missense variants. All variants were situated in highly conserved positions on an important protein domain. L318S and L318F mutations were linked to changes in IRAK4 secondary structures. Eight SNPs were revealed to have a decreasing effect on the stability of IRAK4 via both I-Mutant 2.0 and Mu-Pro tools, while Mu-Pro tool identified a decreasing effect for the G198E SNP. In addition, detrimental effects on the 3D structure of IRAK4 were also discovered for the selected variants. Molecular modeling studies highlighted the detrimental impact of these identified SNP mutant residues on the druggability of the IRAK4 ATP-binding site towards the known target inhibitor, HG-12-6, as compared to the native protein. The loss of important ligand residue-wise contacts, altered protein global flexibility, increased steric clashes, and even electronic penalties at the ligand-binding site interfaces were all suggested to be associated with SNP models for hampering the HG-12-6 affinity towards IRAK4 target protein. This given model lays the foundation for the better prediction of various disorders relevant to IRAK4 malfunction and sheds light on the impact of deleterious IRAK4 variants on IRAK4 inhibitor efficacy.

Structural basis for kinase inhibition in the tripartite E. coli HipBST toxin-antitoxin system

Many bacteria encode multiple toxin-antitoxin (TA) systems targeting separate, but closely related, cellular functions. The toxin of the Escherichia coli hipBA system, HipA, is a kinase that inhibits translation via phosphorylation of glutamyl-tRNA synthetase. Enteropathogenic E. coli O127:H6 encodes the hipBA-like, tripartite TA system; hipBST, in which the HipT toxin specifically targets the tryptophanyl-tRNA synthetase, TrpS. Notably, in the tripartite system, the function as antitoxin has been taken over by the third protein, HipS, but the molecular details of how activity of HipT is inhibited remain poorly understood. Here, we show that HipBST is structurally different from E. coli HipBA and that the unique HipS protein, which is homologous to the N-terminal subdomain of HipA, inhibits the kinase through insertion of a conserved Trp residue into the active site. We also show how auto-phosphorylation at two conserved sites in the kinase toxin serve different roles and affect the ability of HipS to neutralize HipT. Finally, solution structural studies show how phosphorylation affects overall TA complex flexibility.

Antiemetic effects of sclareol, possibly through 5-HT3 and D2 receptor interaction pathways: In-vivo and in-silico studies

BACKGROUND

Emesis is a complex physiological phenomenon that serves as a defense against numerous toxins, stressful situations, adverse medication responses, chemotherapy, and movement. Nevertheless, preventing emesis during chemotherapy or other situations is a significant issue for researchers. Hence, the majority view contends that successfully combining therapy is the best course of action. In-vivo analysis offers a more comprehensive grasp of how compounds behave within a complex biological environment, whereas in-silico evaluation refers to the use of computational models to forecast biological interactions.

OBJECTIVES

The objectives of the present study were to evaluate the effects of Sclareol (SCL) on copper sulphate-induced emetic chicks and to investigate the combined effects of these compounds using a conventional co-treatment approach and in-silico study.

METHODS

SCL (5, 10, and 15 mg/kg) administered orally with or without pre-treatment with anti-emetic drugs (Ondansetron (ODN): 24 mg/kg, Domperidone (DOM): 80 mg/kg, Hyoscine butylbromide (HYS): 100 mg/kg, and Promethazine hydrochloride (PRO): 100 mg/kg) to illustrate the effects and the potential involvement with 5HT3, D2, M3/AChM, H1, or NK1 receptors by SCL. Furthermore, an in-silico analysis was conducted to forecast the role of these receptors in the emetic process.

RESULTS

The results suggest that SCL exerted a dose-dependent anti-emetic effect on the chicks. Pretreatment with SCL-10 significantly minimized the number of retches and lengthened the emesis tendency of the experimental animals. SCL-10 significantly increased the anti-emetic effects of ODN and DOM. However, compared to the ODN-treated group, (SCL-10 + ODN) group considerably (p < 0.0001) extended the latency duration (109.40 ± 1.03 s) and significantly (p < 0.01) decreased the number of retches (20.00 ± 0.70), indicating an anti-emetic effect on the test animals. In in-silico analysis, SCL exhibited promising binding affinities with suggesting receptors.

CONCLUSION

SCL-10 exerted an inhibitory-like effect on emetic chicks, probably through the interaction of the 5HT3 and D2 receptors. Further studies are highly appreciated to validate this study and determine the precise mechanism(s) behind the anti-emetic effects of SCL. We expect that SCL-10 may be utilized as an antiemetic treatment in a single dosage form or that it may function as a synergist with other traditional medicines.

An amino-domino model described by a cross-peptide-bond Ramachandran plot defines amino acid pairs as local structural units

Protein structure, both at the global and local level, dictates function. Proteins fold from chains of amino acids, forming secondary structures, α-helices and β-strands, that, at least for globular proteins, subsequently fold into a three-dimensional structure. Here, we show that a Ramachandran-type plot focusing on the two dihedral angles separated by the peptide bond, and entirely contained within an amino acid pair, defines a local structural unit. We further demonstrate the usefulness of this cross-peptide-bond Ramachandran plot by showing that it captures β-turn conformations in coil regions, that traditional Ramachandran plot outliers fall into occupied regions of our plot, and that thermophilic proteins prefer specific amino acid pair conformations. Further, we demonstrate experimentally that the effect of a point mutation on backbone conformation and protein stability depends on the amino acid pair context, i.e., the identity of the adjacent amino acid, in a manner predictable by our method.

Molecular interaction analysis of anti-IL-8 scFv-10F8-6His against IL-8 monomer through molecular docking and molecular dynamic simulations

Elevated interleukin 8 (IL-8) expression has been linked to unfavorable outcomes in a range of inflammatory conditions, such as rheumatoid arthritis, psoriasis, and cancer. The human monoclonal antibody (HuMab) 10F8 and the hybridoma 35B11-B bind to an epitope on human IL-8, respectively. 10F8 inhibited interaction between IL-8 and neutrophils in eczema and pustulosis palmoplantaris patients while 35B11-B decreased size lesion in rat model. The binding interaction of monoclonal antibodies and IL-8, especially how complementarity-determining region (CDR) loops could bind the N-terminal of IL-8, has not been fully deliberated at molecular-level. Here, we used a combination of molecular docking, heated and long coarse-grained molecular dynamics simulations to identify key residues of established interaction. Based on heated MD simulation, docked pose of complexes generated by ClusPro showed good binding stability throughout of 70 ns simulation. Based on long molecular dynamic simulations, key residues for the binding were identified throughout of 1000 ns simulation. TYR-53, ASP-99, and ARG-100 of heavy chain CDR together with TYR-33 of light chain CDR are among the highest contributing energy residues within the binding interaction. Meanwhile, LYS11 and TYR13 of IL-8 are important for the determination of overall binding energy. Furthermore, the result of decomposition residues analysis is in good agreement with the interaction analysis data. Current study provides a list of important interacting residues and further scrutiny on these residues is essential for future development and design of a new and stable recombinant antibody against IL-8.Communicated by Ramaswamy H. Sarma.

Precision drug design against Acidovorax oryzae: leveraging bioinformatics to combat rice brown stripe disease

Bacterial brown stripe disease caused by Acidovorax oryzae is a major threat to crop yields, and the current reliance on pesticides for control is unsustainable due to environmental pollution and resistance. To address this, bacterial-based ligands have been explored as a potential treatment solution. In this study, we developed a protein-protein interaction (PPI) network for A. oryzae by utilizing shared differentially expressed genes (DEGs) and the STRING database. Using a maximal clique centrality (MCC) approach through CytoHubba and Network Analyzer, we identified hub genes within the PPI network. We then analyzed the genomic data of the top 10 proteins, and further narrowed them down to 2 proteins by utilizing betweenness, closeness, degree, and eigenvector studies. Finally, we used molecular docking to screen 100 compounds against the final two proteins (guaA and metG), and Enfumafungin was selected as a potential treatment for bacterial resistance caused by A. oryzae based on their binding affinity and interaction energy. Our approach demonstrates the potential of utilizing bioinformatics and molecular docking to identify novel drug candidates for precision treatment of bacterial brown stripe disease caused by A. oryzae, paving the way for more targeted and sustainable control strategies. The efficacy of Enfumafungin in inhibiting the growth of A. oryzae strain RS-1 was investigated through both computational and wet lab methods. The models of the protein were built using the Swiss model, and their accuracy was confirmed via a Ramachandran plot. Additionally, Enfumafungin demonstrated potent inhibitory action against the bacterial strain, with an MIC of 100 µg/mL, reducing OD600 values by up to 91%. The effectiveness of Enfumafungin was further evidenced through agar well diffusion assays, which exhibited the highest zone of inhibition at 1.42 cm when the concentration of Enfumafungin was at 100 µg/mL. Moreover, Enfumafungin was also able to effectively reduce the biofilm of A. oryzae RS-1 in a concentration-dependent manner. The swarming motility of A. oryzae RS-1 was also found to be significantly inhibited by Enfumafungin. Further validation through TEM observation revealed that bacterial cells exposed to Enfumafungin displayed mostly red fluorescence, indicating destruction of the bacterial cell membrane.

Asp179 in the class A β-lactamase from Mycobacterium tuberculosis is a conserved yet not essential residue due to epistasis

Conserved residues are often considered essential for function, and substitutions in such residues are expected to have a negative influence on the properties of a protein. However, mutations in a few highly conserved residues of the β-lactamase from Mycobacterium tuberculosis, BlaC, were shown to have no or only limited negative effect on the enzyme. One such mutant, D179N, even conveyed increased ceftazidime resistance upon bacterial cells, while displaying good activity against penicillins. The crystal structures of BlaC D179N in resting state and in complex with sulbactam reveal subtle structural changes in the Ω-loop as compared to the structure of wild-type BlaC. Introducing this mutation in four other β-lactamases, CTX-M-14, KPC-2, NMC-A and TEM-1, resulted in decreased antibiotic resistance for penicillins and meropenem. The results demonstrate that the Asp in position 179 is generally essential for class A β-lactamases but not for BlaC, which can be explained by the importance of the interaction with the side chain of Arg164 that is absent in BlaC. It is concluded that Asp179 though conserved is not essential in BlaC, as a consequence of epistasis.