Molecular Dynamics Simulation Approach Research Articles

Integrated virtual screening and molecular dynamics simulation approaches revealed potential natural inhibitors for DNMT1 as therapeutic solution for triple negative breast cancer

Triple negative breast cancers (TNBC) are clinically heterogeneous but mostly aggressive malignancies devoid of expression of the estrogen, progesterone, and HER2 (ERBB2 or NEU) receptors. It accounts for 15–20% of all cases. Altered epigenetic regulation including DNA hypermethylation by DNA methyltransferase 1 (DNMT1) has been implicated as one of the causes of TNBC tumorigenesis. The antitumor effect of DNMT1 has also been explored in TNBC that currently lacks targeted therapies. However, the actual treatment for TNBC is yet to be discovered. This study is attributed to the identification of novel drug targets against TNBC. A comprehensive docking and simulation analysis was performed to optimize promising new compounds by estimating their binding affinity to the target protein. Molecular dynamics simulation of 500 ns well complemented the binding affinity of the compound and revealed strong stability of predicted compounds at the docked site. Calculation of binding free energies using MMPBSA and MMGBSA validated the strong binding affinity between compound and binding pockets of DNMT1. In a nutshell, our study uncovered that Beta-Mangostin, Gancaonin Z, 5-hydroxysophoranone, Sophoraflavanone L, and Dorsmanin H showed maximum binding affinity with the active sites of DNMT1. Furthermore, all of these compounds depict maximum drug-like properties. Therefore, the proposed compounds can be a potential candidate for patients with TNBC, but, experimental validation is needed to ensure their safety. Communicated by Ramaswamy H. Sarma

Evolutionary relevance of single nucleotide variants within the forebrain exclusive human accelerated enhancer regions

BackgroundHuman accelerated regions (HARs) are short conserved genomic sequences that have acquired significantly more nucleotide substitutions than expected in the human lineage after divergence from chimpanzees. The fast evolution of HARs may reflect their roles in the origin of human-specific traits. A recent study has reported positively-selected single nucleotide variants (SNVs) within brain-exclusive human accelerated enhancers (BE-HAEs) hs1210 (forebrain), hs563 (hindbrain) and hs304 (midbrain/forebrain). By including data from archaic hominins, these SNVs were shown to be Homo sapiens-specific, residing within transcriptional factors binding sites (TFBSs) for SOX2 (hs1210), RUNX1/3 (hs563), and FOS/JUND (hs304). Although these findings suggest that the predicted modifications in TFBSs may have some role in present-day brain structure, work is required to verify the extent to which these changes translate into functional variation.ResultsTo start to fill this gap, we investigate the SOX2 SNV, with both forebrain expression and strong signal of positive selection in humans. We demonstrate that the HMG box of SOX2 binds in vitro with Homo sapiens-specific derived A-allele and ancestral T-allele carrying DNA sites in BE-HAE hs1210. Molecular docking and simulation analysis indicated highly favourable binding of HMG box with derived A-allele containing DNA site when compared to site carrying ancestral T-allele.ConclusionThese results suggest that adoptive changes in TF affinity within BE-HAE hs1210 and other HAR enhancers in the evolutionary history of Homo sapiens might have brought about changes in gene expression patterns and have functional consequences on forebrain formation and evolution.MethodsThe present study employ electrophoretic mobility shift assays (EMSA) and molecular docking and molecular dynamics simulations approaches.

Unraveling the molecular mechanism of novel leukemia mutations on NTRK2 (A203T & R458G) and NTRK3 (E176D & L449F) genes using molecular dynamics simulations approach

