PDB Entry Research Articles

Improving the diffraction quality of heat-shock protein 47 crystals.

Heat-shock protein 47 (HSP47) is a potential target for inhibitors that ameliorate fibrosis by reducing collagen assembly. In an effort to develop a structure-based drug-design system, it was not possible to replicate a previous literature result (PDB entry 4au4) for apo dog HSP47; instead, crystal forms were obtained in which pairs of dog HSP47 molecules interacted through a noncleavable C-terminal His-tag to build up tetramers, all of which had multiple molecules of HSP47 in the asymmetric unit and none of which diffracted as wellas the literature precedent. To overcome these difficulties, a two-pronged approach was followed: (i) the His-tag was moved from the C-terminus to the N-terminus and was made cleavable, and (ii) Adnectin (derived from the tenth domain of human fibronectin type III) crystallization chaperones were developed. Both approaches provided well diffracting crystals, but the latter approach yielded crystal forms with only one or two HSP47 complexes per asymmetric unit, which made model building less onerous.

Crystal structures of NAD(P)H nitroreductases from Klebsiella pneumoniae.

Klebsiella pneumoniae (Kp) is an infectious disease pathogen that poses a significant global health threat due to its potential to cause severe infections andits tendency to exhibit multidrug resistance. Understanding the enzymatic mechanisms of the oxygen-insensitive nitroreductases (Kp-NRs) from Kp is crucial for the development of effective nitrofuran drugs, such as nitrofurantoin, that can be activated as antibiotics. In this paper, three crystal structures of two Kp-NRs (PDB entries 7tmf/7tmg and 8dor) are presented, and an analysis of their crystal structures and their flavin mononucleotide (FMN)-binding mode is provided. The structures with PDB codes 7tmf (Kp-NR1a), 7tmg (Kp-NR1b) and 8dor (Kp-NR2) were determined at resolutions of 1.97, 1.90 and 1.35 Å, respectively. The Kp-NR1a and Kp-NR1b structures adopt an αβ fold, in which four-stranded antiparallel β-sheets are surrounded by five helices. With domain swapping, the β-sheet was expanded with a β-strand from the other molecule of the dimer. The difference between the structures lies in the loop spanning Leu173-Ala185: in Kp-NR1a the loop is disordered, whereas the loop adopts multiple conformations in Kp-NR1b. The FMN interactions within Kp-NR1/NR2 involve hydrogen-bond and π-stacking interactions. Kp-NR2 contains four-stranded antiparallel β-sheets surrounded by eight helices with two short helices and one β-sheet. Structural and sequence alignments show that Kp-NR1a/b and Kp-NR2 are homologs of the Escherichia coli oxygen-insensitive NRs YdjA and NfnB and of Enterobacter cloacae NR, respectively. By homology inference from E. coli, Kp-NR1a/b and Kp-NR2 may detoxify polynitroaromatic compounds and Kp-NR2 may activate nitrofuran drugs to cause bactericidal activity through a ping-pong bi-bi mechanism, respectively.

PPI3D: a web server for searching, analyzing and modeling protein-protein, protein-peptide and protein-nucleic acid interactions.

Structure-resolved protein interactions with other proteins, peptides and nucleic acids are key for understanding molecular mechanisms. The PPI3D web server enables researchers to query preprocessed and clustered structural data, analyze the results and make homology-based inferences for protein interactions. PPI3D offers three interaction exploration modes: (i) all interactions for proteins homologous to the query, (ii) interactions between two proteins or their homologsand (iii) interactions within a specific PDB entry. The server allows interactive analysis of the identified interactions in both summarized and detailed manner. This includes protein annotations, structures, the interface residues and the corresponding contact surface areas. In addition, users can make inferences about residues at the interaction interface for the query protein(s) from the sequence alignments and homology models. The weekly updated PPI3D database includes all the interaction interfaces and binding sites from PDB, clustered based on both protein sequence and structural similarity, yielding non-redundant datasets without loss of alternative interaction modes. Consequently, the PPI3D users avoid being flooded with redundant information, a typical situation for intensely studied proteins. Furthermore, PPI3D provides a possibility to download user-defined sets of interaction interfaces and analyze them locally. The PPI3D web server is available at https://bioinformatics.lt/ppi3d.

