Molecular dynamics simulations suggest ligand's binding to nicotinamidase/pyrazinamidase

Molecular pathway and structural mechanism of human oncochannel TRPV6 inhibition by the phytocannabinoid tetrahydrocannabivarin

The calcium-selective oncochannel TRPV6 is an important driver of cell proliferation in human cancers. Despite increasing interest of pharmacological research in developing synthetic inhibitors of TRPV6, natural compounds acting at this channel have been largely neglected. On the other hand, pharmacokinetics of natural small-molecule antagonists optimized by nature throughout evolution endows these compounds with a medicinal potential to serve as potent and safe next-generation anti-cancer drugs. Here we report the structure of human TRPV6 in complex with tetrahydrocannabivarin (THCV), a natural cannabinoid inhibitor extracted from Cannabis sativa. We use cryo-electron microscopy combined with electrophysiology, calcium imaging, mutagenesis, and molecular dynamics simulations to identify THCV binding sites in the portals that connect the membrane environment surrounding the protein to the central cavity of the channel pore and to characterize the allosteric mechanism of TRPV6 inhibition. We also propose the molecular pathway taken by THCV to reach its binding site. Our study provides a foundation for the development of new TRPV6-targeting drugs.

Binding characteristics of staphylococcal protein A and streptococcal protein G for fragment crystallizable portion of human immunoglobulin G

In the wide array of physiological processes, protein-protein interactions and their binding are the most basal activities for achieving adequate biological metabolism. Among the studies on binding proteins, the examination of interactions between immunoglobulin G (IgG) and natural immunoglobulin-binding ligands, such as staphylococcal protein A (spA) and streptococcal protein G (spG), is essential in the development of pharmaceutical science, biotechnology, and affinity chromatography. The widespread utilization of IgG-spA/spG binding characteristics has allowed researchers to investigate these molecular interactions. However, the detailed binding strength of each ligand and the corresponding binding mechanisms have yet to be fully investigated. In this study, the authors analyzed the binding strengths of IgG-spA and IgG-spG complexes and identified the mechanisms enabling these bindings using molecular dynamics simulation, steered molecular dynamics, and advanced Poisson-Boltzmann Solver simulations. Based on the presented data, the binding strength of the spA ligand was found to significantly exceed that of the spG ligand. To find out which non-covalent interactions or amino acid sites have a dominant role in the tight binding of these ligands, further detailed analyses of electrostatic interactions, hydrophobic bonding, and binding free energies have been performed. In investigating their binding affinity, a relatively independent and different unbinding mechanism was found in each ligand. These distinctly different mechanisms were observed to be highly correlated to the protein secondary and tertiary structures of spA and spG ligands, as explicated from the perspective of hydrogen bonding.

Molecular dynamics simulation and free energy calculation studies of Coagulin L as dipeptidyl peptidase-4 inhibitor

Plant derived product can be used as other alternatives to currently used drugs for controlling chronic diseases like Diabetes mellitus. The potential of Coagulin L (a constituent of Withania coagulans) as dipeptidyl peptidase-4 (DPP-4) inhibitor was evaluated by molecular modelling study. It was observed that amino acid residues such as Glu205, Glu206, Tyr 547, His 740, and Try662 interacts with Coagulin L and Saxagliptin (a known DPP-4 inhibitor). Other nonbonded interactions of Coagulin L and Saxagliptin with DPP-4 binding residues were also found similar. The docking energy of Coagulin L was found to be -7.69 Kcal/mol whereas -8.44 kcal/mol was recorded for Saxagliptin. MD simulation study revealed stable binding throughout 100 ns simulation. RMSD plot of the complex was stabilized in 43 ns and remained stable during entire simulation(100 ns). RMSF plot of DPP-4 Coagulin L interaction showed major fluctuations at residue 246 and 766, however, Arg 125, Glu 205, Ser 209 and His 740 showed no major perturbations. Principal Component Analysis showed that important dynamics of the protein remain unchanged during entire simulation since the non-polar, van der waals, ionic interaction and solvation energy, altogether play important role in the complex stability. The molecular modelling study of DPP-4 with Coagulin L was an effort to establish correlation with traditional practices of Withania coagulans as antidiabetic agent in Indian subcontinent.Communicated by Ramaswamy H. Sarma.

