Molecule Building Examples and Strategies

The following examples illustrate how molecules are built in GaussView. As some of the tutorials build on previous material, be sure to save each file when finished.

1 - Building Pyridine: Using Untyped Atoms

The following is an example of how Pyridine is built in GaussView:

1-1 Start a new file by selecting the Create New Molecule Group option from the File menu in the main window.

1-2 Opens the Ring Fragment palette by double clicking on the Ring Fragment icon . Select benzene from the ring fragments; a benzene molecule appears in the Current Fragment window. (Alternatively, you can select benzene from the drop-down list in the main toolbar.)

1-4 Return to the View window. Click once: benzene displays.

1-5 Double-click on Element on the Builder window to open the Element Fragment palette. Click on nitrogen within the periodic table to select it. Then click on the N Atom button at the bottom of the palette.

Do not select one of the hybridizations. It is usually easier to use untyped atoms.

1-7 In the View window, click on the carbon atom you want to change to nitrogen.

1-8 Select the Delete Atom button on the Builder window.

1-9 On the View window, delete the hydrogen attached to the nitrogen by clicking on it.

1-10 Click on the Rebond button on the Builder window. Then, click on the Clean button.

You will often want to rebond before cleaning.

Pyridine is now complete. Save your work.

2 - Building Phenylpyridine: Setting the Angle Between Rings

This example shows how to build Phenylpyridine based on an existing file, in this case, Pyridine.

2-1 We will be continuing from the previous example, Building Pyridine.

2-2 Double-click on Ring Fragment on the Builder window. Be sure that the benzene ring is selected.

2-3 Click on the hydrogen atom located on the carbon atom opposite the nitrogen atom on the View window. A second ring is attached to the first ring in place of the hydrogen.

2-4 Select the Modify Dihedral button on the Builder window.

2-5 Select atoms between the two rings, as shown in the following figure:

The SmartSlide window opens when four atoms have been selected; make sure that Rotate Groups is selected for both Atoms 1 and 4.

Modifying the Dihedral Angle

2-6 Move the slider to create a dihedral until the rings are perpendicular to each other (or enter 90.0 in the text box).

Click on the OK button.

Phenylpyridine is now complete.

3 - Building Iron Pentacarbonyl: Selecting the Best Templates for a Given Task

This example builds iron pentacarbonyl.

3-1 Start a new file by selecting the New=>Create Molecule Group option from the File menu in the View window. Select the element iron, and then choose the trigonal bipyramidal iron template:

3-2 Return to the View window. Click once: the iron molecule displays.

3-3 Double-click on R-Group Fragment on the Builder window to open the R-Group Fragment palette. Select the linear C-N template, and place this group on each bond to the iron atom. Note that we select this template rather than one of the C-O groups because the former is linear.

3-4 Change the nitrogen atoms to oxygen atoms: select the element oxygen and then the free atom template (using either the choices at the bottom of the Element Fragment palette or the drop-down list on the main screen).

The molecule is now complete.

4 - Building Tyrian Purple Dye: Appending a Structure

The following is an example of how a Tyrian purple dye 6,6'-dibromoindigo is built in GaussView:

4-1 Create a new file. Double-click on Ring Fragment on the Builder window to open the Ring Fragment palette.

4-2 Select the fused 5 and 6 carbon rings on the Builder window (see illustration immediately above), and then click on the View window to place a fragment.

4-3 Select the oxygen element and change the appropriate hydrogen atom in the five-member ring to oxygen by clicking on it.

4-4 Change the carbon atom opposite the oxygen to nitrogen.

4-5 Change the hydrogen atom on the outermost carbon atom in the six-member ring (on the same side as the nitrogen) to bromine in a similar manner.

4-6 Select Clean on the Builder window. Do not rebond. Save the file.

4-7 Append the file you just saved to the display on the View window: copy the file to make it appear in the Current Fragment window, and then click in the View window to paste the current fragment.