Computational Binding Study Hints at Ecdysone 20-Mono-Oxygenase as the Hitherto Unknown Target for Ring C-Seco Limonoid-Type Insecticides.

The insecticidal property of ring C-seco limonoids has been discovered empirically and the target protein identified, but, to date, the molecular mechanism of action has not been described at the atomic scale. We elucidate on computational grounds whether nine C-seco limonoids present sufficiently high affinity to bind specifically with the putative target enzyme of the insects (ecdysone 20-monooxygenase). To this end, 3D models of ligands and the receptor target were generated and their interaction energies estimated by docking simulations. As a proof of concept, the tetrahydro-isoquinolinyl propenamide derivative QHC is the reference ligand bound to aldosterone synthase in the complex with PDB entry 4ZGX. It served as the 3D template for target modeling via homology. QHC was successfully docked back to its crystal pose in a one-digit nanomolar range. The reported experimental binding affinities span over the nanomolar to lower micromolar range. All nine limonoids were found with strong affinities in the range of -9 < ΔG < -13 kcal/mol. The molt hormone ecdysone showed a comparable ΔG energy of -12 kcal/mol, whereas -11 kcal/mol was the back docking result for the liganded crystal 4ZGX. In conclusion, the nine C-seco limonoids were strong binders on theoretical grounds in an activity range between a ten-fold lower to a ten-fold higher concentration level than insecticide ecdysone with its known target receptor. The comparable or even stronger binding hints at ecdysone 20-monooxygenase as their target biomolecule. Our assumption, however, is in need of future experimental confirmation before conclusions with certainty can be drawn about the true molecular mechanism of action for the C-seco limonoids under scrutiny.

Binding interactions and in silico ADME prediction of isoconessimine derivatives as potent acetylcholinesterase inhibitors

In pursuit of new acetylcholinesterase (AChE) inhibitors for treating Alzheimer's disease (AD), a series of ten previously synthesized isoconessimine compounds (7a-7j) was in silico investigated for their binding interactions with AChE and pharmacokinetics based on absorption, distribution, metabolism, and excretion (ADME) properties using molecular docking, ONIOM (Our own N-layered Integrated molecular Orbital and molecular Mechanics) method and SwissADME tools. Docking experiments showed that all compounds bind within the active site gorge of AChE (PDB entry 1C2B), posing its aryloxy-substitutional ethyl group to catalytic site and conessine skeleton to peripheral anionic site. ONIOM interaction energy was used as an ONIOM score to improve docking score, and it ranked 7b as the most potent AChE inhibitor, in agreement with previous experiment. Residues, ASP74, TRP86, GLY122, GLU202, TRP286, GLU292, SER293, ILE294, TYR337, TYR341, and HIS447 were identified as important for the binding of the AChE-isoconessimine complex. The SwissADME investigation suggested that four compounds (7a, 7c, 7d and 7f) agree with the rules of drug-likeness. The steric and electronic effects on the aryloxy-substitutional ethyl group as important factors in the AChE inhibition were also discussed, which brings a better understanding of Alzheimer's disease drug development.

Molecular Docking Study of New Pyrrolyl Benzohydrazide Schiff Bases as M. tuberculosis InhA Enzyme Inhibitors