Background: NTRK1, NTRK2, and NTRK3 are members of the neurotrophic receptor tyrosine kinases (NTRK) family, which encode TrkA, TrkB, and TrkC receptors, respectively. Hematologic cancers are also linked to point mutations in the NTRK gene's kinase domain. Trk fusions are the most common genetic change associated with oncogenic activity in Trk-driven liquid tumors. Thus, point mutations in NTRK genes may also play a role in tumorigenesis. The structural and functional effect of mutations in Trk-B & Trk-C proteins remains unclear. Methods: In this research, Homology (threading-based approach) modeling and the all-atom molecular dynamics simulations approaches are applied to examine the structural and functional behavior of native and mutant Trk-B and Trk-C proteins at the molecular level. Results: The result of this study reveals how the mutations in Trk-B (A203T & R458G) and Trk-C (E176D & L449F) proteins lost their stability and native conformations. The Trk-B mutant A203T became more flexible than the native protein, whereas the R458G mutation became more rigid than the native conformation of the Trk-B protein. Also, the Trk-C mutations (E176D & L449F) become more rigid compared to the native structure. Conclusions: This structural transition may interrupt the function of Trk-B and Trk-C proteins. Observing the impact of NTRK-2/3 gene alterations at the atomic level could aid in discovering a viable treatment for Trk-related leukemias.

Repurposing FIASMAs against Acid Sphingomyelinase for COVID-19: A Computational Molecular Docking and Dynamic Simulation Approach.

Over the past few years, COVID-19 has caused widespread suffering worldwide. There is great research potential in this domain and it is also necessary. The main objective of this study was to identify potential inhibitors against acid sphingomyelinase (ASM) in order to prevent coronavirus infection. Experimental studies revealed that SARS-CoV-2 causes activation of the acid sphingomyelinase/ceramide pathway, which in turn facilitates the viral entry into the cells. The objective was to inhibit acid sphingomyelinase activity in order to prevent the cells from SARS-CoV-2 infection. Previous studies have reported functional inhibitors against ASM (FIASMAs). These inhibitors can be exploited to block the entry of SARS-CoV-2 into the cells. To achieve our objective, a drug library containing 257 functional inhibitors of ASM was constructed. Computational molecular docking was applied to dock the library against the target protein (PDB: 5I81). The potential binding site of the target protein was identified through structural alignment with the known binding pocket of a protein with a similar function. AutoDock Vina was used to carry out the docking steps. The docking results were analyzed and the inhibitors were screened based on their binding affinity scores and ADME properties. Among the 257 functional inhibitors, Dutasteride, Cepharanthine, and Zafirlukast presented the lowest binding affinity scores of -9.7, -9.6, and -9.5 kcal/mol, respectively. Furthermore, computational ADME analysis of these results revealed Cepharanthine and Zafirlukast to have non-toxic properties. To further validate these findings, the top two inhibitors in complex with the target protein were subjected to molecular dynamic simulations at 100 ns. The molecular interactions and stability of these compounds revealed that these inhibitors could be a promising tool for inhibiting SARS-CoV-2 infection.

Thermodiffusion of CO2 in Water by Nonequilibrium Molecular Dynamics Simulations

The components of a fluid mixture may segregate due to the Soret effect, a coupling phenomenon in which mass flux can be induced by a thermal gradient. In this work, we evaluate systematically the thermodiffusion of the CO2-H2O mixture, and the influence of the geothermal gradient on CO2 segregation in deep saline aquifers in CO2 storage. The eHeX method, a nonequilibrium molecular dynamics simulation approach, is judiciously selected to simulate the phenomenon. At 350 K, 400 bar, and CO2 mole fraction of 0.02 (aquifer conditions), CO2 accumulates on the cold side, and the thermal diffusion factor is close to 1 in a number of force fields. The lower the temperature, the higher is the separation and the thermal diffusion factor. In colder regions, water self-association is stronger, whereas the CO2-H2O cross-association and the CO2-CO2 interactions enhance at higher temperatures. Thermodiffusion and gravitational segregation have opposite effects on CO2 segregation. At typical subsurface conditions, the Soret effect is more pronounced than gravity segregation, and CO2 concentrates in the top (colder region). Our work sets the stage to model the effect of electrolytes on CO2 segregation in subsurface aquifers and other areas of interest.

Pharmacoinformatics and molecular dynamics simulation approach to identify anti-diarrheal potentials of Centella asiatica (L.) Urb. against Vibrio cholerae

