Integrating protein structural dynamics and evolutionary analysis with Bio3D

Mechanistic insights into P-glycoprotein ligand transport and inhibition revealed by enhanced molecular dynamics simulations

P-glycoprotein (P-gp) plays a crucial role in cellular detoxification and drug efflux processes, transitioning between inward-facing (IF) open, occluded, and outward-facing (OF) states to facilitate substrate transport. Its role is critical in cancer therapy, where P-gp contributes to the multidrug resistance phenotype. In our study, classical and enhanced molecular dynamics (MD) simulations were conducted to dissect the structural and functional features of the P-gp conformational states. Our advanced MD simulations, including kinetically excited targeted MD (ketMD) and adiabatic biasing MD (ABMD), provided deeper insights into state transition and translocation mechanisms. Our findings suggest that the unkinking of TM4 and TM10 helices is a prerequisite for correctly achieving the outward conformation. Simulations of the IF-occluded conformations, characterized by kinked TM4 and TM10 helices, consistently demonstrated altered communication between the transmembrane domains (TMDs) and nucleotide binding domain 2 (NBD2), suggesting the implication of this interface in inhibiting P-gp's efflux function. A particular emphasis was placed on the unstructured linker segment connecting the NBD1 to TMD2 and its role in the transporter's dynamics. With the linker present, we specifically noticed a potential entrance of cholesterol (CHOL) through the TM4-TM6 portal, shedding light on crucial residues involved in accommodating CHOL. We therefore suggest that this entry mechanism could be employed for some P-gp substrates or inhibitors. Our results provide critical data for understanding P-gp functioning and developing new P-gp inhibitors for establishing more effective strategies against multidrug resistance.

Construction of highly active and stable recombinant nattokinase by engineered bacteria and computational design

Nattokinase (NK) is an enzyme that has been recognized as a new potential thrombolytic drug due to its strong thrombolytic activity. However, it is difficult to maintain the enzyme activity of NK during high temperature environment of industrial production. In this study, we constructed six NK mutants with potential for higher thermostability using a rational protein engineering strategy integrating free energy-based methods and molecular dynamics (MD) simulation. Then, wild-type NK and NK mutants were expressed in Escherichia coli (E. coli), and their thermostability and thrombolytic activity were tested. The results showed that, compared with wild-type NK, the mutants Y256P, Q206L and E156F all had improved thermostability. The optimal mutant Y256P showed a higher melting temperature (Tm) of 77.4 °C, an increase of 4 °C in maximum heat-resistant temperature and an increase of 51.8 % in activity at 37 °C compared with wild-type NK. Moreover, we also explored the mechanism of the increased thermostability of these mutants by analysing the MD trajectories under different simulation temperatures.

Investigation of the Effect of Peptide p5 Targeting CDK5-p25 Hyperactivity on Munc18-1 (P67) Regulating Neuronal Exocytosis Using Molecular Simulations

Munc18-1 is an SM (sec1/munc-like) family protein involved in vesicle fusion and neuronal exocytosis. Munc18-1 is known to regulate the exocytosis process by binding with closed- and open-state conformations of Syntaxin1, a protein belonging to the SNARE family established to be central to the exocytosis process. Our previous work studied peptide p5 as a promising drug candidate for CDK5-p25 complex, an Alzheimer's disease (AD) pathological target. Experimental in vivo and in vitro studies suggest that Munc18-1 promotes p5 to selectively inhibit the CDK5-p25 complex without affecting the endogenous CDK5 activity, a characteristic of remarkable therapeutic implications. In this paper, we identify several binding modes of p5 with Munc18-1 that could potentially affect the Munc18-1 binding with SNARE proteins and lead to off-target effects on neuronal communication using molecular dynamics simulations. Recent studies indicate that disruption of Munc18-1 function not only disrupts neurotransmitter release but also results in neurodegeneration, exhibiting clinical resemblance to other neurodegenerative conditions such as AD, causing diagnostic and treatment challenges. We characterize such interactions between p5 and Munc18-1, define the corresponding pharmacophores, and provide guidance for the in vitro validation of our findings to improve therapeutic efficacy and safety of p5.

Explaining Conformational Diversity in Protein Families through Molecular Motions

Proteins play a central role in biological processes, and understanding their conformational variability is crucial for unraveling their functional mechanisms. Recent advancements in high-throughput technologies have enhanced our knowledge of protein structures, yet predicting their multiple conformational states and motions remains challenging. This study introduces Dimensionality Analysis for protein Conformational Exploration (DANCE) for a systematic and comprehensive description of protein families conformational variability. DANCE accommodates both experimental and predicted structures. It is suitable for analysing anything from single proteins to superfamilies. Employing it, we clustered all experimentally resolved protein structures available in the Protein Data Bank into conformational collections and characterized them as sets of linear motions. The resource facilitates access and exploitation of the multiple states adopted by a protein and its homologs. Beyond descriptive analysis, we assessed classical dimensionality reduction techniques for sampling unseen states on a representative benchmark. This work improves our understanding of how proteins deform to perform their functions and opens ways to a standardised evaluation of methods designed to sample and generate protein conformations.

MDFF_NM: Improved Molecular Dynamics Flexible Fitting into Cryo-EM Density Maps with a Multireplica Normal Mode-Based Search

Molecular Dynamics Flexible Fitting (MDFF) is a widely used tool to refine high-resolution structures into cryo-EM density maps. Despite many successful applications, MDFF is still limited by its high computational cost, overfitting, accuracy, and performance issues due to entrapment within wrong local minima. Modern ensemble-based MDFF tools have generated promising results in the past decade. In line with these studies, we present MDFF_NM, a stochastic hybrid flexible fitting algorithm combining Normal Mode Analysis (NMA) and simulation-based flexible fitting. Initial tests reveal that, besides accelerating the fitting process, MDFF_NM increases the diversity of fitting routes leading to the target, uncovering ensembles of conformations in closer agreement with experimental data. The potential integration of MDFF_NM with other existing methods and integrative modeling pipelines is also discussed.

Deciphering the mechanism of HM43239 inhibiting the mutant F691L resistant to gilteritinib in FMS-like tyrosine kinase 3

FMS-like tyrosine kinase (FLT3) has become the legitimate molecular therapeutic target for acute myeloid leukemia therapy. Though FLT3 inhibitors have impact on disease progression, drug resistance induced by secondary point mutations is the primary mechanism and urgent to overcome. Herein, we sought to decipher the mechanism of HM43239 inhibiting the mutant F691L resistant to gilteritinib in FLT3. A series of molecular modeling studies, including molecular dynamics (MD) simulation, dynamic cross-correlation (DCC) analysis, binding free energy (MM-GBSA) and docking study were explored to elucidate the differential tolerance mechanisms of two inhibitors to the same mutant. The F691L mutation had relatively larger effect on gilteritinib than HM43239, which showed as the changed and fixed conformation, respectively. These observations rationalized that the binding affinity of gilteritinib decreased more than that of HM43239 in the F691L mutant.Communicated by Ramaswamy H. Sarma.

Nanomechanical collective vibration of SARS-CoV-2 spike proteins

The development of effective therapeutics against COVID-19 requires a thorough understanding of the receptor recognition mechanism of the SARS-CoV-2 spike (S) protein. Here the multidomain collective dynamics on the trimer of the spike protein has been analyzed using normal mode analysis (NMA). A common nanomechanical profile was identified in the spike proteins of SARS-CoV-2 and its variants. The profile involves collective vibrations of the receptor-binding domain (RBD) and the N-terminal domain (NTD), which may mediate the physical interaction process. Quantitative analysis of the collective modes suggests a nanomechanical property involving large-scale conformational changes, which explains the difference in receptor binding affinity among different variants. These results support the use of intrinsic global dynamics as a valuable perspective for studying the allosteric and functional mechanisms of the S protein. This approach also provides a low-cost theoretical toolkit for screening potential pathogenic mutations and drug targets.

Computational evaluation and benchmark study of 342 crystallographic holo-structures of SARS-CoV-2 Mpro enzyme

The coronavirus disease 19 pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has led to a global health crisis with millions of confirmed cases and related deaths. The main protease (Mpro) of SARS-CoV-2 is crucial for viral replication and presents an attractive target for drug development. Despite the approval of some drugs, the search for effective treatments continues. In this study, we systematically evaluated 342 holo-crystal structures of Mpro to identify optimal conformations for structure-based virtual screening (SBVS). Our analysis revealed limited structural flexibility among the structures. Three docking programs, AutoDock Vina, rDock, and Glide were employed to assess the efficiency of virtual screening, revealing diverse performances across selected Mpro structures. We found that the structures 5RHE, 7DDC, and 7DPU (PDB Ids) consistently displayed the lowest EF, AUC, and BEDROCK scores. Furthermore, these structures demonstrated the worst pose prediction results in all docking programs. Two structural differences contribute to variations in docking performance: the absence of the S1 subsite in 7DDC and 7DPU, and the presence of a subpocket in the S2 subsite of 7DDC, 7DPU, and 5RHE. These findings underscore the importance of selecting appropriate Mpro conformations for SBVS, providing valuable insights for advancing drug discovery efforts.

