Estrogen receptor alpha (ERα) acts as the transcription factor and the main therapeutic target against breast cancer. One of the compounds that has been shown to act as an ERα is α-mangostin. However, it still has weaknesses due to its low solubility and low potent activity. In this study, α-mangostin was modified by substituting -OH group at C6 using benzoyl derivatives through a step by step in silico study, namely pharmacokinetic prediction (https://preadmet.bmdrc.kr/adme/), pharmacophore modeling (LigandScout 4.1), molecular docking simulation (AutoDock 4.2), molecular dynamics simulation (AMBER 16) and a binding free energy analysis using MM-PBSA method. From the computational studies, three compounds which are derived from α-mangostin (AMB-1 (-9.84 kcal/mol), AMB-2 (-6.80 kcal/mol) and AMB-10 (-12.42 kcal/mol)) have lower binding free energy than α-mangostin (-1.77 kcal/mol), as evidenced by the binding free energy calculation using the MM-PBSA method. They can then be predicted to have potent activities as ERα antagonists.Communicated by Ramaswamy H. Sarma.
birgHPC, a bootable Linux Live CD has been developed to create high-performance clusters for bioinformatics and molecular dynamics studies using any Local Area Network (LAN)-networked computers. birgHPC features automated hardware and slots detection as well as provides a simple job submission interface. The latest versions of GROMACS, NAMD, mpiBLAST and ClustalW-MPI can be run in parallel by simply booting the birgHPC CD or flash drive from the head node, which immediately positions the rest of the PCs on the network as computing nodes. Thus, a temporary, affordable, scalable and high-performance computing environment can be built by non-computing-based researchers using low-cost commodity hardware.
Novel coronavirus disease (COVID-19) has become a pandemic threat to public health. Vaccines and targeted therapeutics to prevent infections and stop virus proliferation are currently lacking. Endoribonuclease Nsp15 plays a vital role in the life cycle, including replication and transcription as well as virulence of the virus. Here, we investigated Vitamin D for its in silico potential inhibition of the binding sites of SARS-CoV-2 endoribonuclease Nsp15. In this study, we selected Remdesivir, Chloroquine, Hydroxychloroquine and Vitamin D to study the potential binding affinity with the putative binding sites of endoribonuclease Nsp15 of COVID-19. The docking study was applied to rationalize the possible interactions of the target compounds with the active site of endoribonuclease Nsp 15. Among the results, Vitamin D was found to have the highest potency with strongest interaction in terms of LBE, lowest RMSD, and lowest inhibition intensity Ki than the other standard compounds. The investigation results of endoribonuclease Nsp15 on the PrankWeb server showed that there are three prospective binding sites with the ligands. The singularity of Vitamin D interaction with the three pockets, particularly in the second pocket, may write down Vitamin D as a potential inhibitor of COVID-19 Nsp15 endoribonuclease binding sites and favour addition of Vitamin D in the treatment plan for COVID-19 alone or in combination with the other used drugs in this purpose, which deserves exploration in further in vitro and in vivo studies.
Terminal moieties of most proteins are long known to be disordered and flexible. To unravel the functional role of these regions on the structural stability and biochemical properties of AT2 lipase, four C-terminal end residues, (Ile-Thr-Arg-Lys) which formed a flexible, short tail-like random-coil segment were targeted for mutation. Swapping of the tail-like region had resulted in an improved crystallizability and anti-aggregation property along with a slight shift of the thermostability profile. The lipolytic activity of mutant (M386) retained by 43 % compared to its wild-type with 18 % of the remaining activity at 45 °C. In silico analysis conducted at 25 and 45 °C was found to be in accordance to the experimental findings in which the RMSD values of M386 were more stable throughout the total trajectory in comparison to its wild-type. Terminal moieties were also observed to exhibit large movement and flexibility as denoted by high RMSF values at both dynamics. Variation in organic solvent stability property was detected in M386 where the lipolytic activity was stimulated in the presence of 25 % (v/v) of DMSO, isopropanol, and diethyl ether. This may be worth due to changes in the surface charge residues at the mutation point which probably involve in protein-solvent interaction.
Genomic surveillance is crucial for tracking emergence and spread of novel variants of pathogens, such as SARS-CoV-2, to inform public health interventions and to enforce control measures. However, in some settings especially in low- and middle- income counties, where sequencing platforms are limited, only certain patients get to be selected for sequencing surveillance. Here, we show that patients with multiple comorbidities potentially harbour SARS-CoV-2 with higher mutation rates and thus deserve more attention for genomic surveillance. The relationship between the patient comorbidities, and type of amino acid mutations was assessed. Correlation analysis showed that there was a significant tendency for mutations to occur within the ORF1a region for patients with higher number of comorbidities. Frequency analysis of the amino acid substitution within ORF1a showed that nsp3 P822L of the PLpro protease was one of the highest occurring mutations. Using molecular dynamics, we simulated that the P822L mutation in PLpro represents a system with lower Root Mean Square Deviation (RMSD) fluctuations, and consistent Radius of gyration (Rg), Solvent Accessible Surface Area (SASA) values-indicate a much stabler protein than the wildtype. The outcome of this study will help determine the relationship between the clinical status of a patient and the mutations of the infecting SARS-CoV-2 virus.
The utilization of organic solvents as reaction media for enzymatic reactions provides numerous industrially attractive advantages. However, an adaptation of enzyme towards organic solvent is unpredictable and not fully understood because of limited information on the organic solvent tolerant enzymes. To understand how the enzyme can adapt to the organic solvent environment, structural and computational approaches were employed. A recombinant elastase from Pseudomonas aeruginosa strain K was an organic solvent tolerant zinc metalloprotease was successfully crystallized and diffracted up to 1.39 Å. Crystal structure of elastase from strain K showed the typical, canonical alpha-beta hydrolase fold consisting of 10-helices (118 residues), 10- β-strands (38 residues) and 142 residues were formed other secondary structure such as loop and coil to whole structure. The elastase from Pseusomonas aeruginosa strain K possess His-140, His-144 and Glu-164 served as a ligand for zinc ion. The conserved catalytic triad was composed of Glu-141, Tyr-155 and His-223. Three-dimensional structure features such as calcium-binding and presence of disulphide-bridge contribute to the stabilizing the elastase structure. Molecular dynamic (MD) simulation of elastase revealed that, amino acid residues located at the surface area and disulphide bridge in Cys-30 to Cys-58 were responsible for enzyme stability in organic solvents.
Herein, we detail an atomic-level investigation of the cutinase enzyme encapsulated within a model metal-organic framework (MOF) platform using quantum mechanics calculations and molecular dynamics simulations. Cutinase, when encapsulated in an isoreticularly expanded MOF-74 (cutinase@IRMOF-74-VI), was proven to maintain its structural stability at temperatures that would otherwise denature the enzyme in its unprotected native state. Hydrogen bonding and salt bridge interactions, most notably involving arginine residues at the surface of the enzyme, were critical for stabilizing cutinase within the pore channels of IRMOF-74-VI. The findings reported support the viability of enzyme encapsulation in a porous material by demonstrating that a model enzyme not only retains its structural integrity but also remains accessible and active under extreme and foreign conditions.
We utilized molecular dynamics (MD) simulations and Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) free energy calculations to investigate the specificity of two oligonucleotide probes, namely probe B and probe D, in detecting single-stranded DNA (ssDNA) within three bacteria families: Enterobacteriaceae, Pasteurellaceae, and Vibrionaceae. Due to the limited understanding of molecular mechanisms in the previous research, we have extended the discussion to focus specifically on investigating the binding process of bacteria-probe DNA duplexes, with an emphasis on analyzing the binding free energy. The role of electrostatic contributions in the specificity between the oligonucleotide probes and the bacterial ssDNAs was investigated and found to be crucial. Our calculations yielded results that were highly consistent with the experimental data. Through our study, we have successfully exhibited the benefits of utilizing in-silico approaches as a powerful virtual-screening tool, particularly in research areas that demand a thorough comprehension of molecular interactions.
The activation of lipases has been postulated to proceed by interfacial activation, temperature switch activation, or aqueous activation. Recently, based on molecular dynamics (MD) simulation experiments, the T1 lipase activation mechanism was proposed to involve aqueous activation in addition to a double-flap mechanism. Because the open conformation structure is still unavailable, it is difficult to validate the proposed theory unambiguously to understand the behavior of the enzyme. In this study, we try to validate the previous reports and uncover the mystery behind the activation process using structural analysis and MD simulations. To investigate the effects of temperature and environmental conditions on the activation process, MD simulations in different solvent environments (water and water-octane interface) and temperatures (20, 50, 70, 80, and 100°C) were performed. Based on the structural analysis of the lipases in the same family of T1 lipase (I.5 lipase family), we proposed that the lid domain comprises α6 and α7 helices connected by a loop, thus forming a helix-loop-helix motif involved in interfacial activation. Throughout the MD simulations experiments, lid displacements were only observed in the water-octane interface, not in the aqueous environment with respect to the temperature effect, suggesting that the activation process is governed by interfacial activation coupled with temperature switch activation. Examining the activation process in detail revealed that the large structural rearrangement of the lid domain was caused by the interaction between the hydrophobic residues of the lid with octane, a nonpolar solvent, and this conformation was found to be thermodynamically favorable.
Molecular dynamics simulation was used to study the dynamic differences between native Aspergillus niger PhyA phytase and a mutant with 20 % greater thermostability. Atomic root mean square deviation, radius of gyration, and number of hydrogen bonds and salt bridges are examined to determine thermostability factors. The results suggest that, among secondary structure elements, loops have the most impact on the thermal stability of A. niger phytase. In addition, the location rather than the number of hydrogen bonds is found to have an important contribution to thermostability. The results also show that salt bridges may have stabilizing or destabilizing effect on the enzyme and influence its thermostability accordingly.
Candidates generated from unsaturated ketone (chalcone) demonstrated as strong, reversible and specific monoamine oxidase-B (MAO-B) inhibitory activity. For the research on MAO-B inhibition, our team has synthesized and evaluated a panel of aldoxime-chalcone ethers (ACE) and hydroxylchalcones (HC). The MAO-B inhibitory activity of several candidates is in the micro- to nanomolar range in these series. The purpose of this research was to develop predictive QSAR models and look into the relation between MAO-B inhibition by aldoxime and hydroxyl-functionalized chalcones. It was shown that the molecular descriptors ETA Shape P, MDEO-12, ETA dBetaP, SpMax1 Bhi and ETA EtaP B are significant in the inhibitory action of the MAO-B target. Using the current 2D QSAR models, potential chalcone-based MAO-B inhibitors might be created. The lead molecules were further analyzed by the detailed molecular dynamics study to establish the stability of the ligand-enzyme complex.Communicated by Ramaswamy H. Sarma.
The utilization of cold active lipases in organic solvents proves an excellent approach for chiral synthesis and modification of fats and oil due to the inherent flexibility of lipases under low water conditions. In order to verify whether this lipase can function as a valuable synthetic catalyst, the mechanism concerning activation of the lid and interacting solvent residues in the presence of organic solvent must be well understood. A new alkaline cold-adapted lipase, AMS8, from Pseudomonas fluorescens was studied for its structural adaptation and flexibility prior to its exposure to non-polar, polar aprotic and protic solvents. Solvents such as ethanol, toluene, DMSO and 2-propanol showed to have good interactions with active sites. Asparagine (Asn) and tyrosine (Tyr) were key residues attracted to solvents because they could form hydrogen bonds. Unlike in other solvents, Phe-18, Tyr-236 and Tyr-318 were predicted to have aromatic-aromatic side-chain interactions with toluene. Non-polar solvent also was found to possess highest energy binding compared to polar solvents. Due to this circumstance, the interaction of toluene and AMS8 lipase was primarily based on hydrophobicity and molecular recognition. The molecular dynamic simulation showed that lid 2 (residues 148-167) was very flexible in toluene and Ca(2+). As a result, lid 2 moves away from the catalytic areas, leaving an opening for better substrate accessibility which promotes protein activation. Only a single lid (lid 2) showed the movement following interactions with toluene, although AMS8 lipase displayed double lids. The secondary conformation of AMS8 lipase that was affected by toluene observed a reduction of helical strands and increased coil structure. Overall, this work shows that cold active lipase, AMS8 exhibits distinguish interfacial activation and stability in the presence of polar and non-polar solvents.
Dengue virus Type 2 (DENV-2) is predominant serotype causing major dengue epidemics. There are a number of studies carried out to find its effective antiviral, however to date, there is still no molecule either from peptide or small molecules released as a drug. The present study aims to identify small molecules inhibitor from National Cancer Institute database through virtual screening. One of the hits, D0713 (IC50 = 62 μM) bearing thioguanine scaffold was derivatised into 21 compounds and evaluated for DENV-2 NS2B/NS3 protease inhibitory activity. Compounds 18 and 21 demonstrated the most potent activity with IC50 of 0.38 μM and 16 μM, respectively. Molecular dynamics and MM/PBSA free energy of binding calculation were conducted to study the interaction mechanism of these compounds with the protease. The free energy of binding of 18 calculated by MM/PBSA is -16.10 kcal/mol compared to the known inhibitor, panduratin A (-11.27 kcal/mol), which corroborates well with the experimental observation. Results from molecular dynamics simulations also showed that both 18 and 21 bind in the active site and stabilised by the formation of hydrogen bonds with Asn174.
Dehalogenases continue to garner interest of the scientific community due to their potential applications in bioremediation of halogen-contaminated environment and in synthesis of various industrially relevant products. Example of such enzymes is DehL, an L-2-haloacid dehalogenase (EC 3.8.1.2) from Rhizobium sp. RC1 that catalyses the specific cleavage of halide ion from L-2-halocarboxylic acids to produce the corresponding D-2-hydroxycarboxylic acids. Recently, the catalytic residues of DehL have been identified and its catalytic mechanism has been fully elucidated. However, the enantiospecificity determinants of the enzyme remain unclear. This information alongside a well-defined catalytic mechanism are required for rational engineering of DehL for substrate enantiospecificity. Therefore, using quantum mechanics/molecular mechanics and molecular mechanics Poisson-Boltzmann surface area calculations, the current study theoretically investigated the molecular basis of DehL enantiospecificity. The study found that R51L mutation cancelled out the dehalogenation activity of DehL towards it natural substrate, L-2-chloropropionate. The M48R mutation, however introduced a new activity towards D-2-chloropropionate, conveying the possibility of inverting the enantiospecificity of DehL from L-to d-enantiomer with a minimum of two simultaneous mutations. The findings presented here will play important role in the rational design of DehL dehalogenase for improving substrate utility.
Congenital myopathy is a broad category of muscular diseases with symptoms appearing at the time of birth. One type of congenital myopathy is Congenital Fiber Type Disproportion (CFTD), a severely debilitating disease. The G48D and G48C mutations in the D-loop and the actin-myosin interface are the two causes of CFTD. These mutations have been shown to significantly affect the structure and function of muscle fibers. To the author's knowledge, the effects of these mutations have not yet been studied. In this work, the power stroke structure of the head domain of myosin and the wild and mutated types of actin were modeled. Then, a MD simulation was run for the modeled structures to study the effects of these mutations on the structure, function, and molecular dynamics of actin. The wild and mutated actins docked with myosin showed differences in hydrogen bonding patterns, free binding energies, and hydrogen bond occupation frequencies. The G48D and G48C mutations significantly impacted the conformation of D-loops because of their larger size compared to Glycine and their ability to interfere with the polarity or hydrophobicity of this neutralized and hydrophobic loop. Therefore, the mutated loops were unable to fit properly into the hydrophobic groove of the adjacent G-actin. The abnormal structure of D-loops seems to result in the abnormal assembly of F-actins, giving rise to the symptoms of CFTD. It was also noted that G48C and G48D did not form hydrogen bonds with myosin in the residue 48 location. Nevertheless, in this case, muscles are unable to contract properly due to muscle atrophy.
Brugia malayi is a filarial nematode, which causes lymphatic filariasis in humans. In 1995, the disease has been identified by the World Health Organization (WHO) as one of the second leading causes of permanent and long-term disability and thus it is targeted for elimination by year 2020. Therefore, accurate filariasis diagnosis is important for management and elimination programs. A recombinant antigen (BmR1) from the Bm17DIII gene product was used for antibody-based filariasis diagnosis in "Brugia Rapid". However, the structure and dynamics of BmR1 protein is yet to be elucidated. Here we study the three dimensional structure and dynamics of BmR1 protein using comparative modeling, threading and ab initio protein structure prediction. The best predicted structure obtained via an ab initio method (Rosetta) was further refined and minimized. A total of 5 ns molecular dynamics simulation were performed to investigate the packing of the protein. Here we also identified three epitopes as potential antibody binding sites from the molecular dynamics average structure. The structure and epitopes obtained from this study can be used to design a binder specific against BmR1, thus aiding future development of antigen-based filariasis diagnostics to complement the current diagnostics.
Dehalogenase E (DehE) is a non-stereospecific enzyme produced by the soil bacterium, Rhizobium sp. RC1. Till now, the catalytic mechanism of DehE remains unclear although several literature concerning its structure and function are available. Since DehE is non-stereospecific, the enzyme was hypothesized to follow a 'direct attack mechanism' for the catalytic breakdown of a haloacid. For a molecular insight, the DehE modelled structure was docked in silico with the substrate 2-chloropropionic acid (2CP) in the active site. The ideal position of DehE residues that allowed a direct attack mechanism was then assessed via molecular dynamics (MD) simulation. It was revealed that the essential catalytic water was hydrogen bonded to the 'water-bearer', Asn114, at a relatively constant distance of ∼2.0 Å after 50 ns. The same water molecule was also closely sited to the catalytic Asp189 at an average distance of ∼2.0 Å, signifying the imperative role of the latter to initiate proton abstraction for water activation. This reaction was crucial to promote a direct attack on the α-carbon of 2CP to eject the halide ion. The water molecule was oriented favourably towards the α-carbon of 2CP at an angle of ∼75°, mirrored by the formation of stable enzyme-substrate orientations throughout the simulation. The data therefore substantiated that the degradation of a haloacid by DehE followed a 'direct attack mechanism'. Hence, this study offers valuable information into future advancements in the engineering of haloacid dehalogenases with improved activity and selectivity, as well as functionality in solvents other than water.
Benzothiazole and its natural or synthetic derivatives have been used as precursors for several pharmacological agents for neuroprotective, anti-bacterial, and anti-allergic activities. The objective of the present study was to evaluate effects of benzothiazole analogs (compounds 1-26) for their immunomodulatory activities. Eight compounds (2, 4, 5, 8-10, 12, and 18) showed potent inhibitory activity on PHA-activated peripheral blood mononuclear cells (PBMCs) with IC50 ranging from 3.7 to 11.9 μM compared to that of the standard drug, prednisolone <1.5 μM. Some compounds (2, 4, 8, and 18) were also found to have potent inhibitory activities on the production of IL-2 on PHA/PMA-stimulated PBMCs with IC50 values ranging between <4.0 and 12.8 μM. The binding interaction of these compounds was performed through silico molecular docking. Compounds 2, 8, 9, and 10 significantly suppressed oxidative burst ROS production in phagocytes with IC50 values between <4.0 and 15.2 μM. The lipopolysaccharide (LPS)-induced nitrites in murine macrophages cell line J774 were found to be inhibited by compounds 4, 8, 9, and 18 at a concentration of 25 μg/mL by 56%, 91%, 58%, and 78%, respectively. Furthermore, compounds 5, 8, 12, and 18 showed significant (P<0.05) suppressive activity on Th-2 cytokine, interleukin 4 (IL-4) with an IC50 range of <4.0 to 40.3 μM. Interestingly compound 4 has shown a selective inhibitory activity on IL-2 and T cell proliferation (naïve T cell proliferation stage) rather than on IL-4 cytokine, while compound 12 displayed an interference with T-cell proliferation and IL-4 generation. Moreover compound 8 and 18 exert non-selective inhibition on both IL-2 and IL-4 cytokines, indicating a better interference with stage leading to humoral immune response and hence possible application in autoimmune diseases.
Both the inactive- and active-state CB1 receptor crystal structures have now been solved, allowing their application in various structure-based drug design methods. One potential method utilizing these crystal structures is the Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) method of predicting relative binding free.
Crystal structures of holo vitamin D receptor (VDR) revealed a canonical conformation in which the ligand is entrapped in a hydrophobic cavity buried in the ligand-binding domain (LBD). The mousetrap model postulates that helix 12 is positioned away from the domain to expose the interior cavity. However, the extended form of helix 12 is likely due to artifacts during crystallization. In this study, we set out to investigate conformational dynamics of apo VDR using molecular dynamics simulation on microsecond timescale. Here we show the neighboring backbones of helix 2-helix 3n and beta strand 2-helix 6 of LBD, instead of the helix 12, undergo large-scale motion, possibly gating the entrance of ligand to the ligand binding domain. Docking analysis to the simulated open structure of VDR with the estimated free energy of -37.0 kJ/mol, would emphasise the role of H2-H3n and S2-H6 in facilitating the entrance of calcitriol to the LBD of VDR.