Secondary structural changes of metmyoglobin and apomyoglobin in anionic and cationic surfactant solutions: Effect of the hydrophobic chain length of the surfactants on the structural changes

Secondary structural changes of metmyoglobin and apomyoglobin were examined in solutions of sodium alkylsulfates with hydrocarbon numbers of 8 and 12, and alkyltrimethylammonium bromides with hydrocarbon numbers of 10, 12, 14, and 16. The relative proportion ofa-helical structure was estimated by the curve-fitting method of circular dichroic spectrum. The helical proportions of metmyoglobin and apomyoglobin were 82 and 63%, respectively. The shorter the hydrocarbon chain the surfactant had, the higher the concentration necessary to disrupt the secondary structures of these proteins. However, the helical proportion had a tendency to decrease down to lower values in solutions of the cationic surfactants with short hydrophobic groups. On the other hand, thea-helical structure of apomyoglobin was disrupted in lower concentrations of each cationic surfactant than that of metmyoglobin, although the disruptions of the same structures in both the proteins occurred in the same concentration range of each anionic surfactant. It appeared likely that the removal of the heme group unstabilized the myoglobin conformation only in the cationic surfactant solutions.     

Spectroscopy and molecular docking study on the interaction of daidzein and genistein with pepsin

The interaction of pepsin with daidzein (Dai) or genistein (Gen) was investigated using spectroscopic techniques under simulated physiological conditions. Dai and Gen can quench the fluorescence of pepsin and the quenching mechanism was a static process. The binding site number n and apparent binding constant K were measured at different temperatures. The thermodynamic parameters ΔΗ, ΔG and ΔS were calculated. The results indicated that van der Waals forces and hydrogen bond formation played major roles in the interaction of Dai or Gen with pepsin. The binding distance between pepsin and Dai or Gen was calculated according to energy transfer theory. The results of synchronous fluorescence spectra showed that the microenvironment and conformation of pepsin were changed. UV absorption and 3D fluorescence spectra showed that the binding interaction disturbed the microenvironment of amino acid residues and induced conformational changes in pepsin. Molecular docking results showed that Dai and Gen entered into the hydrophobic cavity of pepsin and two hydrogen bonds formed between Dai or Gen and pepsin. The results demonstrated that the interaction behavior between Dai and Gen with pepsin was slightly different, which denoted that the 5‐hydroxyl group of Gen, to a certain extent, had an effect on ligand binding to proteins. Copyright © 2016 John Wiley & Sons, Ltd.     

Computational and experimental validation of morin as adenosine deaminase inhibitor

Adenosine deaminase (ADA) is one of the major enzymes involved in purin metabolism, it has a significant role in cell growth and differentiation. Over-activity of ADA has been noticed in some pathology, like malignancy and inflammation and makes it an attractive target for the development of drugs for such diseases. In the present study, ADA inhibitory activity of morin, a bioactive flavonoid, was assessed through computational and biophysical methods. The enzyme kinetics data showed that morin is a competitive inhibitor of ADA. Binding energy calculated from ITC analysis was ?7.11?kcal/mol. Interaction of morin with ADA was also studied using fluorescence quenching method. Molecular docking studies revealed the structural details of the interaction. Molecular dynamics study in explicit solvent was also conducted to assess the structural stability of protein ligand complex.     

Stereoselectivity of phosphotriesterase with paraoxon derivatives: a computational study

The bacterial enzyme phosphotriesterase (PTE) exhibits stereoselectivity toward hydrolysis of chiral substrates with a preference for the S_p enantiomer. In this work, docking analysis and two explicit-solvent molecular dynamics (MD) simulations were performed to characterize and differentiate the structural dynamics of PTE bound to the S_p and R_p paraoxon derivative enantiomers (R_p-1 and S_p-1) hydrolyzed with distinct catalytic efficiencies. Comparative analysis of the molecular trajectories for PTE bound to R_p-1 and S_p-1 suggested that substrate binding induced conformational changes in the loops near the active site. After 100 ns of MD simulation, the Zn β²⁺ metal ion formed hexacoordinated- and tetracoordinated geometries in the S_p-1-PTE and R_p-1-PTE ensembles, respectively. Simulation results further showed that the hydrogen bond between Asp301 and His254 occurred with a higher probability after S_p-1 binding to PTE (47.5%) than that after R_p-1 binding (22.2%). These results provide a qualitative and molecular-level explanation for the 10 orders of magnitude increase in the catalytic efficiency of PTE toward the S_p enantiomer of paraoxon.     

Homology modeling of human serum paraoxonase1 and its molecular interaction studies with aspirin and cefazolin