The Legionella collagen-like protein employs a distinct binding mechanism for the recognition of host glycosaminoglycans

Bacterial adhesion is a fundamental process which enables colonisation of niche environments and is key for infection. However, in Legionella pneumophila, the causative agent of Legionnaires' disease, these processes are not well understood. The Legionella collagen-like protein (Lcl) is an extracellular peripheral membrane protein that recognises sulphated glycosaminoglycans on the surface of eukaryotic cells, but also stimulates bacterial aggregation in response to divalent cations. Here we report the crystal structure of the Lcl C-terminal domain (Lcl-CTD) and present a model for intact Lcl. Our data reveal that Lcl-CTD forms an unusual trimer arrangement with a positively charged external surface and negatively charged solvent exposed internal cavity. Through molecular dynamics simulations, we show how the glycosaminoglycan chondroitin-4-sulphate associates with the Lcl-CTD surface via distinct binding modes. Our findings show that Lcl homologs are present across both the Pseudomonadota and Fibrobacterota-Chlorobiota-Bacteroidota phyla and suggest that Lcl may represent a versatile carbohydrate-binding mechanism.

A slow but steady nanoLuc: R162A mutation results in a decreased, but stable, nanoLuc activity

NanoLuc (NLuc) luciferase has found extensive application in designing a range of biological assays, including gene expression analysis, protein-protein interaction, and protein conformational changes due to its enhanced brightness and small size. However, questions related to its mechanism of interaction with the substrate, furimazine, as well as bioluminescence activity remain elusive. Here, we combined molecular dynamics (MD) simulation and mutational analysis to show that the R162A mutation results in a decreased but stable bioluminescence activity of NLuc in living cells and in vitro. Specifically, we performed multiple, all-atom, explicit solvent MD simulations of the apo and furimazine-docked (holo) NLuc structures revealing differential dynamics of the protein in the absence and presence of the ligand. Further, analysis of trajectories for hydrogen bonds (H-bonds) formed between NLuc and furimazine revealed substantial H-bond interaction between R162 and Q32 residues. Mutation of the two residues in NLuc revealed a decreased but stable activity of the R162A, but not Q32A, mutant NLuc in live cell and in vitro assays performed using lysates prepared from cells expressing the proteins and with the furimazine substrate. In addition to highlighting the role of the R162 residue in NLuc activity, we believe that the mutant NLuc will find wide application in designing in vitro assays requiring extended monitoring of NLuc bioluminescence activity. SIGNIFICANCE: Bioluminescence has been extensively utilized in developing a variety of biological and biomedical assays. In this regard, engineering of brighter bioluminescent proteins, i.e. luciferases, has played a significant role. This is acutely exemplified by the engineering of the NLuc luciferase, which is small in size and displays much enhanced bioluminescence and thermal stability compared to previously available luciferases. While enhanced bioluminescent activity is desirable in a multitude of biological and biomedical assays, it would also be useful to develop variants of the protein that display a prolonged bioluminescence activity. This is specifically relevant in designing assays that require bioluminescence for extended periods, such as in the case of biosensors designed for monitoring slow enzymatic or cellular signaling reactions, without necessitating multiple rounds of luciferase substrate addition or any specialized reagents that result in increased assay costs. In the current manuscript, we report a mutant NLuc that possesses a stable and prolonged bioluminescence activity, albeit lower than the wild-type NLuc, and envisage a wider application of the mutant NLuc in designing biosensors for monitoring slower biological and biomedical events.

Identification of a potential DNA methyltransferase (DNMT) inhibitor

DNA methyltransferases (DNMTs) play an important role in the epigenetic regulation of gene expression through the methylation of DNA. Since hypermethylation and consequent suppression of tumor suppressor genes are associated with cancer development and progression, DNA hypomethylating agents (HMAs) such as DNMT inhibitors have been proposed for cancer therapy. Two nucleoside analogues approved for the treatment of hematological cancers, decitabine and azacytidine, have poor pharmacokinetic properties, and hence there is a critical need for identifying novel HMAs. Virtual screening of a library of ∼40,000 compounds from the ZINC database, followed by molecular docking of 4,000 compounds having potential druggable properties with DNMT1, DNMT3A and DNMT3B were performed. One unique inhibitor (ZINC167686681) was identified that successfully passed through the Lipinski Rule of 5, geometry constraints as well as ADME/Tox filters and having strong binding energy for DNMTs. Further, molecular dynamics simulations of the docked complexes showed detailed structural features critical for its binding with the DNMTs and the stability of their interaction. Our study identified a compound with potential druggable properties and predicted to bind and inhibit DNMTs. Further investigations involving cellular and animal models of ZINC167686681 will help in potentially taking it into clinical trials for the treatment of cancers.Communicated by Ramaswamy H. Sarma.

Identification of FTY720 and COH29 as novel topoisomerase I catalytic inhibitors by experimental and computational studies

The development of novel topoisomerase I (TOP1) inhibitors is crucial for overcoming the drawbacks and limitations of current TOP1 poisons. Here, we identified two potential TOP1 inhibitors, namely, FTY720 (a sphingosine 1-phosphate antagonist) and COH29 (a ribonucleotide reductase inhibitor), through experimental screening of known active compounds. Biological experiments verified that FTY720 and COH29 were nonintercalative TOP1 catalytic inhibitors that did not induce the formation of DNA-TOP1 covalent complexes. Molecular docking revealed that FTY720 and COH29 interacted favorably with TOP1. Molecular dynamics simulations revealed that FTY720 and COH29 could affect the catalytic domain of TOP1, thus resulting in altered DNA-binding cavity size. The alanine scanning and interaction entropy identified Arg536 as a hotspot residue. In addition, the bioinformatics analysis predicted that FTY720 and COH29 could be effective in treating malignant breast tumors. Biological experiments verified their antitumor activities using MCF-7 breast cancer cells. Their combinatory effects with TOP1 poisons were also investigated. Further, FTY720 and COH29 were found to cause less DNA damage compared with TOP1 poisons. The findings provide reliable lead compounds for the development of novel TOP1 catalytic inhibitors and offer new insights into the potential clinical applications of FTY720 and COH29 in targeting TOP1.

Exploring optimal drug targets through subtractive proteomics analysis and pangenomic insights for tailored drug design in tuberculosis

Tuberculosis (TB), caused by Mycobacterium tuberculosis, ranks among the top causes of global human mortality, as reported by the World Health Organization's 2022 TB report. The prevalence of M. tuberculosis strains that are multiple and extensive-drug resistant represents a significant barrier to TB eradication. Fortunately, having many completely sequenced M. tuberculosis genomes available has made it possible to investigate the species pangenome, conduct a pan-phylogenetic investigation, and find potential new drug targets. The 442 complete genome dataset was used to estimate the pangenome of M. tuberculosis. This study involved phylogenomic classification and in-depth analyses. Sequential filters were applied to the conserved core genome containing 2754 proteins. These filters assessed non-human homology, virulence, essentiality, physiochemical properties, and pathway analysis. Through these intensive filtering approaches, promising broad-spectrum therapeutic targets were identified. These targets were docked with FDA-approved compounds readily available on the ZINC database. Selected highly ranked ligands with inhibitory potential include dihydroergotamine and abiraterone acetate. The effectiveness of the ligands has been supported by molecular dynamics simulation of the ligand-protein complexes, instilling optimism that the identified lead compounds may serve as a robust basis for the development of safe and efficient drugs for TB treatment, subject to further lead optimization and subsequent experimental validation.

The roles of ABCB1/P-glycoprotein drug transporters in regulating gut microbes and inflammation: insights from animal models, old and new

Commensal enteric bacteria have evolved systems that enable growth in the ecologic niche of the host gastrointestinal tract. Animals evolved parallel mechanisms to survive the constant exposure to bacteria and their metabolic by-products. We propose that drug transporters encompass a crucial system to managing the gut microbiome. Drug transporters are present in the apical surface of gut epithelia. They detoxify cells from small molecules and toxins (xenobiotics) in the lumen. Here, we review what is known about commensal structure in the absence of the transporter ABCB1/P-glycoprotein in mammalian models. Knockout or low-activity alleles of ABCB1 lead to dysbiosis, Crohn's disease and ulcerative colitis in mammals. However, the exact function of ABCB1 in these contexts remain unclear. We highlight emerging models-the zebrafish Danio rerio and sea urchin Lytechinus pictus-that are poised to help dissect the fundamental mechanisms of ATP-binding cassette (ABC) transporters in the tolerance of commensal and pathogenic communities in the gut. We and others hypothesize that ABCB1 plays a direct role in exporting inflammatory bacterial products from host epithelia. Interdisciplinary work in this research area will lend novel insight to the transporter-mediated pathways that impact microbiome community structure and accelerate the pathogenesis of inflammatory bowel disease when perturbed. This article is part of the theme issue 'Sculpting the microbiome: how host factors determine and respond to microbial colonization'.

