Molecular characterization and structural dynamics of Aquaporin1 from walking catfish in lipid bilayers

Alternative splicing generates a novel ferroportin isoform with a shorter C-terminal and an intact iron- and hepcidin-binding property

Ferroportin (FPN) is a transmembrane protein and is the only known iron exporter that helps in maintaining iron homeostasis in vertebrates. To maintain stable iron equilibrium in the body, ferroportin works in conjunction with a peptide called hepcidin. In this study, we have identified an alternatively spliced novel isoform of the human SLC40A1 gene, which encodes for the FPN protein and is found to be expressed in different tissues. The novel transcript has an alternate last exon and encodes 31-amino acid long peptide sequence that replaces 104 amino acids at C-terminal in the novel transcript. Molecular modelling and molecular dynamics (MD) simulation studies revealed key structural features of the novel isoform (FPN-N). FPN-N was predicted to have 12 transmembrane domains similar to the reported isoform (FPN), despite being much smaller in size. FPN-N was found to interact with hepcidin, a key regulator of ferroportin activity. Also, the iron-binding sites were retained in the novel isoform as revealed by the MD simulation of FPN-N in bilipid membrane. The novel isoform identified in this study may play important role in iron homeostasis. However, further studies are required to characterize the FPN-N isoform and decipher its role inside the cell.

Machine learning and classical MD simulation to identify inhibitors against the P37 envelope protein of monkeypox virus

Monkeypox virus (MPXV) outbreak is a serious public health concern that requires international attention. P37 of MPXV plays a pivotal role in DNA replication and acts as one of the promising targets for antiviral drug design. In this study, we intent to screen potential analogs of existing FDA approved drugs of MPXV against P37 using state-of-the-art machine learning and computational biophysical techniques. AlphaFold2 guided all-atoms molecular dynamics simulations optimized P37 structure is used for molecular docking and binding free energy calculations. Similar to members of Phospholipase-D family , the predicted P37 structure also adopts a β-α-β-α-β sandwich fold, harbouring strongly conserved HxKxxxxD motif. The binding pocket comprises of Tyr48, Lys86, His115, Lys117, Ser130, Asn132, Trp280, Asn240, His325, Lys327 and Tyr346 forming strong hydrogen bonds and dense hydrophobic contacts with the screened analogs and is surrounded by positively charged patches. Loops connecting the two domains and C-terminal region exhibit high degree of flexibility. In some structural ensembles, the partial disorderness in the C-terminal region is presumed to be due to its low confidence score, acquired during structure prediction. Transition from loop to β-strands (244-254 aa) in P37-Cidofovir and its analog complexes advocates the need for further investigations. MD simulations support the accuracy of the molecular docking results, indicating the potential of analogs as potent binders of P37. Taken together, our results provide preferable understanding of molecular recognition and dynamics of ligand-bound states of P37, offering opportunities for development of new antivirals against MPXV. However, the need of in vitro and in vivo assays for confirmation of these results still persists.Communicated by Ramaswamy H. Sarma.

All-atoms molecular dynamics study to screen potent efflux pump inhibitors against KpnE protein of Klebsiella pneumoniae

The Small Multidrug Resistance efflux pump protein KpnE, plays a pivotal role in multi-drug resistance in Klebsiella pneumoniae. Despite well-documented study of its close homolog, EmrE, from Escherichia coli, the mechanism of drug binding to KpnE remains obscure due to the absence of a high-resolution experimental structure. Herein, we exclusively elucidate its structure-function mechanism and report some of the potent inhibitors through drug repurposing. We used molecular dynamics simulation to develop a dimeric structure of KpnE and explore its dynamics in lipid-mimetic bilayers. Our study identified both semi-open and open conformations of KpnE, highlighting its importance in transport process. Electrostatic surface potential map suggests a considerable degree of similarity between KpnE and EmrE at the binding cleft, mostly occupied by negatively charged residues. We identify key amino acids Glu14, Trp63 and Tyr44, indispensable for ligand recognition. Molecular docking and binding free energy calculations recognizes potential inhibitors like acarbose, rutin and labetalol. Further validations are needed to confirm the therapeutic role of these compounds. Altogether, our membrane dynamics study uncovers the crucial charged patches, lipid-binding sites and flexible loop that could potentiate substrate recognition, transport mechanism and pave the way for development of novel inhibitors against K. pneumoniae.Communicated by Ramaswamy H. Sarma.

Toxoplasma gondii aspartic protease 5: structural basis of substrate binding and inhibition mechanism