Background Tuberculosis TB is an ancient chronic infectious disease affecting primarily the lungs of humans and rapidly co-evolved with humans since several million years. It is mainly caused by a group of bacterial species termed Mycobacterium tuberculosis a principal strain of this infection. Each second millions are at risk of acquiring new infection of becoming sick prone to death and a conclusive understanding of pathogenesis is still deficient.Objective To design a series of new pyrrole benzohydrazide schiff bases as M. tuberculosis InhA inhibitors.Methods Docking software Sybyl-X version 2.0 was used for molecular modeling. To dock intended compounds Surflex-Dock procedure of sybyl was castoff. M. tuberculosis enoyl reductase InhA inhibitor target protein 4TZK complexed with ligand 1-cyclohexyl-N-3 5-dichlorophenyl-5-oxopyrrolidine3-carboxamide was taken from Protein Data Bank PDB entry code 4TZK httpwww.rcsb.orgpdb.Results Docking scores of the compounds ranged between 5.47 - 3.32 summarizing the interface force field between enzyme and the ligands.Conclusion Sixteen molecules bearing pyrrole hydrazone moiety were docked showing promising M. tuberculosis inhibition activity. In silico molecular docking analysis exhibited that peptide linkage -CO-NH- and imine -CN- of pyrrole hydrazones were important in binding with the targeted receptor.

Phytochemical investigation and anti-inflammatory potential of Atriplex leucoclada Boiss

BackgroundThe plant kingdom has long been considered a valuable source for therapeutic agents, however, some plant species still untapped and need to be phytochemically and biologically explored. Although several Atriplex species have been investigated in depth, A. leucoclada, a halophytic plant native to Saudi Arabian desert, remains to be explored for its phytochemical content and biological potentials. Herein, the current study investigated the metabolic content and the anti-inflammatory potential of A. leucoclada.MethodsPowdered aerial parts of the plant were defatted with n-hexane then the defatted powder was extracted with 80% methanol. n-Hexane extract (ATH) was analyzed using GC–MS, while the defatted extract (ATD) was subjected to different chromatographic methods to isolate the major phytoconstituents. The structures of the purified compounds were elucidated using different spectroscopic methods including advanced NMR techniques. Anti-inflammatory activity of both extracts against COX-1 and COX-2 enzymes were examined in vitro. Molecular docking of the identified compounds into the active sites of COX-1 and COX-2 enzymes was conducted using pdb entries 6Y3C and 5IKV, respectively.ResultsPhytochemical investigation of ATD extract led to purification and identification of nine compounds. Interestingly, all the compounds, except for 20-hydroxy ecdysone (1), are reported for the first time from A. leucoclada, also luteolin (6) and pallidol (8) are isolated for the first time from genus Atriplex. Inhibitory activity of ATD and ATH extracts against COX-1 and COX-2 enzymes revealed concentration dependent activity of both fractions with IC50 41.22, 14.40 μg/ml for ATD and 16.74 and 5.96 μg/ml for ATH against COX-1 and COX-2, respectively. Both extracts displayed selectivity indices of 2.86 and 2.80, respectively as compared to 2.56 for Ibuprofen indicating a promising selectivity towards COX-2. Molecular docking study supported in vitro testing results, where purified metabolites showed binding affinity scores ranged from -9 to -6.4 and -8.5 to -6.6 kcal/mol for COX-1 and 2, respectively, in addition the binding energies of GC–MS detected compounds ranged from -8.9 to -5.5 and -8.3 to -5.1 kcal/mol for COX-1 and 2, respectively as compared to Ibuprofen (-6.9 and -7.5 kcal/mol, respectively), indicating high binding affinities of most of the compounds. Analysis of the binding orientations revealed variable binding patterns depending on the nature of the compounds. Our study suggested A. leucoclada as a generous source for anti-inflammatory agents.

In silico development of potential InhA inhibitors through 3D-QSAR analysis, virtual screening and molecular dynamics