The α C helix is a central regulator of PKR activation

Protein kinase R (PKR) functions in the eukaryotic innate immune system as a first-line defense against viral infections. PKR binds viral dsRNA, leading to autophosphorylation and activation. In its active state, PKR can phosphorylate its primary substrate, eIF2 α , which blocks initiation of translation in the infected cell. It has been established that PKR activation occurs when the kinase domain dimerizes in a back-to-back configuration. However, the mechanism by which dimerization leads to enzymatic activation is not fully understood. Here, we investigate the structural mechanistic basis and energy landscape for PKR activation, with a focus on the α C helix - a kinase activation and signal integration hub - using all-atom equilibrium and enhanced sampling molecular dynamics simulations. By employing window-exchange umbrella sampling, we compute free energy profiles of activation which show that back-to-back dimerization stabilizes a catalytically competent conformation of PKR. Key hydrophobic residues in the homodimer interface contribute to stabilization of the α C helix in an active conformation and the position of its glutamate residue. Using linear mutual information analysis, we analyze allosteric communication connecting the protomers' N-lobes and the α C helix dimer interface with the α C helix.

Prebiotic levan type fructan from Bacillus subtilis PR-C18 as a potent antibiofilm agent: Structural elucidation and in silico analysis

The global demand for therapeutic prebiotics persuades the quest for novel exopolysaccharides that can retard the growth of pathobionts and healthcare-associated pathogens. In this regard, an exopolysaccharide (3.69 mg/mL) producing strain showing prebiotic and antibiofilm activity was isolated from indigenous pineapple pomace of Tripura and identified as Bacillus subtilis PR-C18. Zymogram analysis revealed EPS PR-C18 was synthesized by levansucrase (∼57 kDa) with a maximal activity of 4.62 U/mg. Chromatography techniques, FTIR, and NMR spectral data revealed the homopolymeric nature of purified EPS with a molecular weight of 3.40 × 104 Da. SEM and rheological study unveiled its microporous structure and shear-thinning effect. Furthermore, EPS PR-C18 showed remarkable emulsification, flocculation, water retention, water solubilization, and antioxidant activity. DSC-TGA data demonstrated its high thermostability and cytotoxicity analysis verified its nontoxic biocompatible nature. In addition, the antibiofilm activity of EPS PR-C18 was validated using molecular docking, molecular simulation, MM-GBSA and PCA studies, which exhibited its strong binding affinity (-20.79 kcal/moL) with PelD, a virulence factor from Pseudomonas aeruginosa. Together, these findings support the future exploitation of EPS PR-C18 as an additive or adjuvant in food and pharmaceutical sectors.

Dissecting the Allosteric Fine-Tuning of Enzyme Catalysis

Fully understanding the mechanism of allosteric regulation in biomolecules requires separating and examining all of the involved factors. In enzyme catalysis, allosteric effector binding shifts the structure and dynamics of the active site, leading to modified energetic (e.g., energy barrier) and dynamical (e.g., diffusion coefficient) factors underlying the catalyzed reaction rate. Such modifications can be subtle and dependent on the type of allosteric effector, representing a fine-tuning of protein function. The microscopic description of allosteric regulation at the level of function-dictating factors has prospective applications in fundamental and pharmaceutical sciences, which is, however, largely missing so far. Here, we characterize the allosteric fine-tuning of enzyme catalysis, using human Pin1 as an example, by performing more than half-millisecond all-atom molecular dynamics simulations. Changes of reaction kinetics and the dictating factors, including the free energy surface along the reaction coordinate and the diffusion coefficient of the reaction dynamics, under various enzyme and allosteric effector binding conditions are examined. Our results suggest equal importance of the energetic and dynamical factors, both of which can be modulated allosterically, and the combined effect determines the final allosteric output. We also reveal the potential dynamic basis for allosteric modulation using an advanced statistical technique to detect function-related conformational dynamics. Methods developed in this work can be applied to other allosteric systems.

In Silico Insights into the Arsenic Binding Mechanism Deploying Application of Computational Biology-Based Toolsets

An assortment of environmental matrices includes arsenic (As) in its different oxidation states, which is often linked to concerns that pose a threat to public health worldwide. The current difficulty lies in addressing toxicological concerns and achieving sustained detoxification of As. Multiple conventional degradation methods are accessible; however, they are indeed labor-intensive, expensive, and reliant on prolonged laboratory evaluations. Molecular interaction and atomic level degradation mechanisms for enzyme-As exploration are, however, underexplored in those approaches. A feasible approach in this case for tackling this accompanying concern of As might be to cope with undertaking multivalent computational methodologies and tools. This work aimed to provide molecular-level insight into the enzyme-aided As degradation mechanism. AutoDock Vina, CABS-flex 2.0, and Desmond high-performance molecular dynamics simulation (MDS) were utilized in the current investigation to simulate multivalent molecular processes on two protein sets: arsenate reductase (ArsC) and laccase (LAC) corresponding arsenate (ART) and arsenite (AST), which served as model ligands to comprehend binding, conformational, and energy attributes. The structural configurations of both proteins exhibited variability in flexibility and structure framework within the range of 3.5-4.5 Å. The LAC-ART complex exhibited the lowest calculated binding affinity, measuring -5.82 ± 0.01 kcal/mol. Meanwhile, active site residues ILE-200 and HIS-206 were demonstrated to engage in H-bonding with the ART ligand. In contrast to ArsC, the ligand binding affinity of this bound complex was considerably greater. Additional validation of docked complexes was carried out by deploying Desmond MDS of 100 ns to capture protein and ligand conformation behavior. The system achieved stability during the 100 ns simulation run, as confirmed by the average P-L RMSD, which was ∼1 Å. As a preliminary test of the enzyme's ability to catalyze As species, corresponding computational insights might be advantageous for bridging gaps and regulatory consideration.

A Network Pharmacology and Molecular Dynamics Simulation-Based Study of Qing Run Hua Jie Decoction in Interstitial Pneumonia Treatment

Objective

This study is dedicated to revealing the potential mechanism of Qin Run Hua Jie (QRHJ) decoction in Interstitial pneumonia (IP) treatment.

Methods

The TCMSP database predicted the chemical components and targets of QRHJ decoction, and the IP-related genes were from the Genecards database. Cytoscape software was used to establish the interaction network. R package clusterProfiler was utilized for Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. The molecular docking analysis of target proteins and the corresponding active pharmaceutical ingredients in the core position of the interaction network was conducted. Then, molecular dynamics (MD) simulations of a potential active substance and its key targets were performed. The binding efficiency of EGFR and luteolin, HIF1A and diosgenin was detected by cellular thermal shift assay (CETSA), and protein expression was measured by Western blot. CCK-8 was used to detect cell activity.

Results

A total of 153 active ingredients, 127 targets and 362 IP-related genes were obtained. KEGG enrichment analysis identified IP-related signaling pathways including HIF-1 signaling pathway and TNF signaling pathway. The two key components luteolin and diosgenin stably bound to the key targets EGFR and HIF1A. Cell experiments further showed that EGFR and luteolin, HIF1A and diosgenin bound to exert anti-fibrotic effects.

Conclusion

As an active ingredient of QRHJ decoction, luteolin and diosgenin may exert therapeutic effect on IP through binding to the key target EGFR and HIF1A. This work initially revealed the key molecular mechanism of QRHJ decoction in IP treatment and offered theoretical evidence.

Identification and mechanistic exploration of structural and conformational dynamics of NF-kB inhibitors: rationale insights from in silico and in vitro studies