Toxoplasma gondii, a worldwide prevalent parasite is responsible for causing toxoplasmosis in almost all warm-blooded animals, including humans. Golgi-resident T. gondii aspartic protease 5 (TgASP5) plays an essential role in the maturation and export of the effector proteins those modulate the host immune system to establish a successful infection. Hence, inhibiting this enzyme can be a possible therapeutic strategy against toxoplasmosis. This is the first report of the detailed structural investigations of the TgASP5 mature enzyme using molecular modeling and an all-atom simulation approach which provide in-depth knowledge of the active site architecture of TgASP5. The analysis of the binding mode of the TEXEL (Toxoplasma EXport Element) substrate to TgASP5 highlighted the importance of the active site residues. Ser505, Ala776 and Tyr689 in the S2 binding pocket are responsible for the specificity towards Arg at the P2 position of TEXEL substrate. The molecular basis of inhibition by the only known inhibitor RRLStatine has been identified, and our results show that it blocks the active site by forming a hydrogen bond with a catalytic aspartate. Besides that, known aspartic protease inhibitors were screened against TgASP5 using docking, MD simulations and MM-PBSA binding energy calculations. The top-ranked inhibitors (SC6, ZY1, QBH) showed higher binding energy than RRLStatine. Understanding the structural basis of substrate recognition and the binding mode of these inhibitors will help to develop potent mechanistic inhibitors against TgASP5. This study will also provide insights into the structural features of pepsin-like aspartic proteases from other apicomplexan parasites for developing antiparasitic agents.Communicated by Ramaswamy H. Sarma.

Spectroscopic insights with molecular docking and molecular dynamic simulation studies of anticancer drug 5-Fluorouracil targeting human pyruvate kinase m2

This study was conducted to test the efficacy of 5-fluorouracil (5-FU) as an anticancer drug against the human pyruvate kinase isozyme M2 (PKM2) using spectroscopic, molecular docking and molecular dynamic simulation studies. PKM2 fluorescence quenching studies in the presence of 5-FU performed at three different temperatures indicates dynamic quenching processes with single-set of binding (n ≈ 1) profile. The biomolecular quenching constants (kq) and the effective binding constants (Kb) obtained are shown to increase with temperature. The calculated enthalpy (ΔH) and entropy changes (ΔS) are estimated to be -118.06 kJ/mol and 146.14 kJ/mol/K respectively, which suggest the possible mode of interaction as electrostatic and hydrogen bonding. Further, these values were used to estimate the free energy changes (ΔG) and that increases with temperature. The negative ΔG values clearly indicates spontaneous binding process that stabilizes the complex formed between 5-FU and PKM2. Far-UV CD spectra of PKM2 in the presence of 5-FU shows decrease in α-helix contents which point towards the destabilization of secondary structure that weakens the biological activity of PKM2. The intrinsic fluorescence study and circular dichroism (CD) spectra showed minor conformational changes of PKM2 in the presence of 5-FU. Additionally, the results obtained from molecular docking and all-atom molecular dynamic simulation study supports the insight of the spectroscopic binding studies, and strengthens the dynamic stability of the complex between 5-FU and PKM2 through H-bonding. This study establishes a paradigm of 5-FU-PKM2 complexation and the efficacy of 5-FU that compromises the biological activity of the targeted PKM2.Communicated by Ramaswamy H. Sarma.

Comparative analysis of SilA-laccase mediated degradation of ciprofloxacin, norfloxacin and ofloxacin and interpretation of the possible catalytic mechanism

Fluoroquinolones (FQs) are the most commonly used antimicrobial drugs and regardless of their advantages in the healthcare sector, the pollution of these antimicrobial drugs in the environment has big concerns about human and environmental health. The presence of these antibiotic drugs even at the lowest concentrations in the environment has resulted in the emergence and spread of antibiotic resistance. Hence, it is necessary to remediate these pollutants from the environment. Previously alkaline laccase (SilA) from Streptomyces ipomoeae has been demonstrated to show degrading potentials against two of the FQs, Ciprofloxacin (CIP) and Norfloxacin (NOR); however, the molecular mechanism was not elucidated in detail. In this study, we have analyzed the possible molecular catalytic mechanism of FQ degrading SilA-laccase for the degradation of the FQs, CIP, NOR and Ofloxacin (OFL) using three-dimensional protein structure modeling, molecular docking and molecular dynamic (MD) studies. The comparative protein sequence analysis revealed the presence of tetrapeptide conserved catalytic motif, His102-X-His104-Gly105. After evaluating the active site of the enzyme in depth using CDD, COACH and S-site tools, we have identified the catalytic triad composed of three conserved amino acid residues, His102, Val103 and Tyr108 with which ligands interacted during the catalysis process. By analyzing the MD trajectories, it is revealed that the highest degradation potential of SilA is for CIP followed by NOR and OFL. Ultimately, this study provides the possible comparative catalytic mechanism for the degradation of CIP, NOR and OFL by the SilA enzyme.Communicated by Ramaswamy H. Sarma.