Applications of Molecular Simulation in the Discovery of Antituberculosis Drugs: A Review

After decades of efforts, tuberculosis has been well controlled in most places. The existing drugs are no longer sufficient for the treatment of drug-resistant Mycobacterium tuberculosis due to significant toxicity and selective pressure, especially for XDR-TB. In order to accelerate the development of high-efficiency, low-toxic antituberculosis drugs, it is particularly important to use Computer Aided Drug Design (CADD) for rational drug design. Here, we systematically reviewed the specific role of molecular simulation in the discovery of new antituberculosis drugs. The purpose of this review is to overview current applications of molecular simulation methods in the discovery of antituberculosis drugs. Furthermore, the unique advantages of molecular simulation was discussed in revealing the mechanism of drug resistance. The comprehensive use of different molecular simulation methods will help reveal the mechanism of drug resistance and improve the efficiency of rational drug design. With the help of molecular simulation methods such as QM/MM method, the mechanisms of biochemical reactions catalyzed by enzymes at atomic level in Mycobacterium tuberculosis has been deeply analyzed. QSAR and virtual screening both accelerate the development of highefficiency, low-toxic potential antituberculosis drugs. Improving the accuracy of existing algorithms and developing more efficient new methods for CADD will always be a hot topic in the future. It is of great value to utilize molecular dynamics simulation to investigate complex systems that cannot be studied in experiments, especially for drug resistance of Mycobacterium tuberculosis.

MD Simulation Investigation on the Binding Process of Smoke-Derived Germination Stimulants to Its Receptor

Karrikins (KARs) are a class of smoke-derived seed germination stimulants with great significance in both agriculture and plant biology. By means of direct binding to the receptor protein KAI2, the compounds can initiate the KAR signal transduction pathway, hence triggering germination of the dormant seeds in the soil. In the research, several molecular dynamics (MD) simulation techniques were properly integrated to investigate the binding process of KAR₁ to KAI2 and reveal the details of the whole binding event. The calculated binding free energy, -7.00 kcal/mol, is in good agreement with the experimental measurement, -6.83 kcal/mol. The obtained PMF profile indicates the existence of three intermediate states in the binding process. The analysis of the simulation trajectories demonstrates that, in the intermediate structures, KAR₁ is stabilized by some hydrophobic residues (Phe26, Phe134, Leu142, Trp153, Phe157, Leu160, Phe194), along with several bridging water molecules, and meanwhile, the significant shifting occurs in the local conformation of the protein as the ligand's binding. A series of the residues (Gln141-Phe157) on the so-called "cap domain" are proposed to be responsible for capturing the ligand at the initial stage of the binding. Besides, the changes of the ligand's poses are also quantitatively characterized by the proper choice of the coordinate system. Our work will contribute to the more penetrating understanding of the ligand binding process and the receptor affinity difference between several members in the KAR family and help design new, more effective germination stimulants.

Computational insights into pH-dependence of structure and dynamics of pyrazinamidase: A comparison of wild type and mutants

The mycobacterial enzyme pyrazinamidase (PZase) is the target of key tuberculosis drug, pyrazinamide. Mutations in PZase cause drug resistance. Herein, three point mutations, W68G, L85P, and V155G, were investigated through over 8 µs of molecular dynamics simulations coupled with essential dynamics and binding pocket analysis at neutral (pH = 7) and acidic (pH = 4) ambient conditions. The 51-71 flap region exhibited drastic displacement leading to enlargement of binding cavity, especially at the lower pH. Accessibility of solvent to the active site of the mutant enzymes was also reduced. The protonation of key surface residues at low pH results in more contribution of these residues to structural stability and integrity of the enzyme and reduced interactions with solvent molecules, which acts as a cage, keeping the enzyme together. The observed results suggest a pattern of structural alterations due to point mutations in PZase, which is consistent with other experimental and theoretical investigations and, can be harnessed for drug design purposes.

In silico approaches for investigating the binding propensity of apigenin and luteolin against class I HDAC isoforms

AIM