Increased expression of target genes that code for proinflammatory chemical mediators results from a series of intracellular cascades triggered by activation of dysregulated NF-κB signaling pathway. Dysfunctional NF-kB signaling amplifies and perpetuates autoimmune responses in inflammatory diseases, including psoriasis. This study aimed to identify therapeutically relevant NF-kB inhibitors and elucidate the mechanistic aspects behind NF-kB inhibition. After virtual screening and molecular docking, five hit NF-kB inhibitors opted, and their therapeutic efficacy was examined using cell-based assays in TNF-α stimulated human keratinocyte cells. To investigate the conformational changes of target protein and inhibitor-protein interaction mechanisms, molecular dynamics (MD) simulations, binding free energy calculations together with principal component (PC) analysis, dynamics cross-correlation matrix analysis (DCCM), free energy landscape (FEL) analysis and quantum mechanical calculations were carried out. Among identified NF-kB inhibitors, myricetin and hesperidin significantly scavenged intracellular ROS and inhibited NF-kB activation. Analysis of the MD simulation trajectories of ligand-protein complexes revealed that myricetin and hesperidin formed energetically stabilized complexes with the target protein and were able to lock NF-kB in a closed conformation. Myricetin and hesperidin binding to the target protein significantly impacted conformational changes and internal dynamics of amino acid residues in protein domains. Tyr57, Glu60, Lys144 and Asp239 residues majorly contributed to locking the NF-kB in a closed conformation. The combinatorial approach employing in silico tools integrated with cell-based approaches substantiated the binding mechanism and NF-kB active site inhibition by the lead molecule myricetin, which can be explored as a viable antipsoriatic drug candidate associated with dysregulated NF-kB.Communicated by Ramaswamy H. Sarma.

Structure-oriented engineering of nitrile hydratase: Reshaping of substrate access tunnel and binding pocket for efficient synthesis of cinnamamide

Cinnamamide and its derivatives are the most common and important building blocks widely present in natural products. Currently, nitrile hydratase (NHase, EC 4.2.1.84) has been widely used in large-scale industrial production of nicotinamide and acrylamide, while its catalytic activity is extremely low or inactive for bulky nitrile substrates such as cinnamonitrile. Therefore, beneficial variant βF37P/L48P/F51N were obtained from PtNHase of Pseudonocardia thermophila JCM3095 by reshaping of substrate access tunnel and binding pocket, which exhibited 14.88-fold improved catalytic efficiency compared to the wild-type PtNHase. Structure analysis, molecular dynamics simulations and dynamical cross-correlation matrix (DCCM) analysis revealed that the introduced mutations enlarged the substrate access tunnel and binding pocket, enhanced overall anti-correlated movements of enzymes, which would promote product release during the dynamic process of catalysis. In a hydration process, the complete conversion of 5 mM cinnamonitrile was achieved by βF37P/L48P/F51N in a 50 mL reaction, with cinnamamide yield of almost 100 % and productivity of 0.736 g L-1 h-1. The study demonstrates the co-evolution of substrate access tunnel and binding pocket is an effective strategy, and provides a valuable reference for future research. Furthermore, NHases have huge potential for catalyzing bulky nitriles to form corresponding amides in large-scale industrial production.

Unveiling local and global conformational changes and allosteric communications in SOD1 systems using molecular dynamics simulation and network analyses

BACKGROUND

Amyotrophic lateral sclerosis (ALS) is a serious neurodegenerative disorder affecting nerve cells in the brain and spinal cord that is caused by mutations in the superoxide dismutase 1 (SOD1) enzyme. ALS-related mutations cause misfolding, dimerisation instability, and increased formation of aggregates. The underlying allosteric mechanisms, however, remain obscure as far as details of their fundamental atomistic structure are concerned. Hence, this gap in knowledge limits the development of novel SOD1 inhibitors and the understanding of how disease-associated mutations in distal sites affect enzyme activity.

METHODS

We combined microsecond-scale based unbiased molecular dynamics (MD) simulation with network analysis to elucidate the local and global conformational changes and allosteric communications in SOD1 Apo (unmetallated form), Holo, Apo_CallA (mutant and unmetallated form), and Holo_CallA (mutant form) systems. To identify hotspot residues involved in SOD1 signalling and allosteric communications, we performed network centrality, community network, and path analyses.

RESULTS

Structural analyses showed that unmetallated SOD1 systems and cysteine mutations displayed large structural variations in the catalytic sites, affecting structural stability. Inter- and intra H-bond analyses identified several important residues crucial for maintaining interfacial stability, structural stability, and enzyme catalysis. Dynamic motion analysis demonstrated more balanced atomic displacement and highly correlated motions in the Holo system. The rationale for structural disparity observed in the disulfide bond formation and R143 configuration in Apo and Holo systems were elucidated using distance and dihedral probability distribution analyses.

CONCLUSION

Our study highlights the efficiency of combining extensive MD simulations with network analyses to unravel the features of protein allostery.

Exploration of human pancreatic alpha-amylase inhibitors from Physalis peruviana for the treatment of type 2 diabetes

Type 2 Diabetes (T2D), a chronic metabolic disorder characterized by persistent hyperglycemia, accounts for ∼90% of all types of diabetes. Pancreatic α-amylase is a potential drug target for preventing postprandial hyperglycemia and inhibiting T2D in humans. Although many synthetic drugs have been identified against pancreatic α-amylase, however, reported several side effects, and plant-derived natural products are less explored against T2D. This study tested 34 flavonoids derived from the plant Physalis peruviana against the human pancreatic α-amylase (HPA) using in silico computational approaches such as molecular docking and molecular dynamics simulation approaches. Schrödinger, a drug discovery package with modules applicable for molecular docking, protein-ligand interaction analysis, molecular dynamics, post-dynamics simulation, and binding free energy calculation, was employed for all computational studies. Four flavonoids, namely, Chlorogenic acid, Withaperuvin F, Withaperuvin H, and Rutin, were picked based on their docking score ranging between -7.03 kcal/mol and -11.35 kcal/mol compared to the docking score -7.3 kcal/mol of reference ligand, i.e. Myricetin. The molecular dynamics analysis suggested that all flavonoids showed considerable stability within the protein's catalytic pocket, except chlorogenic acid, which showed high deviation during the last 15 ns. However, the interactions observed in initial docking and extracted from the simulation trajectory involved > 90% identical residues, indicating the affinity and stability of the docked flavonoids with the protein. Therefore, all four compounds identified in this study are proposed as promising antidiabetic candidates and should be further considered for their in vitro and in vivo validation.Communicated by Ramaswamy H. Sarma.

The Legionella collagen-like protein employs a unique binding mechanism for the recognition of host glycosaminoglycans

Bacterial adhesion is a fundamental process which enables colonisation of niche environments and is key for infection. However, in Legionella pneumophila, the causative agent of Legionnaires' disease, these processes are not well understood. The Legionella collagen-like protein (Lcl) is an extracellular peripheral membrane protein that recognises sulphated glycosaminoglycans (GAGs) on the surface of eukaryotic cells, but also stimulates bacterial aggregation in response to divalent cations. Here we report the crystal structure of the Lcl C-terminal domain (Lcl-CTD) and present a model for intact Lcl. Our data reveal that Lcl-CTD forms an unusual dynamic trimer arrangement with a positively charged external surface and a negatively charged solvent exposed internal cavity. Through Molecular Dynamics (MD) simulations, we show how the GAG chondroitin-4-sulphate associates with the Lcl-CTD surface via unique binding modes. Our findings show that Lcl homologs are present across both the Pseudomonadota and Fibrobacterota-Chlorobiota-Bacteroidota phyla and suggest that Lcl may represent a versatile carbohydrate binding mechanism.

Illuminating the structural basis of human neurokinin 1 receptor (NK1R) antagonism through classical all-atoms molecular dynamics simulations

Advances in structural biology have bestowed insights into the pleiotropic effects of neurokinin 1 receptors (NK1R) in diverse patho-physiological processes, thereby highlighting the potential therapeutic value of antagonists directed against NK1R. Herein, we investigate the mode of antagonist recognition to discern the obscure atomic facets germane for the function and molecular determinants of NK1R. To commence discernment of potent antagonists and the conformational changes in NK1R, induced upon antagonist binding, state-of-the-art classical all-atoms molecular dynamics (MD) simulations in lipid mimetic bilayers have been utilized. MD simulations of structural ensembles reveals the involvement of TM5 and TM6 in tight anchoring of antagonists through a network of interhelical hydrogen-bonds, while, the extracellular loop 2 (ECL2) governs the overall size and nature of the pocket, thereby modulating NK1R. Consistent comparison between experiments and MD simulation results discerns the predominant role of TM3, TM4, and TM6 in lipid-NK1R interaction. Correlation between hydrophobic index and helicity of TM domains elucidates their importance in maintaining the structural stability in addition to regulating NK1R antagonism. Taken together, we anticipate that our computational study marks a comprehensive structural basis of NK1R antagonism in lipid bilayers, which may facilitate designing of new therapeutics against associated diseases targeting human neurokinin receptors.

Identifying therapeutic antibacterial peptides against Vibrio cholerae to inhibit the function of Na(+)-translocating NADH-quinone reductase