Tuberculosis is one of the most ancient infectious diseases known to mankind predating upper Paleolithic era. In the current scenario, treatment of drug resistance tuberculosis is the major challenge as the treatment options are limited, less efficient and more toxic. In our study we have developed an atom based 3D QSAR model, statistically validated sound with R2 > 0.90 and Q2 > 0.72 using reported direct inhibitors of InhA (2018–2022), validated by enzyme inhibition assay. The model was used to screen a library of 3958 molecules taken from Binding DB and candidates molecules with promising predicted activity value (pIC50) > 5) were selected for further analyzed screening by using molecular docking, ADME profiling and molecular dynamic simulations. The lead molecule, ZINC11536150 exhibited good docking score (glideXP = −11.634 kcal/mol) compared to standard triclosan (glideXP = −7.129 kcal/mol kcal/mol) and through molecular dynamics study it was observed that the 2nv6-complex of ZINC11536150 with Mycobacterium tuberculosis InhA (PDB entry: 2NV6) complex remained stable throughout the entire simulation time of 100 ns.Communicated by Ramaswamy H. Sarma

Mitochondrial Acyl Carrier Protein of Leishmania major Displays Features Distinct from the Canonical Type II ACP.

Prokaryotes synthesize fatty acids using a type II synthesis pathway (FAS). In this process, the central player, i.e., the acyl carrier protein (ACP), sequesters the growing acyl chain in its internal hydrophobic cavity. As the acyl chain length increases, the cavity expands in size, which is reflected in the NMR chemical shift perturbations and crystal structures of the acyl-ACP intermediates. A few eukaryotic organelles, such as plastids and mitochondria, also harbor type II fatty acid synthesis machinery. Plastid FAS from spinach and Plasmodium falciparum has been characterized at the molecular level, but the mitochondrial pathway remains unexplored. Here, we report NMR studies of the mitochondrial acyl-acyl carrier protein intermediates of Leishmania major (acyl-LmACP). Our studies show that LmACP experiences remarkably small conformational changes upon acylation, with perturbations confined to helices II and III only. CastP determined that the cavity size of apo-LmACP (PDB entry 5ZWT) is less than that of Escherichia coli ACP (PDB 1T8K). Thus, the small chemical shift perturbations observed in the LmACP intermediates, coupled with CastP results, suggest an unusually small cavity when fully expanded. The faster rate of C8-LmACP chain hydrolysis compared to E. coli ACP (EcACP) also supports these convictions. Structure comparison of LmACP with other type II ACP disclosed unique differences in the helix I and loop I conformations, as well as several residues present there. Numerous hydrophobic residues in helix I and loop I (conserved in all mitochondrial ACPs) are substituted with hydrophilic residues in the bacterial/plastid type II ACP. For instance, Phe and leucine at positions 14 and 34 in LmACP are substituted with a hydrophilic residue and Ala in bacterial/plastid type II ACP. Mutation of Leu 34 to Ala (corresponding residue in EcACP) resulted in a complete loss of structure, underscoring its importance in maintaining the ACP fold. Thus, our NMR studies, combined with insights from the crystal structure, highlight several unique features of LmACP, distinct from the prokaryote and plastid type II ACP. Given the high sequence identity, the features might be conserved in all mitochondrial ACPs.

Biochemical properties of Scedosporium aurantiacum extracellular elastase-like protease

Extracellular elastase-like protease is one of the key virulence proteases of Scedosporium aurantiacum. To date, little is known about this enzyme in terms of genetic information, structure, properties and virulence mechanism due to the difficulties in purification caused by its low secretion amount, high specific activity, uncompleted genome sequencing and annotation. This work investigated the gene, structure and enzymatic properties of this enzyme. The S. aurantiacum elastase-like protease from the fungal culture supernatant was analyzed through tandem mass spectrometry (MS/MS) approach, illustrating its primary structure. Bioinformatics tools were employed to predict the conserved domain and tertiary structure, the enzymatic properties were also studied. It turned out that S. aurantiacum extracellular elastase-like protease demonstrated well hydrolysis towards elastin and bovine achilles tendon collagen, with Vmax of 18.14 μg/s and 17.57 μg/s respectively, better than fish scale gelatin, with the lowest hydrolysis effect on casein. Its activity towards elastin was lower than that of the elastase from porcine pancreas, with values of Kcat/Km of 3.541 (μg/s) and 4.091 (μg/s), respectively. It was an alkaline protease, with optimal pH 8.2 and temperature 37 oC. Zn2+ promoted the enzymatic activity while Ca2+, Mg2+, Na+, elastatinal and PMSF inhibited its activity. Its sequence was similar to Paecilomyces lilacinus secreted serine protease (PDB Entry: c3f7oB_) with multiple conserved fractions each containing more than 7 amino acids, thus suitable for design of PCR primer. This study increased our knowledge on S. aurantiacum extracellular elastase-like protease in terms of structure and enzymatic properties, and may facilitate later studies on protein expression and virulence mechanism.