Aberrant activity of class I histone deacetylases (HDACs) has strong implications for various cancers. Targeting these HDACs with synthetic HDAC inhibitors has shown significant side effects such as atrial fibrillation and QT prolongation emphasizing the need of natural inhibitors as substitutes to synthetic ones.

RESULTS

The binding propensity of the two plant-derived inhibitors apigenin and luteolin towards class I HDAC isoforms was checked using extra-precision molecular docking and implicit solvation MMGBSA. Apigenin showed a superior binding affinity against these isoforms as compared to luteolin. Both inhibitors docked stable to the binding pocket of these HDACs as determined by molecular dynamics simulation study.

CONCLUSION

Apigenin and luteolin may serve as substitutes to synthetic inhibitors for effective HDAC based anticancer therapy.

Unbinding pathway energy of glyphosate from the EPSPs enzyme binding site characterized by Steered Molecular Dynamics and Potential of Mean Force

The quantification of herbicides in the environment, like glyphosate, is extremely important to prevent contamination. Nanobiosensors stands out in the quantization process, because of the high selectivity, sensitivity and short response time of the method. In order to emulate the detection of glyphosate using a specific nanobiossensor through an Atomic Force Microscope (AFM), this work carried out Steered Molecular Dynamics simulations (SMD) in which the herbicide was unbinded from the active site of the enzyme 5- enolpyruvylshikimate 3 phosphate synthase (EPSPS) along three different directions.After the simulations, Potential of Mean Force calculations were carried, from a cumulant expansion of Jarzynski's equation to obtain the profile of free energy of interaction between the herbicide and the active site of the enzyme in the presence of shikimate-3 substrate phosphate (S3P). The set of values for external work, had a Gaussian distribution. The PMF values ranged according to the directions of the unbindong pahway of each simulation, displaying energy values of 10.7, 14.7 and 19.5KJmol-1. The results provide a theoretical support in order to assist the construction of a specific nanobiossensor to quantify the glyphosate herbicide.

Product Release Pathways in Human and Plasmodium falciparum Phosphoribosyltransferase

Atomistic molecular dynamics (MD) simulations coupled with the metadynamics technique were carried out to delineate the product (PPi.2Mg and IMP) release mechanisms from the active site of both human (Hs) and Plasmodium falciparum (Pf) hypoxanthine-guanine-(xanthine) phosphoribosyltransferase (HG(X)PRT). An early movement of PPi.2Mg from its binding site has been observed. The swinging motion of the Asp side chain (D134/D145) in the binding pocket facilitates the detachment of IMP, which triggers the opening of flexible loop II, the gateway to the bulk solvent. In PfHGXPRT, PPi.2Mg and IMP are seen to be released via the same path in all of the biased MD simulations. In HsHGPRT too, the product molecules follow similar routes from the active site; however, an alternate but minor escape route for PPi.2Mg has been observed in the human enzyme. Tyr 104 and Phe 186 in HsHGPRT and Tyr 116 and Phe 197 in PfHGXPRT are the key residues that mediate the release of IMP, whereas the motion of PPi.2Mg away from the reaction center is guided by the negatively charged Asp and Glu and a few positively charged residues (Lys and Arg) that line the product release channels. Mutations of a few key residues present in loop II of Trypanosoma cruzi (Tc) HGPRT have been shown to reduce the catalytic efficiency of the enzyme. Herein, in silico mutation of corresponding residues in loop II of HsHGPRT and PfHGXPRT resulted in partial opening of the flexible loop (loop II), thus exposing the active site to bulk water, which offers a rationale for the reduced catalytic activity of these two mutant enzymes. Investigations of the product release from these HsHGPRT and PfHGXPRT mutants delineate the role of these important residues in the enzymatic turnover.

Ranking of Molecular Biomarker Interaction with Targeted DNA Nucleobases via Full Atomistic Molecular Dynamics