Vibrio cholerae is the bacteria responsible for cholera, which is a significant threat to many nations. Curing and treating this infection requires identification of the critical protein and development of a drug to inhibit its function. In this context, Na(+)-translocating NADH-quinone reductase was considered a potential therapeutic target. A library of antibacterial peptides with residue lengths of 50 was screened using a docking method, and the five most potent peptides were selected on the basis of a weighted score derived from solvent accessible surface area and docking score. To investigate the stability of the protein-peptide complex, a 100-ns molecular dynamics simulation was performed. These peptides targeted the native dimeric binding interface of Na(+)-transporting NADH-quinone reductase. This study evaluated the binding affinity and conformational stability of these peptides with the protein using different post-simulation metrics. A peptide, CCL28, exhibited steady RMSD characteristics; nonetheless, it modified the docked conformation but stabilized in the new conformation. This peptide also demonstrated the best performance in addressing the protein's native binding interface. It demonstrated a binding free energy of -120 kcal/mol with the protein. Principal component analysis (PCA) revealed that the first PC had the lowest conformational variation and the greatest coverage. Eventually, these peptides were also evaluated using steered molecular dynamics, and it was discovered that CCL28 had a greater maximum force than the other five peptides, at 1139.08 kJ/mol/nm. Targeting the native binding interface, we present a CCL28 peptide with a strong potential to block the biological activity of Vibrio cholerae's Na(+)-translocating NADH-quinone reductase.Communicated by Ramaswamy H. Sarma.

Normal mode analysis and comparative study of intrinsic dynamics of alcohol oxidase enzymes from GMC protein family

Glucose-Methanol-Choline (GMC) family enzymes are very important in catalyzing the oxidation of a wide range of structurally diverse substrates. Enzymes that constitute the GMC family, share a common tertiary fold but < 25% sequence identity. Cofactor FAD, FAD binding signature motif, and similar structural scaffold of the active site are common features of oxidoreductase enzymes of the GMC family. Protein functionality mainly depends on protein three-dimensional structures and dynamics. In this study, we used the normal mode analysis method to search the intrinsic dynamics of GMC family enzymes. We have explored the dynamical behavior of enzymes with unique substrate catabolism and active site characteristics from different classes of the GMC family. Analysis of individual enzymes and comparative ensemble analysis of enzymes from different classes has shown conserved dynamic motion at FAD binding sites. The present study revealed that GMC enzymes share a strong dynamic similarity (Bhattacharyya coefficient >90% and root mean squared inner product >52%) despite low sequence identity across the GMC family enzymes. The study predicts that local deformation energy between atoms of the enzyme may be responsible for the catalysis of different substrates. This study may help that intrinsic dynamics can be used to make meaningful classifications of proteins or enzymes from different organisms.Communicated by Ramaswamy H. Sarma.

Machine learning-based modulation of Ca2+-binding affinity in EF-hand proteins and comparative structural insights into site-specific cooperative binding

Ca2+-binding proteins are present in almost all living organisms and different types display different levels of binding affinities for the cation. Here, we report two new scoring schemes enabling the user to estimate and manipulate the calcium binding affinities in EF hand containing proteins. To validate this, we designed a unique EF-hand loop capable of binding calcium with high affinity by altering five residues. The N-terminal domain of Entamoeba histolytica calcium-binding protein1 (NtEhCaBP1) is used for site-directed mutagenesis to incorporate the designed loop sequence into the second EF-hand motif of this protein, referred as Nt-EhCaBP1-EF2 mutant. The binding isotherms calculated using ITC calorimetry showed that Nt-EhCaBP1-EF2 mutant site binds Ca2+ with higher affinity than Wt-Nt-EhCaBP1, by ∼600 times. The crystal structure of the mutant displayed more compact Ca2+-coordination spheres in both of its EF loops than the structure of the wildtype protein. The compact coordination sphere of EF-2 causes the bend in the helix-3, which leads to the formation of unexpected hexamer of NtEhCaBP1-EF2 mutant structure. Further dynamic correlation analysis revealed that the mutation in the second EF loop changed the entire residue network of the monomer, resulting in stronger coordination of Ca2+ even in another EF-hand loop.

Sequence-Dependent Properties of the RNA Duplex

Sequence-dependent properties of the DNA duplex have been accurately described using extensive molecular dynamics simulations. The RNA duplex meanwhile─which is typically represented as a sequence-averaged rigid rod─does not benefit from having equivalent molecular dynamics simulations. In this paper, we present a massive simulation effort using a set of ABC-optimized duplexes from which we derived tetramer-resolution properties of the RNA duplex and a simple mesoscopic model that can represent elastic properties of long RNA duplexes. Despite the extreme chemical similarity between DNA and RNA, the local and global elastic properties of the duplexes are very different. DNA duplexes show a complex and nonelastic pattern of flexibility, for instance, while RNA duplexes behave as an elastic system whose deformations can be represented by simple harmonic potentials. In RNA duplexes (RNA2), not only are intra- and interbase pair parameters (equilibrium and mechanical) different from those in the equivalent DNA duplex sequences (DNA2) but the correlations between movements also differ. Simple statements on the relative flexibility or stability of both polymers are meaningless and should be substituted by a more detailed description depending on the sequence and the type of deformation considered.

Cheminformatics Strategies Unlock Marburg Virus VP35 Inhibitors from Natural Compound Library

The Ebola virus and its close relative, the Marburg virus, both belong to the family Filoviridae and are highly hazardous and contagious viruses. With a mortality rate ranging from 23% to 90%, depending on the specific outbreak, the development of effective antiviral interventions is crucial for reducing fatalities and mitigating the impact of Marburg virus outbreaks. In this investigation, a virtual screening approach was employed to evaluate 2042 natural compounds for their potential interactions with the VP35 protein of the Marburg virus. Average and worst binding energies were calculated for all 20 poses, and compounds that exhibited binding energies <-6 kcal/mol in both criteria were selected for further analysis. Based on binding energies, only six compounds (Estradiol benzoate, INVEGA (paliperidone), Isosilybin, Protopanaxadiol, Permethrin, and Bufalin) were selected for subsequent investigations, focusing on interaction analysis. Among these selected compounds, Estradiol benzoate, INVEGA (paliperidone), and Isosilybin showed strong hydrogen bonds, while the others did not. In this study, the compounds Myricetin, Isosilybin, and Estradiol benzoate were subjected to a molecular dynamics (MD) simulation and free binding energy calculation using MM/GBSA analysis. The reference component Myricetin served as a control. Estradiol benzoate exhibited the most stable and consistent root-mean-square deviation (RMSD) values, whereas Isosilybin showed significant fluctuations in RMSD. The compound Estradiol benzoate exhibited the lowest ΔG binding free energy (-22.89 kcal/mol), surpassing the control compound's binding energy (-9.29 kcal/mol). Overall, this investigation suggested that Estradiol benzoate possesses favorable binding free energies, indicating a potential inhibitory mechanism against the VP35 protein of the Marburg virus. The study proposes that these natural compounds could serve as a therapeutic option for preventing Marburg virus infection. However, experimental validation is required to further corroborate these findings.

ACKR3-arrestin2/3 complexes reveal molecular consequences of GRK-dependent barcoding

Atypical chemokine receptor 3 (ACKR3, also known as CXCR7) is a scavenger receptor that regulates extracellular levels of the chemokine CXCL12 to maintain responsiveness of its partner, the G protein-coupled receptor (GPCR), CXCR4. ACKR3 is notable because it does not couple to G proteins and instead is completely biased towards arrestins. Our previous studies revealed that GRK2 and GRK5 install distinct distributions of phosphates (or "barcodes") on the ACKR3 carboxy terminal tail, but how these unique barcodes drive different cellular outcomes is not understood. It is also not known if arrestin2 (Arr2) and 3 (Arr3) bind to these barcodes in distinct ways. Here we report cryo-electron microscopy structures of Arr2 and Arr3 in complex with ACKR3 phosphorylated by either GRK2 or GRK5. Unexpectedly, the finger loops of Arr2 and 3 directly insert into the detergent/membrane instead of the transmembrane core of ACKR3, in contrast to previously reported "core" GPCR-arrestin complexes. The distance between the phosphorylation barcode and the receptor transmembrane core regulates the interaction mode of arrestin, alternating between a tighter complex for GRK5 sites and heterogenous primarily "tail only" complexes for GRK2 sites. Arr2 and 3 bind at different angles relative to the core of ACKR3, likely due to differences in membrane/micelle anchoring at their C-edge loops. Our structural investigations were facilitated by Fab7, a novel Fab that binds both Arr2 and 3 in their activated states irrespective of receptor or phosphorylation status, rendering it a potentially useful tool to aid structure determination of any native GPCR-arrestin complex. The structures provide unprecedented insight into how different phosphorylation barcodes and arrestin isoforms can globally affect the configuration of receptor-arrestin complexes. These differences may promote unique downstream intracellular interactions and cellular responses. Our structures also suggest that the 100% bias of ACKR3 for arrestins is driven by the ability of arrestins, but not G proteins, to bind GRK-phosphorylated ACKR3 even when excluded from the receptor cytoplasmic binding pocket.

