Development of new antimalarial drugs continues to be of huge importance because of the resistance of malarial parasite towards currently used drugs. Due to the reliance of parasite on glycolysis for ...energy generation, glycolytic enzymes have played important role as potential targets for the development of new drugs. Plasmodium falciparum lactate dehydrogenase (PfLDH) is a key enzyme for energy generation of malarial parasites and is considered to be a potential antimalarial target. Presently, there are nearly 15 crystal structures bound with inhibitors and substrate that are available in the protein data bank (PDB). In the present work, we attempted to consider multiple crystal structures with bound inhibitors showing affinity in the range of 1.4 × 10
2
-1.3 × 10
6
nM efficacy and optimized the pharmacophore based on the energy involved in binding termed as e-pharmacophore mapping. A high throughput virtual screening (HTVS) combined with molecular docking, ADME predictions and molecular dynamics simulation led to the identification of 20 potential compounds which could be further developed as novel inhibitors for PfLDH.
Human (Hu) familial prion diseases are associated with about 40 point mutations of the gene coding for the prion protein (PrP). Most of the variants associated with these mutations are located in the ...globular domain of the protein. We performed 50 ns of molecular dynamics for each of these mutants to investigate their structure in aqueous solution. Overall, 1.6 μs of molecular dynamics data is presented. The calculations are based on the AMBER(parm99) force field, which has been shown to reproduce very accurately the structural features of the HuPrP wild type and a few variants for which experimental structural information is available. The variants present structural determinants different from those of wild-type HuPrP and the protective mutation HuPrP(E219K-129M). These include the loss of salt bridges in α2–α3 regions and the loss of π-stacking interactions in the β2–α2 loop. In addition, in the majority of the mutants, the α3 helix is more flexible and Y169 is more solvent exposed. The presence of similar traits in this large spectrum of mutations hints to a role of these fingerprints in their known disease-causing properties. Overall, the regions most affected by disease-linked mutations in terms of structure and/or flexibility are those involved in the pathogenic conversion to the scrapie form of the protein and in the interaction with cellular partners. These regions thus emerge as optimal targets for antibody- and ligand-binding studies.
Display omitted
Singleton-Merten syndrome, a critical and rare multifactorial disorder that is closely linked to R516Q mutation in MDA5 protein associated with an enhanced interferon response in the affected ...individual. In the present study, we provide conclusive key evidence on R516Q mutation and their connectivity towards sequence-structural basis dysfunction of MDA5 protein. Among the various mutations, we found R516Q is the most pathogenic mutation based on mutational signature Q-A-RE-G-R-GA-R-A-ED-DE-S-ST-Y-TSAV-L-V designed from our work. Further, we derived a distant ortholog for this mutational signature from which we identified 343 intra-residue interactions that fall communally in the position required to maintain the structural and functional integration of protein architecture. This identification served us to understand the critical role of hot spots in residual aggregation that holds a native form of folding conformation in the functional region. In addition, the long-range molecular dynamics simulation demarcated the residual dependencies of conformational transition in distinct regions (
360-370
,
380-410
,
430-480
and
510-520
) occurring upon R516Q mutation. Together, our results emphasise that the dislocation of functional hot spots Pro229, Arg414, Val498, Met510, Ala513, Gly515 and Arg516 in MDA5 protein which is important for interior structural packing and fold arrangements. In a nutshell, our findings are perfectly conceded with other experimental reports and will have potential implications in immune therapeutical advancement for rare singleton-merten syndrome.
RDL receptors are GABA-activated inhibitory Cys-loop receptors found throughout the insect CNS. They are a key target for insecticides. Here, we characterize the GABA binding site in RDL receptors ...using computational and electrophysiological techniques. A homology model of the extracellular domain of RDL was generated and GABA docked into the binding site. Molecular dynamics simulations predicted critical GABA binding interactions with aromatic residues F206, Y254, and Y109 and hydrophilic residues E204, S176, R111, R166, S176, and T251. These residues were mutated, expressed in Xenopus oocytes, and their functions assessed using electrophysiology. The data support the binding mechanism provided by the simulations, which predict that GABA forms many interactions with binding site residues, the most significant of which are cation-π interactions with F206 and Y254, H-bonds with E204, S205, R111, S176, T251, and ionic interactions with R111 and E204. These findings clarify the roles of a range of residues in binding GABA in the RDL receptor, and also show that molecular dynamics simulations are a useful tool to identify specific interactions in Cys-loop receptors.
Prostate cancer (PC) is one of the leading cancers in men, raising a serious health issue worldwide. Due to lack of suitable biomarker, their inhibitors and the platform for testing those inhibitors ...result in poor prognosis of PC. AMP-activated protein kinase (AMPK) is a highly conserved protein kinase found in eukaryotes that is involved in growth and development, and also acts as a therapeutic target for PC. The aim of the present study is to identify novel potent inhibitors of AMPK and propose a simple cellular model system for understanding its biology. Structural modelling and MD simulations were performed to construct and refine the 3D models of Dictyostelium and human AMPK. Binding mechanisms of different drug compounds were studied by performing molecular docking, molecular dynamics and MM-PBSA methods. Two novel drugs were isolated having higher binding affinity over the known drugs and hydrophobic forces that played a key role during protein-ligand interactions. The study also explored the simple cellular model system for drug screening and understanding the biology of a therapeutic target by performing in vitro experiments.
Display omitted
•Ferulic Acid Decarboxylase (FADase) catalyzes the conversion of ferulic acid to 4-vinylguaiacol.•Sequential site directed mutations were introduced on active site of FADase and ...validated with DCCM, PCA and MD simulation.•Docking studies were performed between ferulic acid and different mutants of FADase.•Best mutant-ligand complex was validated through stable energy optimization.•New mutants were reported with higher affinity towards and more hydrogen bonds formation.
Microbial Ferulic Acid Decarboxylase (FADase) catalyses the conversion of ferulic acid to 4-hydroxy-3-methoxystyrene (4-vinylguaiacol) via non-oxidative decarboxylation. In this article, we present a computational, three-dimensional structural and functional analysis of FADase from Enterobacter sp. P × 6-4 (3NX1) which can be used to generate enhanced bindings of substrates. The enzymatic catalytic site and binding sites have been critically evaluated. Sequential site directed mutations on enzyme have also been introduced for formation of a greater number of hydrogen bonds. Four mutants were generated based on our hypothesis. Active sites of mutated FADases have been analyzed with dynamic cross-correlation maps and principle components analysis. All structures were validated and optimized through energy minimization. Docking studies were also carried out between ferulic acid and different mutated enzymes. The protein (wild and mutants) complexes were further validated with molecular simulation. Mutant3 was found to have better affinity towards ferulic acid. Mutant3 also forms a higher number of hydrogen bonds with the substrate to facilitate greater interaction. This current work will help industry to create new and novel mutants to produce vanillin.
Lipoic acid synthase (LIAS) is an iron-sulfur cluster mitochondrial enzyme which catalyzes the final step in the de novo pathway for the biosynthesis of lipoic acid, a potent antioxidant. Recently ...there has been significant interest in its role in metabolic diseases and its deficiency in LIAS expression has been linked to conditions such as diabetes, atherosclerosis and neonatal-onset epilepsy, suggesting a strong inverse correlation between LIAS reduction and disease status. In this study we use a bioinformatics approach to predict its structure, which would be helpful to understanding its role. A homology model for LIAS protein was generated using X-ray crystallographic structure of Thermosynechococcus elongatus BP-1 (PDB ID: 4U0P). The predicted structure has 93% of the residues in the most favour region of Ramachandran plot. The active site of LIAS protein was mapped and docked with S-Adenosyl Methionine (SAM) using GOLD software. The LIAS–SAM complex was further refined using molecular dynamics simulation within the subsite 1 and subsite 3 of the active site. To the best of our knowledge, this is the first study to report a reliable homology model of LIAS protein. This study will facilitate a better understanding mode of action of the enzyme-substrate complex for future studies in designing drugs that can target LIAS protein.
Streptokinase (SK), a plasminogen activator (PA) that converts inactive plasminogen (Pg) to plasmin (Pm), is a protein secreted by groups A, C, and G streptococci (GAS, GCS, and GGS, respectively), ...with high sequence divergence and functional heterogeneity. While roles of some residual changes in altered SK functionality are shown, the underlying structural mechanisms are less known. Herein, using computational approaches, we analyzed the conformational basis for the increased activity of SK from a GGS (SKG132) isolate with four natural residual substitutions (Ile33Phe, Arg45Gln, Asn228Lys, Phe287Ile) compared to the standard GCS (SKC). Using the crystal structure of SK.Pm catalytic complex as main template SKC.μPm catalytic complex was modeled through homology modeling process and validated by several online validation servers. Subsequently, SKG132.μPm structure was constructed by altering the corresponding residual substitutions. Results of three independent MD simulations showed increased RMSF values for SKG132.μPm, indicating the enhanced structural flexibility compared to SKC.μPm, specially in 170 and 250 loops and three regions: R1 (149-161), R2 (182-215) and R3 (224-229). In parallel, the average number of Hydrogen bonds in 170 loop, R2 and R3 (especially for Asn228Lys) of SKG132 compared to that of the SKC was decreased. Accordingly, residue interaction networks (RINs) analyses indicated that Asn228Lys might induce more level of structural flexibility by generation of free Lys256, while Phe287Ile and Ile33Phe enhanced the stabilization of the SKG132.μPm catalytic complex. These results denoted the potential role of the optimal dynamic state and stabilized catalytic complex for increased PA potencies of SK as a thrombolytic drug.
In standard implementations of computational protein design, a positive-design approach is used to predict sequences that will be stable on a given backbone structure. Possible competing states are ...typically not considered, primarily because appropriate structural models are not available. One potential competing state, the domain-swapped dimer, is especially compelling because it is often nearly identical with its monomeric counterpart, differing by just a few mutations in a hinge region. Molecular dynamics (MD) simulations provide a computational method to sample different conformational states of a structure. Here, we tested whether MD simulations could be used as a post-design screening tool to identify sequence mutations leading to domain-swapped dimers. We hypothesized that a successful computationally designed sequence would have backbone structure and dynamics characteristics similar to that of the input structure and that, in contrast, domain-swapped dimers would exhibit increased backbone flexibility and/or altered structure in the hinge-loop region to accommodate the large conformational change required for domain swapping. While attempting to engineer a homodimer from a 51-amino-acid fragment of the monomeric protein engrailed homeodomain (ENH), we had instead generated a domain-swapped dimer (ENH_DsD). MD simulations on these proteins showed increased B-factors derived from MD simulation in the hinge loop of the ENH_DsD domain-swapped dimer relative to monomeric ENH. Two point mutants of ENH_DsD designed to recover the monomeric fold were then tested with an MD simulation protocol. The MD simulations suggested that one of these mutants would adopt the target monomeric structure, which was subsequently confirmed by X-ray crystallography.
Display omitted
•Computational protein design (CPD) calculations that do not consider competing states may lead to off-target folding.•We developed an MD simulation protocol as a post-CPD screening tool.•The MD protocol identifies mutations leading to undesired competing states.•The MD protocol predicts mutations that favor the target fold.•CPD combined with MD screening can greatly improve design success rates.
Hydroxylated polybrominated diphenyl ethers (OH-PBDEs) are a class of toxic environmental pollutants that are persistent, bioaccumulative, and difficult to degrade. Their structure is very similar to ...the thyroid hormone (T4) and uses the body's thyroid transporter (TTR) binding to interfere with the endocrine balance, disrupting the body's normal physiological activity. According to Fourier transform infrared spectroscopy and dynamics simulation of do_dssp module analysis, there are three kinds of OH-PBDEs that can induce TTR secondary structural changes. Fluorescence spectra and UV-Vis spectra show that for the three kinds of OH-PBDEs for TTR, the main methods of quenching are static quenching and non-radiative energy transfer. According to thermodynamic analysis, ΔG < 0, ΔH > 0, and ΔS > 0 combine to show that the hydrophobic interaction is the main driving force of the combination. From the molecular docking analysis, it was found that 4′-hydroxy-2,2′,4,5′- tetrabromodiphenyl ether (4′-OH-BDE49) and 4 hydroxy-2,2′,3,4′,5,6,6′- heptabromodiphenyl ether (4-OH-BDE188) had a cationic-π interaction with TTR, whereas 4 hydroxy-2,2′,3,4,5,5′,6- heptabromodiphenyl ether (4-OH-BDE187) was bonded to TTR by hydrogen bonds to form stable complexes. In this paper, we highlight the consistency of spectroscopic experiments and computer simulations so as to provide a reliable analytical method for the toxicological properties of small molecule contaminants.