Vibrio cholerae, the etiological agent of cholera, causes dehydration and severe diarrhea with the production of cholera toxin. Due to the acquired antibiotic resistance, V. cholerae has drawn attention to the establishment of novel medications to counteract the virulence and viability of the pathogen. Centella asiatica is a medicinal herb native to Bangladesh that has a wide range of medicinal and ethnobotanical applications including anti-bacterial properties. In the present investigation, a total of 25 bioactive phytochemicals of C. asiatica have been screened virtually through molecular docking, ADMET (Absorption, Distribution, Metabolism, Excretion, and Toxicity) analyses, and molecular dynamics simulation. Our results revealed four lead compounds as Viridiflorol (−8.7 Kcal/mol), Luteolin (−8.1 Kcal/mol), Quercetin (−8.0 Kcal/mol) and, Geranyl acetate (−7.1 Kcal/mol) against V. cholerae Toxin co-regulated pilus virulence regulatory protein (ToxT). All the lead compounds have been found to possess favorable pharmacokinetic, pharmacodynamics, and molecular dynamics properties. Toxicity analysis revealed satisfactory results with no major side effects. Molecular dynamics simulation was performed for 100 ns that revealed noteworthy conformational stability and structural compactness for all the lead compounds, especially for Quercetin. Target class prediction unveiled enzymes in most of the cases and some experimental and investigational drugs were found as structurally similar analogs of the lead compounds. These findings could aid in the development of novel therapeutics targeting Cholera disease and we strongly recommend in vitro trials of our experimental findings. Communicated by Ramaswamy H. Sarma

Delineating the impact of pathogenic mutations on the conformational dynamics of HDL’s vital protein ApoA1: a combined computational and molecular dynamic simulation approach

Apolipoprotein A1 (ApoA1), is the important component of high‐density lipoproteins (HDL), that has key role in HDL biogenesis, cholesterol trafficking, and reverse cholesterol transport (RCT). Non-synonymous Single Nucleotide Polymorphisms (nsSNPs) in ApoA1 have been linked to cardiovascular diseases and amyloidosis as they alter the protein’s native structure and function. Therefore in this study, we attempted to understand the molecular pathogenicity profile of nsSNPs of ApoA1 using various computational approaches. We used state-of-the-art computational methods to thoroughly investigate the 295 ApoA1 nsSNPs at sequence and structural levels. Seven nsSNPs (L13R, L84R, L84P, L99P, R173P, L187P, and L238P) out of 295 were classified as the most deleterious and destabilizing. In order to estimate the effect of such destabilizing mutations on the protein conformation, all-atom molecular dynamics simulations (MDS) of ApoA1 wild-type (WT), L99P and R173P for 100 ns, was carried out using GROMACS 5.0.1 package. The MD simulation investigation revealed significant structural alterations in L99P and R173P. In addition, they had changed principal component analysis and electrostatic surface potential, decreased structural compactness, and intramolecular hydrogen bonds, which supported the rationale underpinning ApoA1 dysfunction with such mutations. This work sheds light on ApoA1 dysfunction due to single amino acid alterations, and offers new insight into the molecular basis of ApoA1-related diseases progression. Communicated by Ramaswamy H. Sarma

Molecular docking and dynamics simulation study of quinones and pyrones from Alternaria solani and Alternaria alternata with HSP90: an important therapeutic target of cancer

Although cancer continues to be one of the world’s major causes of death, current cancer drugs have many serious side effects. There remains a need for new anticancer agents to overcome these shortcomings. Alternaria is one of the most widespread fungal genera, many species of which produce several classes of metabolites with potential polypharmacological activities. A few quinones and pyrones from Alternaria spp. have proven to exert cytotoxic effects against certain cancer cell lines, but their molecular mode of action is not known. The current study aimed to investigate the potential mechanisms that underlie the anticancer activity of a few selected quinones and pyrones from Alternaria solani and Alternaria alternata by molecular docking and dynamic simulation approaches. The selected metabolites were screened for their binding affinity to Heat shock protein 90 (HSP90), which is a known anticancer drug target. Molecular docking studies have revealed that Macrosporin, Altersolanol B, Fonsecin, and Neoaltenuene have good binding affinities with the target protein and the stabilities of the formed complexes were evaluated through molecular dynamics simulations. By analyzing the Root Mean Square Distance (RMSD), Root Mean Square Fluctuation (RMSF), and Principal Component Analysis (PCA) plots obtained from molecular dynamics simulations, this study shows that the complexes of all 4 lead molecules with target protein are stable over a 100 ns period. Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) calculations were used to compute the binding free energies. The lead molecules were studied using in-silico analysis to determine their drug-likeness based on their Absorption, Distribution, Metabolism, Excretion and Toxicity (ADMET) and physicochemical properties. The results demonstrate that Macrosporin, Fonsecin, and Neoaltenuene could become promising anticancer molecules that target HSP90. Communicated by Ramaswamy H. Sarma