Mammalian ATG8 proteins maintain autophagosomal membrane integrity through ESCRTs

The canonical autophagy pathway in mammalian cells sequesters diverse cytoplasmic cargo within the double membrane autophagosomes that eventually convert into degradative compartments via fusion with endolysosomal intermediates. Here, we report that autophagosomal membranes show permeability in cells lacking principal ATG8 proteins (mATG8s) and are unable to mature into autolysosomes. Using a combination of methods including a novel in vitro assay to measure membrane sealing, we uncovered a previously unappreciated function of mATG8s to maintain autophagosomal membranes in a sealed state. The mATG8 proteins GABARAP and LC3A bind to key ESCRT-I components contributing, along with other ESCRTs, to the integrity and imperviousness of autophagic membranes. Autophagic organelles in cells lacking mATG8s are permeant, are arrested as amphisomes, and do not progress to functional autolysosomes. Thus, autophagosomal organelles need to be maintained in a sealed state in order to become lytic autolysosomes.

Investigating the Mechanism of Action of Anti-Dengue Compounds as Potential Binders of Zika Virus RNA-Dependent RNA Polymerase

The World Health Organization (WHO) has designated the Zika virus (ZIKV) as a significant risk to the general public's health. Currently, there are no vaccinations or medications available to treat or prevent infection with the Zika virus. Thus, it is urgently required to develop a highly efficient therapeutic molecule. In the presented study, a computationally intensive search was carried out to identify potent compounds that have the potential to bind and block the activity of ZIKV NS5 RNA-dependent RNA polymerase (RdRp). The anti-dengue chemical library was subjected to high-throughput virtual screening and MM/GBSA analysis in order to rate the potential candidates. The top three compounds were then chosen. According to the MM/GBSA analysis, compound 127042987 from the database had the highest binding affinity to the protein with a minimum binding free energy of -77.16 kcal/mole. Compound 127042987 had the most stable RMSD trend and the greatest number of hydrogen bond interactions when these chemical complexes were evaluated further under a 100 ns molecular dynamics simulation. Compound 127042987 displayed the best binding free energy (GBind) of -96.50 kcal/mol, surpassing the native ligand binding energy (-66.17 kcal/mole). Thereafter, an MM/GBSA binding free energy study was conducted to validate the stability of selected chemical complexes. Overall, this study illustrated that compound 127042987 showed preferred binding free energies, suggesting a possible inhibitory mechanism against ZIKV-RdRp. As per this study, it was proposed that compound 127042987 could be used as a therapeutic option to prevent Zika virus infection. These compounds need to be tested in experiments for further validation.

AI in drug discovery and its clinical relevance

The COVID-19 pandemic has emphasized the need for novel drug discovery process. However, the journey from conceptualizing a drug to its eventual implementation in clinical settings is a long, complex, and expensive process, with many potential points of failure. Over the past decade, a vast growth in medical information has coincided with advances in computational hardware (cloud computing, GPUs, and TPUs) and the rise of deep learning. Medical data generated from large molecular screening profiles, personal health or pathology records, and public health organizations could benefit from analysis by Artificial Intelligence (AI) approaches to speed up and prevent failures in the drug discovery pipeline. We present applications of AI at various stages of drug discovery pipelines, including the inherently computational approaches of de novo design and prediction of a drug's likely properties. Open-source databases and AI-based software tools that facilitate drug design are discussed along with their associated problems of molecule representation, data collection, complexity, labeling, and disparities among labels. How contemporary AI methods, such as graph neural networks, reinforcement learning, and generated models, along with structure-based methods, (i.e., molecular dynamics simulations and molecular docking) can contribute to drug discovery applications and analysis of drug responses is also explored. Finally, recent developments and investments in AI-based start-up companies for biotechnology, drug design and their current progress, hopes and promotions are discussed in this article.

Influence of pH on indole-dependent heterodimeric interactions between Anopheles gambiae odorant-binding proteins OBP1 and OBP4

Insect Odorant Binding Proteins (OBPs) constitute important components of their olfactory apparatus, as they are essential for odor recognition. OBPs undergo conformational changes upon pH change, altering their interactions with odorants. Moreover, they can form heterodimers with novel binding characteristics. Anopheles gambiae OBP1 and OBP4 were found capable of forming heterodimers possibly involved in the specific perception of the attractant indole. In order to understand how these OBPs interact in the presence of indole and to investigate the likelihood of a pH-dependent heterodimerization mechanism, the crystal structures of OBP4 at pH 4.6 and 8.5 were determined. Structural comparison to each other and with the OBP4-indole complex (3Q8I, pH 6.85) revealed a flexible N-terminus and conformational changes in the α4-loop-α5 region at acidic pH. Fluorescence competition assays showed a weak binding of indole to OBP4 that become further impaired at acidic pH. Additional Molecular Dynamic and Differential Scanning Calorimetry studies displayed that the influence of pH on OBP4 stability is significant compared to the modest effect of indole. Furthermore, OBP1-OBP4 heterodimeric models were generated at pH 4.5, 6.5, and 8.5, and compared concerning their interface energy and cross-correlated motions in the absence and presence of indole. The results indicate that the increase in pH may induce the stabilization of OBP4 by increasing its helicity, thereby enabling indole binding at neutral pH that further stabilizes the protein and possibly promotes the creation of a binding site for OBP1. A decrease in interface stability and loss of correlated motions upon transition to acidic pH may provoke the heterodimeric dissociation allowing indole release. Finally, we propose a potential OBP1-OBP4 heterodimer formation/disruption mechanism induced by pH change and indole binding.

Prion Protein Octarepeat Domain Forms Transient β-Sheet Structures upon Residue-Specific Binding to Cu(II) and Zn(II) Ions

Misfolding of the cellular prion protein (PrP^C) is associated with the development of fatal neurodegenerative diseases called transmissible spongiform encephalopathies (TSEs). Metal ions appear to play a crucial role in PrP^C misfolding. PrP^C is a combined Cu(II) and Zn(II) metal-binding protein, where the main metal-binding site is located in the octarepeat (OR) region. Thus, the biological function of PrP^C may involve the transport of divalent metal ions across membranes or buffering concentrations of divalent metal ions in the synaptic cleft. Recent studies have shown that an excess of Cu(II) ions can result in PrP^C instability, oligomerization, and/or neuroinflammation. Here, we have used biophysical methods to characterize Cu(II) and Zn(II) binding to the isolated OR region of PrP^C. Circular dichroism (CD) spectroscopy data suggest that the OR domain binds up to four Cu(II) ions or two Zn(II) ions. Binding of the first metal ion results in a structural transition from the polyproline II helix to the β-turn structure, while the binding of additional metal ions induces the formation of β-sheet structures. Fluorescence spectroscopy data indicate that the OR region can bind both Cu(II) and Zn(II) ions at neutral pH, but under acidic conditions, it binds only Cu(II) ions. Molecular dynamics simulations suggest that binding of either metal ion to the OR region results in the formation of β-hairpin structures. As the formation of β-sheet structures can be a first step toward amyloid formation, we propose that high concentrations of either Cu(II) or Zn(II) ions may have a pro-amyloid effect in TSE diseases.

Efficient production of hydroxytyrosol by directed evolution of HpaB in Escherichia coli

Hydroxytyrosol (HT) is an olive-derived phenolic phytochemical that has gained increasing commercial interest due to its natural antioxidant properties. It is widely used in the field of food supplement and medicine. It is reported that 4-hydroxyphenylacetate 3-hydroxylase (EcHpaB) and flavin reductase (EcHpaC) from E. coli BL21(DE3) can successfully express and catalyze the production of HT from tyrosol. In this study, the tyrosol production strain YMG5∗R as chassis cells, and a random mutant library of EcHpaB was established using error-prone PCR to improve the ability of EcHpaB to convert tyrosol to HT. Finally, a highly efficient HT synthetic mutant strainYMG5∗R-HpaB^TLEHC with high transformation efficiency was screened by directed evolution. The YMG5∗R-HpaB^TLEHC strain efficiently converted 50 mM tyrosol, with a yield of hydroxytyrosol reaching 48.2 mM (7.43 g/L) and a space-time yield reached 0.62 g/L·h. Overall, our study demonstrates the successful development of a highly efficient synthetic enzyme mutant for the production of HT, which has the potential to significantly improve the commercial viability of this natural antioxidant.

Conserved Conformational Dynamics Reveal a Key Dynamic Residue in the Gatekeeper Loop of Human Cyclophilins

