Souring control in fluid samples of oil industry using a multiple ligand simultaneous docking (MLSD) strategy

We have used docking techniques in order to propose potential inhibitors to the enzymes adenosine phosphosulfate reductase and adenosine triphosphate sulfurylase that are responsible, among other deleterious effects, for causing souring of oil and gas reservoirs. Three candidates selected through molecular docking revealed new and improved polar and hydrophobic interactions with the above-mentioned enzymes. Microbiological laboratory assays performed subsequently corroborated the results of computer modelling that the three compounds can efficiently control the biogenic sulfide production.     

Towards an Exhaustive Sampling of the Configurational Spaces of the Two Forms of the Peptide Hormone Guanylin

Abstract

The recently introduced Essential Dynamics sampling method is extended such that an exhaustive sampling of the available (backbone) configurational space can be achieved. From an initial Molecular Dynamics simulation an approximated definition of the essential subspace is obtained. This subspace is used to direct subsequent simulations by means of constraint forces. The method is applied to the peptide hormone guanylin, solvated in water, of which the structure was determined recently. The peptide exists in two forms and for both forms, an extensive sampling was produced. The sampling algorithm fills the available space (of the essential coordinates used in the procedure) at a rate that is approximately six to seven times larger than that for traditional Molecular Dynamics. The procedure does not cause any significant perturbation, which is indicated by the fact that free Molecular Dynamics simulations started at several places in the space defined by the Essential Dynamics sample that complete space. Moreover, analyses of the average free Molecular Dynamics step have shown that nowhere except close to the edge of the available space, there are regions where the system shows a drift in a particular direction. This result also shows that in principle, the essential subspace is a constant free energy surface, with well-defined and steep borders, in which the system moves diffusively. In addition, a comparison between two independent essential dynamics sampling runs, of one form of the peptide, shows that the obtained essential subspaces are virtually identical.     

Electron Microscopic and Physico-Chemical Studies of DNA Complexes with Synthetic Oligopeptides: Binding Specificity and DNA Compact Structures

Abstract

Binding to DNA of two synthetic peptides, Val-Thr-Thr-Val-Val-NH-NH-Dns and Thr-Val- Thr-Lys-Val-Gly-Thr-Lsy-Val-Gly-Thr-Val-Val-NH-NH-Dns (where Dns is a residue of 5- dimethylaminonaphthalene-l-sulfonic acid), has been studied by circular dichroism, electron microscopy and fluorescence methods. It has been found that these two peptides can self- associate in aqueous solution as follows from the fact that concentration-dependent changes are observed in the UV absorbance and fluorescence spectra. The two peptides can bind to DNA both in self-associated and monomeric forms. The pentapeptide in the β-associated form binds more strongly to poly(dG) · poly(dC) than to poly[d(A-C)] · poly[d(G-T)] and poly(dA) · poly(dT) whereas the tridecapeptide exhibits an opposite order of preferences binding more strongly to poly[d(A-C)] · poly[d(G-T)] and poly(dA) · poly(dT) than to poly(dG) · poly(dC).

Binding is a cooperative process which is accompanied by the DNA compaction at peptide/DNA base pair ratios greater than l. At the initial stage of the compaction process, the coalescence of DNA segments covered by bound peptide molecules leads to the formation of DNA loops stabilized by the interaction between peptide molecules bound to different DNA segments. Further increase in the peptide/DNA ratio leads to the formation of rod-like structures each consisting of two or more double-stranded DNA segments. The final stage of the compaction process involves folding of fibrillar macromolecular complexes into a globular structure containing only one DNA molecule.     

Study on the interactions between diketo-acid inhibitors and prototype foamy virus integrase-DNA complex via molecular docking and comparative molecular dynamics simulation methods