Virtual Design of Novel Coumarinyl-Substituted Sulfonamide Inhibitors of Carbonic Anhydrase II as Potential Drugs against Glaucoma

Background: We have carried out virtual design of coumarinyl-substituted sulfonamides (CSAM) analogs as inhibitors of human carbonic anhydrase II (hCA II) endowed with favorable predicted pharmacokinetic profiles and potential therapeutic effects against glaucoma.&#x0D; Methods: modifying in situ the x-ray structure of 2-(7-methoxy-2-oxo-2H-chromen-4-yl)-N-(4-sulfamoyl-phenyl)-acetamide (CSAM0)-hCA II complex (PDB entry 3ML2), permitted 3D models of hCA II-CSAMx complexes preparation for a TS of 14 CSAMs the experimental activities of which are available in the literature (IC50exp). Active bound conformation of the CSAM1-14 assessment led to linear correlation between computed enthalpy of hCA II-CSAMx complexes formation in gas phase (ΔΔHMM) with the IC50exp. Moreover, considering the solvation and ligand loss of entropy upon binding led to a superior QSAR model where a better linear correlation is established between calculated GFE (ΔΔGcom) and IC50exp. The successive 3D-pharmacophore (PH4) built from CSAMs active conformations helped to virtually screen CSAM populating a sulfonamides scaffolds virtual combinatorial library (VCL), focused by Lipinski’s rule-of-five to reach novel CSAMs.&#x0D; Results: Enthalpy QSAR model: pIC50exp = -0.068×ΔΔHMM + 7.722, R2 = 0.82; GFE QSAR model: pIC50exp = - 0.061× ΔΔGcom + 7.647, R2 = 0.92 and PH4 model: pIC50exp = 1.095×pIC50pre – 0.680, R2 = 0.87. The VCL of more than 1,500,625 CSAMs was lowered to 865,670 drug likely compounds by the Lipinski’s rule (except the restriction Mw ≤ 500 g/mol). The three-point PH4-based screening identified 81 novel CSAMs with predicted IC50pre reaching 78-times better than CSAM1 potency (IC50exp = 23 nM). Computed pharmacokinetic profile of the new candidates showed enhanced cell membrane permeability and high human oral absorption compared to current anti-glaucoma agents.&#x0D; Conclusions: Combination of QSAR models of the CSAMs’ affinity to the hCA II, pharmacophore model, and ADME profile guided to the sulfonamides inhibitors identification and helped to in silico screen VCL sulfonamides scaffold bearing analogs and allowed emerging of novel more potent compounds with predicted IC50 and favorable pharmacokinetic profiles.

Correspondence on “Structural Insight into the Catalytic Mechanism of the Endoperoxide Synthase FtmOx1”

AbstractIn their recent Angewandte Chemie publication (doi: 10.1002/anie.202112063), Cen, Wang, Zhou et al. reported the crystal structure of a ternary complex of the non‐heme iron endoperoxidase FtmOx1 (PDB entry 7ETK). The biochemical data assessed in this study were from a retracted study (doi: 10.1038/nature15519) by Zhang, Liu, Zhang et al.; no additional biochemical data were included, yet there was no discussion on the source of the biochemical data in the report by Cen, Wang, Zhou et al. Based on this new crystal structure and subsequent QM/MM‐MD calculations, Cen, Wang, Zhou et al. concluded that their work provided evidence supporting the CarC‐like mechanistic model for FtmOx1 catalysis. However, the authors did not accurately describe either the CarC‐like model or the COX‐like model, and they did not address the differences between them. Further, and contrary to their interpretations in the manuscript, the authors’ data are consistent with the COX‐like model once the details of the CarC‐like and COX‐like models have been carefully analyzed.