Human serum paraoxonase1 (HuPON1) belongs to the family of A-esterases (EC.3.1.8.1). It is associated with HDL particle and prevents atherosclerosis by cleaving lipid hydroperoxides and other proatherogenic molecules of oxidized low density lipoproteins (LDL). Since the precise structure of HuPON1 is not yet available, the structure-function relationship between HuPON1 and activators/inhibitors is still unknown. Therefore, a theoretical model of HuPON1 was generated using homology modelling and precise molecular interactions of an activator aspirin and an inhibitor cefazolin with PON1 were studied using Autodock software. The ligand binding residues were found to be similar to the predicted active site residues. Both cefazolin and aspirin were found to dock in the vicinity of the predicted active sites of PON1; cefazolin bound at residues N166, S193 and Y71, while aspirin at residues N309, I310 and L311. Binding region in the PON1 by prediction (3D2GO server) and docking studies provide useful insight into mechanism of substrate and inhibitor binding to the enzyme active site.     

Theoretical studies on the interaction of guanine riboswitch with guanine and its closest analogues

Experimental studies (M. Mandal, B. Boese, J.E. Barrick, W.C. Winkler and R.R. Breaker, Riboswitches control fundamental biochemical pathways in bacillus subtilis and other bacteria, Cell 113 (2003), pp. 577–586) demonstrated that, besides recognising guanine with high specificity, guanine riboswitch could also bind guanine analogues, but the alteration of every functionalised position on the guanine heterocycle could cause a substantial loss of binding affinity. To investigate the nature of guanine riboswitch recognising metabolites, molecular docking and molecular dynamics simulation were carried out on diverse guanine analogues. The calculation results reveal that (1) most guanine analogues could bind to guanine riboswitch at the same binding pocket, with identical orientations and dissimilar binding energies, which is related to the positions of the functional groups; (2) the two tautomers of xanthine adopt different binding modes, and the enol-tautomer shows similar binding mode and affinity of hypoxanthine, which agrees well with the experimental results and (3) the riboswitch could form stable complexes with guanine analogues by hydrogen bonding contacts with U51 and C74. Particularly, U51 plays an important role in stabilising the complexes.     

Multi-spectroscopic and molecular docking studies on the interaction of neotame with calf thymus DNA

Abstract

In this paper, we have studied the in vitro binding of neotame (NTM), an artificial sweetener, with native calf thymus DNA using different methods including spectrophotometric, spectrofluorometric, competition experiment, circular dichroism (CD), and viscosimetric techniques. From the spectrophotometric studies, the binding constant (K_b) of NTM-DNA was calculated to be 2?×?10³ M^?1. The quenching of the intrinsic fluorescence of NTM in the presence of DNA at different temperatures was also used to calculate binding constants (K_b) as well as corresponding number of binding sites (n). Moreover, the obtained results indicated that the quenching mechanism involves static quenching. By comparing the competitive fluorimetric studies with Hoechst 33258, as a known groove probe, and methylene blue, as a known intercalation probe, and iodide quenching experiments it was revealed that NTM strongly binds in the grooves of the DNA helix, which was further confirmed by CD and viscosimetric studies. In addition, a molecular docking method was employed to further investigate the binding interactions between NTM and DNA, and confirm the obtained results.     

Exploring the structural determinants of novel xanthine derivatives as A2B adenosine receptor antagonists: a computational study

Adenosine is a ubiquitous endogenous nucleoside that controls numerous physiological functions via interacting with its specific G-coupled receptors. Activation of adenosine receptors (AdoRs), particularly A_2B AdoRs promotes the release of inflammatory cytokines; reduces vascular permeabilization and induces angiogenesis, thereby making A_2B AdoR becomes a potentially pharmacological target for drug development. Presently, for investigating the structural determinants of 164 xanthine derivatives as A_2B AdoR antagonists, we performed an in silico study integrating with 3D-QSAR, docking and molecular dynamics (MD) simulation. The obtained optimal model shows strong predictability (Q²?=?0.647, R²_ncv?=?0.955, and R²_pred?=?0.848). Additionally, to explore the binding mode of the ligand with A_2B AdoR and to understand their binding mechanism, docking analysis, MD simulations (20?ns), and the calculation of binding free energy were also carried out. Finally, the structural determinants of these xanthine derivatives were identified and a total of 20 novel A_2B AdoR antagonists with improved potency were computationally designed, and their synthetic feasibility and selectivity were also evaluated. The information derived from the present study offers a better appreciation for exploring the interaction mechanism of the ligand with A_2B AdoR, which could be helpful for designing novel potent A_2B AdoR antagonists.