4-8 Rotate and move the second ring so that the two carbon atoms to be bonded are close together and oriented properly with respect to one another.

To affect just one of the two fragments, hold down the Alt key while dragging the mouse. Using a second window is often helpful when performing complex positioning operations. Use the Add View option on the View menu to create one.

4-9 Delete the hydrogen atoms on the two carbon atoms to be bonded.

4-10 Create a double bond between the two carbon atoms: select the Modify Bond icon on the Builder window, and then click on the two carbon atoms to be bonded. The SmartSlide window opens; choose the double bond option from the Bond-Type list, and click OK.

4-11 Click on Clean.

The 6,6'-dibromoindigo molecule is now complete.

Save your work.

5 - Docking Two Structures: Building Bichloro-Diphenyl

This example illustrates a method for building a molecular structure consisting of docked framents.

5-1 Create a new View window, and place a benzene ring in it. Change one of the hydrogen atoms to Chlorine. Repeat this process to create a second chlorobenzene ring.

5-2 Rotate one of the rings so that the chlorine atom is in the opposite position to the one in the first ring. You can accomplish this by holding down the Alt key while holding down the left mouse button and dragging:

5-3 Position the rings so that they overlap and the carbon atoms to be bonded are in the correct position:

5-4 Delete the hydrogen atoms and dangling bonds in the overlap region.

5-5 Create a bond between the two carbon atoms, completing the molecule:

5-6 Click on Clean; the molecule is now complete.

6 - Assigning Atoms to ONIOM Layers in a Small Peptide Chain

Note: This example is devised for its illustrative clarity rather than for scientific interest. Your molecule may look different from the figures below.

6-1 Create a peptide chain containing these amino acids: Phe-Ala-Gly: use the Biological Fragment palette to choose an Amino-Terminal Phenylalanine, a Central Fragment Alanine, and a Carboxyl-Terminal Glycine, adding each in turn to the the same View window.

6-2 Assign all atoms to the Low layer. Choose Select Layer from the Edit menu, and then click Select All. Verify that the current layer is set to Low. The display will look similar to the following figure.

Click Apply.

6-3 Assign the following atoms to the Medium layer: all atoms in Ala, carbons in Ring in Phe, and atoms up to the double bond in Gly. Use any selection method you want; it will usually take two or more operations to assign all these atoms. For example, you can select the central carbon atom in Ala and drag the Expand Selection slider to add more atoms. Adding too few atoms is usually easier than adding too many. Click Apply.

Add additional atoms as necessary (use the center mouse button [or Shift+Click] to select, and then click Apply) until the selection is correct as in the following figure. The atoms in the Medium layer are now displayed as tubes, and the ones in the Low layer are in wireframe in the display:

6-4 Assign the following atoms to the High layer: all atoms in Ala, including heavy atoms bonded to them. One selection method is to click the central carbon in Ala, change the Using field to Distance, and click twice at the right end of Expand Selection. The display will be as in the following figure:

Click Apply. The display is as in the following figure:

Completed ONIOM-Assigned Molecule

Remove the hydrogen atom that was just placed into the High layer by selecting it and then placing it in the Medium layer.

Click Close.

The atom layer assignment is now complete.

Save your work.

7 - Setting Up an ONIOM Job

This example uses the molecule that was built and the layers that were defined in the final example in the Building Tutorial, and it continues on from those steps.

7-1 Select Gaussian Calculation Setup from the Calculate menu.

7-2 On the Method tab, check Multilayer ONIOM Model. Three tabs appear (High, Medium, and Low) in the Method section, which you can use to select the methods for treating each layer of atoms.

7-3 Click on the High tab, and then select an RHF/6-31G job:

7-4 Click on the Medium tab, then select an AM1 job:

7-5 Click on the Low tab, and specify a molecular mechanics calculation using the Amber force field:

7-6 These settings place the following line in the job’s route section:

# ONIOM=(RHF/6-31G:AM1:Amber)

You can now go on to fill out the rest of the dialog box, and the job will be ready to be saved and run.