INTERACTION OF 3L3K PROTEIN WITH LIGANDS USING AUTODOCK
B.V.S ANIKETH (13BCE0031)
B.Tech Computer Science and Engineering
B. ROHITH NAIDU (13BCE0424)
B.Tech Computer Science and Engineering
SOURABH SHARMA (13BCE0791)
B.Tech Computer Science and Engineering
Abstract—Proteins are the workhouse molecules of life which take part in every essential structure and activity in our life. Ligand is a molecule which bonds to any central atom which makes the co-ordination complex. 3L3K is a 3 chain structure with sequence that is from humans which deals with their immune systems. We interact 3L3K with 3 ligands namely 6-deacetylnimbin, Nimbin pp, Nimbin mw which are related to medicinal plants mainly neem. Nimbin and other associated group of compounds were secluded from a methanol extract of neem oil at room temperature. Interaction of the protein and ligand will be performed and an analysis of their interaction will be available through this. AutoDock has shown great results for the application of practical viewing to a range of targets and is the software that will be used. Applications considered in our laboratory have ranged from the advance of enzyme inhibitors, compounds that steady protein demobilization, antitoxins, and nucleic acid-strap compounds. Through the use of basic and preset scripting methods such as virtual screening with AutoDock is at present accessible to a broad user area.
3L3K protein deals with the human system. The immune system is a system of various biological structures and functions within an organism that prevents against disease. To function accurately, an immune system must detect a broad variety of agents, known as pathogens, from viruses to parasitic worms, and differentiate them from the organism's own healthy tissue. Nimbin is the active compound from the Neem leaf extract which is our main ligand. Nimbin has some antiseptic properties and it can also used for manufacturing antiseptics/medicines and nimbin is an alkaloid. Nimbin is the first bitter compound remote from neem oil, more than 166 compounds have been found, which are split into two major divisions: Isoprenoids. The isoprenoids include diterpenoids and triterpenoids contain protomeliacin, limonoid, azadirones and its derivative, gedunin and its derivative, vilasinin type of compounds and Csecomeliacins such as nimbin, salanin and azadirachtin. Neem has anti-bacterial properties and cures that help in fighting against skin infections such as acne, psoriasis, eczema,scabies etc. Neem extracts also help in treating diabetes, AIDS, cancers, heart diseases, herpes, allergies, ulcers, hepatitis and several other diseases. There are quite a few neem uses, and they contain a allay for skin diseases, plus the flexibility to function as an aid to one’s hygiene, and even an insecticide. It can be used as a source of preventive oil that may be involved to safeguard against disease, plus it’s also majorly lauded as a aromic spice. Logically, the whole neem plant, including the fruit, leaves, and bark, has been cranked, dried or ground into fine powder to keep one safe and take care of illnesses. Because the general ability has been shown by the plant to eliminate viruses, fungi, and microorganisms that bring misshaping of the gum and seeds have been processed for medication. Interaction of these medicinal ligands related to neem will be known.
Molecular Formula C30H3606 and
Index Name 2H-Cyclopent naphtha 2, 3-d furan-10-acetic acid.
2.1 6-deacetylnimbin ( C28H34O8)
A limonoid that is nimbin in which the acetyloxy group at position 6 is replaced by a hydroxy group. It has been isolated from Azadirachta indica.
Limonin is a limonoid, and a dash of bitter, white, crystal substance found in citrus and other plants. Limonin is riched in citrus fruits and is often present higher concentrations in seeds, for example as of orange and lemon seeds.
2.2 Nimbin pp ( C30H36O11)
Chemical name is Isonimbinolide. Ploysaccharides from bark have high antitumor and anti-inflammable action. Bark is regarded as a bitter tonic, contracting,useful in fever,thirstness,nauseas, vomiting and many skin diseases. Almost every area of the tree has been worn for treatment of variety of human cures, majorly against bacterial and fungal diseases. Various factors derived from neem contain various biological effects against insects like repellence, phagode.The leaves have been used in the cure of ulcer, eczema.
2.3 Nimbin mw ( C30H36O10)
Nimbin mw is a triterpenoid isolated from Neem. Nimbin is thought to be a major factor for much of the major biological activities of neem oil, and other essentials. It has Hydrogen Bond donor count of 0 and Hydrogen Bond acceptor count of 10. It has rotatable bond count of 8.
3 REQUIREMENTS AND METHODS
3.1 Data set definition and preparation
The crystalline structure of 3L3K is obtained from the RCSB Protein Data Bank.
The 3 ligands namely 6-deacetylnimbin, Nimbin pp, Nimbin mw were obtained from PCB NCBI website. So we have 1 protein and 3 ligands.
3.2 Virtual Screening
Virtual screening (VS) is a computational performance technique used in drug detection to search libraries of minute molecules in order to spot those structures which are most likely to unite to a drug target, typically a protein receptor or enzyme. Virtual screening techniques play vital role in the lead innovation process. Virtual screening is a computer-based technique for identification of capable compounds to fasten to a target molecule of known structure. Given the rapidly escalating number of protein and nucleic acid structures, virtual screening continues to grow as an efficient method for the sighting of new inhibitors and drug molecules.
3.3 Software Required
Autodock and AutoGrid Suite, MGL Tools 1.5.6, Open Babel, Python molecular viewer, Discovery Visual Studio.
Operating system: Microsoft windows 8.1 / Ubuntu 14
Download the ligands from the NCBI website and save them in sdf format. Now they are converted to pdbqt layout in openBabael and save in the working directory.
The receptor is downloaded from the RCBS website by its PDB ID and saved in pdb format in the working directory. Unrequired chains are deleted from it by virtual screening techniques in Discovery Visual Studio software and saved.
3.5 Run Time Setup
Open MGL tools and create a working directory. Then read the pdb file of the receptor. Delete selected atoms of chains that are not required and delete water, add hydrogens and merge non polar. Add gastegier Charge.
Write the PDB and save it. Load the pdbqt file of the receptor, detect its root and we are ready for auto grid. If Structure is in 2D convert it to 3D using open Babel software.
3.6 Auto Grid
Choose the receptor as the macromolecule and chose the ligand and save as pdbqt. Set a covalent map if required. Set a grid box of 60X60X60 and spacing, center if required. Save as .gpf format and go to run and launch auto grid. Wait for the process manager to close and the gpf file is created in our working directory. We are ready for auto dock. Any sort of error would be viewed in the terminal window.
3.7 Auto Dock
Autodock is a suite of software for predicting the optimal bound conformations of ligands to proteins. The preliminary applications of AutoDock were in the analysis of binding modes and catalytic properties of protein and nucleic acid complexes, and a typical study would include results from numerous dozen docking simulations. More recently, however, enhancement in the performance of AutoDock combined with the accessibility of high speed computers and cluster of computers has allowed much larger experiments, where whole compound libraries are screened against pharmaceutically-related targets
Choose macromolecule and set rigid filename and select the .pdbqt macromolecule created during the setup of autogrid. In docking, opt the ligand and click defaults. For search parameters use Genetic Algorithm and set number of GA runs as 10 (higher is suggested for better results but consumes more time) and use remaining defaults. Set docking output as Lamarckian GA (4.2) and save the output as .dpf. Now run AutoDock and click launch and wait for the program manager to close (may take some time). Any sort of error would be displayed in the terminal window.
Once AutoDock completes the calculations, you can then set off to: Analyze > Confirmations > Play, ranked by energy. This will show all binding confirmations in order of energy. At this point, you may analyze the results, connections, and confirmations as needed.
We note down the results after running the interactions between the protein 3L3K and ligands. The values and their results are depicted in tabular fashion below. We have mainly tabulated the binding energy and inhibition constant.
S.no Protein_Ligand Binding Energy Inhibition Constant
1. 3L3K- 6deacetylnimbin -5.02 208.26
2. 3L3K-nimbin pp -3.58 2.38
3. 3L3K-nimbin mw -4.64 399.6
5.1) 6-deacetylnimbin with 3L3K
The above figure shows the values of binding energy, ligand efficiency, inhibation constant, intermol energy, torsional energy and others.
On increase of the value the ligand moves over the protein and the values of the binding energy changes.
5.2) Nimbin pp with 3L3K
The ligand interaction is clearly visible with the yellow portion moving throughout the protein on increase of the running values.
5.3) Nimbin mw with 3L3K
The interaction between the protein 3L3K and the ligands related to medicinal structure have performed and analyzed with their values on their interaction. Auto Dock is a popular non-commercial docking program that docks a ligand to its protein and performs well and fast it is user-friendly for docking ligands with target protein. In this paper we propose an easier user-friendly docking protocol for docking ligands with target protein that utilizes Auto Dock . Auto Dock helps in reducing the cost and time for drug discovery process which otherwise takes many years. Virtual screening and docking studies helped to obtain ligand molecules that can inhibit the important proteins involved. Ligand can react with the protein in all patterns and their values can be depicted in change with their binding energy. If the range of binding energy varies from -3 to -6 then the result would be a clear value that is required. Using all these Auto dock and open babel tools we can infer the binding energy and inhibition constant values of the given ligands and proteins.
(Detection of Ligand interaction points using Cygwin)
We would like to thank Prof. Kannadasan R for supporting our work throughout.
V. Edible Medicinal and Non Medicinal Plants: Volume 8 Flowers
VI. NEEM TODAY AND IN THE NEW MILLENNIUM - edited by Opender Koul andSeema Wahab
VII. Anderson AC. The process of the structure based drug design. Chem Biol 2003;
VIII. Lazarova M. Virtual screening-models, methods and software systems. InternationalScientific Conference of Computer Science 2008;.
IX. Morris GM, Goodsell DS, Halliday RS, Huey R, Hart WE, Belew RK et al. Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function. J Comput Chem 1998;19: 1639–62.
X. Warren G, Andrews C, Capelli A, Clarke B, LaLonde J, Lambert MH et al. A critical assessment of docking programs and scoring functions. J Med Chem 2006;49:5912-31.
...(download the rest of the essay above)