Protein structure modeling with MODELLER

Sequence-based identification of amyloidogenic β-hairpins reveals a prostatic acid phosphatase fragment promoting semen amyloid formation

β-Structure-rich amyloid fibrils are hallmarks of several diseases, including Alzheimer's (AD), Parkinson's (PD), and type 2 diabetes (T2D). While amyloid fibrils typically consist of parallel β-sheets, the anti-parallel β-hairpin is a structural motif accessible to amyloidogenic proteins in their monomeric and oligomeric states. Here, to investigate implications of β-hairpins in amyloid formation, potential β-hairpin-forming amyloidogenic segments in the human proteome were predicted based on sequence similarity with β-hairpins previously observed in Aβ, α-synuclein, and islet amyloid polypeptide, amyloidogenic proteins associated with AD, PD, and T2D, respectively. These three β-hairpins, established upon binding to the engineered binding protein β-wrapin AS10, are characterized by proximity of two sequence segments rich in hydrophobic and aromatic amino acids, with high β-aggregation scores according to the TANGO algorithm. Using these criteria, 2505 potential β-hairpin-forming amyloidogenic segments in 2098 human proteins were identified. Characterization of a test set of eight protein segments showed that seven assembled into Thioflavin T-positive aggregates and four formed β-hairpins in complex with AS10 according to NMR. One of those is a segment of prostatic acid phosphatase (PAP) comprising amino acids 185-208. PAP is naturally cleaved into fragments, including PAP(248-286) which forms functional amyloid in semen. We find that PAP(185-208) strongly decreases the protein concentrations required for fibril formation of PAP(248-286) and of another semen amyloid peptide, SEM1(86-107), indicating that it promotes nucleation of semen amyloids. In conclusion, β-hairpin-forming amyloidogenic protein segments could be identified in the human proteome with potential roles in functional or disease-related amyloid formation.

Virtual screening, molecular docking and dynamics simulation studies to identify potential agonists of orphan receptor GPR78 targeting CNS disorders

G protein-coupled receptors (GPCRs) are important targets in drug discovery because of their roles in physiological and pathological processes. Orphan GPCRs are GPCR proteins for which endogenous ligands have not yet been identified and they present interesting avenues for therapeutic intervention. This study focuses on GPR78, an orphan GPCR that is expressed in the central nervous system and linked to neurological disorders. GPR78 has no reported crystal structure and there is limited research. In this study, we have predicted the three dimensional model of GPR78 and its probable binding pocket. Structure-based virtual screening was carried out using the ChemDiv and Enamine REAL databases, followed by induced-fit docking studies to identify potential lead compounds with favorable interactions. These lead compounds were then embedded into a POPC lipid bilayer for a 200 ns molecular dynamics simulation. Free energy landscapes and MM-PBSA analyses were performed to assess the binding energies and conformational dynamics. The results highlight the dynamic nature of GPR78 in the presence of lead compounds and show favorable binding interactions. This study aims to predict a reliable three dimensional model of GPR78 and identify novel lead compounds through a comprehensive in silico approach. The identification of these potential GPR78 agonists represents a significant step in the development of new therapeutics for neurological disorders, highlighting the therapeutic potential of orphan GPR78 in CNS disorders.

Metal-ion-chelating phenylalanine nanostructures reverse immune dysfunction and sensitize breast tumour to immune checkpoint blockade

An immunosuppressive tumour microenvironment strongly influences response rates in patients receiving immune checkpoint blockade-based cancer immunotherapies, such as programmed death-1 (PD-1) and programmed death-ligand 1 (PD-L1). Here we demonstrate that metal-ion-chelating L-phenylalanine nanostructures synergize with short-term starvation (STS) to remodel the immunosuppressive microenvironment of breast and colorectal tumours. These nanostructures modulate the electrophysiological behaviour of dendritic cells and activate them through the NLRP3 inflammasome and calcium-mediated nuclear factor-κB pathway. STS promotes the cellular uptake of nanostructures through amino acid transporters and plays a key role in dendritic cell maturation and tumour-specific cytotoxic T lymphocyte responses. This study demonstrates the potential role of metal-ion-chelating L-phenylalanine nanostructures in activating immune responses and the effect of STS treatment in improving nanomaterial-mediated cancer immunotherapy.

Computational Modeling Study of the Molecular Basis of dNTP Selectivity in Human Terminal Deoxynucleotidyltransferase

Human terminal deoxynucleotidyl transferase (TdT) can catalyze template-independent DNA synthesis during the V(D)J recombination and DNA repair through nonhomologous end joining. The capacity for template-independent random addition of nucleotides to single-stranded DNA makes this polymerase useful in various molecular biological applications involving sequential stepwise synthesis of oligonucleotides using modified dNTP. Nonetheless, a serious limitation to the applications of this enzyme is strong selectivity of human TdT toward dNTPs in the order dGTP > dTTP ≈ dATP > dCTP. This study involved molecular dynamics to simulate a potential impact of amino acid substitutions on the enzyme's selectivity toward dNTPs. It was found that the formation of stable hydrogen bonds between a nitrogenous base and amino acid residues at positions 395 and 456 is crucial for the preferences for dNTPs. A set of single-substitution and double-substitution mutants at these positions was analyzed by molecular dynamics simulations. The data revealed two TdT mutants-containing either substitution D395N or substitutions D395N+E456N-that possess substantially equalized selectivity toward various dNTPs as compared to the wild-type enzyme. These results will enable rational design of TdT-like enzymes with equalized dNTP selectivity for biotechnological applications.

Outcomes of the EMDataResource cryo-EM Ligand Modeling Challenge

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein-nucleic acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: Escherichia coli beta-galactosidase with inhibitor, SARS-CoV-2 virus RNA-dependent RNA polymerase with covalently bound nucleotide analog and SARS-CoV-2 virus ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. The quality of submitted ligand models and surrounding atoms were analyzed by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics and contact scores. A composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Unveiling Tst3, a Multi-Target Gating Modifier Scorpion α Toxin from Tityus stigmurus Venom of Northeast Brazil: Evaluation and Comparison with Well-Studied Ts3 Toxin of Tityus serrulatus

Studies on the interaction sites of peptide toxins and ion channels typically involve site-directed mutations in toxins. However, natural mutant toxins exist among them, offering insights into how the evolutionary process has conserved crucial sequences for activities and molecular target selection. In this study, we present a comparative investigation using electrophysiological approaches and computational analysis between two alpha toxins from evolutionarily close scorpion species of the genus Tityus, namely, Tst3 and Ts3 from T. stigmurus and T. serrulatus, respectively. These toxins exhibit three natural substitutions near the C-terminal region, which is directly involved in the interaction between alpha toxins and Nav channels. Additionally, we characterized the activity of the Tst3 toxin on Nav1.1-Nav1.7 channels. The three natural changes between the toxins did not alter sensitivity to Nav1.4, maintaining similar intensities regarding their ability to alter opening probabilities, delay fast inactivation, and induce persistent currents. Computational analysis demonstrated a preference for the down conformation of VSD4 and a shift in the conformational equilibrium towards this state. This illustrates that the sequence of these toxins retained the necessary information, even with alterations in the interaction site region. Through electrophysiological and computational analyses, screening of the Tst3 toxin on sodium isoform revealed its classification as a classic α-NaTx with a broad spectrum of activity. It effectively delays fast inactivation across all tested isoforms. Structural analysis of molecular energetics at the interface of the VSD4-Tst3 complex further confirmed this effect.

Introgression of Δ1-pyrroline-5-carboxylate synthetase (PgP5CS) confers enhanced resistance to abiotic stresses in transgenic tobacco