Study on Thermal Conductivity and Mechanical Properties of Cyclotriphosphazene Resin-Forced Epoxy Resin Composites

Cyclotriphosphazenes, a variety of inorganic rings together with a curing ingredient, 4,4′-methylene dianiline (MDA), are mainly used to enhance the thermal conductivity and mechanical characteristics of epoxy resin (DGEBA). Three DGEBA@MDA, HGCP@MDA, and thermosets were produced, and their curing behaviors were investigated. Using a molecular dynamics (MD) approach, the impact of cyclotriphosphazene on the characteristics of DGEBA composites is thoroughly explored in this paper. Results indicated that the glass transition temperatures (Tg) of DGEBA containing HGCP had slightly decreased compared to DGEBA. With the addition of HGCP to DGEBA, epoxy resin (DGEBA@HGCP@MDA) has a high thermal conductivity of 0.215284 W/m·K, with an increase of 116.04% more than pure DGEBA (0.185524 W/m·K). Moreover, the DGEBA@HGCP@MDA composite has high mechanical strength with a specific Young’s modulus of 5.4902 GPa. In order to forecast and analyze certain performances directly associated with the microstructure characteristics of the various cross-linked resin systems and their composite materials, an MD simulation approach will be quite valuable.

Storage and diffusion of CO2 in covalent organic frameworks-A neural network-based molecular dynamics simulation approach.

As a consequence of the accelerated climate change, solutions to capture, store and potentially activate carbon dioxide received increased interest in recent years. Herein, it is demonstrated, that the neural network potential ANI-2x is able to describe nanoporous organic materials at approx. density functional theory accuracy and force field cost, using the example of the recently published two- and three-dimensional covalent organic frameworks HEX-COF1 and 3D-HNU5 and their interaction with CO2 guest molecules. Along with the investigation of the diffusion behaviour, a wide range of properties of interest is analyzed, such as the structure, pore size distribution and host-guest distribution functions. The workflow developed herein facilitates the estimation of the maximum CO2 adsorption capacity and is easily generalizable to other systems. Additionally, this work illustrates, that minimum distance distribution functions can be a highly useful tool in understanding the nature of interactions in host-gas systems at the atomic level.

An Investigation on the Compatibility of Epoxy Resin and Asphalt Based on Molecular Simulation

Epoxy resin and asphalt compatibility have a significant impact on the material qualities of epoxy-modified asphalt. At present, there is a shortage of the study on the compatibility mechanism from the standpoint of nano. This paper evaluated the influence of temperature on the compatibility of epoxy resin (ER) and asphalt using the molecular dynamics (MD) simulation approach to explore the compatibility between ER and asphalt in greater detail. First, a twelve-component molecular model of asphalt was built, as well as a model of the system following cross-linking of epoxy resin and curing agent. And model optimization and model validation were performed. Secondly, the model of ER and asphalt co-blended system was constructed by combining the optimized model with the Amorphous Cell module. Then, MD simulations were done at various temperatures to derive the solubility parameters and intermolecular interaction energy of the mixed system for each system. Finally, the compatibility of ER with asphalt at different temperatures was analyzed based on the calculation results. It was found that the compatibility of ER and asphalt was temperature dependent. As the temperature increased, the solubility parameters of ER and asphalt decreased, and the difference in solubility parameters Δδ between them took the minimum value at 170 °C. As the temperature rises, so does the total potential energy in each system, and the intermolecular interaction energy of ER and asphalt mixed system achieves the maximum value at 170 °C. Currently, the structure of ER and asphalt mixed system is the most stable, and the compatibility between the two is the best.

