Molecular Modeling of Chemicals Products Inhibitors of Growth Struvite Crystal

Struvite (MgNH4PO46H2O) crystals were produced by infection associated with urea generating organisms.The aim of this study is to examine the interactions between the enzyme urease and two inhibitors, the first is an inhibitor monoatomic: Aluminum and the second is a polyatomic: Citrate by the methods of molecular modeling: molecular mechanics, molecular dynamics (MM + , AMBER) and molecular docking (FleX). Supersaturated solutions induce crystallization by nucleation and subsequent crystal growth .The mechanisms for the formation of calcium phosphate urinary stones are still not understood. Chemicals product has been studied extensively as inhibitors and has been observed in the attachment of crystals to in vitro study. As a complement we have using an electron microscope Hitachi TM1000, we examined specimens of crystals struvite. The various figures show a set of grains of sizes of the order of 20 μm. The majority of these particles present regular forms. This suggests the crystal growing. This result to an alteration in the expression of these faces and the development of a characteristic architectural struvite morphology. Similar changes were observed in the presence of identical concentrations of citrate acid, and Alluminuium, emphasizing the unique interaction of phosphocitrate with the struvite crystal.


Introduction
Infection lithiasis are composed primarily of magnesium ammonium phosphate hexahydrate (MgNH 4 PO 4 6H 2 O) but may, in addition, contain calcium phosphate in the form of carbonate apatite (Ca 10 [PO 4 ]6CO 3 ) [1,2].Infection of the urinary tract with urease-producing organisms is required for the formation of infection stones in humans.Although the species Enterobacteriaceae comprises the majority of urease-producing pathogens, a variety of grampositive and gramnegative bacteria and some yeasts and mycoplasma species have the capacity to synthesize urease [3].An elevated urinary pH reduces the solubility of magnesium ammonium phosphate and favors precipitation of Struvite crystals.Higher intake of phosphate (from Proteins) and magnesium based food and lower intake of water gives rise to the PO 4 3-and Mg 2+ ions in the supersaturated urine, which leads to the conditions of formation of Struvite [4].Probably, the negative charged points of the external side of the cellular structures could reduce the metastability field of struvite and calcite, acting as heterogeneous nuclei of crystallization [5].For example, during digestion, are enzymes that accelerate the decomposition and food processing.An enzyme catalyzes only one reaction in general, because it is able to bind efficiently than a single substrate.This specificity results from the fact that when two molecules meet, their association may be stabilized by the establishment of a large number of weak bonds (hydrogen bonding, ionic bonding, Van Der Waals forces) these bonds are about 100 times weaker than a covalent bond.Modeling and simulation of struvite growth, incorporating solution chemistry and thermodynamics, growth kinetic and process description of the recovery system.An ensemble of experimental data is combined with the dynamic model to estimate struvite growth kinetics [6].
The active site in catalysis refers to the portion of the catalyst that will interact with (s) substrate (s) to form (s) product (s).The function of enzymes is related to the presence within their structure of a specific site called the active site which has the shape of a cavity or a groove.The molecules on which an enzyme acts are defined as substrates of the enzymatic reaction.Each enzyme "recognizes" one or more specific substrate molecules as a principle of complementarily-type lock and key.
An inhibitor is a compound (usually a molecule) whose action is to inhibit (that is to say slow down or stop) a chemical reaction, that is to say that he acts more or less on the speed of reaction (kinetic factor).Can be described as opposite of a catalyst.The study of the effect of inhibitor can:  To refine the catalytic mechanism of an enzymatic reaction,  To learn more about the specificity of the enzyme,  To obtain physical and chemical data about the active site.Urease (discovered in 1876 by Musculus) is an enzyme that catalyzes the conversion of urea into carbon dioxide and ammonia according to the following: Scheme 1.