Experimentally Observed Conformational Changes in Antibodies Due to Binding and Paratope-epitope Asymmetries

Understanding binding related changes in antibody conformations is important for epitope prediction and antibody refinement. The increase of available data in the PDB allowed a more detailed investigation of the conformational landscape for free and bound antibodies. A dataset containing a total of 835 unique PDB entries of antibodies that were crystallized in complex with their antigen and in a free state was constructed. It was examined for binding related conformation changes. We present further evidence supporting the theory of a pre-existing-equilibrium in experimental data. Multiple sequence alignments did not show binding induced tendencies in the solvent accessibility of residues in any specific position. Evaluating the changes in solvent accessibility per residue revealed a certain binding induced increase for several amino acids. Antibody-antigen interaction statistics were established and quantify a significant directional asymmetry between many interacting antibody and antigen residue pairs, especially a richness in tyrosine in the antibody epitope compared to its paratope. This asymmetry could potentially facilitate an increase in the success rate of computationally guided antibody refinement.

An In Silico Molecular Modelling-Based Prediction of Potential Keap1 Inhibitors from Hemidesmus indicus (L.) R.Br. against Oxidative-Stress-Induced Diseases.

The present study investigated the antioxidant potential of aqueous methanolic extracts of Hemidesmus indicus (L.) R.Br., followed by a pharmacoinformatics-based screening of novel Keap1 protein inhibitors. Initially, the antioxidant potential of this plant extract was assessed via antioxidant assays (DPPH, ABTS radical scavenging, and FRAP). Furthermore, 69 phytocompounds in total were derived from this plant using the IMPPAT database, and their three-dimensional structures were obtained from the PubChem database. The chosen 69 phytocompounds were docked against the Kelch-Neh2 complex protein (PDB entry ID: 2flu, resolution 1.50 Å) along with the standard drug (CPUY192018). H. indicus (L.) R.Br. extract (100 µg × mL-1) showed 85 ± 2.917%, 78.783 ± 0.24% of DPPH, ABTS radicals scavenging activity, and 161 ± 4 μg × mol (Fe (II)) g-1 ferric ion reducing power. The three top-scored hits, namely Hemidescine (-11.30 Kcal × mol-1), Beta-Amyrin (-10.00 Kcal × mol-1), and Quercetin (-9.80 Kcal × mol-1), were selected based on their binding affinities. MD simulation studies showed that all the protein-ligand complexes (Keap1-HEM, Keap1-BET, and Keap1-QUE) were highly stable during the entire simulation period, compared with the standard CPUY192018-Keap1 complex. Based on these findings, the three top-scored phytocompounds may be used as significant and safe Keap1 inhibitors, and could potentially be used for the treatment of oxidative-stress-induced health complications.

Crystal structure of a three-tetrad, parallel, K+-stabilized human telomeric G-quadruplex at 1.35 Å resolution.

The crystal structure of the G-rich human telomeric DNA Tel22 has been determined at 1.35 Å resolution in space group P6. Tel22 forms a non-canonical DNA structure called the G-quadruplex. The space group and unit-cell parameters are comparable to those in the crystal structures with PDB codes 6ip3 (1.40 Å resolution) and 1kf1 (2.15 Å resolution). The G-quadruplexes are highly similar in all of the structures. However, this structure of Tel22 displays clear density for polyethylene glycol and two potassium ions, which are located outside the ion channel in the G-quadruplex and play an important role in stabilizing the crystal contacts. In addition, 111 water molecules were identified (compared with 79 and 68 in PDB entries 6ip3 and 1kf1, respectively) that participate in intricate and extensive networks providing high stability to the G-quadruplex.