Human immunodeficiency virus type 1 (HIV-1) integrase (IN) is an important drug target for anti-acquired immune deficiency disease (AIDS) treatment and diketo-acid (DKA) inhibitors are potent and selective inhibitors of HIV-1 IN. Due to lack of three-dimensional structures including detail interactions between HIV-1 IN and its substrate viral DNA, the drug design and screening platform remains incompleteness and deficient. In addition, the action mechanism of DKA inhibitors with HIV-1 IN is not well understood. In view of the high homology between the structure of prototype foamy virus (PFV) IN and that of HIV-1 IN, we used PFV IN as a surrogate model for HIV-1 IN to investigate the inhibitory mechanism of raltegravir (RLV) and the binding modes with a series of DKA inhibitors. Firstly, molecular dynamics simulations of PFV IN, IN-RLV, IN-DNA, and IN-DNA-RLV systems were performed for 10?ns each. The interactions and inhibitory mechanism of RLV to PFV IN were explored through overall dynamics behaviors, catalytic loop conformation distribution, and hydrogen bond network analysis. The results show that the coordinated interactions of RLV with IN and viral DNA slightly reduce the flexibility of catalytic loop region of IN, and remarkably restrict the mobility of the CA end of viral DNA, which may lead to the partial loss of the inhibitory activity of IN. Then, we docked a series of DKA inhibitors into PFV IN-DNA receptor and obtained the IN-DNA-inhibitor complexes. The docking results between PFV IN-DNA and DKA inhibitors agree well with the corresponding complex of HIV-1 IN, which proves the dependability of PFV IN-DNA used for the anti-AIDS drug screening. Our study may help to make clear some theoretical questions and to design anti-AIDS drug based on the structure of IN.     

Design of inhibitors of thymidylate kinase from Variola virus as new selective drugs against smallpox: part II

Abstract

Acknowledging the importance of studies toward the development of measures against terrorism and bioterrorism, this study aims to contribute to the design of new prototypes of potential drugs against smallpox. Based on a former study, nine synthetic feasible prototypes of selective inhibitors for thymidylate kinase from Variola virus (VarTMPK) were designed and submitted to molecular docking, molecular dynamics simulations and binding energy calculations. The compounds are simplifications of two more complex scaffolds, with a guanine connected to an amide or alcohol through a spacer containing ether and/or amide groups, formerly suggested as promising for the design of selective inhibitors of VarTMPK. Our study showed that, despite the structural simplifications, the compounds presented effective energy values in interactions with VarTMPK and HssTMPK and that the guanine could be replaced by a simpler imidazole ring linked to a –NH₂ group, without compromising the affinity for VarTMPK. It was also observed that a positive charge in the imidazole ring is important for the selectivity toward VarTMPK and that an amide group in the spacer does not contribute to selectivity. Finally, prototype 3 was pointed as the most promising to be synthesized and experimentally evaluated.

Communicated by Ramaswamy H. Sarma     

KAISO inhibition: an atomic insight

In today’s world, the pursuit of a novel anti-cancer agent remains top priority because of the fact that the global burden of this malady is continuously increasing. Our work is no different from others in searching for new therapeutic solutions. To achieve this, we are looking into Epigenetics, the phenomenon governed by hypermethylation and hypomethylation of tumor suppressor genes and oncogenes, respectively. Our target for this study is an important intermediary methyl-CpG binding protein named kaiso. In our study, we have used the X-ray crystallographic structure of Kaiso for virtual screening and molecular dynamics simulations to study the binding modes of possible inhibitors. The C2H2 domain comprising LYS539 was used for screening the inter bio screen Database having 48,531 natural compounds. Our approach of using computer-aided drug designing methods helped us to remove the execrable compounds and narrowed our focus on a selected few for molecular simulation studies. The top ranked compound (chem. ID 28127) exhibited the highest binding affinity and was also found to be stable throughout the 20 ns timeframe. This compound is therefore a good starting point for developing strong inhibitors.     

In silico 3D structure modeling and inhibitor binding studies of human male germ cell-associated kinase

Human male germ cell-associated kinase (hMAK) is an androgen-inducible gene in prostate epithelial cells, and it acts as a coactivator of androgen receptor signaling in prostate cancer. The 3D structure of the hMAK kinase was modeled based on the crystal structure of CDK2 kinase using comparative modeling methods, and the ATP-binding site was characterized. We have collected five inhibitors of hMAK from the literature and docked into the ATP-binding site of the kinase domain. Solvated interaction energies (SIE) of inhibitor binding are calculated from the molecular dynamics simulations trajectories of protein–inhibitor complexes. The contribution from each active site residue in hMAK toward inhibitor binding revealed the nature and extent of interactions between inhibitors and individual residues. The main chain atoms of Met79 invariably form hydrogen bonds with all five inhibitors. The amino acids Leu7, Val15, and Leu129 stabilize the inhibitors via CH–pi interactions. The Asp140 in the active site and Glu77 in hinge region show characteristic hydrogen bonding interactions with inhibitors. From SIE, the residue-wise interactions revealed the nature of non-bonding contacts and modifications required to increase the inhibitor activity. Our work provides 3D model structure of hMAK and molecular basis for the mechanisms of hMAK inhibition at atomic level that aid in designing new potent inhibitors.     