DNA-based sensors can detect disease biomarkers, including acetone and ethanol for diabetes and H2S for cardiovascular diseases. Before experimenting on thousands of potential DNA segments, we conduct full atomistic steered molecular dynamics (SMD) simulations to screen the interactions between different DNA sequences with targeted molecules to rank the nucleobase sensing performance. We study and rank the strength of interaction between four single DNA nucleotides (Adenine (A), Guanine (G), Cytosine (C), and Thymine (T)) on single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA) with acetone, ethanol, H2S and HCl. By sampling forward and reverse interaction paths, we compute the free-energy profiles of eight systems for the four targeted molecules. We find that dsDNA react differently than ssDNA to the targeted molecules, requiring more energy to move the molecule close to DNA as indicated by the potential of mean force (PMF). Comparing the PMF values of different systems, we obtain a relative ranking of DNA base for the detection of each molecule. Via the same procedure, we could generate a library of DNA sequences for the detection of a wide range of chemicals. A DNA sensor array built with selected sequences differentiating many disease biomarkers can be used in disease diagnosis and monitoring.

Molecular simulation assisted identification of Ca(2+) binding residues in TMEM16A

Calcium-activated chloride channels (CaCCs) play vital roles in a variety of physiological processes. Transmembrane protein 16A (TMEM16A) has been confirmed as the molecular counterpart of CaCCs which greatly pushes the molecular insights of CaCCs forward. However, the detailed mechanism of Ca(2+) binding and activating the channel is still obscure. Here, we utilized a combination of computational and electrophysiological approaches to discern the molecular mechanism by which Ca(2+) regulates the gating of TMEM16A channels. The simulation results show that the first intracellular loop serves as a Ca(2+) binding site including D439, E444 and E447. The experimental results indicate that a novel residue, E447, plays key role in Ca(2+) binding. Compared with WT TMEM16A, E447Y produces a 30-fold increase in EC50 of Ca(2+) activation and leads to a 100-fold increase in Ca(2+) concentrations that is needed to fully activate the channel. The following steered molecular dynamic (SMD) simulation data suggests that the mutations at 447 reduce the Ca(2+) dissociation energy. Our results indicated that both the electrical property and the size of the side-chain at residue 447 have significant effects on Ca(2+) dependent gating of TMEM16A.

Exploration of the chlorpyrifos escape pathway from acylpeptide hydrolases using steered molecular dynamics simulations

Acylpeptide hydrolases (APH) catalyze the removal of an N-acylated amino acid from blocked peptides. APH is significantly more sensitive than acetylcholinesterase, a target of Alzheimer's disease, to inhibition by organophosphorus (OP) compounds. Thus, OP compounds can be used as a tool to probe the physiological functions of APH. Here, we report the results of a computational study of molecular dynamics simulations of APH bound to the OP compounds and an exploration of the chlorpyrifos escape pathway using steered molecular dynamics (SMD) simulations. In addition, we apply SMD simulations to identify potential escape routes of chlorpyrifos from hydrolase hydrophobic cavities in the APH-inhibitor complex. Two previously proposed APH pathways were reliably identified by CAVER 3.0, with the estimated relative importance of P1 > P2 for its size. We identify the major pathway, P2, using SMD simulations, and Arg526, Glu88, Gly86, and Asn65 are identified as important residues for the ligand leaving via P2. These results may help in the design of APH-targeting drugs with improved efficacy, as well as in understanding APH selectivity of the inhibitor binding in the prolyl oligopeptidase family.

Insight into the interactive residues between two domains of human somatic Angiotensin-converting enzyme and Angiotensin II by MM-PBSA calculation and steered molecular dynamics simulation

Angiotensin-converting enzyme (ACE), a membrane-bound zinc metallopeptidase, catalyzes the formation of Angiotensin-II (AngII) and the deactivation of bradykinin in the renin-angiotensin-aldosterone and kallikrein-kinin systems. As a hydrolysis product of ACE, AngII is regarded as an inhibitor and displays stronger competitive inhibition in the C-domain than the N-domain of ACE. However, the AngII binding differences between the two domains and the mechanisms behind AngII dissociation from the C-domain are rarely explored. In this work, molecular docking, Molecular Mechanics/Poisson-Boltzmann Surface Area calculation, and steered molecular dynamics (SMD) are applied to explore the structures and interactions in the binding or unbinding of AngII with the two domains of human somatic ACE. Calculated free energy values suggest that the C-domain-AngII complex is more stable than the N-domain-AngII complex, consistent with available experimental data. SMD simulation results imply that electrostatic interaction is dominant in the dissociation of AngII from the C-domain. Moreover, Gln106, Asp121, Glu123, and Tyr213 may be the key residues in the unbinding pathway of AngII. The simulation results in our work provide insights into the interactions between the two domains of ACE and its natural peptide inhibitor AngII at a molecular level. Moreover, the results provide theoretical clues for the design of new inhibitors.