De Novo Computational Design of a Lipase with Hydrolysis Activity towards Middle-Chained Fatty Acid Esters.

Innovations in biocatalysts provide great prospects for intolerant environments or novel reactions. Due to the limited catalytic capacity and the long-term and labor-intensive characteristics of mining enzymes with the desired functions, de novo enzyme design was developed to obtain industrial application candidates in a rapid and convenient way. Here, based on the catalytic mechanisms and the known structures of proteins, we proposed a computational protein design strategy combining de novo enzyme design and laboratory-directed evolution. Starting with the theozyme constructed using a quantum-mechanical approach, the theoretical enzyme-skeleton combinations were assembled and optimized via the Rosetta "inside-out" protocol. A small number of designed sequences were experimentally screened using SDS-PAGE, mass spectrometry and a qualitative activity assay in which the designed enzyme 1a8uD1 exhibited a measurable hydrolysis activity of 24.25 ± 0.57 U/g towards p-nitrophenyl octanoate. To improve the activity of the designed enzyme, molecular dynamics simulations and the RosettaDesign application were utilized to further optimize the substrate binding mode and amino acid sequence, thus keeping the residues of theozyme intact. The redesigned lipase 1a8uD1-M8 displayed enhanced hydrolysis activity towards p-nitrophenyl octanoate-3.34 times higher than that of 1a8uD1. Meanwhile, the natural skeleton protein (PDB entry 1a8u) did not display any hydrolysis activity, confirming that the hydrolysis abilities of the designed 1a8uD1 and the redesigned 1a8uD1-M8 were devised from scratch. More importantly, the designed 1a8uD1-M8 was also able to hydrolyze the natural middle-chained substrate (glycerol trioctanoate), for which the activity was 27.67 ± 0.69 U/g. This study indicates that the strategy employed here has great potential to generate novel enzymes exhibiting the desired reactions.

Conformational transition of the Ixodes ricinus salivary serpin Iripin-4.

Iripin-4, one of the many salivary serpins from Ixodes ricinus ticks with an as-yet unexplained function, crystallized in two different structural conformations, namely the native partially relaxed state and the cleaved serpin. The native structure was solved at a resolution of 2.3 Å and the structure of the cleaved conformation was solved at 2.0 Å resolution. Furthermore, structural changes were observed when the reactive-centre loop transitioned from the native conformation to the cleaved conformation. In addition to this finding, it was confirmed that Glu341 represents a primary substrate-recognition site for the inhibitory mechanism. The presence of glutamate instead of the typical arginine in the P1 recognition site of all structurally characterized I. ricinus serpins (PDB entries 7b2t, 7pmu and 7ahp), except for the tyrosine in the P1 site of Iripin-2 (formerly IRS-2; PDB entry 3nda), would explain the absence of inhibition of the tested proteases that cleave their substrate after arginine. Further research on Iripin-4 should focus on functional analysis of this interesting serpin.

Unified access to up-to-date residue-level annotations from UniProtKB and other biological databases for PDB data

More than 61,000 proteins have up-to-date correspondence between their amino acid sequence (UniProtKB) and their 3D structures (PDB), enabled by the Structure Integration with Function, Taxonomy and Sequences (SIFTS) resource. SIFTS incorporates residue-level annotations from many other biological resources. SIFTS data is available in various formats like XML, CSV and TSV format or also accessible via the PDBe REST API but always maintained separately from the structure data (PDBx/mmCIF file) in the PDB archive. Here, we extended the wwPDB PDBx/mmCIF data dictionary with additional categories to accommodate SIFTS data and added the UniProtKB, Pfam, SCOP2, and CATH residue-level annotations directly into the PDBx/mmCIF files from the PDB archive. With the integrated UniProtKB annotations, these files now provide consistent numbering of residues in different PDB entries allowing easy comparison of structure models. The extended dictionary yields a more consistent, standardised metadata description without altering the core PDB information. This development enables up-to-date cross-reference information at the residue level resulting in better data interoperability, supporting improved data analysis and visualisation.