Identification of novel C-15 fluoro isosteviol derivatives for GABA-AT inhibition by in silico investigations

The treatment of epilepsy is associated with the inhibition of γ-aminobutyric acid-aminotransferase (GABA-AT), which suppresses the concentration of a key neurotransmitter GABA. Isosteviol, a natural bioactive molecule, has been reported to possess an anticonvulsant property. In this work, we first reported a series of C-15 fluoro isosteviol analogs which are bearing different functional groups at C-16 to investigate the interactions with GABA-AT by applying molecular docking and molecular dynamic simulation approach. The results revealed that all fluoro isosteviol analogs displayed a greater binding affinity than references vigabatrin, an FDA-approved GABA-AT inactivator, and CPP-115, which has Orphan Drug Designation status, and positioned at the same binding site as references. Furthermore, molecular dynamic (MD) simulation studies on minimum (A1), maximum (E1) binding energy score of fluoro isosteviol analogs, and isosteviol (G1) revealed their stable complex formation in terms of RMSD, RMSF, RG, and hydrogen bond formation. All analogs were found to have drug-like nature, non-toxic, >80% absorption, and the majority tend to penetrate brain-blood-barrier (BBB). The investigations found in this study can help in the development of isosteviol derivatives as drugs for the treatment of epilepsy. The two-dimensional (2D) ligand structures were drawn using ChembioDraw Ultra 14.0. Molecular docking with Autodock4 and molecular dynamic simulation with GROMACS version 2020.1 were performed. The CHARMM27 all-atom force field was applied for writing the topology. Biovia Discovery Studio DS2021 was used for viewing and analyzing the protein-ligand complexes. The data generated from molecular dynamic simulation trajectories were plotted using the Origin® 8 software. The Open Babel software was utilized for extracting SMILEs files of all the fluoro isosteviol analogs. The drug-likeness and ADMET of the molecules were evaluated by SwissADME and ADMETlab 2.0 web tools.

Investigating the dynamic behavior of the nano-bubble in two-phase systems of argon and water: A molecular dynamics simulation approach

In this work, the dynamic behavior of the nano-bubble and the factors affecting its stability in argon two-phase systems in the absence and presence of water have been investigated using molecular dynamics simulations. According to the simulation results, addition of water to argon has led to the following changes in the nano-bubble: faster formation as well as increase in its stability, lifespan, and size. In addition, in the finite time stability of the nano-bubble, its radius of curvature slightly shift. These changes lead to oscillate motions, which can be considered as the nano-bubble breathing. The breathing movements cause more flexibility of the vapor nano-bubble, which makes it more stable. Moreover, the obtained results indicate that in the presence of water the contact angle increases.Importantly, we observed that the three-phase contact line is not only unstable during the stability of the nano-bubble, but also its values are slightly changing. Therefore, we showed what factors could be related to changes and fluctuations in the three-phase contact line.

Thermoresponsive Conductivity of Graphene-Based Fibers.

Smart materials are versatile material systems which exhibit a measurable response to external stimuli. Recently, smart material systems have been developed which incorporate graphene in order to share on its various advantageous properties, such as mechanical strength, electrical conductivity, and thermal conductivity as well as to achieve unique stimuli-dependent responses. Here, a graphene fiber-based smart material that exhibits reversible electrical conductivity switching at a relatively low temperature (60°C), is reported. Using molecular dynamics (MD) simulation and density functional theory-based non-equilibrium Green's function (DFT-NEGF) approach, it is revealed that this thermo-response behavior is due to the change in configuration of amphiphilic triblock dispersant molecules occurring in the graphene fiber during heating or cooling. These conformational changes alter the total number of graphene-graphene contacts within the composite material system, and thus the electrical conductivity as well. Additionally, this graphene fiber fabrication approach uses a scalable, facile, water-based method, that makes it easy to modify material composition ratios. In all, this work represents an important step forward to enable complete functional tuning of graphene-basedsmart materials at the nanoscale while increasing commercialization viability.