Lipid and cholesterols modulate the dynamics of SARS-CoV-2 viral ion channel ORF3a and its pathogenic variants

SARS-CoV-2 accessory protein, ORF3a is a putative ion channel which immensely contributes to viral pathogenicity by modulating host immune responses and virus-host interactions. Relatively high expression of ORF3a in diseased individuals and implication with inflammasome activation, apoptosis and autophagy inhibition, ratifies as an effective target for developing vaccines and therapeutics. Herein, we present the elusive dynamics of ORF3a-dimeric state using all-atoms molecular dynamics (MD) simulations at μ-seconds scale in a heterogeneous lipid-mimetic system in multiple replicates. Additionally, we also explore the effect of non-synonymous pathogenic mutations on ORF3a ion channel activity and viral pathogenicity in different SARS-CoV-2 variants using various structure-based protein stability (ΔΔG) tools and computational saturation mutagenesis. Our study ascertains the role of phosphatidylcholines and cholesterol in modulating the structure of ORF3a, which perturbs the size and flexibility of the polar cavity that allows permeation of large cations. Discrete trend in ion channel pore radius and area per lipid arises the premise that presence of lipids might also affect the overall conformation of ORF3a. MD structural-ensembles, in some replicates rationalize the crucial role of TM2 in maintaining the native structure of ORF3a. We also infer that loss of structural stability primarily grounds for pathogenicity in more than half of the pathogenic variants of ORF3a. Overall, the effect of mutation on alteration of ion permeability of ORF3a, proposed in this study brings mechanistic insights into variant consequences on viral membrane proteins of SARS-CoV-2, which can be utilized for the development of novel therapeutics to treat COVID-19 and other coronavirus diseases.

Effect of AQP4 and its palmitoylation on the permeability of exogenous reactive oxygen species: Insights from computational study

Aquaporin 4 (AQP4) facilitates the transport of reactive oxygen species (ROS). Both cancer cells and the ionizing radiation microenvironment can induce posttranslational modifications (PTMs) in AQP4, which may affect its permeability to ROS. Because this ROS diffusion process is rapid, microscopic, and instantaneous within and outside cells, conventional experimental methods are inadequate for elucidating the molecular mechanisms involved. In this study, computational methods were employed to investigate the permeability of exogenous ROS mediated by radiation in AQP4 at a molecular scale. We constructed a simulation system incorporating AQP4 and AQP4-Cysp13 in a complex lipid environment with ROS. Long-timescale molecular dynamics simulations were conducted to assess the structural stability of both AQP4 and AQP4-Cysp13. Free energy calculations were utilized to determine the ROS transport capability of the two AQP4 proteins. Computational electrophysiology and channel structural analysis quantitatively evaluated changes in ROS transport capacity under various radiation-induced transmembrane voltage microenvironments. Our findings demonstrate the distinct transport capabilities of AQP4 channels for water molecules and various types of ROS and reveal a decrease in transport efficiency when AQP4 undergoes palmitoylation modification. In addition, we have simulated the radiation-induced alteration of cell membrane voltage, which significantly affected the ROS transport capacity. We propose that this research will enhance the understanding of the molecular mechanisms governing the transport of exogenous ROS by AQP4 and elucidate the influence of palmitoylation on ROS transport. This study will also help clarify how different structural features of AQP4 affect the transport of exogenous ROS mediated by radiotherapy, thereby providing a theoretical molecular basis for the development of new treatment strategies that combine with radiotherapy.

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.

Genome-Wide Identification and Gene Expression Analysis of Sweet Cherry Aquaporins (Prunus avium L.) under Abiotic Stresses

Aquaporins (AQPs) are integral transmembrane proteins well known as channels involved in the mobilization of water, small uncharged molecules and gases. In this work, the main objective was to carry out a comprehensive study of AQP encoding genes in Prunus avium (cv. Mazzard F12/1) on a genome-wide scale and describe their transcriptional behaviors in organs and in response to different abiotic stresses. A total of 28 non-redundant AQP genes were identified in Prunus spp. Genomes, which were phylogenetically grouped into five subfamilies (seven PIPs, eight NIPs, eight TIPs, three SIPs and two XIPs). Bioinformatic analyses revealed a high synteny and remarkable conservation of structural features among orthologs of different Prunus genomes. Several cis-acting regulatory elements (CREs) related to stress regulation were detected (ARE, WRE3, WUN, STRE, LTR, MBS, DRE, AT-rich and TC-rich). The above could be accounting for the expression variations associated with plant organs and, especially, each abiotic stress analyzed. Gene expressions of different PruavAQPs were shown to be preferentially associated with different stresses. PruavXIP2;1 and PruavXIP1;1 were up-regulated in roots at 6 h and 72 h of hypoxia, and in PruavXIP2;1 a slight induction of expression was also detected in leaves. Drought treatment strongly down-regulated PruavTIP4;1 but only in roots. Salt stress exhibited little or no variation in roots, except for PruavNIP4;1 and PruavNIP7;1, which showed remarkable gene repression and induction, respectively. Interestingly, PruavNIP4;1, the AQP most expressed in cherry roots subjected to cold temperatures, also showed this pattern in roots under high salinity. Similarly, PruavNIP4;2 consistently was up-regulated at 72 h of heat and drought treatments. From our evidence is possible to propose candidate genes for the development of molecular markers for selection processes in breeding programs for rootstocks and/or varieties of cherry.

