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     

Molecular insights into Aβ42 protofibril destabilization with a fluorinated compound D744: A molecular dynamics simulation study

The aggregation of amyloid β‐peptide (Aβ₄₂) into toxic oligomers, fibrils, has been identified as a key process in Alzheimer's disease (AD) progression. The role of halogen‐substituted compounds have been highlighted in the disassembly of Aβ protofibril. However, the underlying inhibitory mechanism of Aβ₄₂ protofibril destabilization remains elusive. In this regard, a combined molecular docking and molecular dynamics (MD) simulations were performed to elucidate the inhibitory mechanism of a fluorinated compound, D744 , which has been reported previously for potential in vitro and in vivo inhibitory activity against Aβ₄₂ aggregation and reduction in the Aβ‐induced cytotoxicity. The molecular docking analysis highlights that D744 binds and interacts with chain A of the protofibril structure with hydrophobic contacts and orthogonal multipolar interaction. MD simulations reveal destabilization of the protofibril structure in the presence of D744 due to the decrease in β‐sheet content and a concomitant increase of coil and bend structures, increase in the interchain D23‐K28 salt bridge distance, decrease in the number of backbone hydrogen bonds, increase in the average distance between Cα atoms, and decrease in the binding affinity between chains A and B of the protofibril structure. The binding free‐energy analysis between D744 and the protofibril structure with Molecular Mechanics Poisson‐Boltzmann Surface Area (MM‐PBSA) reveal that residues Leu17, Val18, Phe19, Phe20, Ala21, Glu22, Asp23, Leu34, Val36, Gly37, and Gly38 of chain A of the protofibril structure contribute maximum towards binding free energy (ΔG _binding  = −44.87 kcal/mol). The insights into the underlying inhibitory mechanism of small molecules that show potential in vitro anti‐aggregation activity against Aβ₄₂ will be beneficial for the current and future AD therapeutic studies.     

Comparative MD analysis of the stability of transthyretin providing insight into the fibrillation mechanism

Proteins can misfold and aggregate, which is believed to be the cause of a variety of diseases, affecting very diverse organs in the body. Many questions about the nature of aggregation and the proteins that are involved in these events are still left unanswered. One of the proteins that is known to form amyloids is transthyretin (TTR), the secondary transporter of thyroxine, and transporter of retinol-binding protein. Several experimental results have helped to explain this aberrant behavior of TTR; however, structural insights of the amyloidgenic process are still lacking. Therefore, we have used all-atom MD simulation and free energy calculations to study the initial phase of this process. We have calculated the free energy changes of the initial tetramer dissociation under different conditions and in the presence of thyroxine. We show that tetramer formation is indeed only thermodynamically favorable in neutral pH conditions. We find that binding of two thyroxine molecules stabilizes the complex, and that this occurs with negative cooperativity. In addition to the energetic calculations, we have also investigated the dominant motions of the TTR and found that only the dimeric form of the protein could undergo the initial fibril formation.     

Crystal structures of the luciferase and green fluorescent protein from Renilla reniformis

Due to its ability to emit light, the luciferase from Renilla reniformis (RLuc) is widely employed in molecular biology as a reporter gene in cell culture experiments and small animal imaging. To accomplish this bioluminescence, the 37-kDa enzyme catalyzes the degradation of its substrate coelenterazine in the presence of molecular oxygen, resulting in the product coelenteramide, carbon dioxide, and the desired photon of light. We successfully crystallized a stabilized variant of this important protein (RLuc8) and herein present the first structures for any coelenterazine-using luciferase. These structures are based on high-resolution data measured to 1.4 Å and demonstrate a classic α/β-hydrolase fold. We also present data of a coelenteramide-bound luciferase and reason that this structure represents a secondary conformational form following shift of the product out of the primary active site. During the course of this work, the structure of the luciferase's accessory green fluorescent protein (RrGFP) was also determined and shown to be highly similar to that of Aequorea victoria GFP.     

Study on the inhibitory mechanism and binding mode of the hydroxycoumarin compound NSC158393 to HIV‐1 integrase by molecular modeling

