Maestro is a part of the Schrödinger program suite.
However, in variance to most of the other programs, it is free to
academic users (but then, you cannot fill in missing residues or
predict pKa values of ligands).
Ensure that you download it from https://www.schrodinger.com/freemaestro.
The current version is Maestro 2020-1.
The files are located in /temp4/bio/Schrodinger/Maestro-2020-1 and
it is run by maestro.
However, currently it seems to work only if you have own your own
installation (available also for Mac and PC laptops), i.e. you need to run the current installation as user bio.
and then run the program with
You get a number of warnings, but the program seems to work
Maestro is particularly strong for the setup of a new protein and
this is currently our prime use of the program.
Below, I describe a tentative scheme for the setup of a new protein
structure, using a combination of Maestro and our local programs.
- Keep thorough notes on all you done with the file (e.g.
flipped or deleted groups).
- Click on the Prep Wiz (protein preparation wizard).
- Open the desired PDB file by giving the ID or browsing the
Coloured groups are HET compounds or problems.
- Uncheck "Delete waters beyond 5 Å" and then hit Preprocess.
It takes ~1 minute.
- Review the "Problems" "Current overlapping atoms" box.
Too close H atoms will probably solve later on.
- Review the "Problems" "Alternate positions" box.
By "Next position" you can see the other position
By "commit" you select one configuration.
This part is probably better run or at least checked by
changepdb, command occ
- Go to the "Review and Modify" side and click "Analyze
Possibly delete subunits and Het groups by clicking on the and
the on Delete.
- Save a pdb file with export command.
- Go to the "Refine" side and click on "Optimize"
This may take several hours as all proton positions are
- While running, run propKa, note the pKa of the His residues,
and check if any of the Asp and Glu have a pKa above 7 or any
Lys, Arg, and Tyr have pKa below 7.
- Run changepdb, command hi to analyse the His residues.
- Run changepdb, command bc to find all buried charges.
Find all ionic pairs and note what buried charges are not
involved in ionic pairs.
- Check all these with rasmol or similar. There should be no
buried charges, not involved in ionic pairs, except close to
- With all this information, decide the protonation of all His,
Asp, Glu, Lys, and Arg residues.
- When Maestro is finished, click on "Interactive Optimizer" and
then on "Analyze Network"
Again, it may take many minutes.
- When finished, you may go through all Asn and Gln group to see
which may be flipped.
- Likewise, you should check for each His residue the
protonation state, considering 6 possible states.
With the information from changepdb, decide the actual
protonation state and check the lock box for all these.
Also enforce possible unusual protonation state of other
- Finally, click "Optimize".
This may again take hours.
- Save the pdb file.
Note that HID/E/P, ASH, GLH, LYN, and ARN need to be set by hand
- Load the file into amber with tleap, as described in the
Add protons to a ligand
- Import the pdb file (File - Import Structures)
- Mark each double / conjugated bond
- Right click and select Increase bond order (C atoms are then protonated)
- Repeat 2-3 for all bonds.
- Selet all (Ctrl-A)
- Edit - Add Hydrogens
- File - Export structure