Special Issue in Honor of Professor Nick Hadjiliadis’ RetirementView this Special Issue
Computational Modeling of the Mechanism of Urease
In order to elucidate aspects of the mechanism of the hydrolytic enzyme urease, theoretical calculations were undertaken on a model of the active site, using density functional theory. The bridging oxygen donor that has been found in the crystal structures was determined to be a hydroxide ion. The initial coordination of urea at the active site occurs most likely through the urea oxygen to the nickel ion with the lowest coordination number. This coordination can be made without much gain in energy. The calculations also showed that weak coordination of one of the urea amine nitrogen atoms to the second nickel atom is energetically feasible. Furthermore, a proposed mechanism including a tetrahedral intermediate generated by hydrolytic attack on the urea carbon by the bridging hydroxide was modeled, and the tetrahedral intermediate was found to be energetically unfavorable relative to terminal coordination of the substrate (urea).
The hydrolytic enzyme urease is responsible for the catalytic decomposition of urea to volatile ammonia and carbon dioxide . The enzyme releases ammonia and carbamate, which in turn spontaneously generate the products. The enzyme was studied early  and has generated interest for several reasons. It has been suggested to play a role in bacteria-induced ulcers [3, 4], and its activity has also been found to have implications in agriculture through the volatilization of urea—a commonly used fertilizer—that is generated by the enzyme . Urease was also the first enzyme to be found to be dependent on nickel for its function , which has made it an interesting target site for bioinorganic model chemists .
The protein structure of urease from Klebsiella aerogenes was first solved in 1995 in , and since then several other structures of the enzyme, with or without bound inhibitors have been determined [9–12], including structures from Bacillus pasteurii [10–12] and Helicobacter pylori . The active site contains two nickel ions with an interatomic distance of about 3.5 (Figure 1). The ions are bridged by a carbamylated lysine and an oxygen donor. In addition to the bridges, one of the nickel ions (Ni1) is coordinated by two histidines and a water molecule. The coordination of Ni2 is similar to the one of Ni1 and includes two histidine residues, a water molecule and a terminally bound aspartate.
A number of proposals have been made regarding possible reaction mechanisms. Consensus has been reached regarding the initial coordination of urea to the active site, which has been suggested to occur through the urea oxygen attacking the vacant coordination site on Ni1 (Figure 2), but there are divergent proposals regarding the subsequent steps. The initially suggested mechanism involves the attack on the urea carbon by a hydroxide that is terminally bound to Ni2. This leads to an intermediate that bridges the two metals and can release ammonia to form products [14, 15]. Based on the structure of an inhibitor complex, Benini et al. [12, 16] suggested a mechanism involving a secondary coordination of one of the urea nitrogen atoms to Ni2, which positions the substrate for an attack from the bridging hydroxide (water). The bridging coordination of the urea substrate was supported by Pearson et al. , but these investigators have suggested that the nucleophile is not the bridging hydroxide, but rather a water/hydroxide coordinated to Ni2.
A few computational studies on the mechanism of urease have been published during the past decade. These include molecular mechanics studies by Zimmer [18, 19] and Smyj  and DFT studies by Merz et al. [21, 22]. The latter density functional study was carried out at the B3LYP level and modeled the entire direct coordination sphere of the dinickel site with the coordinated amino acid residues truncated so that the full functional group (imidazole (His), carbamate (carbamylated Lys) and carboxylate (Asp)) were represented. This study aimed to discern between the two proposed mechanisms that involve a bridging substrate [16, 17], but it also evaluated the original suggestion that urea is terminally coordinated. For the mechanism involving a bridging intermediate, the density functional study indicated that the bridging hydroxide is the actual nucleophile, but it was found that the rate determining transition states of these mechanisms and that involving a terminally coordinated urea molecule were very similar and that it was not possible to discriminate between these mechanisms on the basis of the calculations. In more recent molecular dynamics and quantum mechanical simulations, Estiu and Merz [23, 24] have detected initial coordination modes of urea that involve either terminal coordination to Ni1 or a bridging coordination with coordination of the urea carbonyl moiety to Ni1 and simultaneous hydrogen-bonding interaction with one of the amine nitrogens with the bridging hydroxide, depending on the protonation state of surrounding amino acid residues.
In this paper, we wish to present DFT-based calculations that have been carried out in order to study the urease mechanism. The structure of the active site after initial urea coordination has been modeled. The theory involving an attack by the bridging oxygen donor has been tested. Activation data from the literature, which list energies of activation to circa 50 kJ/mol [25, 26] were compared with the energy differences between the starting structure and identified intermediates.
2. Results and Discussion
Because of the computational expense of the DFT method, a survey was done to find the smallest dependable model on which to base the rest of the study. Three model systems were investigated. The smallest system contained only the atoms within three to four bonds from the metal ions. The second contained all atoms of the coordinating amino acid residues up to and including the -carbons. The third model system also included some nonbonded residues in close vicinity. As starting structure, the native urease crystal structure from Bacillus pasteurii (PDB code 2UBP) was used . In the model survey, the structure from the protein was extracted and then geometrically minimized. A possible intermediate structure with urea bond to Ni1 was also studied. The result showed that all the structures gave very similar results, both geometrically and in the difference of energy between the two studied models. This result indicated that the smallest of the models was suitable for the study. In the chosen model, only metal-bound ligands were included. The histidine residues were modeled as imidazoles, the aspartate residue was capped as a methyl group at the -carbon while the carbamylated lysine was terminated similarly at the -carbon (Figure 1). The remaining ligands were left unchanged.
2.1. Determining the Resting Structure
Since the information provided by X-ray crystallography, especially protein crystallography, may be insufficient to accurately determine the presence and location of hydrogen atoms, different depictions of the nature of the bridging oxygen donor have been published [28–30]. As a part of preparing a suitable starting structure for the study, three different bridging ligands were considered, namely, O2- (1a), OH- (1), and H2O (1b). In order to avoid hydrogen bonding to the bridging ligand, which could by itself drastically change the energy of the structure, the aspartate ligand was turned away from the nickel center by rotating it 135° counterclockwise around the Ni–O bond. To allow comparison between the different structures throughout the study, the energy of the substrates (urea and water) and further on the intermediate compounds were added to the resulting energy of each calculated structure. In calculating the contributing energies, the energy of a “free” proton was needed. Due to the many histidines conveniently located around the active site, the energy of a proton was determined by the energy difference of a protonated histidine and a neutral histidine, which were both calculated independent of the urease model. The result from the calculations can be seen in Table 1. In the comparison of the different possible oxygen donors, the structure with the lowest energy (by a comfortable margin) was the hydroxide-bridged complex (structure 1). The energies of structure 1a and 1b were 101 and 33 kJ/mol higher than for structure 1, respectively. The identity of the bridging hydroxide ligand is in agreement with the computational results obtained by Suárez et al. . If the bridging ligand was water and the aspartate side chain was allowed to point towards the bridging ligand, a proton transfer from the water to the aspartate occurred during the geometry optimization. This observed proton transfer is an indirect confirmation of the instability of the model complex 1b. The overall structure of the energetically favoured model complex 1 is very similar to the published X-ray structure. The interatomic distance is slightly shorter (3.48 versus 3.6 in 2UBP), and the site also differs in the orientation of the rotated aspartate.
2.2. Coordination of Urea
In order to determine the mode of interaction as urea coordinates to the active site of urease, a docking experiment between the optimized model complex 1 and urea was carried out. A series of optimizations were performed moving urea from a distant position outside the active site and stepwise closer until a binding interaction was attained. Initially, urea was placed in the plane of the nickel ions and the carbamylate oxygen atoms, approximately 5.5 from the site. In order not to favor coordination to a specific nickel, the distance between the carbamylate carbon atom and the urea oxygen was chosen as the progressive distance. In the first optimization, this distance was 8 . After a number of iterations, which did not always result in the most favorable structure, the optimization was stopped, and the distance was shortened. This process was carried out repeatedly until urea was bound. The qualitative result was that urea stayed in a relatively centered position during the approach but as it came closer to the active site, it migrated towards Ni1, where it finally settled in a binding interaction with nickel. The urea coordinated through its oxygen atom in trans position to the carbamylated lysine on Ni1. Another effect of the stepwise optimization was that the water-Ni1 bond distance increased until finally water was released. The water molecule was omitted from the structure, which was reoptimized to structure 2 (Figure 3). The final urea to nickel distance is 2.12 . Compared to structure 1, the internuclear distance is approximately the same (3.47 versus 3.48 ). No major rearrangements are in other words needed in the binding of urea.
To test the viability of complex 2 as the first bound interaction between the urease active site and the substrate, some other possible structures were tested. The first one was coordination by urea trans to one of the histidines on Ni1, replacing the bound water (structure 2a). A second alternative model was achieved by replacing the water molecule on Ni2 by coordinated urea (structure 2b). The resulting energies are shown in Table 2. The two alternative structures are both disfavored by more than 30 kJ/mol relative to 2; this is a significant energy difference comparable to H‡ for the whole reaction.
2.3. Tetrahedral Intermediate
In order to assess the possibility of a tetrahedral intermediate in the catalytic reaction of urease , a number of probable intermediates along the proposed pathway were studied. The choice of intermediates was based on what could be gathered from the suggested pathway and chemical intuition. Starting from 2, the search for the most favorable configuration for complex 3, in which one of the urea nitrogen atoms coordinates to Ni2, was done using a similar approach as in the two earlier cases. The distance between the closest urea nitrogen atom and Ni2 was fixed and shortened stepwise. During the approach, the Ni2-water distance grew longer until the water molecule was finally omitted and the structure with Ni2 and one of the urea nitrogen atoms within binding distance was optimized. The optimized structure for complex 3 is shown in Figure 4. One of the nitrogen atoms has approached Ni2 to a bonding interaction. This replaces the water molecule on Ni2.
Starting from 3, a nucleophilic attack on the bridging substrate was modeled by shortening the distance between the urea carbon atom and the bridging oxygen (Figure 5). The direct approach was again determined by stepwise optimizations. During the approach, the O–H bond in the bridging hydroxyl group was lengthened, and the proton was removed in the final optimized structure (4). Complex 5 was obtained by protonating the uncoordinated urea nitrogen atom in order to facilitate the release of ammonia. This lengthened the carbon nitrogen distance, but the intermediate was kept without completely removing the ammonia molecule. In complex 6, involving the dissociation of ammonia to yield carbamate, the ammonia was completely removed (deleted) from the structure, and the structure of the resultant complex was optimized. The remaining structures were based on the resting state, 1. The differences between complexes 1, 7, and 8 are the molecules that are added to get the comparable energy. In complex 7, which involves the dissociation of carbamate, which gets protonated, and recoordination of water, the released carbamic acid and the previously released ammonia are added in order to obtain a comparable energy. The released carbamic acid is then assumed to spontaneously break down to carbon dioxide and ammonia, which are added in the optimization of complex 8.
The energies of the proposed intermediates and internuclear distances are summarized in Table 3. A graphical representation of the modelled reaction pathway is shown in Figure 5. There was a moderate increase in energy of about 40 kJ/mol before and after binding urea. Rather surprisingly, the formation of the ensuing coordination of one of the nitrogen atoms to Ni2 (structure 3) did not lead to any larger energy change. The bridging coordination of urea was only disfavored over the terminally bound state by approximately 0.9 kJ/mol. However, it should be noted that the transition state is not studied and may be significantly higher. The Ni–Ni distance is again very similar, 3.47 . The effect on the urea molecule was also rather small. The C–N distance was 1.37 for the nitrogen coordinated to nickel, while the free nitrogen carbon bond was found to be approximately 1.35 . The distances in 2 were both about 1.36 . The geometry around the coordinated nitrogen had also shifted from a nearly planar structure, which is typical for urea, to a clear sp3 hybridized configuration. The optimized Ni–N(urea) distance is rather long (2.51 ), suggesting a weak interaction. The Ni1–O(urea) bond distance was more or less unchanged.
The formation of 4 was found to require a fair amount of energy, as a large movement within the dinuclear site was needed to form the tetrahedral intermediate. The bridging oxygen donor between the two metals moved away from the nickel atoms and was positioned within bonding distance to the urea carbon. The Ni–O(phenolate) distance in structure 3 was 2.01 (Ni1) and 2.04 (Ni2) and in structure 4, 2.07 (Ni1) and 2.16 (Ni2), respectively. In 4, the partial double bond character of the C–N bonds in urea has disappeared and the hybridization of the carbon is clearly sp3. The C–N bonds are longer, 1.46 and 1.51 (see Figure 6) versus 1.37 and 1.35 in structure 3. The nickel atoms have moved away from each other, and the internuclear distance is now 3.66 . The coordination of the bridging hydroxide has pulled urea closer to the dinuclear site, and the bonded Ni–N(urea) distance is shortened from 2.51 in 3 to 2.09 in 4. The energy difference to go from structure 3 to 4 has been calculated to 111 kJ/mol. Considering that a transition state will lie higher in energy, it may be argued that the modeled reaction pathway is not likely to be a part of the mechanism of urease, which was earlier found to have an activation energy of 50 kJ/mol.
The protonation of the uncoordinated urea nitrogen atom, to start the dissociation of ammonia, relieved some of the energy in structure 4. The energy for the resultant structure 5 is about 11 kJ/mol lower than in the previous structure. The protonated nitrogen moves away from the product carbon and was in the optimized state 1.70 from the carbon, which is about 0.20 further than in structure 4. The geometry around the carbon has again started to become more planar. This has led to a strain that is relayed to the nickel ions, which have moved even further apart. The Ni–Ni distance was determined to 3.80 , which is a fairly long distance for this site. The move towards a completely planar product continued in structure 6, in which ammonia was formally removed. This again pushes the nickel ions apart to 3.92 and adds energy to the structure. The energy for structure 6 is 142 kJ/mol higher than in the resting state. This energy is relieved in complex 7 and 8 when the product is released and water is recoordinated.
As seen in Figure 5, the energy needed for the suggested pathway exceeds the empirically determined activation energy of 50 kJ/mol. Even if one assumes a large uncertainty in the calculated energies, they indicate that the tetrahedral intermediate is not a part of the urease mechanism. The activation energy for the calculated mechanism is expected to exceed 150 kJ/mol.
The binding and hydrolysis of urea at the active site of urease has been modeled. The resting state of the enzyme (active site) has been calculated. The calculations indicate that the bridging oxygen donor is a hydroxide ion; this structure is 101 and 32.9 kJ/mol more stable than the corresponding oxo- or water complexes. In agreement with previous inorganic model studies [31, 32] and proposed mechanisms [15–17, 21], the calculation further suggested that urea initially binds through its oxygen atom to Ni1 in the active site. The coordination was directed to the trans position relative to the carbamylated lysine, and the water ligand of the resting state was released.
A study of the proposed mechanism involving a tetrahedral intermediate based on the bridging hydroxy group (mechanism A, Figure 2) was carried out. In the study, four additional possible intermediates were studied. It was found that it is energetically possible to coordinate one of the urea nitrogen atoms to Ni2, but further transformation including the formation of a tetrahedral intermediate based on the bridging hydroxyl group is energetically unfavourable. Published empirical data gives activation energies of about 50 kJ/mol while the calculations indicate that 150 kJ/mol would be needed to reach the tetrahedral intermediate. Further studies involving computational modelling of the alternative mechanism involving nucleophilic attack by a terminally bound hydroxide on terminally bound urea (mechanism B, Figure 2) will be undertaken, and detection of transition states for the different mechanisms will be investigated.
The density functional calculations were performed with the Amsterdam Density Functional (ADF) program, version 2003.01 [33–37], using an uncontracted triple- STO basis set with frozen cores and an added polarization function. The implementation of the local density approximation (LDA) uses the standard Slater exchange term  and the correlation term due to Vosko, Wilk, and Nusair . Geometries were optimized at the LDA level using analytical energy gradients within a spin-restricted formalism . Total binding energies were calculated at the LDA geometries using Becke’s 1988  and Perdew’s 1986  gradient-corrected functionals for exchange and correlation, respectively. In all calculations a COSMO type solvent correction was applied by assuming a solvent dielectric constant of 4, mimicking internal protein conditions. Default convergence criteria were employed throughout. The calculations were carried out on the Lund University supercomputer facility LUNARC.
This paper is dedicated to Professor Nick. Hadjiliadis in recognition of his contributions to bioinorganic and inorganic chemistry. This work has been supported by a grant from the Swedish Research Council (VR). The authors thank Professor Ulf Ryde for helpful discussions, and are also indebted to the directors of LUNARC for the use of this supercomputing facility.
P. A. Karplus, M. A. Pearson, and R. P. Hausinger, “70 years of crystalline urease: what have we learned?” Accounts of Chemical Research, vol. 30, no. 8, pp. 330–337, 1997.View at: Google Scholar
J. B. Sumner, “The isolation and crystallization of the enzyme urease,” The Journal of Biological Chemistry, vol. 69, pp. 435–441, 1926.View at: Google Scholar
H. L. T. Mobley, M. D. Island, and R. P. Hausinger, “Molecular biology of microbial ureases,” Microbiological Reviews, vol. 59, no. 3, pp. 451–480, 1995.View at: Google Scholar
K. Stingl and H. De Reuse, “Staying alive overdosed: how does Helicobacter pylori control urease activity,” International Journal of Medical Microbiology, vol. 295, no. 5, pp. 307–315, 2005.View at: Publisher Site | Google Scholar
J. M. Bremner, “Recent research on problems in the use of urea as a nitrogen fertilizer,” Fertilizer Research, vol. 42, no. 1–3, pp. 321–329, 1995.View at: Publisher Site | Google Scholar
N. E. Dixon, C. Gazzola, R. L. Blakeley, and B. Zerner, “Jack Bean urease (EC 184.108.40.206). A metalloenzyme. A simple biological role for nickel?” Journal of the American Chemical Society, vol. 97, no. 14, pp. 4131–4133, 1975.View at: Publisher Site | Google Scholar
D. E. Fenton, “Urease—a muse for the coordination chemist,” Proceedings of the Indian National Science Academy Part A, vol. 70, no. 2, pp. 311–328, 2004.View at: Google Scholar
E. Jabri, M. B. Carr, R. P. Hausinger, and P. A. Karplus, “The crystal structure of urease from Klebsiella aerogenes,” Science, vol. 268, no. 5213, pp. 998–1004, 1995.View at: Google Scholar
M. A. Pearson, L. O. Michel, R. P. Hausinger, and P. A. Karplus, “Structures of Cys319 variants and acetohydroxamate-inhibited Klebsiella aerogenes urease,” Biochemistry, vol. 36, no. 26, pp. 8164–8172, 1997.View at: Publisher Site | Google Scholar
S. Benini, W. R. Rypniewski, K. S. Wilson, S. Ciurli, and S. Mangani, “The complex of Bacillus pasteurii urease with β-mercaptoethanol from X-Ray data at 1.65-Å resolution,” The Journal of Biological Inorganic Chemistry, vol. 3, no. 3, pp. 268–273, 1998.View at: Publisher Site | Google Scholar
S. Benini, W. R. Rypniewski, K. S. Wilson, S. Miletti, S. Ciurli, and S. Mangani, “A new proposal for urease mechanism based on the crystal structures of the native and inhibited enzyme from Bacillus pasteurii: why urea hydrolysis costs two nickels,” Structure, vol. 7, no. 2, pp. 205–216, 1999.View at: Publisher Site | Google Scholar
S. Benini, W. R. Rypniewski, K. S. Wilson, S. Ciurli, and S. Mangani, “Structure-based rationalization of urease inhibition by phosphate: novel insights into the enzyme mechanism,” The Journal of Biological Inorganic Chemistry, vol. 6, no. 8, pp. 778–790, 2001.View at: Publisher Site | Google Scholar
N.-C. Ha, S.-T. Oh, J. Y. Sung, K. A. Cha, M. H. Lee, and B.-H. Oh, “Supramolecular assembly and acid resistance of Helicobacter pylori urease,” Nature Structural Biology, vol. 8, no. 6, pp. 505–509, 2001.View at: Publisher Site | Google Scholar
N. E. Dixon, P. W. Riddles, C. Gazzola, R. L. Blakeley, and B. Zerner, “Jack bean urease (EC 220.127.116.11). V. On the mechanism of action of urease on urea, formamide, acetamide, N-methylurea, and related compounds,” Canadian Journal of Biochemistry, vol. 58, no. 12, pp. 1335–1344, 1980.View at: Google Scholar
S. J. Lippard, “At last—the crystal structure of urease,” Science, vol. 268, no. 5213, pp. 996–997, 1995.View at: Google Scholar
F. Musiani, E. Arnofi, R. Casadio, and S. Ciurli, “Structure-based computational study of the catalytic and inhibition mechanisms of urease,” The Journal of Biological Inorganic Chemistry, vol. 3, no. 8, pp. 300–314, 2001.View at: Google Scholar
M. A. Pearson, I.-S. Park, R. A. Schaller, L. O. Michel, P. A. Karplus, and R. P. Hausinger, “Kinetic and structural characterization of urease active site variants,” Biochemistry, vol. 39, no. 29, pp. 8575–8584, 2000.View at: Publisher Site | Google Scholar
M. Zimmer, “Molecular mechanics evaluation of the proposed mechanisms for the degradation of urea by urease,” Journal of Biomolecular Structure and Dynamics, vol. 17, no. 5, pp. 787–797, 2000.View at: Google Scholar
M. Zimmer, “Are classical molecular mechanics calculations still useful in bioinorganic simulations?” Coordination Chemistry Reviews, vol. 253, no. 5-6, pp. 817–826, 2009.View at: Publisher Site | Google Scholar
R. P. Smyj, “A conformational analysis study of a nickel(II) enzyme: urease,” Journal of Molecular Structure: THEOCHEM, vol. 391, no. 3, pp. 207–223, 1997.View at: Publisher Site | Google Scholar
D. Suárez, N. Díaz, and K. M. Merz Jr., “Ureases: quantum chemical calculations on cluster models,” Journal of the American Chemical Society, vol. 125, no. 50, pp. 15324–15337, 2003.View at: Publisher Site | Google Scholar
G. Estiu and K. M. Merz Jr., “Enzymatic catalysis of urea decomposition: elimination or hydrolysis?” Journal of the American Chemical Society, vol. 126, no. 38, pp. 11832–11842, 2004.View at: Publisher Site | Google Scholar
G. Estiu and K. M. Merz Jr., “Catalyzed decomposition of urea. Molecular dynamics simulations of the binding of urea to urease,” Biochemistry, vol. 45, no. 14, pp. 4429–4443, 2006.View at: Publisher Site | Google Scholar
G. Estiu and K. M. Merz Jr., “Competitive hydrolytic and elimination mechanisms in the urease catalyzed decomposition of urea,” Journal of Physical Chemistry B, vol. 111, no. 34, pp. 10263–10274, 2007.View at: Publisher Site | Google Scholar
R. E. Muck, “Urease activity in bovine feces,” Journal of Dairy Science, vol. 65, no. 11, pp. 2157–2163, 1982.View at: Google Scholar
C. C. Moyo, D. E. Kissel, and M. L. Cabrera, “Temperature effects on soil urease activity,” Soil Biology and Biochemistry, vol. 21, no. 7, pp. 935–938, 1989.View at: Google Scholar
G. Schaftenaar and J. H. Noordik, “Molden: a pre- and post-processing program for molecular and electronic structures,” Journal of Computer-Aided Molecular Design, vol. 14, no. 2, pp. 123–134, 2000.View at: Publisher Site | Google Scholar
A. M. Barrios and S. J. Lippard, “Amide hydrolysis effected by a hydroxo-bridged dinickel(II) complex: insights into the mechanism of urease,” Journal of the American Chemical Society, vol. 121, no. 50, pp. 11751–11757, 1999.View at: Publisher Site | Google Scholar
D. E. Fenton, “Structural diversity in oligonuclear nickel(II) complexes of unsymmetrical compartmental ligands,” Inorganic Chemistry Communications, vol. 5, no. 7, pp. 537–547, 2002.View at: Publisher Site | Google Scholar
B. Hommerich, H. Schwöppe, D. Volkmer, and B. Krebs, “Model complexes for ureases: a dinickel(II) complex with a novel asymmetric ligand and comparative kinetic studies on catalytically active zinc, cobalt, and nickel complexes,” Zeitschrift für anorganische und allgemeine Chemie, vol. 625, no. 1, pp. 75–82, 1999.View at: Google Scholar
H. Carlsson, M. Haukka, and E. Nordlander, “Hydrolytically active tetranuclear nickel complexes with structural resemblance to the active site of urease,” Inorganic Chemistry, vol. 41, no. 20, pp. 4981–4983, 2002.View at: Publisher Site | Google Scholar
H. Carlsson, J.-M. Haukka, A. Bousseksou, J.-M. Latour, and E. Nordlander, “Nickel complexes of carboxylate-containing polydentate ligands as models for the active site of urease,” Inorganic Chemistry, vol. 43, no. 26, pp. 8252–8262, 2004.View at: Publisher Site | Google Scholar
E. J. Baerends, D. E. Ellis, and P. Ros, “Self-consistent molecular Hartree-Fock-Slater calculations I. The computational procedure,” Chemical Physics, vol. 2, no. 1, pp. 41–51, 1973.View at: Google Scholar
G. te Velde and E. J. Baerends, “Numerical integration for polyatomic systems,” Journal of Computational Physics, vol. 99, no. 1, pp. 84–98, 1992.View at: Publisher Site | Google Scholar
C. F. Guerra, J. G. Snijders, G. te Velde, and E. J. Baerends, “Towards an order-N DFT method,” Theoretical Chemistry Accounts, vol. 99, no. 6, pp. 391–403, 1998.View at: Google Scholar
L. Versluis and T. Ziegler, “The determination of molecular structures by density functional theory. The evaluation of analytical energy gradients by numerical integration,” The Journal of Chemical Physics, vol. 88, no. 1, pp. 322–328, 1988.View at: Google Scholar
G. te Velde, F. M. Bickelhaupt, S. J. A. van Ginsbergen et al., “Chemistry with ADF,” Journal of Computational Chemistry, vol. 22, no. 9, pp. 931–967, 2001.View at: Publisher Site | Google Scholar
J. C. Slater, “Statistical exchange-correlation in the self-consistent field,” Advances in Quantum Chemistry, vol. 6, no. 1, pp. 1–92, 1972.View at: Publisher Site | Google Scholar
S. H. Vosko, L. Wilk, and M. Nusair, “Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis,” Canadian Journal of Physics, vol. 58, no. 8, pp. 1200–1211, 1980.View at: Google Scholar
A. D. Becke, “Density-functional exchange-energy approximation with correct asymptotic behavior,” Physical Review A, vol. 38, no. 6, pp. 3098–3100, 1988.View at: Publisher Site | Google Scholar
J. P. Perdew, “Density-functional approximation for the correlation energy of the inhomogeneous electron gas,” Physical Review B, vol. 33, no. 12, pp. 8822–8824, 1986.View at: Publisher Site | Google Scholar