Week 6

Modeling of small molecule-protein interactions

 

 

I. .Basic Operations of SYBYL
II. Basic Docking
III. Building an Inhibitor
IV. Advanced Docking

 

I. BASIC OPERATIONS OF SYBYL

The SYBYL package is optimized for operations with small molecules and their interactions with biopolymers. The program operation is similar to the InsightII and is menu-driven using an intuitive command structure. We will first practice some basic operations that allow you to manipulate biopolymers and their complexes.

Read in a molecule

File

Read

Select 1AQ7.pdb (This file contains the structure of an inhibitor bound to trypsin). (Read about this structure here)

Click OK

A dialog box appears asking "Center the molecule"

Click OK

The molecule is now displayed on the screen. You can control the size and orientation with the mouse buttons identically to how you do it in InsightII.

Use the "center-right" mouse buttons to fit the entire molecule to the screen size

We will now center the rotation axis around the inhibitor.

Build/Edit

Center

Click on the "Sets" button

In the new window, Click on the UNK_ATOMS

Then click OK in the two windows.

We will now 'undisplay' the water molecules.

View

Undisplay Atoms

Click on the "Sets" button

In the new window, Click on the Water

Then click OK in the two windows.

Now we will differentially color the enzyme and the inhibitor

View

Color

Atoms

Click on the "Sets" button

In the new window, Click on the UNK_ATOMS

Then click OK in the two windows.

A new window appears. Select a color (Red is good)

View

Color

Atoms

Click on the "Sets" button

In the new window, Click on the UNK_ATOMS

Click OK in this window.

In the selection window, click on the 'INVERT' Button

Click OK in this window.

A new window appears. Select a color (Green is good)

Now we will display the different parts of the enzyme active site: the residues responsible for substrate binding specificity and catalytic activity

View

Color

Atoms

Click on the "Substructures" button

In the new window, scroll down to Asp 189 and select it with your left mouse button

Click OK in this window.

Then click OK in the two windows.

A new window appears. Select a color (Yellow is good)

We have just highlighted the residue that recognizes the positively charged side chains in trypsin substrates and inhibitors. Now we will examine the "catalytic triad"-i.e., the residues that perform the actual bond cleavage in all serine protease.

View

Color

Atoms

Click on the "Substructures" button

In the new window, select His57, Asp102 and Ser195

Click OK in this window.

Then click OK in the two windows.

A new window appears. Select a color (White is good)

You can now use the mouse buttons to move the molecule and examine the interaction of the inhibitor with the enzyme-its binding site and why it functions as an inhibitor.

Now we will learn basic docking techniques. Before doing that we will have to clear the screen.

Build/Edit

Zap/Delete Molecule

Go to top


II. BASIC DOCKING

Read in a molecule

File

Read

Select 1AQ7.pdb (This file contains the structure of an inhibitor bound to trypsin)

Click OK

A dialog box appears asking "Center the molecule"

Click OK

Center the rotation axis around the inhibitor.

Build/Edit

Center

Click on the "Sets" button

In the new window, Click on the UNK_ATOMS

Then click OK in the two windows.

 

In order to work with one molecule in the absence of another, the molecules must be in different work areas. The trypsin molecules and the inhibitor are stored in the same file in the database, so before we can dock the inhibitor, we have to extract it.

In order to extract a molecule, we have to have somewhere to put it. SYBYL has an option that allows you to divide the screen into several work areas, each of which can hold different molecules. We will start by dividing the screen into halves

Dividing the screen

On the left side of the SYBYL main window there is a series of buttons. The third one down looks like this: .

 

Click the button and then click the "Half" button, followed by the Q (quit) button.

Copy the cofactor to another molecule area

Build/Edit

Extract

Click on the "Sets" button

In the new window, Click on the UNK_ATOMS

Then click OK in the two windows.

A new window appears.

Select M2:<empty> as the molecule area

Click OK

Now we will delete the inhibitor from the original complex.

Build/Edit

Delete

Click on the "Sets" button

In the new window, Click on UNK_ATOMS.

Then click OK in the two windows.

Now we will name our inhibitor so that it will be easier to keep track of.

Build/Edit

Name molecule

Select M2

Click OK

Type inhibitor in the name block

Click OK

The docking process evaluates the suitability of an individual position/orientation by evaluating the interaction energies. As you already know, hydrogen bonds and charge interactions are important in molecular interactions. Thus, we must prepare the molecules to have the appropriate charges and numbers of hydrogens.

The first step in adding hydrogens is to ensure that SYBYL understands what kind of bonds are connecting the atoms in the molecule of interest.

Build/Edit

Modify

Bond

Select 'Type' as the element to modify.

Click OK

Select M2:inhibitor as the molecule

Then press the Bond Types button