Human immunodeficiency virus type 1 integrase (IN) is an essential enzyme in the life cycle of this virus and also an important target for the study of anti‐HIV drugs. In this work, the binding modes of the wild type IN core domain and the two mutants, that is, W132G and C130S, with the 4‐hydroxycoumarin compound NSC158393 were evaluated by using the “relaxed complex” molecular docking approach combined with molecular dynamics (MD) simulations. Based on the monomer MD simulations, both of the two substitutions affect not only the stability of the 128–136 peptides, but also the flexibility of the functional 140s loop. In principle, NSC158393 binds the 128–136 peptides of IN; however, the specific binding modes for the three systems are various. According to the binding mode of NSC158393 with WT, NSC158393 can effectively interfere with the stability of the IN dimer by causing a steric hindrance around the monomer interface. Additionally, through the comparative analysis of the MD trajectories of the wild type IN and the IN‐NSC158393 complex, we found that NSC15893 may also exert its inhibitory function by diminishing the mobility of the function loop of IN. Three key binding residues, that is, W131, K136, and G134, were discovered by energy decomposition calculated with the Molecular Mechanics Generalized Born Surface Area method. Characterized by the largest binding affinity, W131 is likely to be indispensable for the ligand binding. All the above results are consistent with experiment data, providing us some helpful information for understanding the mechanism of the coumarin‐based inhibitors. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 700–709, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Binding mode of chitin and TLR2 via molecular docking and dynamics simulation

Innate immunity is an important part of immune system, providing immediate defence for the host against various infections through phagocytes. Toll-like receptors (TLRs) are major proteins expressed on the cell membrane known as pattern recognition receptors (PRR) that recognise non-self molecules (pathogen-associated molecular patterns (PAMPs)). Because TLRs have been implicated in many inflammatory diseases and cancer, TLRs targeted therapeutics have drawn great attention in clinical application in wide range of conditions. Many of them are undergoing evaluation in clinical trials. Chitin is the second most abundant polysaccharide detected in many insects and fungi. Studies have shown that chitin, as major PAMPs in host-infection, can activate TLR2-dependent innate immunity pathway. Therefore, chitin has potential use as an important agonist or antagonist to control key processes in innate immunity. However, no direct evidence has shown that chitin is the direct target of TLR2. This study first demonstrates a binding model of chitin and TLR2 and then confirmed its stability by molecular dynamic simulation and MM/PBSA (molecular mechanics/Poisson?Boltzmann surface area) calculations. The binding between chitin and TLR2 was taken place inside the binding pocket. Two hydrogen bonds were formed between chitin and TLR2, including Ser320 and Lys321. The van der Waals interaction has the major contribution in stabilising the binding of the chitin molecule with the protein. This study also suggests six hot-spots for specific binding of chitin in the binding site of TLR2, namely, Phe296, Phe299, Leu302, Thr309, Ser320 and Val322. Molecular dynamics simulation demonstrates that the complex of chitin and TLR2 is very stable with a total binding affinity of ?27.2 kcal/mol from MM/PBSA calculation.     

Structural insights into ligand-binding pocket formation in Nurr1 by molecular dynamics simulations

Abstract

The nuclear receptor Nurr1 (NR4A2) has been identified as a potential target for the treatment of Parkinson’s disease. In contrast to most other nuclear receptors, the X-ray crystal structure of the Nurr1 ligand-binding domain (LBD) lacks any ligand-binding pocket (LBP). However, NMR spectroscopy measurements have revealed that the known Nurr1 agonist docosahexaenoic acid (DHA) binds to a region within the LBD that corresponds to the classical NR ligand-binding pocket (LBP). In order to investigate the structural dynamics of the Nurr1 LBD and to study potential LBP formation, the conformational space of the receptor was sampled using a molecular dynamics (MD) simulation. Docking of DHA into 50,000 LBD structures extracted from the simulation revealed the existence of a transient LBP that is capable to fully harbor the compound. The location of the identified pocket overlaps with the ligand-binding site suggested by NMR experiments. Structural analysis of the protein-ligand complex showed that only modest structural rearrangements within the Nurr1 LBD are required for LBP formation. These findings may support structure-based drug discovery campaigns for the development of receptor-specific agonists.     

In-Silico molecular docking and simulation studies on novel chalcone and flavone hybrid derivatives with 1, 2, 3-triazole linkage as vital inhibitors of Plasmodium falciparum dihydroorotate dehydrogenase