Cyclophilins are ubiquitous human enzymes that catalyze peptidyl-prolyl cis-trans isomerization in protein substrates. Of the 17 unique isoforms, five closely related isoforms (CypA-E) are found in various environments and participate in diverse cellular processes, yet all have similar structures and the same core catalytic function. The question is what key residues are behind the conserved function of these enzymes. Here, conformational dynamics are compared across these isoforms to detect conserved dynamics essential for the catalytic activity of cyclophilins. A set of key dynamic residues, defined by the most dynamically conserved positions, are identified in the gatekeeper 2 region. The highly conserved glycine (Gly80) in this region is predicted to underlie the local flexibility, which is further tested by molecular dynamics simulations performed on mutants (G80A) of CypE and CypA. The mutation leads to decreased flexibility of CypE and CypA during substrate binding but increased flexibility during catalysis. Dynamical changes occur in the mutated region and a distal loop downstream of the mutation site in sequence. Examinations of the mutational effect on catalysis show that both mutated CypE and CypA exhibit shifted binding free energies of the substrate under distinct isomer conformations. The results suggest a loss of function in the mutated CypE and CypA. These catalytic changes by the mutation are likely independent of the substrate sequence, at least in CypA. Our work presents a method to identify function-related key residues in proteins.

Functional Protein Dynamics Directly from Sequences

The sequence correlations within a protein multiple sequence alignment are routinely being used to predict contacts within its structure, but here we point out that these data can also be used to predict a protein's dynamics directly. The elastic network protein dynamics models rely directly upon the contacts, and the normal modes of motion are obtained from the decomposition of the inverse of the contact map. To make the direct connection between sequence and dynamics, it is necessary to apply coarse-graining to the structure at the level of one point per amino acid, which has often been done, and protein coarse-grained dynamics from elastic network models has been highly successful, particularly in representing the large-scale motions of proteins that usually relate closely to their functions. The interesting implication of this is that it is not necessary to know the structure itself to obtain its dynamics and instead to use the sequence information directly to obtain the dynamics.

Non-active Site Residue in Loop L4 Alters Substrate Capture and Product Release in d-Arginine Dehydrogenase

Numerous studies demonstrate that enzymes undergo multiple conformational changes during catalysis. The malleability of enzymes forms the basis for allosteric regulation: residues located far from the active site can exert long-range dynamical effects on the active site residues to modulate catalysis. The structure of Pseudomonas aeruginosa d-arginine dehydrogenase (PaDADH) shows four loops (L1, L2, L3, and L4) that span the substrate and the FAD-binding domains. Loop L4 comprises residues 329-336, spanning over the flavin cofactor. The I335 residue on loop L4 is ∼10 Å away from the active site and ∼3.8 Å from N(1)-C(2)═O atoms of the flavin. In this study, we used molecular dynamics and biochemical techniques to investigate the effect of the mutation of I335 to histidine on the catalytic function of PaDADH. Molecular dynamics showed that the conformational dynamics of PaDADH are shifted to a more closed conformation in the I335H variant. In agreement with an enzyme that samples more in a closed conformation, the kinetic data of the I335H variant showed a 40-fold decrease in the rate constant of substrate association (k₁), a 340-fold reduction in the rate constant of substrate dissociation from the enzyme-substrate complex (k₂), and a 24-fold decrease in the rate constant of product release (k₅), compared to that of the wild-type. Surprisingly, the kinetic data are consistent with the mutation having a negligible effect on the reactivity of the flavin. Altogether, the data indicate that the residue at position 335 has a long-range dynamical effect on the catalytic function in PaDADH.

Revealing intrinsic changes of DNA induced by spore photoproduct lesion through computer simulation

In bacterial endospores, a cross-linked thymine dimer, 5-thyminyl-5,6-dihydrothymine, commonly referred to as the spore photoproduct (SP), is found as the dominant DNA photo lesion under UV radiation. During spore germination, SP is faithfully repaired by the spore photoproduct lyase (SPL) for normal DNA replication to resume. Despite this general mechanism, the exact way in which SP modifies the duplex DNA structure so that the damaged site can be recognized by SPL to initiate the repair process is still unclear. A previous X-ray crystallographic study, which used a reverse transcriptase as a DNA host template, captured a protein-bound duplex oligonucleotide containing two SP lesions; the study showed shortened hydrogen bonds between the AT base pairs involved in the lesions and widened minor grooves near the damaged sites. However, it remains to be determined whether the results accurately reflect the conformation of SP-containing DNA (SP-DNA) in its fully hydrated pre-repair form. To uncover the intrinsic changes in DNA conformation caused by SP lesions, we performed molecular dynamics (MD) simulations of SP-DNA duplexes in aqueous solution, using the nucleic acid portion of the previously determined crystal structure as a template. After MD relaxation, our simulated SP-DNAs showed weakened hydrogen bonds at the damaged sites compared to those in the undamaged DNA. Our analyses of the MD trajectories revealed a range of local and global structural distortions of DNA induced by SP. Specifically, the SP region displays a greater tendency to adopt an A-like-DNA conformation, and curvature analysis revealed an increase in the global bending compared to the canonical B-DNA. Although these SP-induced DNA conformational changes are relatively minor, they may provide a sufficient structural basis for SP to be recognized by SPL during the lesion repair process.

Structural evolution of an amphibian-specific globin: A computational evolutionary biochemistry approach

Studies on the globin family are continuously revealing insights into the mechanisms of gene and protein evolution. The rise of a new globin gene type in Pelobatoidea and Neobatrachia (Amphibia:Anura) from an α-globin precursor provides the opportunity to investigate the genetic and physical mechanisms underlying the origin of new protein structural and functional properties. This amphibian-specific globin (globin A/GbA) discovered in the heart of Rana catesbeiana is a monomer. As the ancestral oligomeric state of α-globins is a homodimer, we inferred that the ancestral state was lost somewhere in the GbA lineage. Here, we combined computational molecular evolution with structural bioinformatics to determine the extent to which the loss of the homodimeric state is pervasive in the GbA clade. We also characterized the loci of GbA genes in Bufo bufo. We found two GbA clades in Neobatrachia. One was deleted in Ranidae, but retained and expanded to yield a new globin cluster in Bufonidae species. Loss of the ancestral oligomeric state seems to be pervasive in the GbA clade. However, a taxonomic sampling that includes more Pelobatoidea, as well as early Neobatrachia, lineages would be necessary to determine the oligomeric state of the last common ancestor of all GbA. The evidence presented here points out a possible loss of oligomerization in Pelobatoidea GbA as a result of amino acid substitutions that weaken the homodimeric state. In contrast, the loss of oligomerization in both Neobatrachia GbA clades was linked to independent deletions that disrupted many packing contacts at the homodimer interface.

Molecular surveillance of Kelch-13 gene in Plasmodium falciparum field isolates from Mayurbhanj District, Odisha, India, and in silico artemisinin-Kelch-13 protein interaction study

The global malaria control and elimination program faces major threats due to the emergence and transmission of the anti-malarial drug-resistant strain of Plasmodium falciparum. Monitoring of artemisinin (ART) resistance marker Kelch-13 in the malaria-endemic region is essential in mitigating the disease's morbidity and mortality. The current study aimed to generate baseline information for further surveillance in the future. The current research was designed and conducted from July 2019 to June 2021 to monitor Pfkelch13 mutation at the molecular level in the eastern region of India. We also conducted an in silico study to understand the drug-protein interactions between ART and the protein crystal of PfKelch13 (KELCH) with PDB id:4ZGC. The kelch-13 gene was amplified by nested polymerase chain reaction (PCR) and sequenced through the Sanger sequencing method. Reference 3D7 clone (PF3D7_1343700) was used to align and probe all the sequences. The sequence analysis showed the absence of validated or associated mutation in the Kelch-13 propeller domain. The absence of natural selection in drug resistance was confirmed by the Tajima test. Further, in silico interaction studies between the drug ART and the Kelch propeller domain of P. falciparum were evaluated by structure predictions, molecular docking, molecular dynamics (MD) simulations, and estimations of binding free energies for the KELCH-ART complex. The results were compared with the apoprotein (KELCH-APO). The study confirmed the favorable binding of ART with the Kelch-13 propeller domain.

Evolutionary conservation of protein dynamics: insights from all-atom molecular dynamics simulations of 'peptidase' domain of Spt16