A QM/MM study of the catalytic mechanism of nicotinamidase

Nicotinamidase (Pnc1) is a member of Zn-dependent amidohydrolases that hydrolyzes nicotinamide (NAM) to nicotinic acid (NA), which is a key step in the salvage pathway of NAD(+) biosynthesis. In this paper, the catalytic mechanism of Pnc1 has been investigated by using a combined quantum-mechanical/molecular-mechanical (QM/MM) approach based on the recently obtained crystal structure of Pnc1. The reaction pathway, the detail of each elementary step, the energetics of the whole catalytic cycle, and the roles of key residues and Zn-binding site are illuminated. Our calculation results indicate that the catalytic water molecule comes from the bulk solvent, which is then deprotonated by residue D8. D8 functions as a proton transfer station between C167 and NAM, while the activated C167 serves as the nucleophile. The residue K122 only plays a role in stabilizing intermediates and transition states. The oxyanion hole formed by the amide backbone nitrogen atoms of A163 and C167 has the function to stabilize the hydroxyl anion of nicotinamide. The Zn-binding site rather than a single Zn(2+) ion acts as a Lewis acid to influence the reaction. Two elementary steps, the activation of C167 in the deamination process and the decomposition of catalytic water in the hydrolysis process, correspond to the large energy barriers of 25.7 and 28.1 kcal mol(-1), respectively, meaning that both of them contribute a lot to the overall reaction barrier. Our results may provide useful information for the design of novel and efficient Pnc1 inhibitors and related biocatalytic applications.

A steered molecular dynamics mediated hit discovery for histone deacetylases

The inhibitors of class I histone deacetylases (HDACIs) have gained significant interest in cancer therapeutics. Virtual high throughput screening (vHTS) is one of the popular approaches used in the identification of novel scaffolds of HDACIs. However, an accurate description of ligand-protein flexibilities in the vHTS remains challenging. In this work, we implement an integrated approach, which combines the vHTS with the 'state-of-the-art' steered molecular dynamics (SMD). This approach serves as an efficient tool to identify potential hits and characterize their binding potencies against the class I HDACs in a flexible solvent environment. A hybrid pharmacophore-based and structure-based vHTS method identifies the hits with more favourable physico-chemical features against the class I HDACs. Our pharmacophore-based screening enhanced the quality of the vHTS outcomes. Further, the molecular interactions between the hits and the HDACs are investigated using the SMD-driven force profiles, which in turn resulted in filtering the hits with higher binding potencies against the HDACs. Our results, therefore, reveal that vHTS and SMD can be a complementary and effective analytical tool for accelerating the hit identification phase in structure-based drug design.

The evolutionary portrait of metazoan NAD salvage

Nicotinamide Adenine Dinucleotide (NAD) levels are essential for cellular homeostasis and survival. Main sources of intracellular NAD are the salvage pathways from nicotinamide, where Nicotinamide phosphoribosyltransferases (NAMPTs) and Nicotinamidases (PNCs) have a key role. NAMPTs and PNCs are important in aging, infection and disease conditions such as diabetes and cancer. These enzymes have been considered redundant since either one or the other exists in each individual genome. The co-occurrence of NAMPT and PNC was only recently detected in invertebrates though no structural or functional characterization exists for them. Here, using expression and evolutionary analysis combined with homology modeling and protein-ligand docking, we show that both genes are expressed simultaneously in key species of major invertebrate branches and emphasize sequence and structural conservation patterns in metazoan NAMPT and PNC homologues. The results anticipate that NAMPTs and PNCs are simultaneously active, raising the possibility that NAD salvage pathways are not redundant as both are maintained to fulfill the requirement for NAD production in some species.

Biochemical and mutational analysis of a novel nicotinamidase from Oceanobacillus iheyensis HTE831