The structural motifs of chalcones, flavones, and triazoles with varied substitutions have been studied for the antimalarial activity. In this study, 25 novel derivatives of chalcone and flavone hybrid derivatives with 1, 2, 3-triazole linkage are docked with Plasmodium falciparum dihydroorotate dehydrogenase to establish their inhibitory activity against Plasmodium falciparum. The best binding conformation of the ligands at the catalytic site of dihydroorotate dehydrogenase are selected to characterize the best bound ligand using the best consensus score and the number of hydrogen bond interactions. The ligand namely (2E)-3-(4-{[1-(3-chloro-4-fluorophenyl)-1H-1, 2, 3-triazol-4-yl]methoxy}-3-methoxyphenyl-1-(2-hydroxy-4,6-dimethoxyphenyl)prop-2-en-1-one, is one the among the five best docked ligands, which interacts with the protein through nine hydrogen bonds and with a consensus score of five. To refine and confirm the docking study results, the stability of complexes is verified using Molecular Dynamics Simulations, Molecular Mechanics /Poisson–Boltzmann Surface Area free binding energy analysis, and per residue contribution for the binding energy. The study implies that the best docked Plasmodium falciparum dihydroorotate dehydrogenase–ligand complex is having high negative binding energy, most stable, compact, and rigid with nine hydrogen bonds. The study provides insight for the optimization of chalcone and flavone hybrids with 1, 2, 3-triazole linkage as potent inhibitors.     

Human dipeptidyl peptidase III: insights into ligand binding from a combined experimental and computational approach

Human dipeptidyl peptidase III (DPP III) is a zinc-exopeptidase with implied roles in protein catabolism, pain modulation, and defense against oxidative stress. To understand the mode of ligand binding into its active site, we performed molecular modeling, site-directed mutagenesis, and biochemical analyses. Using the recently determined crystal structure of the human DPP III we built complexes between both, the wild-type (WT) protein and its mutant H568N with the preferred substrate Arg-Arg-2-naphthylamide (RRNA) and a competitive inhibitor Tyr-Phe-hydroxamate (Tyr-Phe-NHOH). The mutation of the conserved His568, structurally equivalent to catalytically important His231 in thermolysin, to Asn, resulted in a 1300-fold decrease of k(cat) for RRNA hydrolysis and in significantly lowered affinity for the inhibitor. Molecular dynamics simulations revealed the key protein-ligand interactions as well as the ligand-induced reorganization of the binding site and its partial closure. Simultaneously, the non-catalytic domain was observed to stretch and the opening at the wide side of the inter-domain cleft became enhanced. The driving force for these changes was the formation of the hydrogen bond between Asp372 and the bound ligand. The structural and dynamical differences, found for the ligand binding to the WT enzyme and the H568N mutant, and the calculated binding free energies, agree well with the measured affinities. On the basis of the obtained results we suggest a possible reaction mechanism. In addition, this work provides a foundation for further site-directed mutagenesis experiments, as well as for modeling the reaction itself.     

The role of S477N mutation in the molecular behavior of SARS-CoV-2 spike protein: An in-silico perspective

The attachment of SARA-CoV-2 happens between ACE2 and the receptor binding domain (RBD) on the spike protein. Mutations in this domain can affect the binding affinity of the spike protein for ACE2. S477N, one of the most common mutations reported in the recent variants, is located in the RBD. Today's computational approaches in biology, especially during the SARS-CoV-2 pandemic, assist researchers in predicting a protein's behavior in contact with other proteins in more detail. In this study, we investigated the interactions of the S477N-hACE2 in silico to find the impact of this mutation on its binding affinity for ACE2 and immunity responses using dynamics simulation, protein–protein docking, and immunoinformatics methods. Our computational analysis revealed an increased binding affinity of N477 for ACE2. Four new hydrogen and hydrophobic bonds in the mutant RBD-ACE2 were formed (with S19 and Q24 of ACE2), which do not exist in the wild type. Also, the protein spike structure in this mutation was associated with an increase in stabilization and a decrease in its fluctuations at the atomic level. N477 mutation can be considered as the cause of increased escape from the immune system through MHC-II.     

Computational investigation of the HIV-1 Rev multimerization using molecular dynamics simulations and binding free energy calculations

The HIV Rev protein mediates the nuclear export of viral mRNA, and is thereby essential for the production of late viral proteins in the replication cycle. Rev forms a large organized multimeric protein-protein complex for proper functioning. Recently, the three-dimensional structures of a Rev dimer and tetramer have been resolved and provide the basis for a thorough structural analysis of the binding interaction. Here, molecular dynamics (MD) and binding free energy calculations were performed to elucidate the forces thriving dimerization and higher order multimerization of the Rev protein. It is found that despite the structural differences between each crystal structure, both display a similar behavior according to our calculations. Our analysis based on a molecular mechanics-generalized Born surface area (MM/GBSA) and a configurational entropy approach demonstrates that the higher order multimerization site is much weaker than the dimerization site. In addition, a quantitative hot spot analysis combined with a mutational analysis reveals the most contributing amino acid residues for protein interactions in agreement with experimental results. Additional residues were found in each interface, which are important for the protein interaction. The investigation of the thermodynamics of the Rev multimerization interactions performed here could be a further step in the development of novel antiretrovirals using structure based drug design. Moreover, the variability of the angle between each Rev monomer as measured during the MD simulations suggests a role of the Rev protein in allowing flexibility of the arginine rich domain (ARM) to accommodate RNA binding.     

Insight into the key active sites of afChiA1 based on molecular dynamics simulations and free energy calculations

In this study, the binding of the enzyme chitinase A1 (afChiA1) from the plant-type Aspergillus fumigatus with four potent inhibitors, allosamidin (ASM), acetazolamide (AZM), 8-chloro-theophylline (CTP) and kinetin (KIT) is investigated by molecular docking, molecular dynamics simulation and binding free energy calculation. The results reveal that the electrostatic interactions play an important role in the stabilisation of the binding of afChiA1 with inhibitors. Based on the binding energy of afChiA1-ligands, the key residues (Gln37 and Trp312) in the active binding pocket of the complex systems are confirmed by molecular mechanics/Poisson–Boltzmann surface area method, and the active inhibitors, ASM and AZM, both could form strong interaction with Gln37 and Trp312, and the non-active ligands, CTP and KIT, could not interact with these two residues, which is consistent with the result of experimental report. Then, it is identified that Gln37 and Trp312 should be one of the important active site residues of afChiA1.     

New molecular insights into CETP structure and function: a review

Cholesteryl ester transfer protein (CETP) is important clinically and is the current target for new drug development. Its structure and mechanism of action has not been well understood. We have combined current new structural and functional methods to compare with relevant prior data. These analyses have led us to propose several steps in CETP's function at the molecular level, in the context of its interactions with lipoproteins, e.g., sensing, penetration, docking, selectivity, ternary complex formation, lipid transfer, and HDL dissociation. These new molecular insights improve our understanding of CETP's mechanisms of action.     

Taxifolin prevents postprandial hyperglycemia by regulating the activity of α-amylase: Evidence from an in vivo and in silico studies

There has been a dramatic increase in the prevalence of diabetes mellitus (DM) and its associated complications globally. The postprandial stage of DM involves prompt elevation in the levels of blood glucose and α-amylase, a carbohydrate-metabolizing enzyme is mainly involved in the regulation of postprandial hyperglycemia. This study was designed to assess the ability of a well-known flavonoid, taxifolin (TFN), against postprandial hyperglycemia and its inhibitory effects on α-amylase activity through the assessment of therapeutic potentials of TFN in an alloxan-induced diabetic animal model. The binding potential TFN with an α-amylase receptor was also investigated through molecular dynamics (MD) simulation and docking of to compare the binding affinities and energies of TFN and standard drug acarbose (ACB) with target enzyme. TFN significantly improved the postprandial hyperglycemia, lipid profile, and serum levels of α-amylase, lipase, and C-reactive protein in a dose-dependent manner when compared with that of either DM-induced and ACB-treated alloxan-induced diabetic rats. Moreover, TFN also enhanced the anti-oxidant status and normal functioning of the liver in alloxan-induced diabetic rats more efficiently as compared to that of ACB-treated alloxan-induced diabetic rats. Therapeutic potentials of TFN were also verified by MD simulation and docking results, which exhibited that the binding energy and affinity of TFN to bind with receptor was significantly higher as compared to that of ACB. Hence, the results of this study signify that TFN might be a potent inhibitor of α-amylase that has the potential to regulate the postprandial hyperglycemia along with its anti-inflammatory and anti-oxidant properties during the treatment of DM.     

Structural insights into a low-spin myoglobin variant with bis-histidine coordination from molecular modeling

Rational design of functional enzymes is a powerful strategy to gain deep insights into more complex native enzymes, such as nitric oxide reductase (NOR). Recently, we engineered a functional model of NOR by creating a two His and one Glu (2‐His‐1‐Glu) non‐heme iron center in sperm whale myoglobin (swMb L29E, F43H, H64, called Fe_BMb(‐His)). It was found that Fe_BMb(‐His) adopts a low‐spin state with bis‐His coordination in the absence of metal ions binding to the designed metal center. However, no structural information was available for the variant in this special spin state. We herein performed molecular modeling of Fe_BMb(‐His) and compared with the X‐ray structure of its copper bound derivative, Cu(II)‐CN^?‐Fe_BMb(‐His), resolved recently at a high resolution (1.65 Å) (PDB entry 3MN0). The simulated structure shows that mutation of Leu to Glu at position 29 in the hydrophobic heme pocket alters the folding behavior of Mb. The hydrogen bond between Glu29 and His64 further plays a role in stabilizing the bis‐His (His64/His93) coordination structure. This study offers an excellent example of using molecular modeling to gain insights in rational design of both structural and functional proteins. Proteins 2011. © 2010 Wiley‐Liss, Inc.