Protein function is encoded in its sequence, manifested in its three-dimensional structure, and facilitated by its dynamics. Studies have suggested that protein structures with higher sequence similarity could have more similar patterns of dynamics. However, such studies of protein dynamics within and across protein families typically rely on coarse-grained models, or approximate metrics like crystallographic B-factors. This study uses µs scale molecular dynamics (MD) simulations to explore the conservation of dynamics among homologs of ∼50 kDa N-terminal module of Spt16 (Spt16N). Spt16N from Saccharomyces cerevisiae (Sc-Spt16N) and three of its homologs with 30-40% sequence identities were available in the PDB. To make our data-set more comprehensive, the crystal structure of an additional homolog (62% sequence identity with Sc-Spt16N) was solved at 1.7 Å resolution. Cumulative MD simulations of 6 µs were carried out on these Spt16N structures and on two additional protein structures with varying degrees of similarity to it. The simulations revealed that correlation in patterns of backbone fluctuations vary linearly with sequence identity. This trend could not be inferred using crystallographic B-factors. Further, normal mode analysis suggested a similar pattern of inter-domain (inter-lobe) motions not only among Spt16N homologs, but also in the M24 peptidase structure. On the other hand, MD simulation results highlighted conserved motions that were found unique for Spt16N protein, this along with electrostatics trends shed light on functional aspects of Spt16N.Communicated by Ramaswamy H. Sarma.

Bioactive plantaricins as potent anti-cancer drug candidates: double docking, molecular dynamics simulation and in vitro cytotoxicity analysis

The medical community is desperate for a reliable source of medications to alleviate the severity of conventional cancer treatments and prevent secondary microbial infections in oncological patients. In this regard, plantaricins from lactic acid bacteria were explored as prospective drug candidates against known anti-cancer drug targets. Three plantaricins, JLA-9, GZ1-27 and BN, have a binding affinity of -8.8, -8.6 and -7.2 kcal/mol, respectively, with squalene synthase (SQS), a key molecule in lung cancer metastasis. All three plantaricins displayed analogous binding patterns as SQS inhibitors and generated hydrogen and hydrophobic interactions with ARG 47, ARG 188, PHE24, LEU183 and PRO292. Structural stability of docked complexes was validated using molecular dynamics simulation derived parameters such as RMSD, RMSF and radius of gyration. Based on MD simulation results, conformational changes and stabilities of docked SQS/plantaricin complexes with respect to the time frame were evaluated using machine learning (logistic regression algorithm). Double docking with SQS/matrix metalloproteinase MMP1 and PCA analysis revealed the potential of plantaricin JLA-9 as a multi-substrate inhibitor. Further, plantaricin JLA-9 induced a significant cytotoxic response against the lung carcinoma cell line (A549) in a dose and time dependent manner with inhibition concentration (IC50) of 0.082 µg/ml after 48 h. However, plantaricin JLA-9 did not induce cytotoxicity in normal lung cells (L-132), as the IC50 value was not obtained even at a higher dose of 0.8 µg/ml. In silico pharmacokinetic (ADMET) profile implies that plantaricin JLA-9 could be developed as new age anti-cancer therapeutic with a preference for parenteral administration.Communicated by Ramaswamy H. Sarma.

Effect of mutation of NS2B cofactor residues on Dengue 2 NS2B-NS3 protease complex - an insight to viral replication

Dengue fever is an endemic virus-borne disease that causes many severe ailments, including dengue hemorrhagic fever and dengue shock syndrome. NS2B-NS3 protease is present in all four strains of the dengue virus. NS2B-NS3 is a non-structural protein that performs three distinct functions: protease activity, helicase activity, and nucleoside triphosphatase activity. NS2B-NS3 pro-complex plays a crucial role in viral replication, and NS2B interacts with NS3 protease at a flat active site with an amino acid of the N-terminal region. NS2B acts as a cofactor for NS3 protease. In the current study, the conserved residues of NS2B were identified. Dengue virus-2 NS2B was mutated at the identified conserved amino acid region to investigate the role of NS2B on activation of NS3 pro. Molecular dynamics simulations were performed to investigate the mutated complex's changes in stability, conformation, and free energy. The EAG mutant complex exhibited more unstable conformation, less hydrogen bond formation, and high binding energy than wild type. This result suggests a vital role of E63, A65, G69 mutation in NS2B for the interruption of activation of the NS3.Communicated by Ramaswamy H. Sarma.

The Binding of Different Substrate Molecules at the Docking Site and the Active Site of γ-Secretase Can Trigger Toxic Events in Sporadic and Familial Alzheimer's Disease

Pathogenic changes in γ-secretase activity, along with its response to different drugs, can be affected by changes in the saturation of γ-secretase with its substrate. We analyze the saturation of γ-secretase with its substrate using multiscale molecular dynamics studies. We found that an increase in the saturation of γ-secretase with its substrate could result in the parallel binding of different substrate molecules at the docking site and the active site. The C-terminal domain of the substrate bound at the docking site can interact with the most dynamic presenilin sites at the cytosolic end of the active site tunnel. Such interactions can inhibit the ongoing catalytic activity and increase the production of the longer, more hydrophobic, and more toxic Aβ proteins. Similar disruptions in dynamic presenilin structures can be observed with different drugs and disease-causing mutations. Both, C99-βCTF-APP substrate and its different Aβ products, can support the toxic aggregation. The aggregation depends on the substrate N-terminal domain. Thus, the C99-βCTF-APP substrate and β-secretase path can be more toxic than the C83-αCTF-APP substrate and α-secretase path. Nicastrin can control the toxic aggregation in the closed conformation. The binding of the C99-βCTF-APP substrate to γ-secretase can be controlled by substrate channeling between the nicastrin and β-secretase. We conclude that the presented two-substrate mechanism could explain the pathogenic changes in γ-secretase activity and Aβ metabolism in different sporadic and familial cases of Alzheimer's disease. Future drug-development efforts should target different cellular mechanisms that regulate the optimal balance between γ-secretase activity and amyloid metabolism.

MDexciteR: Enhanced Sampling Molecular Dynamics by Excited Normal Modes or Principal Components Obtained from Experiments

Molecular dynamics with excited normal modes (MDeNM) is an enhanced sampling method for exploring conformational changes in proteins with minimal biases. The excitation corresponds to injecting kinetic energy along normal modes describing intrinsic collective motions. Herein, we developed a new automated open-source implementation, MDexciteR (https://github.com/mcosta27/MDexciteR), enabling the integration of MDeNM with two commonly used simulation programs with GPU support. Second, we generalized the method to include the excitation of principal components calculated from experimental ensembles. Finally, we evaluated whether the use of coarse-grained normal modes calculated with elastic network representations preserved the performance and accuracy of the method. The advantages and limitations of these new approaches are discussed based on results obtained for three different protein test cases: two globular and a protein/membrane system.

Comparative study of binding pocket structure and dynamics in cardiac and skeletal myosin

The development of small molecule myosin modulators has seen an increased effort in recent years due to their possible use in the treatment of cardiac and skeletal myopathies. Omecamtiv mecarbil (OM) is the first-in-class cardiac myotrope and the first to enter clinical trials. Its selectivity toward slow/beta-cardiac myosin lies at the heart of its function; however, little is known about the underlying reasons for selectivity to this isoform as opposed to other closely related ones such as fast-type skeletal myosins. In this work, we compared the structure and dynamics of the OM binding site in cardiac and in fasttype IIa skeletal myosin to identify possible reasons for OM selectivity. We found that the different shape, size, and composition of the binding pocket in skeletal myosin directly affects the binding mode and related affinity of OM, which is potentially a result of weaker interactions and less optimal molecular recognition. Moreover, we identified a side pocket adjacent to the OM binding site that shows increased accessibility in skeletal myosin compared with the cardiac isoform. These findings could pave the way to the development of skeletal-selective compounds that can target this region of the protein and potentially be used to treat congenital myopathies where muscle weakness is related to myosin loss of function.

Interdomain dynamics in human Replication Protein A regulates kinetics and thermodynamics of its binding to ssDNA

Human Replication Protein A (hRPA) is a multidomain protein that interacts with ssDNA intermediates to provide the latter much-needed stability during DNA metabolism and maintain genomic integrity. Although the ssDNA organization with hRPA was studied recently through experimental means, characterizing the underlying mechanism at the atomic level remains challenging because of the dynamic domain architecture of hRPA and poorly understood heterogeneity of ssDNA-protein interactions. Here, we used a computational framework, precisely tailored to capture protein-ssDNA interactions, and investigated the binding of hRPA with a 60 nt ssDNA. Two distinct binding mechanisms are realized based on the hRPA domain flexibility. For a rigid domain architecture of hRPA, ssDNA binds sequentially with hRPA domains, resulting in slow association kinetics. The binding pathway involves the formation of stable and distinct intermediate states. On contrary, for a flexible domain architecture of hRPA, ssDNA binds synergistically to the A and B domains followed by the rest of hRPA. The domain dynamics in hRPA alleviates the free energy cost of domain orientation necessary for specific binding with ssDNA, leading to fast association kinetics along a downhill binding free energy landscape. An ensemble of free energetically degenerate intermediate states is encountered that makes it arduous to characterize them structurally. An excellent match between our results with the available experimental observations provides new insights into the rich dynamics of hRPA binding to ssDNA and in general paves the way to investigate intricate details of ssDNA-protein interactions, crucial for cellular functioning.