Modelling family 2 cystatins and their interaction with papain

Cystatins are extensively studied cysteine protease inhibitors, found in wide range of organisms with highly conserved structural folds. S-type of cystatins is well known for their abundance in saliva, high selectivity and poorer activity towards host cysteine proteases in comparison to their immediate ancestor cystatin C. Despite more than 90% sequence similarity, the members of this group show highly dissimilar binding affinity towards papain. Cystatin M/E is a potent inhibitor of legumain and papain like cysteine proteases and recognized for its involvement in skin barrier formation and potential role as a tumor suppressor gene. However, the structures of these proteins and their complexes with papain or legumain are still unknown. In the present study, we have employed computational methods to get insight into the interactions between papain and cystatins. Three-dimensional structures of the cystatins are generated by homology modelling, refined with molecular dynamics simulation, validated through numerous web servers and finally complexed with papain using ZDOCK algorithm in Discovery Studio. A high degree of shape complementarity is observed within the complexes, stabilized by numerous hydrogen bonds (HB) and hydrophobic interactions. Using interaction energy, HB and solvent accessible surface area analyses, we have identified a series of key residues that may be involved in papain–cystatin interaction. Differential approaches of cystatins towards papain are also noticed which are possibly responsible for diverse inhibitory activity within the group. These findings will improve our understanding of fundamental inhibitory mechanisms of cystatin and provide clues for further research.     

Transmembrane helix 6 observed at the interface of β2AR homodimers in blind docking studies

Peptide– and protein–protein dockings were carried out on β₂-adrenergic receptor (β₂AR) to confirm the presence of transmembrane helix 6 (TM6) at the interface region between two β₂AR monomers, thereby its possible role in dimerization as suggested in numerous experimental and computational studies. Initially, a portion of TM6 was modeled as a peptide consisting of 23 residues and blindly docked to β₂AR monomer using a rigid body approach. Interestingly, all highest score conformations preferred to be near TM5 and TM6 regions of the receptor. Furthermore, longer peptides generated from a whole TM region were blindly docked to β₂AR using the same rigid body approach. This yielded a total of seven docked peptides, each derived from one TM helix. Most interestingly, for each peptide, TM6 was among the most preferred binding site region in the receptor. Besides the peptide dockings, two β₂AR monomers were blindly docked to each other using a full rigid-body search of docking orientations, which yielded a total of 16,000 dimer conformations. Each dimer was then filtered according to a fitness value based on the membrane topology. Among 149 complexes that met the topology requirements, 102 conformers were composed of two monomers oriented in opposite directions, whereas in the remaining 47, the monomers were arranged in parallel. Lastly, all 149 conformers were clustered based on a root mean-squared distance value of 6 Å. In agreement with the peptide results, the clustering yielded the largest population of conformers with the highest Z-score value having TM6 at the interface region.     

A comprehensive in silico analysis of huntingtin and its interactome

A polyglutamine expansion of the N-terminal region of huntingtin (Htt) causes Huntington’s disease, a severe neurodegenerative disorder. Htt huge multidomain structure, the presence of disordered regions, and the lack of sequence homologs of known structure, so far prevented structural studies of Htt, making the study of its structure-function relationships very difficult. In this work, the presence and location of five Htt ordered domains (named from Hunt1 to Hunt5) has been detected and the structure of these domains has been predicted for the first time using a combined threading/ab initio modeling approach. This work has led to the identification of a previously undetected HEAT repeats region in the Hunt3 domain. Furthermore, a putative function has been assigned to four out of the five domains. Hunt1 and Hunt5, displaying structural similarity with the regulatory subunit A of protein phosphatase 2A, are predicted to play a role in regulating the phosphorylation status of cellular proteins. Hunt2 and Hunt3 are predicted to be homologs of two yeast importins and to mediate vescicles transport and protein trafficking. Finally, a comprehensive analysis of the Htt interactome has been carried out and is discussed to provide a global picture of the Htt’s structure–function relationships.