Trp-(Pro)3-Trp peptide (step by step)

You can check here the GAMESS manual


We propose an example of hole transfer mechanism in order to explain how calculate the electron coupling by e-coupling server. In the model Trp-(Pro)3-Trp a hole is transfered from one radical Trp to the other Trp.


    1. Prepare the system in its red/red form according with the hole transfer mechanism.

    2. Run an optimization of the molecule. In this example the molecule has charge 0, multiplicity 1 and we optimized at Hartree Fock level of theory (HF) and basis set 6-31g*.

    3. Note. Choose an appropiate level of theory and basis set accurate to your system.

      The number of optimization steps depends on the system. Not enough optimization steps could give a spin contaminated wave function. The Spin operator <S2> measures the spin contamination and estimates the quality of the wave function.

      S = (M-1)/2 (M = multiplicity)

      <S2> = S*(S+1) (0, ¾, 2... etc)

      For example, in a doublet state the spin operator without contamination should have a value of 0.75. If the computed value differs largely from this number ( >0.05), consider reoptimizing the wave function.


    1. Include the following keywords in the input:

      NPRINT=3 into the CONTROL group to print the Overlap matrix

      IEDEN=1 in the ELDENS block to compute the electron density

      IEDINT=1 in theELDENS block to print electron density integrals

      MORB=0 in the ELDENS block to compute the total electron density in the molecular orbitals

      IEMINT=1 IEMOM=1 in the ELMOM group to calculate and print the dipole matrix (GMH)

    2. Run a single point calculation from the optimized structure, using the same theory level and basis set.

      Check the full GAMESS input file for the Trp-Trp peptide example.


    From the output file identify donor and acceptor atom numbers and localize orbitals involved in the electronic coupling using a visualization program like Molden or JSMOL (directly from e-coupling server).

    1. Identify atom numbers from donor and acceptor

      Load the structure (.log) and click button label >> atom number

      Load the structure (.log) and click right button in the mouse style >> label >> with atom number. In this example, first Trp (donor) is from atom 1 to 15 and the second (acceptor) from 16 to 30.

    2. Visualize orbitals involved in the reaction.

      From the loaded structure press Density mode >> Orbital >> and click in the window Molden Orbital Select to visualize the orbitals. The reaction in the example is hole transfer and for this reason, we visualize the HOMOs.

      Click the buttons to visualize different HOMO/LUMO orbitals. Detail information about each orbital appears in the left corner.

    3. Select orbitals involved in the electronic coupling

      Select the number of more stable degenerated HOMOs (in case of hole transfer) or LUMOs (electron transfer) placed in both donor (Trp5) and acceptor (Trp1).

      Consider degenerated orbitals when the difference is lower than 0.005 eV (or decide the value using your own criteria and depending on your system)

      In the example, HOMO is placed on the donor (Trp5) and the HOMO-1 is placed on the acceptor (Trp1). The orbital HOMO-2 is placed on the acceptor (Trp 1) but is not degenerated (comparing with the HOMO-1) and we do not consider it in the calculation. We use the same reasoning to reject HOMO-3.

      As a conclusion, in the example we take one orbital from donor (HOMO) and one for the acceptor(HOMO-1)