Communicated by Ramaswamy H. Sarma     

Effects of surfactant, salt and solvent on the structure and activity of adenosine deaminase: molecular dynamic and spectrophotometric studies

Effects of sodium dodecyl sulfate, dodecyltrimethylammonium bromide, sodium chloride, sodium sulfate, methanol and ethanol, on the structure and activity of adenosine deaminase (ADA) were investigated by UV-Vis, circular dichroism spectrophotometry and molecular dynamics (MDs) studies. Relative activity, experimental and computational helix content, total accessible surface area (ASA) and exposed charged surface area (ECSA) were obtained. The relative activity of ADA in the absence and the presence of denaturants were compared with structural results. It was shown that an increase in the surface area and a decrease in the amount of helicity are associated with a decrease in the activity of ADA.     

Experimental and computational characterization on the binding of two fluoroquinolones to bovine hemoglobin

Ciprofloxacin (CPFX) and enrofloxacin (ENFX) are 2 representatives of widely used fluoroquinolones (FQs) with many human and veterinary applications. The residues of FQs in the environment are potentially harmful. Recently, great concern has been paid to their persistence and fate in the environment because of the potential adverse effects on humans and ecosystem functions. In the present study, we examined the interactions of bovine hemoglobin (BHb) with these 2 FQs by means of multiple spectroscopic and molecular docking methods under physiological conditions. The experimental results revealed that both FQs could bind with BHb to form complexes mainly through electrostatic interactions. And CPFX posed more of an affinity threat to BHb than ENFX. On the basis of molecular docking, both FQs could bind into the central cavity of BHb and interact with the residue Trp 37, resulting in the remarkable fluorescence quenching of protein. Additionally, as shown by the synchronous fluorescence, UV‐visible absorption and circular dichroism data, both CPFX and ENFX could lead to the conformational and microenvironmental changes of BHb, which may affect its physiological functions. The work is beneficial for understanding the biological toxicity of FQs in vivo.     

Investigation of binding properties of two ethidium derivatives with serum albumins: spectral and computational approach

The interaction mechanisms of two ethidium derivatives, 3,8-dibenzoylamino-5-ethyl-6-phenylphenantridinium chloride (E2) and 3,8-diphenylacetylamino-5-ethyl-6-phenylphenantridinium chloride (E3) with serum albumins (BSA and HSA) have been investigated by a combined experimental and computational approach. Fluorescence quenching and UV–vis results revealed that the interaction of derivatives with albumins resulted in formation of ground-state complexes and the obtained Stern–Volmer quenching constants designate the presence of a static component in the quenching mechanisms. Thermodynamic parameters (ΔH and ΔS values) point out the ionic interactions play the major role in E2-BSA, E2-HSA and E3-HSA complexes. The van der Waals interactions are dominant forces in E3-BSA complex. Moreover, the obtained results in this study were supported with computational analyzes which have same tendency.     

Novel urushiol derivatives as HDAC8 inhibitors: rational design,virtual screening,molecular docking and molecular dynamics studies

Three series of novel urushiol derivatives were designed by introducing a hydroxamic acid moiety into the tail of an alkyl side chain and substituents with differing electronic properties or steric bulk onto the benzene ring and alkyl side chain. The compounds’ binding affinity toward HDAC8 was screened by Glide docking. The highest-scoring compounds were processed further with molecular docking, MD simulations, and binding free energy studies to analyze the binding modes and mechanisms. Ten compounds had Glide scores of ?8.2 to ?10.2, which revealed that introducing hydroxy, carbonyl, amino, or methyl ether groups into the alkyl side chain or addition of –F, –Cl, sulfonamide, benzamido, amino, or hydroxy substituents on the benzene ring could significantly increase binding affinity. Molecular docking studies revealed that zinc ion coordination, hydrogen bonding, and hydrophobic interactions contributed to the high calculated binding affinities of these compounds toward HDAC8. MD simulations and binding free energy studies showed that all complexes possessed good stability, as characterized by low RMSDs, low RMSFs of residues, moderate hydrogen bonding and zinc ion coordination and low values of binding free energies. Hie147, Tyr121, Phe175, Hip110, Phe119, Tyr273, Lys21, Gly118, Gln230, Leu122, Gly269, and Gly107 contributed favorably to the binding; and Van der Waals and electrostatic interactions provided major contributions to the stability of these complexes. These results show the potential of urushiol derivatives as HDAC8 binding lead compounds, which have great therapeutic potential in the treatment of various malignancies, neurological disorders, and human parasitic diseases.     

NMR and molecular dynamics studies of the interaction of melatonin with calmodulin