Nicotinamidases catalyze the hydrolysis of nicotinamide to nicotinic acid and ammonia, an important reaction in the NAD(+) salvage pathway. This paper reports a new nicotinamidase from the deep-sea extremely halotolerant and alkaliphilic Oceanobacillus iheyensis HTE831 (OiNIC). The enzyme was active towards nicotinamide and several analogues, including the prodrug pyrazinamide. The enzyme was more nicotinamidase (kcat/Km  = 43.5 mM(-1)s(-1)) than pyrazinamidase (kcat/Km  = 3.2 mM(-1)s(-1)). Mutational analysis was carried out on seven critical amino acids, confirming for the first time the importance of Cys133 and Phe68 residues for increasing pyrazinamidase activity 2.9- and 2.5-fold, respectively. In addition, the change in the fourth residue involved in the ion metal binding (Glu65) was detrimental to pyrazinamidase activity, decreasing it 6-fold. This residue was also involved in a new distinct structural motif DAHXXXDXXHPE described in this paper for Firmicutes nicotinamidases. Phylogenetic analysis revealed that OiNIC is the first nicotinamidase described for the order Bacillales.

Drug resistance mechanism of PncA in Mycobacterium tuberculosis

Tuberculosis continues to be a global health threat. Pyrazinamide (PZA) is an important first-line drug in multidrug-resistant tuberculosis treatment. The emergence of strains resistant to PZA represents an important public health problem, as both first- and second-line treatment regimens include PZA. It becomes toxic to Mycobacterium tuberculosis when converted to pyrazinoic acid by the bacterial pyrazinamidase (PncA) enzyme. Resistance to PZA is caused mainly by the loss of enzyme activity by mutation, the mechanism of resistance is not completely understood. In our studies, we analysed three mutations (D8G, S104R and C138Y) of PncA which are involved in resistance towards PZA. Binding pocket analysis solvent accessibility analysis, molecular docking and interaction analysis were performed to understand the interaction behaviour of mutant enzymes with PZA. Molecular dynamics simulations were conducted to understand the three-dimensional (3D) conformational behaviour of native and mutants PncA. Our analysis clearly indicates that the mutation (D8G, S104R and C138Y) in PncA is responsible for rigid binding cavity which in turn abolishes conversion of PZA to its active form and is the sole reason for PZA resistance. Excessive hydrogen bonding between PZA binding cavity residues and their neighbouring residues are the reason of rigid binding cavity during simulation. We present an exhaustive analysis of the binding site flexibility and its 3D conformations that may serve as new starting points for structure-based drug design and helps the researchers to design new inhibitors with consideration of rigid criterion of binding residues due to mutation of this essential target. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:11.

How does (E)-2-(acetamidomethylene)succinate bind to its hydrolase? From the binding process to the final result

The binding of (E)-2-(acetamidomethylene)succinate (E-2AMS) to E-2AMS hydrolase is crucial for biological function of the enzyme and the last step reaction of vitamin B(6) biological degradation. In the present study, several molecular simulation methods, including molecular docking, conventional molecular dynamics (MD), steered MD (SMD), and free energy calculation methods, were properly integrated to investigate the detailed binding process of E-2AMS to its hydrolase and to assign the optimal enzyme-substrate complex conformation. It was demonstrated that the substrate binding conformation with trans-form amide bond is energetically preferred conformation, in which E-2AMS's pose not only ensures hydrogen bond formation of its amide oxygen atom with the vicinal oxyanion hole but also provides probability of the hydrophobic interaction between its methyl moiety and the related enzyme's hydrophobic cavity. Several key residues, Arg146, Arg167, Tyr168, Arg179, and Tyr259, orientate the E-2AMS's pose and stabilize its conformation in the active site via the hydrogen bond interaction with E-2AMS. Sequentially, the binding process of E-2AMS to E-2AMS hydrolase was studied by SMD simulation, which shows the surprising conformational reversal of E-2AMS. Several important intermediate structures and some significant residues were identified in the simulation. It is stressed that Arg146 and Arg167 are two pivotal residues responsible for the conformational reversal of E-2AMS in the binding or unbinding. Our research has shed light onto the full binding process of the substrate to E-2AMS hydrolase, which could provide more penetrating insight into the interaction of E-2AMS with the enzyme and would help in the further exploration on the catalysis mechanism.