Structural Insights into Pseudomonas aeruginosa Exotoxin A-Elongation Factor 2 Interactions: A Molecular Dynamics Study.

Exotoxin A (ETA) is an extracellular secreted toxin and a single-chain polypeptide with A and B fragments that is produced by Pseudomonas aeruginosa. It catalyzes the ADP-ribosylation of a post-translationally modified histidine (diphthamide) on eukaryotic elongation factor 2 (eEF2), which results in the inactivation of the latter and the inhibition of protein biosynthesis. Studies show that the imidazole ring of diphthamide plays an important role in the ADP-ribosylation catalyzed by the toxin. In this work, we employ different in silico molecular dynamics (MD) simulation approaches to understand the role of diphthamide versus unmodified histidine in eEF2 on the interaction with ETA. Crystal structures of the eEF2-ETA complexes with three different ligands NAD+, ADP-ribose, and βTAD were selected and compared in the diphthamide and histidine containing systems. The study shows that NAD+ bound to ETA remains very stable in comparison with other ligands, enabling the transfer of ADP-ribose to the N3 atom of the diphthamide imidazole ring in eEF2 during ribosylation. We also show that unmodified histidine in eEF2 has a negative impact on ETA binding and is not a suitable target for the attachment of ADP-ribose. Analyzing of radius of gyration and COM distances for NAD+, βTAD, and ADP-ribose complexes revealed that unmodified His affects the structure and destabilizes the complex with all different ligands throughout the MD simulations.

Design and simulation of a caprylic acid enzymatically modified phosphatidylcholine micelle using a coarse-grained molecular dynamics simulations approach

Computationally simulated micelle models provide useful structural information on the molecular and biological sciences. One strategy to study the self-aggregation process of surfactant molecules that make up a micelle is through molecular dynamics (MD) simulations. In this study, a theoretical approach with a coarse-grained MD simulation (CG-MD) was employed to evaluate the critical micellar concentration (CMC), the micellization process, building a tridimensional (3D) model system of a micelle using data from the experimentally enzymatically modified phospholipids (PL) by phospholipase A1 (PA1). This required enzymatic interesterification of soybean phosphatidylcholine (PC) with caprylic acid, along with purification and characterization by chromatographic techniques to measure the esterified fatty acids and the corresponding PL composition. The number of molecules used in the CG-MD simulation system was determined from the experimental CMC data which was 0.025%. The molecular composition of the system is: 1 C 18:2, 2 C 8:0/8:0, 3 C 8:0/18:3n-9, 4 C 8:0/18:0, 5 C8:0/18:2n-6, 6 C8:0/18:1n-9, and 7 C 8:0/16:0. According to our theoretical results, the micelle model is structurally stable with an average Rg of 3.64 ± 0.10 Å, and might have an elliptical form with a radius of 24.6 Å. Regarding CMC value there was a relationship between the experimental data of the modified PLs and the theoretical analysis by GC-MD, which suggest that the enzymatic modification of PLs does not affect their self-aggregation properties. Finally, the micellar system obtained in the current research can be used as a simple and useful model to design optimal biocompatible nanoemulsions as possible vehicles for bioactive small molecules.Communicated by Ramaswamy H. Sarma

Discovery of AcrAB-TolC pump inhibitors: Virtual screening and molecular dynamics simulation approach

AcrAB-TolC tripartite efflux pump, which belongs to the RND superfamily, is a main multi-drug efflux system of Escherichia coli (E. coli) because of the broad resistance on various antibiotics. With the discovering of efflux pump inhibitors (EPIs), a combination between these and antibiotics is one of the most promising therapies. Therefore, building a virtual screening model with prediction capacities for the efflux pump inhibitory activities of candidates from DrugBank and ZINC15 dataset, is one of the key goals of this project. Based on the database of 170 diverse chemical structures collected from 28 research journals, two 2D-QSAR models and a 3D-pharmacophore model have been performed. On the AcrB protein (PDB 4DX7), two binding sites have been discovered that match to the hydrophobic trap in the distal pocket and the switch loop in the proximal pocket. After virtual screening processes, twenty candidate AcrAB-TolC inhibitors have been subjected to molecular dynamics simulations, binding free energy calculations and ADMET predictions. The results indicate that three compounds namely DB09233, DB02581, and DB15224 are potential inhibitors with ΔGbind of −42.30 ± 4.58, −40.76 ± 7.30 and −31.06 ± 7.63 kcal.mol−1, respectively. Communicated by Ramaswamy H. Sarma