Pineal hormone melatonin (N-acetyl-5-methoxytryptamine) is thought to modulate the calcium/calmodulin signaling pathway either by changing intracellular Ca(2+) concentration via activation of its G-protein-coupled membrane receptors, or through a direct interaction with calmodulin (CaM). The present work studies the direct interaction of melatonin with intact calcium-saturated CaM both experimentally, by fluorescence and nuclear magnetic resonance spectroscopies, and theoretically, by molecular dynamics simulations. The analysis of the experimental data shows that the interaction is calcium-dependent. The affinity, as obtained from monitoring (15)N and (1)H chemical shift changes for a melatonin titration, is weak (in the millimolar range) and comparable for the N- and C-terminal domains. Partial replacement of diamagnetic Ca(2+) by paramagnetic Tb(3+) allowed the measurement of interdomain NMR pseudocontact shifts and residual dipolar couplings, indicating that each domain movement in the complex is not correlated with the other one. Molecular dynamics simulations allow us to follow the dynamics of melatonin in the binding pocket of CaM. Overall, this study provides an example of how a combination of experimental and theoretical approaches can shed light on a weakly interacting system of biological and pharmacological significance.     

Intermolecular interaction of fosinopril with bovine serum albumin (BSA): The multi‐spectroscopic and computational investigation

The intermolecular interaction of fosinopril, an angiotensin converting enzyme inhibitor with bovine serum albumin (BSA), has been investigated in physiological buffer (pH 7.4) by multi‐spectroscopic methods and molecular docking technique. The results obtained from fluorescence and UV absorption spectroscopy revealed that the fluorescence quenching mechanism of BSA induced by fosinopril was mediated by the combined dynamic and static quenching, and the static quenching was dominant in this system. The binding constant, K_b, value was found to lie between 2.69 × 10³ and 9.55 × 10³ M^?1 at experimental temperatures (293, 298, 303, and 308 K), implying the low or intermediate binding affinity between fosinopril and BSA. Competitive binding experiments with site markers (phenylbutazone and diazepam) suggested that fosinopril preferentially bound to the site I in sub‐domain IIA on BSA, as evidenced by molecular docking analysis. The negative sign for enthalpy change (ΔH⁰) and entropy change (ΔS⁰) indicated that van der Waals force and hydrogen bonds played important roles in the fosinopril‐BSA interaction, and 8‐anilino‐1‐naphthalenesulfonate binding assay experiments offered evidence of the involvements of hydrophobic interactions. Moreover, spectroscopic results (synchronous fluorescence, 3‐dimensional fluorescence, and Fourier transform infrared spectroscopy) indicated a slight conformational change in BSA upon fosinopril interaction.     

Insights into the interaction of ulipristal acetate and human serum albumin using multi-spectroscopic methods,molecular docking,and dynamic simulation

Interaction between ulipristal acetate (UPA) and human serum albumin (HSA) was investigated in simulated physiological environment using multi-spectroscopic and computational methods. Fluorescence experiments showed that the quenching mechanism was static quenching, which was confirmed by the time-resolved fluorescence. Binding constants (K_a) were found to be 1?×?10⁵ L mol^?1, and fluorescence data showed one binding site. Thermodynamic constants suggested the binding process was mainly controlled by electrostatic interactions. Results from the competition experiments indicated that UPA bound to site I of HSA. Fourier transform infrared spectra, circular dichroism spectra, synchronous fluorescence spectra, and 3D fluorescence indicated that UPA can induce conformation change in the HSA. The content of α-helix and β-sheet increased, while β-turn decreased. Hydrophobicity around the tryptophan residues declined, whereas its polarity increased. Molecular docking results were consistent with the experimental results. Results suggested that UPA located at the hydrophobic cavity site I of HSA, and hydrophobic force played the key role in the binding process. Moreover, molecular dynamics simulation was performed to determine the stability of free HSA and HSA-UPA system. Results indicated that UPA can stabilize HSA to a certain degree and enhance the flexibility of residues around site I.

Communicated by Ramaswamy H. Sarma     

Exploring the interaction between tyrphostin 9 and human serum albumin using biophysical and computational methods

Abstract