COVID-19Base v3: Update of the knowledgebase for drugs and biomedical entities linked to COVID-19.

COVID-19 has taken a huge toll on our lives over the last 3 years. Global initiatives put forward by all stakeholders are still in place to combat this pandemic and help us learn lessons for future ones. While the vaccine rollout was not able to curb the spread of the disease for all strains, the research community is still trying to develop effective therapeutics for COVID-19. Although Paxlovid and remdesivir have been approved by the FDA against COVID-19, they are not free of side effects. Therefore, the search for a therapeutic solution with high efficacy continues in the research community. To support this effort, in this latest version (v3) of COVID-19Base, we have summarized the biomedical entities linked to COVID-19 that have been highlighted in the scientific literature after the vaccine rollout. Eight different topic-specific dictionaries, i.e., gene, miRNA, lncRNA, PDB entries, disease, alternative medicines registered under clinical trials, drugs, and the side effects of drugs, were used to build this knowledgebase. We have introduced a BLSTM-based deep-learning model to predict the drug-disease associations that outperforms the existing model for the same purpose proposed in the earlier version of COVID-19Base. For the very first time, we have incorporated disease-gene, disease-miRNA, disease-lncRNA, and drug-PDB associations covering the largest number of biomedical entities related to COVID-19. We have provided examples of and insights into different biomedical entities covered in COVID-19Base to support the research community by incorporating all of these entities under a single platform to provide evidence-based support from the literature. COVID-19Base v3 can be accessed from: https://covidbase-v3.vercel.app/. The GitHub repository for the source code and data dictionaries is available to the community from: https://github.com/91Abdullah/covidbasev3.0.

Computer-assisted Design of Novel Polyketide Synthase 13 of Mycobacterium tuberculosis Inhibitors Using Molecular Modeling and Virtual Screening

Aims: Polyketide synthase 13 (Pks13) is an essential enzyme in the synthesis of mycolic acids biosynthesis pathway of Mycobacterium tuberculosis (Mtb). Therefore, Pks13 is a promising drug target for tuberculosis treatment. Here we report the in silico design and evaluation of novel Pks13 inhibitors made of benzofuran derivatives with favorable predicted pharmacokinetic profiles.&#x0D; Methodology: A 3D model of Pks13-TAMx complexes was prepared for a training set of 18 TAMs with experimentally determined inhibitory potencies (half-maximal inhibitory concentrations ) by using in situ modifications of the crystal structure of the TAM1-Pks13 complex (PDB entry 5V3X). A linear QSAR model was built, correlating computed gas phase enthalpies of formation ( ) of Pks13-TAMx complexes with the in order to find active conformations of the 18 TAMs. Furthermore, taking into account the implicit solvent effect and entropy changes upon ligand binding, a superior QSAR model was brought forth, correlating computed complexation of Gibbs’ free energies . Using the active conformations of the training set TAMs, we built a pharmacophore model (PH4) which was used to virtually screen novel analogs included in a virtual combinatorial library (VCL) of compounds containing benzofuran scaffolds. The PH4 model screened the VCL, which was formerly filtered by Lipinski’s rule-of-five, in order to identify new benzofuran analogs.&#x0D; Results: Gas phase QSAR model: , ; superior aqueous phase QSAR model , and PH4 pharmacophore model: , . The PH4-based screening retained 109 new TAM analogues. Finally, the predicted pharmacokinetic profiles of these new analogues were compared to current orally administered antituberculosis drugs, and the former were found to be almost 92 times more active than TAM2 ( .&#x0D; Conclusion: This computational approach, which combines molecular mechanics and the Poisson–Boltzmann (PB) implicit solvation theory, the pharmacophore model, the analysis of Pks13-TAMs interaction energies, the in silico screening of VCL compounds, and the inference of ADME properties resulted in a set of new suggested Pks13 inhibitors.