Results and Discussion
1) Preparing the model Downloading urease was made from the database Brookhaven Protein Data Bank-PDB-(passcode 1UBP) is co-crystallized with the inhibitor beta-mercaptoethanol.The three dimensional structure of urease was obtained by X-ray diffraction with a resolution 1.65 Å.
Note that urease crystallized as a trimer (Figure 1), each monomer consists of 6041 atoms and more than 25 heteroatoms inhibitor.Is formed of three chains A, B and C.
2) Classification: Hydrolase Weight of the structure: 86,655.36It has simplified model of the enzyme by removing water molecules and the two chains B and C that is to say we caught a monomer (Figure 2), this chain has the amino acids forming the active site.
3) Note: When the hydrolysis is complete, a molecule of urea (i.e. 60 g) generates two molecules of ammonia (i.e.34 g). 5 kg of urea can then produce 2.83 kg of ammonia a) The hydrolysis of urea by urease liberates ammonia which alkalinized the middle and turned the red indicator from yellow to red.b) Urease is an enzyme that is present in many plants and soil.It catalyzes the hydrolysis of urea into ammonium ion and carbamate, which decomposes into carbon dioxide and ammonia.
c) The active site (Scheme 2) contains two atoms of nickel (II) are linked by a carbamate bridge, in addition, two imidazole nitrogen atoms are bound to each nickel atom and a carboxylate group and a water molecule.d) Urea is a colorless crystalline compound, the product of excretion of ammonia from the degradation of proteins in living things.Urea, also known as carbonic diamide, has the formula H 2 N-CO-NH 2 (Scheme 3, Figure 2).

Purpose
This work is to study the interactions between the enzyme urease and two inhibitors, the first is an inhibitor monoatomic: Aluminum and the second inhibitor is a polyatomic: Citrate by the methods of molecular modeling:

Preparing the Model
 Downloading urease was made from the database Brookhaven Protein Data Bank-PDB-(passcode 1UBP) is co-crystallized with the inhibitor beta-mercaptoethanol.The three dimensional structure of urease was obtained by X-ray diffraction with a resolution 1.65 Å.  Note that urease crystallized as a trimer (Figure 1), each monomer consists of 6041 atoms and more than 25 heteroatoms inhibitor is formed of three chains A, B and C Classification: Hydrolase Weight of the structure: 86,655.36. It has simplified model of the enzyme by removing water molecules and the two chains B and C that is to say we caught a monomer (Figure 2), this chain has the amino acids forming the active site.

Inhibitors used in this work are:
A: An inhibitor monoatomic: Aluminum Scheme 4 B: An inhibitor polyatomic: Citrate Scheme 5 We note that the structures of inhibitors have your donor groups as free electron pairs which will foster interactions with active site amino acids (Aluminium door and a lone pair of electron pairs Citrate is free.

Minimization of Energy
The energy minimization is done in order to optimize a structure, that is to say, to eliminate all unfavourable steric interactions (remove bad contacts) and thus we obtain a starting model close to the structure real.

Optimization of Urease
The geometry optimization of urease was carried out using the Amber force field 99 implemented in software Hyperchem 7.5.
First we use the steepest descent method with a gradient of 0.01 and a maximum of 23,655 cycles, with this method the convergence is not reached then we spent the gradient Fleetcher Reeves with the same options, even with this gradient is no has no convergence Revenue Is So we go to the gradient Polak Ribiere Optimizer and after 1093 cycles are reached convergence.

Leveraging Inhibitors
Using the force field MM + 7.5 is made on Hyperchem optimization polyatomic inhibitor (citrate).
Note: The inhibitor monoatomic (Al) does not need an optimization since it is an atom (there is no connection so there is no energy).
In Table 1 classify the different results for the optimization of Citrate.

Calculation of Charges of Polyatomic Inhibitor
For further analysis we did a load calculation for citrate to give each atom a partial charge-mail, for that and after optimization using the AM1 method to calculate Hyperchem on various charges on each atom, the results are shown in Table 2.

Molecular Dynamics
Molecular dynamics attempts to simulate the motion of atoms and molecules over time to a temperature above 0 K.It is to study the trajectory of a molecule by applying the laws of classical Newtonian mechanics and simulates the intramolecular motions that can then view real-time.

Dynamics of Urease
From the optimized enzyme as starting model we will make a two-step molecular dynamics: 1) Heating: Heated system (enzyme) with a starting temperature 0 K to 300 K for 10 ps with a 1 fs integration step for 100 steps.
2) Production or DM: At 300 K there is an equilibration, the velocities are adjusted to keep the temperature constant (there is an exchange between kinetic energy and potential energy).Then there is production of conformation.The simulation time is 80 ps of MD (Figure 3).