CHARMM-GUI Membrane Builder: Past, Current, and Future Developments and Applications

Molecular dynamics simulations of membranes and membrane proteins serve as computational microscopes, revealing coordinated events at the membrane interface. As G protein-coupled receptors, ion channels, transporters, and membrane-bound enzymes are important drug targets, understanding their drug binding and action mechanisms in a realistic membrane becomes critical. Advances in materials science and physical chemistry further demand an atomistic understanding of lipid domains and interactions between materials and membranes. Despite a wide range of membrane simulation studies, generating a complex membrane assembly remains challenging. Here, we review the capability of CHARMM-GUI Membrane Builder in the context of emerging research demands, as well as the application examples from the CHARMM-GUI user community, including membrane biophysics, membrane protein drug-binding and dynamics, protein-lipid interactions, and nano-bio interface. We also provide our perspective on future Membrane Builder development.

Molecular dynamic simulation and functional analysis of pathogenic PTEN mutations in glioblastoma

PTEN, a dual-phosphatase and scaffold protein, is one of the most commonly mutated tumour suppressor gene across various cancer types in human. The aim of this study therefore was to investigate the stability, structural and functional effects, and pathogenicity of 12 missense PTEN mutations (R15S, E18G, G36R, N49I, Y68H, I101T, C105F, D109N, V133I, C136Y, R173C and N276S) found by next generation sequencing of the PTEN gene in tissue samples obtained from glioblastoma patients. Computational tools and molecular dynamic simulation programs were used to identify the deleterious effects of these mutations. Furthermore, PTEN mRNA and protein expression levels were evaluated by qRT-PCR, Western Blot, and immunohistochemistry staining methods. Various computational tools predicted strong deleterious effects for the G36R, C105F, C136Y and N276S mutations. Molecular dynamic simulation revealed a significant decrease in protein stability for the Y68H and N276S mutations when compared with the wild type protein; whereas, C105F, D109N, V133I and R173C showed partial stability reduction. Significant residual fluctuations were observed in the R15S, N49I and C136Y mutations and radius of gyration graphs revealed the most compact structure for D109N and least for C136Y. In summary, our study is the first one to show the presence of PTEN E18G, N49I, D109N and N276S mutations in glioblastoma patients; where, D109N is neutral and N276S is a damaging and disease-associated mutation.Communicated by Ramaswamy H. Sarma.

Structural insights into the RNA interaction with Yam bean Mosaic virus (coat protein) from Pachyrhizus erosus using bioinformatics approach

Plants are constantly threatened by a virus infection, i.e., Potyviruses, the second largest genus of plant viruses which results in several million-dollar losses in various essential crops globally. Yam bean (Pachyrhizus erosus) is considered to be one of the essential tuberous legume crops holding a great potential source of starch. Yam Bean Mosaic Virus (YBMV) of Potyvirus group belonging to the family potyviridae affects Yam bean and several angiosperms both in the tropical and sub-tropical regions causing large economical losses in crops. In this study, we attempted to understand the sequence-structure relationship and mode of RNA binding mechanism in YBMV CP using in silico integrative modeling and all-atoms molecular dynamics (MD) simulations. The assembly of coat protein (CP) subunits from YBMV and the plausible mode of RNA binding were compared with the experimental structure of CP from Watermelon mosaic virus potyvirus (5ODV). The transmembrane helix region is present in the YBMV CP sequence ranging from 76 to 91 amino acids. Like the close structural-homolog, 24 CPs monomeric sub-units formed YBMV a conserved fold. Our computational study showed that ARG124, ARG155, and TYR151 orient towards the inner side of the virion, while, THR122, GLN125, SER92, ASP94 reside towards the outer side of the virion. Despite sharing very low sequence similarity with CPs from other plant viruses, the strongly conserved residues Ser, Arg, and Asp within the RNA binding pocket of YBMV CP indicate the presence of a highly conserved RNA binding site in CPs from different families. Using several bioinformatics tools and comprehensive analysis from MD simulation, our study has provided novel insights into the RNA binding mechanism in YBMV CP. Thus, we anticipate that our findings from this study will be useful for the development of new therapeutic agents against the pathogen, paving the way for researchers to better control this destructive plant virus.