Translocation of cell-penetrating peptides using the weighted ensemble method within molecular dynamics simulations approach.

We report spontaneous translocation of a single cell-penetrating peptide (CPP) across a model membrane using the weighted ensemble (WE) method within the molecular dynamics (MD) simulations approach. Furthermore, we show that a water pore plays a crucial role in translocating the peptide and a tilt angle of the peptide also affects the translocation. Last, using umbrella sampling, we obtain the potential of mean force (PMF) along the translocation path and show that a single CPP's translocation is energetically unfavorable.

Bioremediation potential of laccase for catalysis of glyphosate, isoproturon, lignin, and parathion: Molecular docking, dynamics, and simulation.

The present study aimed to investigate the catalytic degradation produced by laccase in the detoxification of glyphosate, isoproturon, lignin polymer, and parathion. We explored laccase-glyphosate, laccase-lignin polymer, laccase-isoproturon, and laccase-parathion using molecular docking (MD) and molecular dynamics simulation (MDS) approaches. The results suggest that laccase interacts well with glyphosate, lignin polymer, isoproturon, and parathion during biodegradation. We calculated the root mean square deviations (RMSD) of laccase-glyphosate, laccase-lignin polymer, laccase-isoproturon, and laccase-parathion as 0.24±0.02, 0.59±0.32, 0.43±0.07, and 0.43±0.06nm, respectively. In an aqueous solution, the stability of laccase with glyphosate, lignin polymer, isoproturon, and parathion is mediated through the formation of hydrophobic interactions, hydrogen bonds, and van der Waals interactions. The presence of xenobiotic toxic compounds in the active site changed the conformation of laccase. MDS of the laccase-substrate complexes confirmed their stability during catalytic degradation. Laccase assay results confirmed that the degradation of syringol, dihydroconiferyl alcohol, guaiacol, parathion, isoproturon, and glyphosate were 100%, 99.31%, 95.69%, 60.96%, 54.51%, and 48.34% within 2h, respectively. Taken together, we describe a novel method to understand the molecular-level biodegradation of xenobiotic compounds through laccase and its potential application in contaminant removal.

Arginine transportation mechanism through cationic amino acid transporter 1: insights from molecular dynamics studies

Metabolic and signaling mechanisms in mammalian cells are facilitated by the transportation of L-arginine (Arg) across the plasma membrane through cationic amino acid transporter (CAT) proteins. Due to a lack of argininosuccinate synthase (ASS) activity in various tumor cells such as acute myeloid leukemia, acute lymphocytic leukemia, and chronic lymphocytic leukemia, these tumor entities are arginine-auxotrophic and therefore depend on the uptake of the amino acid arginine. Cationic amino acid transporter-1 (CAT-1) is the leading arginine importer expressed in the aforementioned tumor entities. Hence, in the present study, to investigate the transportation mechanism of arginine in CAT-1, we performed molecular dynamics (MD) simulation methods on the modeled human CAT-1. The MM-PBSA approach was conducted to determine the critical residues interacting with arginine within the corresponding binding site of CAT-1. In addition, we found out that the water molecules have the leading role in forming the transportation channel within CAT-1. The conductive structure of CAT-1 was formed only when the water molecules were continuously distributed across the channel. Steered molecular dynamics (SMD) simulation approach showed various energy barriers against arginine transportation through CAT-1, especially while crossing the bottlenecks of the related channel. These findings at the molecular level might shed light on identifying the crucial amino acids in the binding of arginine to eukaryotic CATs and also provide fundamental insights into the arginine transportation mechanisms through CAT-1. Understanding the transportation mechanism of arginine is essential to developing CAT-1 blockers, which can be potential medications for some types of cancers.Communicated by Ramaswamy H. Sarma