Molecular Dynamics Inhibitors
The molecular dynamics inhibitors are showing in the f ollowing figures: Figure 4, Figure 5 and Figure 6.plexes is based on the recognition of three-dimensional structure of a ligand by a receptor site and controls the activity of many molecules.
The docking is to determine the mode of interaction of a complex of two or more molecules, seeking directions in space and conformations favorable for ligand binding to a receptor .The next step is the positioning of the inhibitors into the active site of urease (training complex).And to realize the molecular docking software was used Hex 6.1.
Note: The aluminum does not load.

Molecular Docking 3.6. Molecular Dynamics of Complex
The docking is to predict the structure of a complex formed by two molecules, the formation of these com-There is formation of two complexes, Figures 7 and 8

Experimental
We chose the classical model for the study of phosphate crystallization because of its simplicity and satisfactory reproducibility.This model includes the study of crystallization without inhibitor and with it, in order to assess  TAI-represents the number of crystal phosphate after addition of inhibitor.We applied the previous protocol, which was to study the inhibition of phosphate by products chemicals such as citric acid, Alluminium at different pH.
Citric acid is a tricarboxylic acid of molecular weight of 192 Da.Multiple studies have shown that citrate reduces the crystallization of calcium salts in vitro [7,8].This effect is mainly due to the complication of calcium citrate, which lowers the amount of calcium available to react with ions or oxalates phosphates.This results in a decrease in urinary supersaturation of calcium oxalate.Furthermore, recent studies have shown that citrate acts as an inhibitor also the calcium oxalate, as it is capable, at least in vitro, to bind to the crystal surface, d 'reduce the size and shape change [9,10].
The isocitrate is an isomer of citrate, in small amounts excreted in urine.In a model of crystallization of calcium oxalate in vitro at 37˚C with no seeding, and Daudon M. Hennequin C. [10] found that isocitrate, 10 -3 Mol/L, is able to inhibit 42% of the crystal growth  concentration 0.25 of calcium oxalate.Such as citrate, isocitrate acts through its ability to complex the calcium ion and also by interaction with the surface of calcium crystals.This effect is proportionally greater than that observed with citrate.

Study of the Inhibitor Citrate Ion at pH = 6.5
The results of monitoring changes in the size of brushite crystals by polarized light optical microscope, at pH 6.5 and in the presence of citrate as an inhibitor, are presented as following Figure 9.
We find that inhibition of brushite is the total concentration of 6 mmol/l citrate.The FTIR spectrum confirms this inhibition in the presence of citrate solution at pH = 6.5 (Figure 10).

Study of the Inhibitor Citrate Ion at pH = 8
The results of monitoring changes in the size of struvite crystals by optical microscope, at pH 8 and in the presence of citrate as an inhibitor, are mentioned in the following graph: Figure 11.
Citrate appears to be a good inhibitor of low molecular   weight.Its effect is enhanced with increasing pH.It is more effective than acid in the urine alkaline urine and exerts an inhibitory effect on the most important calcium phosphates as calcium oxalate.Citrate ions are able to significantly reduce struvite crystallization in vitro.

Inhibition in the Presence of Aluminium Ion
We tested the inhibitory effect of aluminum (III) in the Copyright © 2011 SciRes.AMI same conditions as the previous experiments.The experiments were conducted at 37˚C in solutions with a range of concentrations ranging from 0.01 mmol/l to 0.1 mmol/l.Inhibition of growth and aggregation of brushite in the presence of aluminium is remarkable.It is significant even for lower concentrations.This reflects the strong effect of aluminium on the phases of the formation of brushite at pH 6.5, which is in perfect agreement with historical data.Indeed, the latter show that Al (III) acts much more on calcium phosphates, especially on the nucleation and crystal growth [11].
The FTIR spectrum confirm the total inhibition of the brushite at a concentration of 0.01 aluminium mMmol / l. Figure 13.

Study of the Inhibitor Al (III) at pH = 8
Monitoring the changes in the size of struvite crystals by optical microscope in the presence of alluminum at pH = 8, the following results: Figure 14.
The FTIR spectrum shows that after the addition of an aluminum concentration of 0.05 mmol/l and pH= 8, the frequency bands confirm the presence of struvite and disappearance of those characterizing the carbapatite Figure 15.
Inhibition of the crystallization of brushite is remarkable in the presence of alluminum for all concentrations, or for growth and for crystal aggregation at different pH.The results are in perfect agreement with those in the literature that mentions the formation of complexes between metal ions Al (III) and citrate ions present in the crystallization solution.The metal complexes resulting from this association strongly inhibit the growth of crys-