Δ1-pyrroline-5-carboxylate synthetase (P5CS) is one of the key regulatory enzymes involved in the proline biosynthetic pathway. Proline acts as an osmoprotectant, molecular chaperone, antioxidant, and regulator of redox homeostasis. The accumulation of proline during stress is believed to confer tolerance in plants. In this study, we cloned the complete CDS of the P5CS from pearl millet (Pennisetum glaucum (L.) R.Br. and transformed into tobacco. Three transgenic tobacco plants with single-copy insertion were analyzed for drought and heat stress tolerance. No difference was observed between transgenic and wild-type (WT) plants when both were grown in normal conditions. However, under heat and drought, transgenic plants have been found to have higher chlorophyll, relative water, and proline content, and lower malondialdehyde (MDA) levels than WT plants. The photosynthetic parameters (stomatal conductance, intracellular CO2 concentration, and transpiration rate) were also observed to be high in transgenic plants under abiotic stress conditions. qRT-PCR analysis revealed that the expression of the transgene in drought and heat conditions was 2-10 and 2-7.5 fold higher than in normal conditions, respectively. Surprisingly, only P5CS was increased under heat stress conditions, indicating the possibility of feedback inhibition. Our results demonstrate the positive role of PgP5CS in enhancing abiotic stress tolerance in tobacco, suggesting its possible use to increase abiotic stress-tolerance in crops for sustained yield under adverse climatic conditions.

Systematic evaluation of lecithin:cholesterol acyltransferase binding sites in apolipoproteins via peptide based nanodiscs: regulatory role of charged residues at positions 4 and 7

Lecithin:cholesterol acyltransferase (LCAT) exhibits α-activity on high-density and β-activity on low-density lipoproteins. However, the molecular determinants governing LCAT activation by different apolipoproteins remain elusive. Uncovering these determinants would offer the opportunity to design and explore advanced therapies against dyslipidemias. Here, we have conducted coarse-grained and all-atom molecular dynamics simulations of LCAT with nanodiscs made with α-helical amphiphilic peptides either derived from apolipoproteins A1 and E (apoA1 and apoE) or apoA1 mimetic peptide 22A that was optimized to activate LCAT. This study aims to explore what drives the binding of peptides to our previously identified interaction site in LCAT. We hypothesized that this approach could be used to screen for binding sites of LCAT in different apolipoproteins and would provide insights to differently localized LCAT activities. Our screening approach was able to discriminate apoA1 helixes 4, 6, and 7 as key contributors to the interaction with LCAT supporting the previous research data. The simulations provided detailed molecular determinants driving the interaction with LCAT: the formation of hydrogen bonds or salt bridges between peptides E4 or D4 and LCAT S236 or K238 residues. Additionally, salt bridging between R7 and D73 was observed, depending on the availability of R7. Expanding our investigation to diverse plasma proteins, we detected novel LCAT binding helixes in apoL1, apoB100, and serum amyloid A. Our findings suggest that the same binding determinants, involving E4 or D4 -S236 and R7-D73 interactions, influence LCAT β-activity on low-density lipoproteins, where apoE and or apoB100 are hypothesized to interact with LCAT.

Characterization of a putative orexin receptor in Ciona intestinalis sheds light on the evolution of the orexin/hypocretin system in chordates

Tunicates are evolutionary model organisms bridging the gap between vertebrates and invertebrates. A genomic sequence in Ciona intestinalis (CiOX) shows high similarity to vertebrate orexin receptors and protostome allatotropin receptors (ATR). Here, molecular phylogeny suggested that CiOX is divergent from ATRs and human orexin receptors (hOX1/2). However, CiOX appears closer to hOX1/2 than to ATR both in terms of sequence percent identity and in its modelled binding cavity, as suggested by molecular modelling. CiOX was heterologously expressed in a recombinant HEK293 cell system. Human orexins weakly but concentration-dependently activated its Gq signalling (Ca2+ elevation), and the responses were inhibited by the non-selective orexin receptor antagonists TCS 1102 and almorexant, but only weakly by the OX1-selective antagonist SB-334867. Furthermore, the 5-/6-carboxytetramethylrhodamine (TAMRA)-labelled human orexin-A was able to bind to CiOX. Database mining was used to predict a potential endogenous C. intestinalis orexin peptide (Ci-orexin-A). Ci-orexin-A was able to displace TAMRA-orexin-A, but not to induce any calcium response at the CiOX. Consequently, we suggested that the orexin signalling system is conserved in Ciona intestinalis, although the relevant peptide-receptor interaction was not fully elucidated.

An interactive visualization tool for educational outreach in protein contact map overlap analysis

Recent advancements in contact map-based protein three-dimensional (3D) structure prediction have been driven by the evolution of deep learning algorithms. However, the gap in accessible software tools for novices in this domain remains a significant challenge. This study introduces GoFold, a novel, standalone graphical user interface (GUI) designed for beginners to perform contact map overlap (CMO) problems for better template selection. Unlike existing tools that cater more to research needs or assume foundational knowledge, GoFold offers an intuitive, user-friendly platform with comprehensive tutorials. It stands out in its ability to visually represent the CMO problem, allowing users to input proteins in various formats and explore the CMO problem. The educational value of GoFold is demonstrated through benchmarking against the state-of-the-art contact map overlap method, map_align, using two datasets: PSICOV and CAMEO. GoFold exhibits superior performance in terms of TM-score and Z-score metrics across diverse qualities of contact maps and target difficulties. Notably, GoFold runs efficiently on personal computers without any third-party dependencies, thereby making it accessible to the general public for promoting citizen science. The tool is freely available for download for macOS, Linux, and Windows.

APE-Gen2.0: Expanding Rapid Class I Peptide-Major Histocompatibility Complex Modeling to Post-Translational Modifications and Noncanonical Peptide Geometries

The recognition of peptides bound to class I major histocompatibility complex (MHC-I) receptors by T-cell receptors (TCRs) is a determinant of triggering the adaptive immune response. While the exact molecular features that drive the TCR recognition are still unknown, studies have suggested that the geometry of the joint peptide-MHC (pMHC) structure plays an important role. As such, there is a definite need for methods and tools that accurately predict the structure of the peptide bound to the MHC-I receptor. In the past few years, many pMHC structural modeling tools have emerged that provide high-quality modeled structures in the general case. However, there are numerous instances of non-canonical cases in the immunopeptidome that the majority of pMHC modeling tools do not attend to, most notably, peptides that exhibit non-standard amino acids and post-translational modifications (PTMs) or peptides that assume non-canonical geometries in the MHC binding cleft. Such chemical and structural properties have been shown to be present in neoantigens; therefore, accurate structural modeling of these instances can be vital for cancer immunotherapy. To this end, we have developed APE-Gen2.0, a tool that improves upon its predecessor and other pMHC modeling tools, both in terms of modeling accuracy and the available modeling range of non-canonical peptide cases. Some of the improvements include (i) the ability to model peptides that have different types of PTMs such as phosphorylation, nitration, and citrullination; (ii) a new and improved anchor identification routine in order to identify and model peptides that exhibit a non-canonical anchor conformation; and (iii) a web server that provides a platform for easy and accessible pMHC modeling. We further show that structures predicted by APE-Gen2.0 can be used to assess the effects that PTMs have in binding affinity in a more accurate manner than just using solely the sequence of the peptide. APE-Gen2.0 is freely available at https://apegen.kavrakilab.org.

The effect of phosphorylation on the conformational dynamics and allostery of the association of death-associated protein kinase with calmodulin

Protein phosphorylation plays an important role in the signal transduction and is capable of regulation of cell activity. The death-associated protein kinase 1 (DAPK1), as a Ser/Thr kinase, interacts with calmodulin (CaM) to regulate apoptotic and autophagic signaling. Autophosphorylation of DAPK1 at Ser308 located at the autoregulatory domain (ARD) blocks CaM binding and inhibits kinase catalytic activity. However, the mechanism underlying the influence of Ser308 phosphorylation (pS308) on the DAPK1 activity remains unclear. Here, we performed multiple, microsecond length molecular dynamics (MD) simulations, the molecular mechanics generalized Born/surface area (MM-GBSA) binding free energy calculations, principal component analysis, and dynamic cross-correlation analysis to unravel the conformational dynamics and allostery of the DAPK1 - CaM interaction triggered by the pS308 at the ARD. MD simulations showed that pS308 affected the conformational stability of the DAPK1 - CaM complex. Further energetic and structural exploration revealed that pS308 weakened the association of the phosphorylated DAPK1 to CaM, which lowered the susceptibility of DAPK1 to be activated by CaM. This result can provide mechanistic insights into the molecular underpinning through which the DAPK1 kinase activity is modulated by the auto-phosphorylation.Communicated by Ramaswamy H. Sarma.

Cinacalcet Reverses Short QT Interval in Familial Hypocalciuric Hypercalcemia Type 1

CONTEXT

Familial hypocalciuric hypercalcemia type 1 (FHH-1) defines an autosomal dominant disease, related to mutations in the CASR gene, with mild hypercalcemia in most cases. Cases of FHH-1 with a short QT interval have not been reported to date.

OBJECTIVE

Three family members presented with FHH-1 and short QT interval (<360 ms), a condition that could lead to cardiac arrhythmias, and the effects of cinacalcet, an allosteric modulator of the CaSR, in rectifying the abnormal sensitivity of the mutant CaSR and in correcting the short QT interval were determined.

METHODS

CASR mutational analysis was performed by next-generation sequencing and functional consequences of the identified CaSR variant (p.Ile555Thr), and effects of cinacalcet were assessed in HEK293 cells expressing wild-type and variant CaSRs. A cinacalcet test consisting of administration of 30 mg cinacalcet (8 Am) followed by hourly measurement of serum calcium, phosphate, and parathyroid hormone during 8 hours and an electrocardiogram was performed.

RESULTS

The CaSR variant (p.Ile555Thr) was confirmed in all 3 FHH-1 patients and was shown to be associated with a loss of function that was ameliorated by cinacalcet. Cinacalcet decreased parathyroid hormone by >50% within two hours, and decreases in serum calcium and increases in serum phosphate occurred within 8 hours, with rectification of the QT interval, which remained normal after 3 months of cinacalcet treatment.

CONCLUSION

Our results indicate that FHH-1 patients should be assessed for a short QT interval and a cinacalcet test used to select patients who are likely to benefit from this treatment.

Biochemical and cellular studies of three human 3-phosphoglycerate dehydrogenase variants responsible for pathological reduced L-serine levels

In the brain, the non-essential amino acid L-serine is produced through the phosphorylated pathway (PP) starting from the glycolytic intermediate 3-phosphoglycerate: among the different roles played by this amino acid, it can be converted into D-serine and glycine, the two main co-agonists of NMDA receptors. In humans, the enzymes of the PP, namely phosphoglycerate dehydrogenase (hPHGDH, which catalyzes the first and rate-limiting step of this pathway), 3-phosphoserine aminotransferase, and 3-phosphoserine phosphatase are likely organized in the cytosol as a metabolic assembly (a "serinosome"). The hPHGDH deficiency is a pathological condition biochemically characterized by reduced levels of L-serine in plasma and cerebrospinal fluid and clinically identified by severe neurological impairment. Here, three single-point variants responsible for hPHGDH deficiency and Neu-Laxova syndrome have been studied. Their biochemical characterization shows that V261M, V425M, and V490M substitutions alter either the kinetic (both maximal activity and Km for 3-phosphoglycerate in the physiological direction) and the structural properties (secondary, tertiary, and quaternary structure, favoring aggregation) of hPHGDH. All the three variants have been successfully ectopically expressed in U251 cells, thus the pathological effect is not due to hindered expression level. At the cellular level, mistargeting and aggregation phenomena have been observed in cells transiently expressing the pathological protein variants, as well as a reduced L-serine cellular level. Previous studies demonstrated that the pharmacological supplementation of L-serine in hPHGDH deficiencies could ameliorate some of the related symptoms: our results now suggest the use of additional and alternative therapeutic approaches.

Mechanistic insights into glycoside 3-oxidases involved in C-glycoside metabolism in soil microorganisms

C-glycosides are natural products with important biological activities but are recalcitrant to degradation. Glycoside 3-oxidases (G3Oxs) are recently identified bacterial flavo-oxidases from the glucose-methanol-coline (GMC) superfamily that catalyze the oxidation of C-glycosides with the concomitant reduction of O2 to H2O2. This oxidation is followed by C-C acid/base-assisted bond cleavage in two-step C-deglycosylation pathways. Soil and gut microorganisms have different oxidative enzymes, but the details of their catalytic mechanisms are largely unknown. Here, we report that PsG3Ox oxidizes at 50,000-fold higher specificity (kcat/Km) the glucose moiety of mangiferin to 3-keto-mangiferin than free D-glucose to 2-keto-glucose. Analysis of PsG3Ox X-ray crystal structures and PsG3Ox in complex with glucose and mangiferin, combined with mutagenesis and molecular dynamics simulations, reveal distinctive features in the topology surrounding the active site that favor catalytically competent conformational states suitable for recognition, stabilization, and oxidation of the glucose moiety of mangiferin. Furthermore, their distinction to pyranose 2-oxidases (P2Oxs) involved in wood decay and recycling is discussed from an evolutionary, structural, and functional viewpoint.

Enhancing alphafold-multimer-based protein complex structure prediction with MULTICOM in CASP15

To enhance the AlphaFold-Multimer-based protein complex structure prediction, we developed a quaternary structure prediction system (MULTICOM) to improve the input fed to AlphaFold-Multimer and evaluate and refine its outputs. MULTICOM samples diverse multiple sequence alignments (MSAs) and templates for AlphaFold-Multimer to generate structural predictions by using both traditional sequence alignments and Foldseek-based structure alignments, ranks structural predictions through multiple complementary metrics, and refines the structural predictions via a Foldseek structure alignment-based refinement method. The MULTICOM system with different implementations was blindly tested in the assembly structure prediction in the 15th Critical Assessment of Techniques for Protein Structure Prediction (CASP15) in 2022 as both server and human predictors. MULTICOM_qa ranked 3rd among 26 CASP15 server predictors and MULTICOM_human ranked 7th among 87 CASP15 server and human predictors. The average TM-score of the first predictions submitted by MULTICOM_qa for CASP15 assembly targets is ~0.76, 5.3% higher than ~0.72 of the standard AlphaFold-Multimer. The average TM-score of the best of top 5 predictions submitted by MULTICOM_qa is ~0.80, about 8% higher than ~0.74 of the standard AlphaFold-Multimer. Moreover, the Foldseek Structure Alignment-based Multimer structure Generation (FSAMG) method outperforms the widely used sequence alignment-based multimer structure generation.

Dietary γ-mangostin triggers immunogenic cell death and activates cGAS signaling in acute myeloid leukemia

Immunogenic cell death (ICD), one of cell-death types through release of damage-associated molecular patterns from dying tumor cells, activates tumor-specific immune response and elicits anti-tumor immunity by traditional radiotherapy and chemotherapy. However, whether natural products could induce ICD in leukemia is not elucidated. Here, we report dietary γ-mangostin eradicates murine primary leukemic cells and prolongs the survival of leukemic mice. As well, it restrains primary leukemic cells and CD34+ leukemic progenitor cells from leukemia patients. Strikingly, γ-mangostin attenuates leukemic cells by inducing ICD as characterized by expression of HSP90B1, ANXA1 and IL1B. Additionally, γ-mangostin accelerates cytoplasmic chromatin fragments generation, promoting DNA damage response, and enhances cGAS activation, leading to up-regulation of chemokines. Meanwhile, it induces HDAC4 degradation and acetylated histone H3 accumulation, which promotes chemokines transcription. Ultimately, CD8+ T cell is activated and recruited by γ-mangostin-induced chemokines in the microenvironment. Our study identifies γ-mangostin triggers ICD and activates cGAS signaling through DNA damage response and epigenetic modification. Therefore, dietary γ-mangostin would act as a potential agent to provoke anti-tumor immunity in the prevention and treatment of leukemia.

Binding of piperine to mycobacterial RNA polymerase improves the efficacy of rifampicin activity against Mycobacterium leprae and nontuberculous mycobacteria

Piperine (PPN) is a known inhibitor of efflux pumps in Mycobacterium tuberculosis and in vitro synergism with rifampicin (RIF) has been proven. The current study evaluates the activity of PPN and synergism with RIF in rapidly and slowly growing nontuberculous mycobacteria (NTM). Also, to propose a possible mechanism of interaction of PPN with M. leprae (Mlp) RNA polymerase (RNAp). Minimal inhibitory concentration and drug combination assay was determined by resazurin microtiter assay and resazurin drug combination assay, respectively. In silico evaluation of PPN binding was performed by molecular docking and molecular dynamics (MD). PPN showed higher antimicrobial activity against rapidly growing NTM (32-128 mg/L) rather than for slowly growing NTM (≥ 256 mg/L). Further, 77.8% of NTM tested exhibited FICI ≤ 0.5 when exposed to PPN and RIF combination, regardless of growth speed. Docking and MD simulations showed a possible PPN binding site at the interface between β and β' subunits of RNAp, in close proximity to the trigger-helix and bridge-helix elements. MD results indicated that PPN binding hindered the mobility of these elements, which are essential for RNA transcription. We hypothesize that PPN binding might affect mycobacterial RNAp activity, and, possibly, RIF activity and that this mechanism is partially responsible for synergic behaviors with RIF reported in vitro. Communicated by Ramaswamy H. Sarma.

The Wild-Type tRNA Adenosine Deaminase Enzyme TadA Is Capable of Sequence-Specific DNA Base Editing

Base editors are genome editing tools that enable site-specific base conversions through the chemical modification of nucleobases in DNA. Adenine base editors (ABEs) convert A ⋅ T to G ⋅ C base pairs in DNA by using an adenosine deaminase enzyme to modify target adenosines to inosine intermediates. Due to the lack of a naturally occurring adenosine deaminase that can modify DNA, ABEs were evolved from a tRNA-deaminating enzyme, TadA. Previous experiments with an ABE comprising a wild-type (wt) TadA showed no detectable activity on DNA, and directed evolution was therefore required to enable this enzyme to accept DNA as a substrate. Here we show that wtTadA can perform base editing in DNA in both bacterial and mammalian cells, with a strict sequence motif requirement of TAC. We leveraged this discovery to optimize a reporter assay to detect base editing levels as low as 0.01 %. Finally, we used this assay along with molecular dynamics simulations of full ABE:DNA complexes to better understand how the sequence recognition of mutant TadA variants change as they accumulate mutations to better edit DNA substrates.

Computational Peptide Design Cotargeting Glucagon and Glucagon-like Peptide-1 Receptors

Peptides are sustainable alternatives to conventional therapeutics for G protein-coupled receptor (GPCR) linked disorders, promising biocompatible and tailorable next-generation therapeutics for metabolic disorders including type-2 diabetes, as agonists of the glucagon receptor (GCGR) and the glucagon-like peptide-1 receptor (GLP-1R). However, single agonist peptides activating GLP-1R to stimulate insulin secretion also suppress obesity-linked glucagon release. Hence, bioactive peptides cotargeting GCGR and GLP-1R may remediate the blood glucose and fatty acid metabolism imbalance, tackling both diabetes and obesity to supersede current monoagonist therapy. Here, we design and model optimized peptide sequences starting from peptide sequences derived from earlier phage-displayed library screening, identifying those with predicted molecular binding profiles for dual agonism of GCGR and GLP-1R. We derive design rules from extensive molecular dynamics simulations based on peptide-receptor binding. Our newly designed coagonist peptide exhibits improved predicted coupled binding affinity for GCGR and GLP-1R relative to endogenous ligands and could in the future be tested experimentally, which may provide superior glycemic and weight loss control.

Targeting of multiple tumor-associated antigens by individual T cell receptors during successful cancer immunotherapy

The T cells of the immune system can target tumors and clear solid cancers following tumor-infiltrating lymphocyte (TIL) therapy. We used combinatorial peptide libraries and a proteomic database to reveal the antigen specificities of persistent cancer-specific T cell receptors (TCRs) following successful TIL therapy for stage IV malignant melanoma. Remarkably, individual TCRs could target multiple different tumor types via the HLA A∗02:01-restricted epitopes EAAGIGILTV, LLLGIGILVL, and NLSALGIFST from Melan A, BST2, and IMP2, respectively. Atomic structures of a TCR bound to all three antigens revealed the importance of the shared x-x-x-A/G-I/L-G-I-x-x-x recognition motif. Multi-epitope targeting allows individual T cells to attack cancer in several ways simultaneously. Such "multipronged" T cells exhibited superior recognition of cancer cells compared with conventional T cell recognition of individual epitopes, making them attractive candidates for the development of future immunotherapies.

Structure based virtual screening and molecular dynamics of natural anti-biofilm compounds against SagS response regulator/sensor kinase in Pseudomonas aeruginosa

SagS sensor regulator plays a vital role in biofilm development of Pseudomonas aeruginosa which subsequently makes the cells more tolerant to various antimicrobials. The multidrug resistance (MDR) issue has risen substantially in recent years and is considered a global threat. Therefore, alternative compounds should be unearthed immediately to address the issues related to P. aeruginosa drug resistance for which SagS could be a candidate. The present study is an attempt to screen natural anti-biofilm compounds as the potent inhibitors of SagS. Twenty natural anti-biofilm/quorum sensing inhibiting compounds were retrieved from various literatures with significant inhibitory effects against P. aeruginosa biofilm from in-vitro experiments which were screened using various pharmacokinetic parameters. The screened and three standard drugs were docked against SagS-HisKA using AutoDock 4.2 tool, which were further analysed by MD simulations to understand the binding mode of compounds and dynamic behaviour of the complexes. Two potential anti-biofilm natural compounds, pinocembrin with binding affinity (-7.19 kcal/mol), vestitol (-7.18 kcal/mol) and the standard drug ceftazidime (-8.89 kcal/mol) were selected based on filtered parameters and better binding affinity. The trajectory analysis of MD simulations reflected Pinocembrin in stabilizing the system compared to ceftazidime. The existing reports state that the natural products represent promising source of therapy with least or almost nil adverse effect compared to synthetic drugs which is well collated with our in-silico findings. This investigation can save both time and cost required for in-vitro and in-vivo analysis for designing of a novel anti-biofilm agent against P. aeruginosa biofilm-associated infections.Communicated by Ramaswamy H. Sarma.

Molecular simulations of DEAH-box helicases reveal control of domain flexibility by ligands: RNA, ATP, ADP, and G-patch proteins

DEAH-box helicases use the energy from ATP hydrolysis to translocate along RNA strands. They are composed of tandem RecA-like domains and a C-terminal domain connected by flexible linkers, and the activity of several DEAH-box helicases is regulated by cofactors called G-patch proteins. We used all-atom molecular dynamics simulations of the helicases Prp43, Prp22, and DHX15 in various liganded states to investigate how RNA, ADP, ATP, or G-patch proteins influence their conformational dynamics. The simulations suggest that apo helicases are highly flexible, whereas binding of RNA renders the helicases more rigid. ATP and ADP control the stability of the RecA1-RecA2 interface, but they have only a smaller effect on domain flexibility in absence of a RecA1-RecA2 interface. Binding of a G-patch protein to DHX15 imposes a more structured conformational ensemble, characterized by more defined relative domain arrangements and by an increased conformational stability of the RNA tunnel. However, the effect of the G-patch protein on domain dynamics is far more subtle as compared to the effects of RNA or ATP/ADP. The simulations characterize DEAH-box helicase as dynamic machines whose conformational ensembles are strongly defined by the presence of RNA, ATP, or ADP and only fine-tuned by the presence of G-patch proteins.

Structural insights into RNA bridging between HIV-1 Vif and antiviral factor APOBEC3G

Great effort has been devoted to discovering the basis of A3G-Vif interaction, the key event of HIV's counteraction mechanism to evade antiviral innate immune response. Here we show reconstitution of the A3G-Vif complex and subsequent A3G ubiquitination in vitro and report the cryo-EM structure of the A3G-Vif complex at 2.8 Å resolution using solubility-enhanced variants of A3G and Vif. We present an atomic model of the A3G-Vif interface, which assembles via known amino acid determinants. This assembly is not achieved by protein-protein interaction alone, but also involves RNA. The cryo-EM structure and in vitro ubiquitination assays identify an adenine/guanine base preference for the interaction and a unique Vif-ribose contact. This establishes the biological significance of an RNA ligand. Further assessment of interactions between A3G, Vif, and RNA ligands show that the A3G-Vif assembly and subsequent ubiquitination can be controlled by amino acid mutations at the interface or by polynucleotide modification, suggesting that a specific chemical moiety would be a promising pharmacophore to inhibit the A3G-Vif interaction.

High-throughput computational investigation of protein electrostatics and cavity for SAM-dependent methyltransferases

S-adenosyl methionine (SAM)-dependent methyl transferases (MTases) are a ubiquitous class of enzymes catalyzing dozens of essential life processes. Despite targeting a large space of substrates with diverse intrinsic reactivity, SAM MTases have similar catalytic efficiency. While understanding of MTase mechanism has grown tremendously through the integration of structural characterization, kinetic assays, and multiscale simulations, it remains elusive how these enzymes have evolved to fit the diverse chemical needs of their respective substrates. In this work, we performed a high-throughput molecular modeling analysis of 91 SAM MTases to better understand how their properties (i.e., electric field [EF] strength and active site volumes) help achieve similar catalytic efficiency toward substrates of different reactivity. We found that EF strengths have largely adjusted to make the target atom a better methyl acceptor. For MTases that target RNA/DNA and histone proteins, our results suggest that EF strength accommodates formal hybridization state and variation in cavity volume trends with diversity of substrate classes. Metal ions in SAM MTases contribute negatively to EF strength for methyl donation and enzyme scaffolds tend to offset these contributions.

Molecular Evolutionary Analyses of the RNA-Dependent RNA Polymerase (RdRp) Region and VP1 Gene in Human Norovirus Genotypes GII.P6-GII.6 and GII.P7-GII.6

To understand the evolution of GII.P6-GII.6 and GII.P7-GII.6 strains, the prevalent human norovirus genotypes, we analysed both the RdRp region and VP1 gene in globally collected strains using authentic bioinformatics technologies. A common ancestor of the P6- and P7-type RdRp region emerged approximately 50 years ago and a common ancestor of the P6- and P7-type VP1 gene emerged approximately 110 years ago. Subsequently, the RdRp region and VP1 gene evolved. Moreover, the evolutionary rates were significantly faster for the P6-type RdRp region and VP1 gene than for the P7-type RdRp region and VP1 genes. Large genetic divergence was observed in the P7-type RdRp region and VP1 gene compared with the P6-type RdRp region and VP1 gene. The phylodynamics of the RdRp region and VP1 gene fluctuated after the year 2000. Positive selection sites in VP1 proteins were located in the antigenicity-related protruding 2 domain, and these sites overlapped with conformational epitopes. These results suggest that the GII.6 VP1 gene and VP1 proteins evolved uniquely due to recombination between the P6- and P7-type RdRp regions in the HuNoV GII.P6-GII.6 and GII.P7-GII.6 virus strains.

Exploring Pyrrolo-Fused Heterocycles as Promising Anticancer Agents: An Integrated Synthetic, Biological, and Computational Approach

Five new series of pyrrolo-fused heterocycles were designed through a scaffold hybridization strategy as analogs of the well-known microtubule inhibitor phenstatin. Compounds were synthesized using the 1,3-dipolar cycloaddition of cycloimmonium N-ylides to ethyl propiolate as a key step. Selected compounds were then evaluated for anticancer activity and ability to inhibit tubulin polymerization in vitro. Notably, pyrrolo[1,2-a]quinoline 10a was active on most tested cell lines, performing better than control phenstatin in several cases, most notably on renal cancer cell line A498 (GI₅₀ 27 nM), while inhibiting tubulin polymerization in vitro. In addition, this compound was predicted to have a promising ADMET profile. The molecular details of the interaction between compound 10a and tubulin were investigated through in silico docking experiments, followed by molecular dynamics simulations and configurational entropy calculations. Of note, we found that some of the initially predicted interactions from docking experiments were not stable during molecular dynamics simulations, but that configurational entropy loss was similar in all three cases. Our results suggest that for compound 10a, docking experiments alone are not sufficient for the adequate description of interaction details in terms of target binding, which makes subsequent scaffold optimization more difficult and ultimately hinders drug design. Taken together, these results could help shape novel potent antiproliferative compounds with pyrrolo-fused heterocyclic cores, especially from an in silico methodological perspective.

Molecular docking in organic, inorganic, and hybrid systems: a tutorial review

Molecular docking simulation is a very popular and well-established computational approach and has been extensively used to understand molecular interactions between a natural organic molecule (ideally taken as a receptor) such as an enzyme, protein, DNA, RNA and a natural or synthetic organic/inorganic molecule (considered as a ligand). But the implementation of docking ideas to synthetic organic, inorganic, or hybrid systems is very limited with respect to their use as a receptor despite their huge popularity in different experimental systems. In this context, molecular docking can be an efficient computational tool for understanding the role of intermolecular interactions in hybrid systems that can help in designing materials on mesoscale for different applications. The current review focuses on the implementation of the docking method in organic, inorganic, and hybrid systems along with examples from different case studies. We describe different resources, including databases and tools required in the docking study and applications. The concept of docking techniques, types of docking models, and the role of different intermolecular interactions involved in the docking process to understand the binding mechanisms are explained. Finally, the challenges and limitations of dockings are also discussed in this review.

Graphical abstract

Enhancing AlphaFold-Multimer-based Protein Complex Structure Prediction with MULTICOM in CASP15

AlphaFold-Multimer has emerged as the state-of-the-art tool for predicting the quaternary structure of protein complexes (assemblies or multimers) since its release in 2021. To further enhance the AlphaFold-Multimer-based complex structure prediction, we developed a new quaternary structure prediction system (MULTICOM) to improve the input fed to AlphaFold-Multimer and evaluate and refine the outputs generated by AlphaFold2-Multimer. Specifically, MULTICOM samples diverse multiple sequence alignments (MSAs) and templates for AlphaFold-Multimer to generate structural models by using both traditional sequence alignments and new Foldseek-based structure alignments, ranks structural models through multiple complementary metrics, and refines the structural models via a Foldseek structure alignment-based refinement method. The MULTICOM system with different implementations was blindly tested in the assembly structure prediction in the 15th Critical Assessment of Techniques for Protein Structure Prediction (CASP15) in 2022 as both server and human predictors. Our server (MULTICOM_qa) ranked 3rd among 26 CASP15 server predictors and our human predictor (MULTICOM_human) ranked 7th among 87 CASP15 server and human predictors. The average TM-score of the first models predicted by MULTICOM_qa for CASP15 assembly targets is ~0.76, 5.3% higher than ~0.72 of the standard AlphaFold-Multimer. The average TM-score of the best of top 5 models predicted by MULTICOM_qa is ~0.80, about 8% higher than ~0.74 of the standard AlphaFold-Multimer. Moreover, the novel Foldseek Structure Alignment-based Model Generation (FSAMG) method based on AlphaFold-Multimer outperforms the widely used sequence alignment-based model generation. The source code of MULTICOM is available at: https://github.com/BioinfoMachineLearning/MULTICOM3.

Native doublet microtubules from Tetrahymena thermophila reveal the importance of outer junction proteins

Cilia are ubiquitous eukaryotic organelles responsible for cellular motility and sensory functions. The ciliary axoneme is a microtubule-based cytoskeleton consisting of two central singlets and nine outer doublet microtubules. Cryo-electron microscopy-based studies have revealed a complex network inside the lumen of both tubules composed of microtubule-inner proteins (MIPs). However, the functions of most MIPs remain unknown. Here, we present single-particle cryo-EM-based analyses of the Tetrahymena thermophila native doublet microtubule and identify 42 MIPs. These data shed light on the evolutionarily conserved and diversified roles of MIPs. In addition, we identified MIPs potentially responsible for the assembly and stability of the doublet outer junction. Knockout of the evolutionarily conserved outer junction component CFAP77 moderately diminishes Tetrahymena swimming speed and beat frequency, indicating the important role of CFAP77 and outer junction stability in cilia beating generation and/or regulation.

Molecular analyses of the C-terminal CRAF variants associated with cardiomyopathy reveal their opposing impacts on the active conformation of the kinase domain

The germline mutations in the C-terminus of CRAF kinase, particularly L603, and S612T/L613V, are associated with congenital heart disorders, for example, dilated cardiomyopathy (DCM) and hypertrophic cardiomyopathy (HCM). The experimental data suggest that genetic alternation at position 603 impairs, while those at positions 612/613 enhance the CRAF kinase activity. However, the underlying mechanistic details by which these mutations increase or decrease kinase activity remain elusive. Therefore, we applied molecular dynamic simulation to investigate the impacts of these point mutations on the conformation of the CRAF kinase domain. The results revealed that the substitution of Leucine 603 for proline transits the kinase domain to a state that exhibits the molecular hallmarks of an inactive kinase, for example, a closed activation loop, 'αC-helix out' conformation and a distorted regulatory hydrophobic spine. However, two HCM-associated variants (S612T and L613V) show features of an active conformation, such as an open activation loop conformation, 'αC-helix in', the assembly of the hydrophobic spine, and more surface-exposed catalytic residues of phosphoryl transfer reaction. Overall, our study provides a mechanistic basis for the contradictory effects of the CRAF variants associated with HCM and DCM.

Tyrosine phosphorylation controlled poly(A) polymerase I activity regulates general stress response in bacteria

RNA 3'-end polyadenylation that marks transcripts for degradation is implicated in general stress response in Escherichia coli Yet, the mechanism and regulation of poly(A) polymerase I (PAPI) in stress response are obscure. We show that pcnB (that encodes PAPI)-null mutation widely stabilises stress response mRNAs and imparts cellular tolerance to multiple stresses, whereas PAPI ectopic expression renders cells stress-sensitive. We demonstrate that there is a substantial loss of PAPI activity on stress exposure that functionally phenocopies pcnB-null mutation stabilising target mRNAs. We identify PAPI tyrosine phosphorylation at the 202 residue (Y202) that is enormously enhanced on stress exposure. This phosphorylation inhibits PAPI polyadenylation activity under stress. Consequentially, PAPI phosphodeficient mutation (tyrosine 202 to phenylalanine, Y202F) fails to stimulate mRNA expression rendering cells stress-sensitive. Bacterial tyrosine kinase Wzc phosphorylates PAPI-Y202 residue, and that wzc-null mutation renders cells stress-sensitive. Accordingly, wzc-null mutation has no effect on stress sensitivity in the presence of pcnB-null or pcnB-Y202F mutation. We also establish that PAPI phosphorylation-dependent stress tolerance mechanism is distinct and operates downstream of the primary stress regulator RpoS.

Emergence of a proton exchange-based isomerization and lactonization mechanism in the plant coumarin synthase COSY

Plants contain rapidly evolving specialized enzymes that support the biosynthesis of functionally diverse natural products. In coumarin biosynthesis, a BAHD acyltransferase-family enzyme COSY was recently discovered to accelerate coumarin formation as the only known BAHD enzyme to catalyze an intramolecular acyl transfer reaction. Here we investigate the structural and mechanistic basis for COSY's coumarin synthase activity. Our structural analyses reveal an unconventional active-site configuration adapted to COSY's specialized activity. Through mutagenesis studies and deuterium exchange experiments, we identify a unique proton exchange mechanism at the α-carbon of the o-hydroxylated trans-hydroxycinnamoyl-CoA substrates during the catalytic cycle of COSY. Quantum mechanical cluster modeling and molecular dynamics further support this key mechanism for lowering the activation energy of the rate-limiting trans-to-cis isomerization step in coumarin production. This study unveils an unconventional catalytic mechanism mediated by a BAHD-family enzyme, and sheds light on COSY's evolutionary origin and its recruitment to coumarin biosynthesis in eudicots.

In silico investigation of the interactions of certain drugs proposed for the treatment of Covid-19 with the paraoxonase-1

Coronavirus disease 2019 (Covid-19) has caused one of the biggest pandemics of modern times, infected over 240 million people and killed over 4.9 million people, and continues to do so. Although many drugs are widely recommended in the treatment of this disease, the interactions of these drugs with an anti-atherosclerotic enzyme, paraoxonase-1 (PON1), are not well known. In our study, we investigated the interactions of 18 different drugs, which are claimed to be effective against covid-19, with the PON1 enzyme and its genetics variants L55M and Q192R with molecular docking, molecular dynamics simulation and free energy calculation method MM/PBSA. We found that ruxolitinib, dexamethasone, colchicine; dexamethasone, sitagliptin, baricitinib and galidesivir, ruxolitinib, hydroxychloroquine were the most effective compounds in binding PON1-w, PON1L55M and PON1Q192R respectively. Mainly, sitagliptin, galidesivir and hydroxychloroquine have attracted attention by showing very high affinity (<-300 kJ/mol) according to the MM/PBSA method. We concluded that the drug interactions should be considered and more attention should be paid in the use of these drugs.Communicated by Ramaswamy H. Sarma.

Regulatory mechanisms of the dynein-2 motility by post-translational modification revealed by MD simulation

Intraflagellar transport for ciliary assembly and maintenance is driven by dynein and kinesins specific to the cilia. It has been shown that anterograde and retrograde transports run on different regions of the doublet microtubule, i.e., separate train tracks. However, little is known about the regulatory mechanism of this selective process. Since the doublet microtubule is known to display specific post-translational modifications of tubulins, i.e., "tubulin code", for molecular motor regulations, we investigated the motility of ciliary specific dynein-2 under different post-translational modification by coarse-grained molecular dynamics. Our setup allows us to simulate the landing behaviors of dynein-2 on un-modified, detyrosinated, poly-glutamylated and poly-glycylated microtubules in silico. Our study revealed that poly-glutamylation can play an inhibitory effect on dynein-2 motility. Our result indicates that poly-glutamylation of the B-tubule of the doublet microtubule can be used as an efficient means to target retrograde intraflagellar transport onto the A-tubule.

Comparative Dynamic Features of Apo and Bound MDM2 Protein Reveal the Mechanism of Inhibitor Recognition for Anti-Cancer Activity

BACKGROUND

Mouse Double Minute 2 Homolog (MDM2) oncogenic protein is the principal cellular antagonist of the p53 tumor suppressor gene. Restoration of p53 activity by inhibiting the MDM2-P53 interactions at the molecular level has become the cornerstone of cancer research due to its promising anticancer effects. Natural medicinal products possess various chemical structures and represent an essential source for drug discovery. α-Mangostin (AM) and gambogic acid (G250) are plant-derived compounds that showed inhibitory effects on MDM2-P53 interactions in vitro and in vivo.

METHODS

Despite the many clinical studies which performed deeper insight about the molecular understanding of the structural mechanisms exhibited by α-Mangostin and Gambogic acid-binding to MDM2 remains critical. In this study, comparative molecular dynamics simulations were performed for each Apo and bound p53 and MDM2 proteins to shed light on the MDM2-p53 interactions and get a better understanding of the inhibition mechanisms.

RESULTS

Results revealed atomistic interaction of AM and G250 within the MDM2-p53 interaction cleft. Both compounds mediate the interaction between the α-helix motifs of the p53 amino-terminal domain, which caused a significant separation between orthogonally opposed residues, specifically Lys8 and Gly47 residues of the p53 and MDM2, respectively. Contrasting changes in magnitudes were observed in per-residue fluctuation on AM and G250 (~0.04 nm and ~2.3 nm, respectively). The Radius of gyration (~0.03 nm and 0.04 nm, respectively), C-alpha deviations (~0.06 nm and 0.1 nm, respectively). The phenolic group of AM was found to establish hydrogen interactions with Glu28 and His96 residues of MDM2. The trioxahexacyclo-ring of G250 also forms hydrogen bond interactions with Lys51 and Leu26 residues of MDM2.

CONCLUSION

Utilizing the information provided on the inhibitory binding mode adopted by each compound in this study may further assist in the tailored designs for cancer therapeutics.

A proteome-scale analysis of vertebrate protein amino acid occurrence: Thermoadaptation shows a correlation with protein solvation but less so with dynamics

Despite differences in behaviors and living conditions, vertebrate organisms share the great majority of proteins, often with subtle differences in amino acid sequence. Here, we present a simple way to analyze the difference in amino acid occurrence by comparing highly homologous proteins on a subproteome level between several vertebrate model organisms. Specifically, we use this method to identify a pattern of amino acid conservation as well as a shift in amino acid occurrence between homeotherms (warm-blooded species) and poikilotherms (cold-blooded species). Importantly, this general analysis and a specific example further establish a broad correlation, if not likely connection between the thermal adaptation of protein sequences and two of their physical features: on average a change in their protein dynamics and, even more strongly, in their solvation. For poikilotherms, such as frog and fish, the lower body temperature is expected to increase the protein-protein interaction due to a decrease in protein internal dynamics. In order to counteract the tendency for enhanced binding caused by low temperatures, poikilotherms enhance the solvation of their proteins by favoring polar amino acids. This feature appears to dominate over possible changes in dynamics for some proteins. The results suggest that a general trend for amino acid choice is part of the mechanism for thermoadaptation of vertebrate organisms at the molecular level.

Insights into pralsetinib resistance to the non-gatekeeper RET kinase G810C mutation through molecular dynamics simulations

OBJECTIVE

RET (rearranged during transfection) kinase, as a transmembrane receptor tyrosine kinase, is a therapeutic target for several human cancer such as non-small cell lung cancer (NSCLC) and thyroid cancer. Pralsetinib is a recently approved drug for the treatment of RET-driven NSCLC and thyroid cancers. A single point mutation G810C at the C-lobe of the RET kinase causes pralsetinib resistance to this non-gatekeeper variant. However, the detailed mechanism remains poorly understood.

METHODS

Here, multiple microsecond molecular dynamics (MD) simulations, molecular mechanics/generalized born surface area (MM/GBSA) binding free energy calculations, and community network analysis were performed to reveal the mechanism of pralsetinib resistance to the RET G810C mutant.

RESULTS

The simulations showed that the G810C mutation had a minor effect on the overall conformational dynamics of the RET kinase domain. Energetic analysis suggested that the G810C mutation reduced the binding affinity of pralsetinib to the mutant. Per-residue energy contribution and structural analyses revealed that the hydrogen bonding interactions between pralsetinib and the hinge residues Glu805 and Ala807 were disrupted in the G810C mutant, which were responsible for the decreased binding affinity of pralsetinib to the mutant.

CONCLUSIONS

The obtained results may provide understanding of the mechanism of pralsetinib resistance to the non-gatekeeper RET G810C mutant.

Molecular modeling and dynamics simulation of alcohol dehydrogenase enzyme from high efficacy cellulosic ethanol-producing yeast mutant strain Pichia kudriavzevii BGY1-γm

One of the major constraints limiting the use of abundantly available lignocellulosic biomass as potential feedstock for alcohol industry is the lack of C₆/C₅ co-sugar fermenting yeast. The present study explores a mutant yeast Pichia kudriavzevii BGY1-γm as a potential strain for bioconversion of glucose/xylose sugars of green biomass into ethanol under batch fermentation. The mutant strain having higher alcohol dehydrogenase activity (11.31%) showed significantly higher ethanol concentration during co-fermentation of glucose/xylose sugars (14.2%) as compared to the native strain. Based on 99% sequence similarity of ADH encoding gene from the mutant with the gene sequences from other yeast strains, the ADH enzyme was identified as ADH-1 type. The study reveals first three-dimensional model of ADH-1 utilizing glucose/xylose sugars from P. kudriavzevii BGY1-γm (PkADH mutant). The refined and validated model of PkADH mutant was used for molecular docking against the substrate (acetaldehyde) and product (ethanol). Molecular docking results showed that substrate and product exhibited a binding affinity of -4.55 and -4.5 kcal/mol with PkADH mutant. Acetaldehyde and ethanol interacted at the active site of PkADH mutant via hydrogen bonds (Ser42, His69 and Asp163) and hydrophobic interactions (Cys40, Ser42, His69, Cys95, Trp123 and Asp163) to form the stable protein-ligand complex. Molecular dynamics analysis revealed that PkADH-mutant acetaldehyde and PkADH-mutant ethanol complexes were more stable than PkADH mutant. MMPBSA binding energy confirmed that binding of substrate and product results in the formation of a lower energy stable protein-ligand complex.Communicated by Ramaswamy H. Sarma.

Identification of the structural features of quinazoline derivatives as EGFR inhibitors using 3D-QSAR modeling, molecular docking, molecular dynamics simulations and free energy calculations

Epidermal growth factor receptor (EGFR) is a promising target for the treatment of different types of malignant tumors. Therefore, a combined molecular modeling study was performed on a series of quinazoline derivatives as EGFR inhibitors. The optimum ligand-based CoMFA and CoMSIA models showed reliable and satisfactory predictability (with R²_cv=0.681, R²_ncv=0.844, R²_pred=0.8702 and R²_cv=0.643, R²_ncv=0.874, R²_pred=0.6423). The derived contour maps provide structural features to improve inhibitory activity. Furthermore, the contour maps, molecular docking, and molecular dynamics (MD) simulations have good consistency, illustrating that the derived models are reliable. In addition, MD simulations and binding free energy calculations were also carried out to understand the conformational fluctuations at the binding pocket of the receptor. The results indicate that hydrogen bond, hydrophobic and electrostatic interactions play significant roles on activity and selectivity. Furthermore, amino acids Val31, Lys50, Thr95, Leu149 and Asp160 are considered as essential residues to participate in the ligand-receptor interactions. Overall, this work would offer reliable theoretical basis for future structural modification, design and synthesis of novel EGFR inhibitors with good potency.Communicated by Ramaswamy H. Sarma.

Molecular dynamics simulations reveal membrane lipid interactions of the full-length lymphocyte specific kinase (Lck)

The membrane-bound lymphocyte-specific protein-tyrosine kinase (Lck) triggers T cell antigen receptor signalling to initiate adaptive immune responses. Despite many structure-function studies, the mode of action of Lck and the potential role of plasma membrane lipids in regulating Lck's activity remains elusive. Advances in molecular dynamics simulations of membrane proteins in complex lipid bilayers have opened a new perspective in gathering such information. Here, we have modelled the full-length Lck open and closed conformations  using data available from different crystalographic studies and simulated its interaction with the inner leaflet of the T cell plasma membrane. In both conformations, we found that the unstructured unique domain and the structured domains including the kinase interacted with the membrane with a preference for PIP lipids. Interestingly, our simulations suggest that the Lck-SH2 domain interacts with lipids differently in the open and closed Lck conformations, demonstrating that lipid interaction can potentially regulate Lck's conformation and in turn modulate T cell signalling. Additionally, the Lck-SH2 and kinase domain residues that significantly contacted PIP lipids are found to be conserved among the Src family of kinases, thereby potentially representing similar PIP interactions within the family.

Data on structural analysis of cholesterol binding and sterol selectivity by ABCG5/G8

ATP-Binding cassette subfamily G (ABCG) sterol transporters maintain whole body endogenous and exogenous sterol homeostasis. A substantial portion of exogenous sterols are undigestible phytosterols (plant sterols), which can introduce complications when accumulated. ABCG5/G8 is the main protein functioning to remove ingested plant sterols providing protection from their toxic effects, although, the structural features behind substrate binding in ABCG5/G8 remain poorly resolved. Within this data article, we present extended preceding in the determination of the cholesterol-bound crystal structure and the sterol docking analysis. The crystal structure was deposited in the Protein Data Bank with the accession number of 8CUB, whereas the diffraction images were deposited at the SBGrid Data Bank. This dataset follows the research article entitled as "Structural analysis of cholesterol binding and sterol selectivity by ABCG5/G8" (doi: 10.1016/j.jmb.2022.167795).

Phosphatidylinositol 4,5-bisphosphate (PIP2) facilitates norepinephrine transporter dimerization and modulates substrate efflux

The plasmalemmal norepinephrine transporter (NET) regulates cardiovascular sympathetic activity by clearing extracellular norepinephrine in the synaptic cleft. Here, we investigate the subunit stoichiometry and function of NET using single-molecule fluorescence microscopy and flux assays. In particular, we show the effect of phosphatidylinositol 4,5-bisphosphate (PIP2) on NET oligomerization and efflux. NET forms monomers (~60%) and dimers (~40%) at the plasma membrane. PIP2 depletion results in a decrease in the average oligomeric state and decreases NET-mediated substrate efflux while not affecting substrate uptake. Mutation of the putative PIP2 binding residues R121, K334, and R440 to alanines does not affect NET dimerization but results in decreased substrate efflux that is not altered upon PIP2 depletion; this indicates that PIP2 interactions with these residues affect NET-mediated efflux. A dysregulation of norepinephrine and PIP2 signaling have both been implicated in neuropsychiatric and cardiovascular diseases. This study provides evidence that PIP2 directly regulates NET organization and function.

The Third Extracellular Loop of Mammalian Odorant Receptors Is Involved in Ligand Binding

Mammals recognize chemicals in the air via G protein-coupled odorant receptors (ORs). In addition to their orthosteric binding site, other segments of these receptors modulate ligand recognition. Focusing on human hOR1A1, which is considered prototypical of class II ORs, we used a combination of molecular modeling, site-directed mutagenesis, and in vitro functional assays. We showed that the third extracellular loop of ORs (ECL3) contributes to ligand recognition and receptor activation. Indeed, site-directed mutations in ECL3 showed differential effects on the potency and efficacy of both carvones, citronellol, and 2-nonanone.

Protein Data Bank: A Comprehensive Review of 3D Structure Holdings and Worldwide Utilization by Researchers, Educators, and Students

The Research Collaboratory for Structural Bioinformatics Protein Data Bank (RCSB PDB), funded by the United States National Science Foundation, National Institutes of Health, and Department of Energy, supports structural biologists and Protein Data Bank (PDB) data users around the world. The RCSB PDB, a founding member of the Worldwide Protein Data Bank (wwPDB) partnership, serves as the US data center for the global PDB archive housing experimentally-determined three-dimensional (3D) structure data for biological macromolecules. As the wwPDB-designated Archive Keeper, RCSB PDB is also responsible for the security of PDB data and weekly update of the archive. RCSB PDB serves tens of thousands of data depositors (using macromolecular crystallography, nuclear magnetic resonance spectroscopy, electron microscopy, and micro-electron diffraction) annually working on all permanently inhabited continents. RCSB PDB makes PDB data available from its research-focused web portal at no charge and without usage restrictions to many millions of PDB data consumers around the globe. It also provides educators, students, and the general public with an introduction to the PDB and related training materials through its outreach and education-focused web portal. This review article describes growth of the PDB, examines evolution of experimental methods for structure determination viewed through the lens of the PDB archive, and provides a detailed accounting of PDB archival holdings and their utilization by researchers, educators, and students worldwide.

Antibacterial properties and in silico modeling perspective of nano ZnO transported oxytetracycline-Zn2+ complex [ZnOTc]+ against oxytetracycline-resistant Aeromonas hydrophila

Emergence of antibiotics resistance has threatening consequences not only for human health but also for animal health issues in agriculture. Several animal pathogenic bacteria have developed antibiotic resistance and managing same has tremendous cost repercussions and may lead to total harvest loss. Hence in the present study, efforts are made to revitalize an old antibiotic molecule, oxytetracycline (OTc), through nanodelivery approaches using zinc oxide nanoparticles (nZnO) to confront OTc resistant fish pathogenic bacteria Aeromonas hydrophila. OTc was impregnated in nZnO through in situ precipitation method to develop OTc loaded ZnO nanoparticles (OTc@nZnO) with average size of 99.42 nm. Spectroscopic investigation of same revealed complexation of Zn²⁺ with amide and aromatic carbonyl moieties of OTc [ZnOTc]⁺. The complex performed better against A. hydrophila with 7-15 mm inhibition zone as compared to nil for bare OTc at same dose. OTc also showed MIC of 150 µg ml^-1 and for OTc@nZnO it was 7.02 µg ml^-1 with faster killing rate (k, -0.95). In silico docking simulation suggest that [ZnOTc]⁺ had low binding affinity (LBE > -5.00 kcal mol^-1) toward TetR(E) and TetA(E) proteins of A. hydrophila as compared to OTc (LBE < -8.00 kcal mol^-1). This study postulates that [ZnOTc]⁺ released from OTc@nZnO can escape TetR(E) and TetA(E) resistance proteins and bind at 30S ribosomal subunit with high affinity (<-11.00 kcal mol^-1) to exert antibacterial properties. In the recent scenario of recurrent antimicrobial resistance, the develop antibiotic-nanocomposites could come out as potential solution, however further study is required for its feasibility for use in animal health care.

In silico study on the Hepatitis E virus RNA Helicase and its inhibition by silvestrol, rocaglamide and other flavagline compounds

Hepatitis E Virus (HEV) follows waterborne or zoonotic/foodborne transmission. Genotype 3 HEV infections are worldwide spread, especially in swine populations, representing an emerging threat for human health, both for farm workers and pork meat consumers. Unfortunately, HEV in vitro culture and analysis are still difficult, resulting in a poor understanding of its biology and hampering the implementation of counteracting strategies. Indeed, HEV encodes for only one non-structural multifunctional and multidomain protein (ORF1), which might be a good candidate for anti-HEV drugging strategies. In this context, an in silico molecular modelling approach that consisted in homology modelling to derive the 3D model target, docking study to simulate the binding event, and molecular dynamics to check complex stability over time was used. This workflow succeeded to describe ORF1 RNA Helicase domain from a molecular standpoint allowing the identification of potential inhibitory compounds among natural plant-based flavagline-related molecules such as silvestrol, rocaglamide and derivatives thereof. In the context of scouting potential anti-viral compounds and relying on the outcomes presented, further dedicated investigations on silvestrol, rocaglamide and a promising oxidized derivative have been suggested. For the sake of data reproducibility, the 3D model of HEV RNA Helicase has been made publicly available.