Tyrphostin 9 (Tyr 9) is a potent platelet-derived growth factor receptor (PDGFR) inhibitor, which induces apoptosis in various cancer cell types. The binding of Tyr 9 to the major transport protein, human serum albumin (HSA) was investigated using several spectroscopic techniques and molecular docking method. Fluorescence quenching titration results showed progressive decrease in the protein fluorescence with increasing drug concentrations. A decreasing trend of the Stern-Volmer constant, K _sv with increasing temperature characterized the drug-induced quenching as static quenching, thus pointed towards the formation of Tyr 9–HSA complex. The binding constant of Tyr 9–HSA interaction was found to lie within the range 3.48–1.69?×?10⁵ M^?1 at three different temperatures, i.e. 15 °C, 25 °C and 35?°C, respectively and suggested intermediate binding affinity between Tyr 9 and HSA. The drug–HSA complex seems to be stabilized by hydrophobic forces, van der Waals forces and hydrogen bonds, as suggested from the thermodynamic data as well as molecular docking results. The far-UV and the near-UV CD spectral results showed slight alteration in the secondary and tertiary structures, respectively, of the protein upon Tyr 9 binding. Interaction of Tyr 9 with HSA also produced microenvironmental perturbations around protein fluorophores, as evident from the three-dimensional fluorescence spectral results but increased protein’s thermal stability. Both competitive drug binding results and molecular docking analysis suggested Sudlow’s Site I of HSA as the preferred Tyr 9 binding site.

Communicated by Ramaswamy H. Sarma     

Characterizing the noncovalent binding behavior of tartrazine to lysozyme: A combined spectroscopic and computational analysis

Tartrazine is a stable water‐soluble azo dye widely used as a food additive, which could pose potential threats to humans and the environment. In this paper, we evaluated the response mechanism between tartrazine and lysozyme under simulated conditions by means of biophysical methods, including multiple spectroscopic techniques, isothermal titration calorimetry (ITC), and molecular docking studies. From the multispectroscopic analysis, we found that tartrazine could effectively quench the intrinsic fluorescence of lysozyme to form a complex and lead to the conformational and microenvironmental changes of the enzyme. The ITC measurements suggested that the electrostatic forces played a major role in the binding of tartrazine to lysozyme with two binding sites. Finally, the molecular docking indicated that tartrazine had specific interactions with the residues of Trp108. The study provides an important insight within the binding mechanism of tartrazine to lysozyme in vitro.     

A newly synthesized 6-methyl-7H,8H,9H-[1,2,4]triazolo[4,3-b][1,2,4]triazepin-8-one for potential inhibitor of adenosine A1 receptor: a combined experimental and computational studies

Abstract

6-Methyl-7H,8H,9H-[1,2,4]triazolo[4,3-b][1,2,4]triazepin-8-onehas been synthesized, characterized by spectroscopic techniques (FT-IR, ¹H and ¹³C NMR) and finally the structure was confirmed by single crystal X-ray diffraction studies. In the title molecule, C₆H₇N₅O, the 7-membered ring adopts a bowl-like conformation. In the crystal, the molecules form stacks along the c-axis direction through offset π-stacking interactions between the 5-membered rings and C–H···N hydrogen bonds. The stacks are associated via a combination of N–H···N, C–H···O and C–H···N hydrogen bonds. Further, the Hirshfeld surface analysis reveals the nature of molecular interactions and the fingerprint plot provides information about the percentage contribution from each individual molecular contact to the surface. In addition, due to its biological interest the target molecule adenosine A1 receptor was found based on Structural Activity Relationship (SAR) analysis and, further, subjected into molecular docking and molecular dynamics analysis to understand the binding interaction and stability of the molecule in adenosine A1 receptor system. Furthermore, the Density Functional Theory (DFT) calculations were carried for free compound and the compound in active site (single point DFT), to know the internal stability.

Communicated by Ramaswamy H. Sarma     

Comparative study of the interaction mechanism of astilbin,isoastilbin, and neoastilbin with CYP3A4

The interactions of human CYP3A4 with three selected isomer flavonoids, such as astilbin, isoastilbin and neoastilbin, were clarified using spectral analysis, molecular docking, and molecular dynamics simulation. During binding with the three flavonoids, the intrinsic fluorescence of CYP3A4 was statically quenched in static mode with nonradiative energy conversion. The fluorescence and ultraviolet/visible (UV/vis) data revealed that the three flavonoids had a moderate and stronger binding affinity with CYP3A4 due to the order of the K_a1 and K_a2 values ranging from 10⁴ to 10⁵ L·mol⁻¹. In addition, astilbin had the highest affinity with CYP3A4, then isoastilbin and neoastilbin, at the three experimental temperatures. Multispectral analysis confirmed that binding of the three flavonoids resulted in clear changes in the secondary structure of CYP3A4. It was found from fluorescence, UV/vis and molecular docking analyses that these three flavonoids strongly bound to CYP3A4 by means of hydrogen bonds and van der Waals forces. The key amino acids around the binding site were also elucidated. Furthermore, the stabilities of the three CYP3A4 complexes were evaluated using molecular dynamics simulation.