G protein-coupled receptors (GPCRs) belong to the largest family of protein targets comprising over 800 members in which at least 500 members are the therapeutic targets. Among the GPCRs, G protein-coupled estrogen receptor-1 (GPER-1) has shown to have the ability in estrogen signaling. As GPER-1 plays a critical role in several physiological responses, GPER-1 has been considered as a potential therapeutic target to treat estrogen-based cancers and other non-communicable diseases. However, the progress in the understanding of GPER-1 structure and function is relatively slow due to the availability of a only a few selective GPER-1 modulators. As with many GPCRs, the X-ray crystal structure of GPER-1 is yet to be resolved and thus has led the researchers to search for new GPER-1 modulators using homology models of GPER-1. In this review, we aim to summarize various approaches used in the generation of GPER-1 homology model and their applications that have resulted in new GPER-1 ligands.
The disease Tuberculosis (TB) is caused by a bacterium called Mycobacterium tuberculosis (Mtb). The bacterial cell-wall consists of peptidoglycan layer maintains the cellular integrity and cell viability. The main problem resides in the cell cycle of Mycobacterium tuberculosis in its quiescent form which is not targeted by any drugs hence there is an immediate need for new antibiotics to target the cell wall. The current study deals with the dTDP-4-dehydrorahmnose reductase (RmlD) which is the final enzyme in the series of cell-wall proteins of Mtb. The RmlD is a part of Carbohydrate biosynthesis has been considered as a good drug target for the novel class of antibiotics. Our study begins with the protein structure prediction, Homology studies were conducted using the Phyre2 web server. The structure is then refined and subjected to molecular dynamics simulations for 50 ns using GROMACS. The clustering analysis has been carried out and generated 41 clusters with 2 Å as the cut-off. Blind docking virtual screening was performed against RmlD protein using the Super Natural-II database with AutoDock4.0. its results helped to screen top ligands based on best binding energies. In both dockings, there are some common residues in which the ligands are interacting and forming the Hydrogen bonds such as Asp-105, Val-158, Thr-160, Gly-161, Arg-224, Arg-256. The ligand-567 giving the best results by being in the top-3 of all the clusters in both blind docking as well as the active-site docking. Hence ligand-567 can be a potential inhibitor of RmlD which can further inhibit the cell-wall synthesis of Mycobacterium tuberculosis.Communicated by Ramaswamy H. Sarma.
Dengue virus (DV) infection is one of the main public health concerns, affecting approximately 390 million people worldwide, as reported by the World Health Organization. Yet, there is no antiviral treatment for DV infection. Therefore, the development of potent and nontoxic anti-DV, as a complement for the existing treatment strategies, is urgently needed. Herein, we investigate a series of small peptides inhibitors of DV antiviral activity targeting the entry process as the promising strategy to block DV infection. The peptides were designed based on our previously reported peptide sequence, DN58opt (TWWCFYFCRRHHPFWFFYRHN), to identify minimal effective inhibitory sequence through molecular docking and dynamics studies. The in silico designed peptides were synthesized using conventional Fmoc solid-phase peptide synthesis chemistry, purified by RP-HPLC and characterized using LCMS. Later, they were screened for their antiviral activity. One of the peptides, AC 001, was able to reduce about 40% of DV plaque formation. This observation correlates well with the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) analysis - AC 001 showed the most favorable binding affinity through 60 ns simulations. Pairwise residue decomposition analysis has revealed four key residues that contributed to the binding of these peptides into the DV2 E protein pocket. This work identifies the minimal peptide sequence required to inhibit DV replication and explains the behavior observed on an atomic level using computational study.Communicated by Ramaswamy H. Sarma.
The high dependency and surplus use of agrochemical products have liberated enormous quantities of toxic halogenated pollutants into the environment and threaten the well-being of humankind. Herein, this study performed molecular docking, molecular dynamic (MD) simulations, molecular mechanics-Poisson Boltzmann Surface Area (MM-PBSA) calculations on the DehH2 from Bacillus thuringiensis, to identify the order of which the enzyme degrades different substrates, haloacids, haloacetate and chlorpyrifos. The study discovered that the DehH2 favored the degradation of haloacids and haloacetates (-3.3 - 4.6 kcal/mol) and formed three hydrogen bonds with Asp125, Arg201 and Lys202. Despite the inconclusive molecular docking result, chlorpyrifos was consistently shown to be the least favored substrate of the DehH2 in MD simulations and MM-PBSA calculations. Results of MD simulations revealed the DehH2-haloacid- (RMSD 0.15 - 0.25 nm) and DehH2-haloacetates (RMSF 0.05 - 0.25 nm) were more stable, with the DehH2-L-2CP complex being the most stable while the least was the DehH2-chlorpyrifos (RMSD 0.295 nm; RMSF 0.05 - 0.59 nm). The Molecular Mechanics Poisson-Boltzmann Surface Area calculations showed the DehH2-L-2CP complex (-24.27 kcal/mol) having the lowest binding energy followed by DehH2-MCA (-22.78 kcal/mol), DehH2-D-2CP (-21.82 kcal/mol), DehH2-3CP (-21.11 kcal/mol), DehH2-2,2-DCP (-18.34 kcal/mol), DehH2-2,3-DCP (-8.34 kcal/mol), DehH2-TCA (-7.62 kcal/mol), while chlorpyrifos was unable to spontaneously bind to DehH2 (+127.16 kcal/mol). In a nutshell, the findings of this study offer valuable insights into the rational tailoring of the DehH2 for expanding its substrate specificity and catalytic activity in the near future.Communicated by Ramaswamy H. Sarma.
Cancer ranks in second place among the cause of death worldwide. Cancer progress in multiple stages of carcinogenesis and metastasis programs through complex pathways. Sex hormones and their receptors are the major factors in promoting cancer progression. Among them, G protein-coupled estrogen receptor-1 (GPER) has shown to mediate cellular signaling pathways and cancer cell proliferation. However, the lack of GPER protein structure limited the search for new modulators. In this study, we curated an extensive database of natural products to discover new potential GPER modulators. We used a combination of virtual screening techniques to generate a homology model of GPER and subsequently used that for the screening of 30,926 natural products from a public database to identify potential active modulators of GPER. The best hits were further screened through the ADMET filter and confirmed by docking analysis. Moreover, molecular dynamics simulations of best hits were also carried out to assess the stability of the ligand-GPER complex. This study predicted several potential GPER modulators with novel scaffolds that could be further investigated and used as the core for the development of novel GPER modulators.Communicated by Ramaswamy H. Sarma.
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.
We previously reported on a mutant lipase KV1 (Mut-LipKV1) from Acinetobacter haemolyticus which optimal pH was raised from 8.0 to 11.0 after triple substitutions of surface aspartic acid (Asp) with lysine (Lys). Herein, this study further examined the Mut-LipKV1 by molecular docking, molecular dynamics (MD) simulations and molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) calculations to explore the structural requirements that participated in the effective binding of tributyrin and its catalytic triad (Ser165, Asp259 and His289) and identify detailed changes that occurred post mutation. Mut-LipKV1 bound favorably with tributyrin (-4.1 kcal/mol) and formed a single hydrogen bond with His289, at pH 9.0. Despite the incongruent docking analysis data, results of MD simulations showed configurations of both the tributyrin-Mut-LipKV1 (RMSD 0.3 nm; RMSF 0.05 - 0.3 nm) and the tributyrin-wildtype lipase KV1 (tributyrin-LipKV1) complexes (RMSD 0.35 nm; RMSF 0.05 - 0.4 nm) being comparably stable at pH 8.0. MM-PBSA analysis indicated that van der Waals interactions made the most contribution during the molecular binding process, with the Mut-LipKV1-tributyrin complex (-44.04 kcal/mol) showing relatively lower binding energy than LipKV1-tributyrin (-43.83 kcal/mol), at pH 12.0. All tributyrin-Mut-LipKV1 complexes displayed improved binding free energies over a broader pH range from 8.0 - 12.0, as compared to LipKV1-tributyrin. Future empirical works are thus, important to validate the improved alkaline-stability of Mut-LipKV1. In a nutshell, our research offered a considerable insight for further improving the alkaline tolerance of lipases.Communicated by Ramaswamy H. Sarma.
Interaction behaviour of an anticancer drug, saracatinib (SCB) with human serum albumin (HSA), the major carrier protein in human blood circulation was investigated using fluorescence and absorption spectroscopy as well as computational methods. Analysis of the fluorescence quenching data along with absorption results confirmed the complex formation between SCB and HSA, based on the inverse correlation of the Stern-Volmer constant (KSV) with temperature and hyperchromic effect in the absorption spectra. Moderate binding affinity between SCB and HSA was evident from the binding constant, Ka value (1.08-0.74 × 104 M-1), while the SCB-HSA complexation was anticipated to be stabilized by hydrophobic and van der Waals interactions along with hydrogen bonds, as revealed from the thermodynamic data (ΔS = + 29.40 J mol-1 K-1 and ΔH = - 13.90 kJ mol-1). Addition of SCB to HSA significantly defended the thermal denaturation of the protein, though it perturbed the surrounding medium around Tyr and Trp residues. Site marker displacement results elucidated Sudlow's site I, positioned in subdomain IIA of HSA as the preferred binding site of SCB, which was well supported by molecular docking. Molecular dynamics simulation results suggested the stability of the SCB-HSA complex.Communicated by Ramaswamy H. Sarma.
Literature has shown that oil palm leaves (OPL) can be transformed into nanocellulose (NC) by fungal lignocellulosic enzymes, particularly those produced by the Trichoderma species. However, mechanism of β-glucosidase and xylanase selectivity to degrade lignin, hemicellulose and cellulose in OPL for NC production remains relatively vague. The study aimed to comprehend this aspect by an in silico approach of molecular docking, molecular dynamics (MD) simulation and Molecular-mechanics Poisson-Boltzmann surface area (MM-PBSA) analysis, to compare interactions between the β-glucosidase- and xylanase from Trichoderma asperellum UC1 in complex with each substrate. Molecular docking of the enzyme-substrate complex showed residues Glu165-Asp226-Glu423 and Arg155-Glu210-Ser160 being the likely catalytic residues of β-glucosidase and xylanase, respectively. The binding affinity of β-glucosidase for the substrates are as follows: cellulose (-8.1 kcal mol-1) > lignin (-7.9 kcal mol-1) > hemicellulose (-7.8 kcal mol-1), whereas, xylanase showed a corresponding preference for; hemicellulose (-6.7 kcal mol-1) > cellulose (-5.8 kcal mol-1) > lignin (-5.7 kcal mol-1). Selectivity of both enzymes was reiterated by MD simulations where interactions between β-glucosidase-cellulose and xylanase-hemicellulose were the strongest. Notably low free-binding energy (ΔGbind) of β-glucosidase and xylanase in complex with cellulose (-207.23 +/- 47.13 kJ/mol) and hemicellulose (-131.48 +/- 24.57 kJ/mol) were observed, respectively. The findings thus successfully identified the cellulose component selectivity of the polymer-acting β-glucosidase and xylanase of T. asperellum UC1.Communicated by Ramaswamy H. Sarma.
Dehalogenases are of high interest due to their potential applications in bioremediation and in synthesis of various industrial products. DehL is an L-2-haloacid dehalogenase (EC 3.8.1.2) that catalyses the cleavage of halide ion from L-2-halocarboxylic acid to produce D-2-hydroxycarboxylic acid. Although DehL utilises the same substrates as the other L-2-haloacid dehalogenases, its deduced amino acid sequence is substantially different (<25%) from those of the rest L-2-haloacid dehalogenases. To date, the 3D structure of DehL is not available. This limits the detailed understanding of the enzyme's reaction mechanism. The present work predicted the first homology-based model of DehL and defined its active site. The monomeric unit of the DehL constitutes α/β structure that is organised into two distinct structural domains: main and subdomains. Despite the sequence disparity between the DehL and other L-2-haloacid dehalogenases, its structural model share similar fold as the experimentally solved L-DEX and DehlB structures. The findings of the present work will play a crucial role in elucidating the molecular details of the DehL functional mechanism.
Multiple spectroscopic techniques, such as fluorescence, absorption, and circular dichroism along with in silico studies were used to characterize the binding of a potent inhibitor molecule, CCG1423 to the major transport protein, human serum albumin (HSA). Fluorescence and absorption spectroscopic results confirmed CCG1423-HSA complex formation. A strong binding affinity stabilized the CCG1423-HSA complex, as evident from the values of the binding constant (Ka = 1.35 × 106-5.43 × 105 M-1). The KSV values for CCG1423-HSA system were inversely correlated with temperature, suggesting the involvement of static quenching mechanism. Thermodynamic data anticipated that CCG1423-HSA complexation was mainly driven by hydrophobic and van der Waals forces as well as hydrogen bonds. In silico analysis also supported these results. Three-dimensional fluorescence and circular dichroism spectral analysis suggested microenvironmental perturbations around protein fluorophores and structural (secondary and tertiary) changes in the protein upon CCG1423 binding. CCG1423 binding to HSA also showed some protection against thermal denaturation. Site-specific marker-induced displacement results revealed CCG1423 binding to Sudlow's site I of HSA, which was also confirmed by the computational results. A few common ions were also found to interfere with the CCG1423-HSA interaction.
A series of dimeric naphthoquinones containing natural 2-hydroxy-1-4-naphthoquinone moiety was designed, synthesized, and evaluated against neuraminidase of H5N1 virus. p-hydroxy derivatives showed higher inhibition when compared to p-halogenated compounds. Molecular docking studies conducted with H5N1 neuraminidase clearly demonstrated different binding modes of the most active compound onto the open and closed conformations of loop 150. The results thus provide not only evidences of a novel scaffold evaluated as inhibitor, but also a rational explanation involving molecular modeling and the role of loop 150 in the binding.
Aromatherapy with essential oils (EOs) has been linked to improvement of cognitive function in patients with dementia. In order to act systemically, active EO components must be absorbed through the skin, enter the systemic circulation, and cross the blood brain barrier (BBB). Thus, the aim of this work was to develop quantitative structure activity relationships (QSARs), to predict skin and blood barrier penetrative abilities of 119 terpenoids from EOs used in aromatherapy. The first model was based on experimentally measured skin permeability for 162 molecules, and the second model on BBB permeability for 138 molecules. Each molecule was encoded with 63 calculated molecular descriptors and an artificial neural network was used to correlate molecular descriptors to permeabilities. Developed QSAR models confirm that EOs components penetrate through the skin and across the BBB. Some well-known descriptors, such as log P (lipophilicity), molecular size and shape, dominated the QSAR model for BBB permeability. Compounds with the highest predicted BBB penetration were hydrocarbon terpenes with the smallest molecular size and highest lipophilicity. Thus, molecular size is a limiting factor for penetration. Compounds with the highest skin permeability have slightly higher molecular size, high lipophilicity and low polarity. Our work shows that a major disadvantage of novel multitarget compounds developed for the treatment of Alzheimer's disease is the size of molecules, which cause problems in their delivery to the brain. Therefore, there is a need for smaller compounds, which possess more desirable physicochemical properties and pharmacokinetics, in addition to targeted biological effects.Communicated by Ramaswamy H. Sarma.
Molecular interaction of the 3,4-methylenedioxy-β-nitrostyrene (MNS), an inhibitor of platelet aggregation with the main transport protein, albumin from human serum (HSA) was explored using absorption, fluorescence and circular dichroism (CD) spectroscopy in combination with in silico analyses. The MNS-HSA complexation was corroborated from the fluorescence and absorption spectral results. Implication of static quenching mechanism for MNS-HSA system was predicted from the Stern-Volmer constant, KSV-temperature relationship as well as the bimolecular quenching rate constant, kq values. Stabilization of the complex was affirmed by the value of the binding constant (Ka = 0.56-1.48 × 104 M-1). Thermodynamic data revealed that the MNS-HSA association was spontaneously driven mainly through hydrophobic interactions along with van der Waal's interaction and H-bonds. These results were well supported by in silico interpretations. Far-UV and near-UV CD spectral results manifested small variations in the protein's secondary and tertiary structures, respectively, while three-dimensional fluorescence spectra displayed microenvironmental fluctuations around protein's fluorophores, upon MNS binding. Significant improvement in the protein's thermostability was evident from the temperature-stability results of MNS-bound HSA. Binding locus of MNS, as identified by competitive drug displacement findings as well as in silico analysis, was found to be located in subdomain IIA (Sudlow's site I) of the protein.Communicated by Ramaswamy H. Sarma.
Interests in Acinetobacter haemolyticus lipases are showing an increasing trend concomitant with growth of the enzyme industry and the widening search for novel enzymes and applications. Here, we present a structural model that reveals the key catalytic residues of lipase KV1 from A. haemolyticus. Homology modeling of the lipase structure was based on the structure of a carboxylesterase from the archaeon Archaeoglobus fulgidus as the template, which has a sequence that is 58% identical to that of lipase KV1. The lipase KV1 model is comprised of a single compact domain consisting of seven parallel and one anti-parallel β-strand surrounded by nine α-helices. Three structurally conserved active-site residues, Ser165, Asp259, and His289, and a tunnel through which substrates access the binding site were identified. Docking of the substrates tributyrin and palmitic acid into the pH 8 modeled lipase KV1 active sites revealed an aromatic platform responsible for the substrate recognition and preference toward tributyrin. The resulting binding modes from the docking simulation correlated well with the experimentally determined hydrolysis pattern, for which pH 8 and tributyrin being the optimum pH and preferred substrate. The results reported herein provide useful insights into future structure-based tailoring of lipase KV1 to modulate its catalytic activity.
G protein-coupled receptors (GPCRs) are proteins of pharmaceutical importance, with over 30% of all drugs in clinical use targeting them. Increasing numbers of X-ray crystal (XRC) structures of GPCRs offer a wealth of data relating to ligand binding. For the β-adrenoceptors (β-ARs), XRC structures are available for human β2- and turkey β1-subtypes, in complexes with a range of ligands. While these structures provide insight into the origins of ligand structure-activity relationships (SARs), questions remain. The ligands in all published complexed XRC structures lack extensive substitution, with no obvious way the ligand-binding site can accommodate β1-AR-selective antagonists with extended side-chains para- to the common aryloxypropanolamine pharmacophore. Using standard computational docking tools with such ligands generally returns poses that fail to explain known SARs. Application of our Active Site Pressurization (ASP) modelling method to β-AR XRC structures and homology models however, reveals a dynamic area in the ligand-binding pocket that, through minor changes in amino acid side chain orientations, opens a fissure between transmembrane (TM) helices H4 and H5, exposing intra-membrane space. This fissure, which we term the 'keyhole', is ideally located to accommodate extended moieties present in many high-affinity β1-AR-selective ligands; allowing the rest of the ligand structure to adopt a canonical pose in the orthosteric binding site. We propose the keyhole may be a feature of both β1- and β2-ARs, but that subtle structural differences exist between the two, contributing to subtype-selectivity. This has consequences for the rational design of future generations of subtype-selective ligands for these therapeutically important targets.
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.
Medicinal plants as rich sources of bioactive compounds are now being explored for drug development against COVID-19. 19 medicinal plants known to exhibit antiviral and anti-inflammatory effects were manually curated, procuring a library of 521 metabolites; this was virtually screened against NSP9, including some other viral and host targets and were evaluated for polypharmacological indications. Leads were identified via rigorous scoring thresholds and ADMET filtering. MM-GBSA calculation was deployed to select NSP9-Lead complexes and the complexes were evaluated for their stability and protein-ligand communication via MD simulation. We identified 5 phytochemical leads for NSP9, 23 for Furin, 18 for ORF3a, and 19 for IL-6. Ochnaflavone and Licoflavone B, obtained from Lonicera japonica (Japanese Honeysuckle) and Glycyrrhiza glabra (Licorice), respectively, were identified to have the highest potential polypharmacological properties for the aforementioned targets and may act on multiple pathways simultaneously to inhibit viral entry, replication, and disease progression. Additionally, MD simulation supports the robust stability of Ochnaflavone and Licoflavone B against NSP9 at the active sites via hydrophobic interactions, H-bonding, and H-bonding facilitated by water. This study promotes the initiation of further experimental analysis of natural product-based anti-COVID-19 therapeutics.