Conclusions
The inhibition obtained with aluminium is complete at  6.5.We obtain the complete disappearance of the brushite crystals in a solution of citrate 3 mmol/l.but, its effect on struvite is lower at pH 8.
Comparing the two complexes formed (their energy), we can say that the complex polyatomic inhibitor is more stable than the complex of the inhibitor because the interactions between the inhibitor polyatomic (Citrate) which comprises several pairs donors (electron pairs) and receiver (active site) of the enzyme is much more than the monatomic inhibitor (Al) which is one lone pair.
Using an electron microscope Hitachi TM1000, we examined specimens of a renal lithiasis.The various figures show a set of grains of sizes of the order of 20 µm.
The majority of these particles present regular forms.This suggests the crystal growing.
Frequently a certain number of differently directed crystalline grains from to each other and joined to the interface by grain boundaries.
Qualitatively, we can say that the crystals present an axis of rotation of order 2, from where the possibility that crystals belong to the monoclinic system.
Certain crystals of size 20 µm are in the center of a concretion which seems of amorphous structure on all the scales and the crucial role of which is to stop the growth of the initial crystal (Figure 16).

Figure 1 .
Figure 1.The shape of the trimer Urease.

Scheme 2
Scheme 2 Active site of Urease.

Figure 2 .
Figure 2. The simplified model of Urease.

Figure 3 .
Figure 3. Change in potential energy only 1UBP Function of time.

Figure 6 .
Figure 6.Variation of potential energy of the inhibited citrate only with time.represent the two complexes.For both complexes were found following energies: Complex 1UBP-Al: ET = −74.5 kcal/mol Al Complex 1UBP: ET = −74.5 kcal/mol

Figure 8 .
Figure 8. Complex 1UBP-Citrate.the inhibiting capacity of any chemical species used.Two solutions of following composition were mixed: A: 11.02 g/l Na 2 SO 4 •10H 2 O, 1.46 g/l MgSO 4 •7H 2 O, 4.64 g/l NH 4 Cl, 12.13 g/l KCl et 0.24 g/l Ca 2+ and B: 2.65 g/l NaH 2 PO 4 •2H 2 O, 18.82 g/l Na 2 HPO 4 •12H 2 O, 13.05 g/l NaCl, 1 g/l Na 3 C 6 H 5 O 7 .2H 2 O et 0.05 g/l 2 2 4 C O  .The solution in 2 2 4 C O  is prepared from oxalic acid 0.05 g.The precipitation of the solid phase of phosphates from artificial urine at different initial pH values (pH = 6.5; 7 and 8) was the object of our investigation.Artificial urine is prepared by mixing and stirring two equal volumes of 50 ml of solutions A and B at constant temperature (37˚C) in capped vessels to give final artificial urine at pH = 7.The pH of solution B was adjusted to required value by adding either HCl or NaOH as appropriate.Mixture agitation was maintained to prevent sedimentation.The crystal size development was monitored by polarized microscopy at different time intervals by proceeding as follows: Sample drops were examined every five minutes by polarising optical microscopy.Crystals

Figure 9 .
Figure 9. Evolution of the brushite crystal size over time and in the presence of citrate at pH 6.5.

Figure 11 .
Figure 11.Evolution of crystal size versus time in the presence of citrate, pH = 8.

4. 2 . 1 .
Study of the Inhibitor Al (III) at pH = 6.5The results of monitoring changes in the size of brushite crystals by optical microscope in the presence of aluminium as an inhibitor and pH = 6.5 are presented in following graph: Figure 12.

Figure 12 .
Figure 12.Evolution of the brushite crystal size versus time and in the presence of Al (III) at pH 6.5.

Figure 14 .
Figure 14.Evolution of the size of struvite crystals over time at pH = 8.

Copyright © 2011 Scheme 7 .
Scheme 6.The inhibitor is found in the receiver of the enzyme after complex formation 1UBP-Citrate.

Figure 16 .
Figure 16.The various figures show a set of grains of sizes of the order of 20 µm.