In the next window, select uk (unknown

Click OK

Click OK in the Bond Expression window

The bond between the oxygen and the benzene moiety is highlighted in green, and a new window appears with a series of bond types 1 (=single), 2 (=double), etc.

Select 1 and click OK

The window disappears indicating that SYBYL understands all the other bond types in the molecule.

Now we will actually add the hydrogens.

Build/Edit

Add

Hydrogens

Select M2:inhibitor in the molecule area window

Click OK

Hydrogens are now added to the inhibitor.

Now we must assign the partial charges to the inhibitor atoms.

Compute

Charges

Gasteiger-Hò ckel

Select M2:inhibitor in the molecule area window

Click OK

A dialog box appears asking whether you want to change formal charges before computing charges

Select NO

Before beginning the docking process, color the molecule as we had done before.

View

Color

Atoms

Click on the "Sets" button

In the new window, Click on the UNK_ATOMS

Click OK in this window.

In the selection window, click on the 'INVERT' Button

Click OK in this window.

A new window appears. Select a color (Green is good)

View

Color

Atoms

Click on the "Substructures" button

In the new window, scroll down to Asp 189 and select it with your left mouse button

Click OK in this window.

Then click OK in the two windows.

A new window appears. Select a color (Yellow is good)

View

Color

Atoms

Click on the "Substructures" button

In the new window, select His57, Asp102 and Ser195

Click OK in this window.

Then click OK in the two windows.

A new window appears. Select a color (White is good)

There is a problem with the default settings of SYBYL. Before we can start the Docking procedure, we have to fix it.

Options

Tailor

Set

Select Dock in the Options window-Click OK

Select MAX_LATT_POINTS in the Tailor Option window-Click OK

A new window appears:

Change the value from 200000 to 400000-Click OK

Click END in the two Option windows

Now we are finally ready to begin docking!

To facilitate the docking, change the screen back to a single window:

Click the button and then click the "Full" button, followed by the Q (quit) button.

Compute

Dock

Select M2:inhibitor as the Ligand Molecule

Click OK

Select M1: as the Site Molecule

Click OK

Three new items appear on the screen.

Ordinarily, mouse movements move both the ligand and the target together. To have the mouse control only the ligand, press the F10 key on the keyboard

Now the translations and rotations only affect the ligand's movements.

Slowly move the ligand around and attempt to bring it into an energy minimum.

When you are done, exit back to the main menus by placing the cursor in the terminal window and pressing the ENTER key.

To return the appropriate mouse controls, click the button on the left row of buttons and select Global followed by the Q (quit) button.

To be sure the program behaves properly, we must 'reset' it. Do this by clicking the 'Reset' button on the left row of buttons and then click the 'Everything' button followed by the Q (quit) button.

Now we will learn how to build our own inhibitor and dock it into trypsin. Before doing that we will have to clear the screen.

Build/Edit

Zap/Delete Molecule

Go to top

III. BUILDING YOUR OWN INHIBITOR

We will be building a benzene-based inhibitor of trypsin from benzene and arginine. (Read about the properties of this inhibitor here)

Build/Edit

Get Fragment

Scroll down and select Benzene from the option list

Click OK

Build/Edit

Get Fragment

Scroll down and select Arginine from the option list

Click OK

Select M2 as the molecule area

Click OK

We will now edit the arginine molecule to obtain the guanidine groups.

To move the arginine out where it can be operated on, click F10, and use the mouse to move the group out of the way.

To make life easier in terms of knowing which atoms to delete, we will label them.

View

Label

Atom name

Click on M2 in the Atom Expression window

Click the Atom Types button

In the next window, select C,O,N

Click OK

Click OK

Now we will create the guanidium group

Build/Edit

Delete

Bond

Click on the N2 nitrogen and the C5 carbon in the molecule display window.

Click OK

Now we will delete the "remainder" of the arginine residue.

Build/Edit

Delete

Substructure

Click on the Arg side chain in the molecule display window

Click OK

Now we will prepare for the 'chemical' joining of the benzene and the guanidium group.

Build/Edit

Delete

Atom

Click the proton on the N2 nitrogen in the molecule display window.

Click OK

Now join the benzene and guanidinium moieties.

Build/Edit

Other Tools

Join molecules

When the select Target window appears, click on one of the benzene protons

When the select Target window appears, click on the N2 nitrogen of the guanidinium group

Repeat this procedure placing the guanidinium group on the benzene at a para position with respect to the first guanidinium group.

The final molecule will look like this:

Now we will save this inhibitor

FILE

Save As

Select Benzene

Call the file inhibitor1

Click OK

Build/Edit

Zap/Delete Molecule

Go to top

IV. ADVANCED DOCKING

Now we will dock our inhibitor into the trypsin molecule.

File

Read

Select 1AQ7.pdb (This file contains the structure of an inhibitor bound to trypsin)

Click OK

A dialog box appears asking "Center the molecule"

Click OK

Build/Edit

Delete

Atom

Click on the "Sets" button

In the new window, Click on the Water

Then click OK in the two windows.

Now we will add H to the protein

Biopolymer

Add Hydrogens

Click on the All button

Click OK

In Option window (which hydrogens to add), select All

Click OK

Now we will load the protein's partial charges

Biopolymer

Load Charges

Click on the All button

Click OK

In Option window accept the Koll_Uni default

Click OK

To facilitate building of the new inhibitor complex, we will superimpose the new inhibitor on the old one and then delete the old one and begin the docking procedure.

View

Undisplay Atoms

Click the Sets button

Click OK

Click the Invert button

Click OK

Now read in the new inhibitor

File

Read

Select inhibitor1

Click OK

To help distinguish one from the other color the new inhibitor cyan

View

Color

Atoms

Click on M2

Click the ALL button

Click OK

Choose cyan and click OK

To help you in aligning the correct atoms, label the molecules

View

Label

Atom name

Click on M1 in the Atom Expression window

Click the Atom Types button

In the next window, select C,O,N

Click OK

Click OK

Repeat this procedure using M2 as the atom option.

Try using to align the one of guanidinium groups of the new inhibitor the guanidinium group o the old inhibitor. Use the F10 key to toggle mouse control between the new inhibitor and both molecules to aid this process.

Now we will re-display the protein.

View

Display Atoms

Click on M1

Click on the All button

Click OK

Now delete the old inhibitor

Build/Edit

Delete

Atom

Be sure M1 is selected

Click on Set

Select UNK_ATOMS

Click on OK in this and previous window.

Now Dock your inhibitor using the Docking procedure above.

Go to top