Category Archives: Crystallography
Quick Post on preparing Gaussian input files from PDB files.
If you’re modeling biological systems chances are that, more often than not, you start by retrieving a PDB file. The Protein Data Bank is a repository for all things biochemistry – from oligo-peptides to full DNA sequences with over 140,000 available files encoding the corresponding structure obtained by various experimental means ranging from X-Ray diffraction, NMR and more recently, Cryo Electron Microscopy (CEM).
The PDB file encodes the Cartesian coordinates for each atom present in the structure as well as their in the same way molecular dynamics codes -like AMBER or GROMACS- code the parameters for a force field; this makes the PDB a natural input file for MD.
There are however some considerations to have in mind for when you need to use these coordinates in electronic structure calculations. Personally I give it a pass with OpenBabel to add (or possibly just re-add) all Hydrogen atoms with the following instruction:
$>obabel -ipdb filename.pdb -ogjf -Ofilename.gjf -h
Alternatively, you can select a pH value, say 7.5 with:
$>obabel -ipdb filename.pdb -ogjf -Ofilename.gjf -h -p7.5
You may also use the GUI if by any chance you’re working in Windows:
This sends all H atoms to the end of the atoms list. Usually for us the next step is to optimize their positions with a partial optimization at a low level of theory for which you need to use the ReadOptimize ReadOpt or RdOpt in the route section and then add the atom list at the end of the input file:
Finally, visual inspection of your input structure is always helpful to find any meaningful errors, remember that PDB files come from experimental measurements which are not free of problems.
As usual thanks for reading, commenting, and sharing.
So many things have happened since I last updated this blog but I will come to write on them when appropriate. Right now I’d like to share an invitation by Prof. Ponnadurai Ramasami from the University of Mauritius to the upcoming Virtual Conference on Computational Chemistry from the 1st to the 31st of August. Deadlines can be consulted here and the most important is the abstract submission on June 30th. This conference is part of the official celebrations of the International Year of Crystallography so talks involving experimental determination of electron densities will be well suited.
I participated in the latest edition and I must say it was a very enriching opportunity to learn from so many other researchers from across the world without leaving my desk. I already know what my talk will be about, now that we are so close to finish and submit a paper on the absence of reactivity for an anti-aromatic set of molecules. (I think I’ll call it “The reactivity of molecules that never existed [but that maybe should have.].) All talks are sent either as pdf, powerpoint presentations, youtube videos, etc. and Q&A are done over e-mail.
So this is a calling to other computational chemists out there who want to participate in this virtual conference. Kudos to Prof. Ponnadurai Ramasami and lets hope we can crystallize his visit to Mexico during 2014 the International Year of Crystallography (pun intended) and here’s to me going sometime to Mauritius!
What a happy coincidence -if indeed it was- that #RealTimeChem week happened to coincide with the sixtieth anniversary of the three seminal papers published in Nature on this day back in 1953, one of which was co-authored by J. Watson and F. Crick; of course I mean the publication for the first time of the structure of deoxyribose nucleic acid, or DNA, as we now call it.
You can get the original Nature papers from 1953 here at: http://www.nature.com/nature/dna50/archive.html (costs may apply)
Molecular Structure of Nucleic Acids: A Structure for Deoxyribose Nucleic Acid 737
J. D. WATSON & F. H. C. CRICK
Molecular Structure of Nucleic Acids: Molecular Structure of Deoxypentose Nucleic Acids 738
M. H. F. WILKINS, A. R. STOKES & H. R. WILSON
Molecular Configuration in Sodium Thymonucleate 740
ROSALIND E. FRANKLIN & R. G. GOSLING
Nature’s podcast released two episodes (called ‘pastcast’) to celebrate DNA’s structure’s birthday, one of them is an interview with Dr. Raymond Gosling who in 1953 worked under Dr. Rosalind Franklin at King’s College London in diffractometry of biological molecules. If you haven’t listened to them you can get them here at nature.com/nature/podcasts. Of course, the history around the discovery of DNA’s structure is not without controversy and it has been long argued that the work of Franklin and Gossling didn’t get all deserved credit from Watson and Crick. In their paper W&C acknowledge the contribution of the general nature of DNA from the unpublished results by Franklin’s laboratory but that is as far as they went, they didn’t even mention photo 51 which Crick saw at Wilkins laboratory, who in turn got it from Gossling at Franklin’s suggestion. Still, no one can deny that the helical structure with which we are now familiar is their work, and more importantly the discovery of the specific pairing, which according to Gossling was a stroke of genious that probably couldn’t have happened in his own group, but without Franklin’s diffraction and Gossling’s crystallization there was little they could do. Details about the process used to crystallize DNA can be heard in the aforementioned podcast, along with an inspiring tale of hard work by Dr. Gossling. Go now and listen to it, its truly inspiring.
For me it was not the story of a helix, that I was familiar with; it was the story of the specific pairing of two hélices
– Dr. Raymond Gosling
Above, the iconic Photo 51 taken by Franklin and Gossling (have you ever noticed how most scientists refer to Franklin just as Rosalind but no one refers to Watson as James? Gender bias has a role in this tale too) To a trained crystallographer, the helical symmetry is evident from the diffraction pattern but going from Photo 51 to the representation below was the subject of hard work too.
There are million of pages written during the last 60 years about DNA’s structure and its role in the chemistry of life; the nature of the pairing and the selectivity of base pairs through hydrogen bond interactions, an interaction found ubiquitously in nature; water itself is a liquid due to the intermolecular hydrogen-bonds, which reminds us about the delicate balance of forces in biochemistry making life a delicate matter. But I digress. Millions of pages have been written and I’m no position of adding a meaningful sentence to them; however, it is a fascinating tale that has shaped the course of mankind, just think of the Human Genome Project and all the possibilities both positive and negative! DNA and its discovery tale will continue to amaze us and inspire us, just like in 2011 it inspired the Genetech company to set a Guiness World Record with the largest human DNA helix.
Happy birthday, DNA!