Molecular Modeling and Synthesis of Ethyl Benzyl Carbamates as Possible Ixodicide Activity

Carbamates are molecules that have different types of biological activities and provide a particular chemical control against ticks. The new structures of the proposed compounds were optimized and synthetized respectively, through a molecular model using the methods:PM3, HF and DFT applying the B3LYP functional, with the basis 6-31+G(d) and 6-311+G(d,p), BVP86 and PBEPBE with 6-31+G(d) and the vibrational frequencies computed. These calculated frequencies were compared with the experimental ones to determine the most accurate level of theory for the prediction of vibrational frequencies of the compounds. The best results were obtained through HF/631+G(d). Additionally, we report a modification to obtain this type of compounds, and based on the amino-dehalogenation of ethyl chloroformate, different benzyl ethyl carbamates were synthesized modifying the base molecule. The performances obtained were compared to others already reported. The methodology used allowed us to synthesize new carbamates using benzylamine derivatives through a modification on the basic catalysis of the addition-elimination reaction.


Introduction
The most important tick that parasitize cattle around the world is Rhipicephalus https://doi.org/10.4236/cc.2019.71001Computational Chemistry Due to the great damages that the tick (Rhipicephalus microplus) brings to the cattle industry, millions of losses for wide distribution of the parasite and the resistance that mites have developed to all kinds of commercial ixodicide products on a national and international level [6].The research and development of new anti-mites play an important role to fight and eliminate the ectoparasite.Today two derivative compounds from carbamic acid, have been identified as very effective for such purpose: Figure 1.Nevertheless, it is necessary to determine the stability of the compounds in a long term, as well as a more accurate pharmaceutical formula, its environmental impact, toxicity, among other factors [7].Therefore, even there is a high possibility that the previous studies on carbamates could be useful for fighting the ectoparasite in a near future, it could not be used in case of being toxic or unsteady.Moreover, it is known that in average the resistance mechanism of the ticks is modified every 7 years.For those reasons, a new compound with effective results against ticks should be studied, nevertheless to synthetize new products and prove its utility seems expensive.On the other hand, such studies could be made in a theoretical way, using computational chemistry; this is a process which requires time but is addressed to the prediction of a great number of compounds that in case these turn out to be useful, they could have high probabilities of being approved for clinical use [8].
The carbamates which were developed and patented by our group, are 4-chlorophenylethyl carbamate (LQM996) and 4-bromophenylethyl carbamate (LQM919), which affected negatively to the biological parameters and reproduction of R. microplus by diminishing the oviposition, avoiding completely the viability of eggs and damaging the reproduction system of female.Likewise, in vitro experiments showed that carbamates interact very weakly with AChE and that its ixodicide activity is unrelated to such enzyme [7] [9].The most outstanding fact is the efficiency in vivo shown by carbamates, being 99.9% for LQM 996 and 98.3% for LQM919.All of this was proved in a triple resistant strain to organophosphates pyrethroids and amidines.
The mechanism of the active carbamates like ixodicides is to disable the AChE enzyme.Nevertheless, the compounds synthetized by our group, do not act upon such enzyme.This result opens up the possibility to study the potential use of the derivatives from carbamic acid upon other targets.As a particular interest, the outstanding fact that the compounds LQM996 and LQM919 showed an unexpected result in vivo, causing the death of practically all the R. microplus larvae and affecting the reproduction system of the ticks as well as their eggs consistency.As this effect was not observed for the in vitro experiments, it could be caused by the action of carbamates upon a place that affects the reproductive system.Considering that sooner or later, parasite's nature will create a resistance, the intention of our group is to synthetize new possible alternatives by the time that the first resistant signs are detected, to rote the actives against the parasite.DOI: 10.4236/cc.2019.71001 4 Computational Chemistry The aim is to validate a theoretical method, this one should be capable to describe the theoretical and experimental frequencies (IR studies) of ethyl carbamates, and to compare the experimental angles obtained through X-ray diffraction with those found in the computational methods.Also, the descriptors such as Atomic Charges (ChelpG), Molecular Volume (Mv) and electronic descriptors correlated with E HOMO and E LUMO values, such as molecular orbital energy gap (ΔE), chemical hardness (η), chemical softness (σ), chemical potential (μ), electrophilicity index (ω), electron affinity (A), ionization potential (I) and electronegativity (χ) were computed to corroborate which base is the more suitable for these compounds.The most suitable method will be used in the future for a QSAR study, to generate new molecules which could possess a high possibility to be useful and new efficient ixodicides.
A comparison between experimental and theoretical data vibrations was done.Theoretical frequencies were corrected through an empiric scale factor specific for each method [16] due to the error caused by the harmonic proximity used during the calculation [12].
Following that, theoretical and experimental bond angles and bond distances of the LQM 9006 compound were compared using crystallographic data and geometric optimizations of the same basis mentioned before.Based on these comparisons, the absolute percentage error (% E Abs ) and average error was calculated for each compound and method.The method that better describes compounds geometry and frecuencies was determined based on the minor error percentage (Validation).
In order to compare theoretical methods, single point calculations were performed with atom geometry that comes from the crystallographic data obtained by X Ray diffraction of LQM 9006 compound, using the 5 basis functions mentioned above.

LD50 Prediction Calculations
For LD50 calculations, the QSAR Toolbox software was used.All molecular structures of all compounds were modeled using the Toolbox modeler tool, and all structures were double checked comparing each SMILES code.The profiling was carried out using: OECD HPV chemical categories, Substance type and The column chromatography was carried out with Merck 938S silica gel.
Benzylamine or its derivative was dissolved in acetone and anhydrous K 2 CO 3 were added and maintained with agitation between −10˚C to 0˚C.To this mixture, a solution of ethyl chlorophormate was added.The reaction was monitored through the TLC method, after ending the reaction, the solvent was eliminated, and it was recrystallized by n-hexane/acetone 95:5 which made a pure and corresponding product as a result.

Characterization
Infrared (IR) spectra was carried out in a ThermoFisher Scientific, Nicolet iS10; with Smart Orbit accessory, diamond optical window, in solid phase with a range from 3000 to 700 cm −1 .
Proton NMR Spectrum (300 MHz) and 13 C (75 MHz) were obtained with a Varian equipment, in d-DMSO and CDCl 3 , using TMS as an internal reference.
X-Ray crystallography experiment of the compound LQM9006 was performed with an Enraf-Nonius CCD diffractometer, the radiation and wave length applied was MoKα con λ = 0.71073 Ǻ, hemisphere scanning type, 293.0 K for temperature, 2θ range for recollection of data from 8.24 to 54.90˚, 6439 collected reflections*, 3488 (Rint = 0.0200) independent reflections, the applied program for the solution of the structure is SHELXS-97 [18], refining method: Full-matrix Least-Squares on F. Mass spectrometry analysis was performed on JEOL Accu-TOF JMS-T100LC.Ionization DART+and PEG-600 was used as internal standard.

Molecular Modeling
The spectral data (vibrational frequencies) of the LQM compounds aresumma-rized in Tables 10-14: bond type that develops each vibrational band, experimental frequency (ν, cm −1 ), theoretical frequencies and their respective scaled (ν CORR, cm −1 ) as well the absolute percentage error (%EAbs) for each theory level.According to the previous information: the average error, standard deviation, relative standard deviation percentage and "L" parameter which shows the spread of individual errors regarding the average (liability level 95%), were calculated.The correction factors were used [16] [17] from the computational chemistry comparison and benchmark database.
Table 3 contains the average error calculated for each theory level, being HF/6-31+G(d) the basis which gave better results and a minor variation.
Comparison of bond angles and lengths: The LQM 9006 carbamate, showed in Figure 2, was the only compound that could be characterized through an X ray diffraction.The atomic resolution data obtained were used to compare the accuracy of theoretical optimized structures, the data are presented in Table 4 and Table 5.
The results for the quantum chemical descriptors are summarized in the Table 6.

LD50 Prediction Calculations
The LD50 prediction results are summarized in Table 7 and the statistical parameters are on Table 8.The model equation found for prediction according to the parameters previously mentioned on section 2. On Table 9 is showed the test set used to validate the QSAR prediction.19 reported compounds were selected and calculated using the equation mentioned.

Vibrational Frequencies Analysis
It is important to mention that although low levels of theory were used, the HF/6-31+G(d) represents a good commitment between the computational cost and the efficient performance of harmonic frequencies; since the electronic correlation methods do not offer significantly greater performance, the times of computational procedure are longer and they even present heavier errors [15].
For the next discussion, a 3% is arbitrary defined as the maximum objective error, the HF method has the lowest error (80.00% of computed frequencies under 3% error), followed by B3LYP/6-31+G(d), B3LYP/6-311+G(d,p) and PBEPBE/6-31+G(d) with 68.5%, 65.71% and 62.86% respectively.And the less accurate is PM3 with 48.57%.Therefore, this last method, it is not considered as a useful one and there will not be further discussion about it and presented in Tables 10-14.
On the other hand, the HF/6-31+G(d) level produced very accurate results, having only troubles for the calculation of stretching N-H bond, overestimating its value considerably.Though it happens the same with stretching C = O bond, the generated error is much lower compared with stretching N-H bond, leaving only the corresponding vibrations ν N-H.This method/basis set provided the minor error from the used sets, being this equal to 1.88% ± 0.21%.Also, HF has an RMS error of 12.81 ± 0.21 cm −1 , which is lightly superior to B3LYP/6-311++ G(d), which presents an RMS error of 12.10 ± 0.28 cm −1 .This last set also has trouble with νN-H, overestimating this even more than HF/6-31+G(d).On the other hand, although it does not present complications with ν C = O, it notably underestimates νN-C, this triggers the average error to be higher.The 65.71% of compared vibrations are minor to a 3%.This points out that although B3LYP/6-311+G(d) has a RMS error lightly minor from HF/6-31+G(d), the ab initio method results on being better for the calculation of vibrational frequencies of studied compounds, by generating values with a lower percentage of error and require less time machine in an average of 06:27:29 h less than B3LYP/6-31+G(d).
According to the observations of the B3LYP/6-311+G (d,p) level, the obtained results point out it has the same complications than on B3LYP/6-31+G(d) underestimating even more the vibrations of same bonds.Although it is a bigger basis, most of the calculated frequencies has a lower accuracy than the B3LYP/6-31+G(d).This triggers an average error to a lower efficiency.Therefore, B3LYP/6-311+G(d,p) is not the best set for the calculation of vibrational frequencies of synthetized products.The obtained errors are widely similar (B3LYP) or minor (HF) to the ones reported previously [14] [15] [16].

Length and Angle Bond Analysis
By comparing atomic bonds, a notable finding was that the lengths of all bonds X-H (where X = C, N) show error from 11% to 17%, as all methods used, overestimate the lengths of such bonds between 0.08 Å and 0.15 Å. Due to the preceding point, we decided to take all the other bonds left.
The HF method in general, predicts with a better accuracy the bonds C-C and C-N, while the DFT method using the functional B3LYP, calculates better the C-F bond and two of the C-O bonds.Even though the average error of B3LYP are minor than the HF, the differences are not meaningful [0.16% with 6-311+G(d,p) basis].Considering that the HF did not overcome the 3% error and that the computational cost regards to B3LYP/6-311+G(d,p) is much lower, it is determined that the HF/6-31+G(d) is the best theory level to predict atomic bonds of LQM 9006 compound, having an average error of 1.14%.
In regard of angles, the 5 sets provide really good results, with average error underneath 1%.The HF method predicts most of the angles with better accuracy having the lowest average error: 0.53% ± 0.13%.Thus, the HF/6-31+G(d) basis is the most suitable method to predict atomic angles of LQM 9006.
The dihedral angles are not evaluated as they lack a conformational analysis of molecule, therefore a conformation of lower energy cannot be determined in the geometric minimizations.Nevertheless, since the imaginary frequencies were not observed, is assumed that the obtained geometries correspond to a local minimum and are useful for the comparison of vibrational frequencies, angles and bond lengths [12].
With the ChelpG data obtained the absolute percentage error (% E Abs ) was calculated for each theory level method, coming afterwards standard deviation (Std.Dev.), relative standard deviation (%RSD) and Confidence Interval (CI α = 0.05) were calculated as well, Table 15.By comparing the values of E HOMO and E LUMO given by the different theory levels, the following similarity can be noticed: the basis pbepbe/6-31+G(d) and bvp86/6-31+G(d) predict the energy values with a maximum difference of 0.005 eV.The same happens between B3LYP/6-31+G(d) and B3LYP6-311+G(d,p).Thus, these similarities extend to the derived electronic descriptors calculated.The HF/6-31+G(d) method is the only basis that show a positive value in some of the E LUMO calculated therefore, the predicted data by this method is the one that varies the most.
The ionization potential (I) is a measure of the tendency of a molecule to donate electrons to an electron accepting specie.In the other hand the electron affinity (A) of a molecule is its tendency to accept electrons.The LQM molecules have an overall higher Ionization potential than Electron affinity and its tendency is: I 9010 > I 9006 > I 9005 > I 9007 > I 9011 The molecular orbital energy gap (ΔE) represents the stability of a compound.
A larger ΔE indicates that the compound is more stable, and it's highly related to the chemical hardness (η) and chemical softness (σ): a hard molecule is said to have a large ΔE value, and a soft molecule has small ΔE value.
The tendency in stability is the next: pbepbe and bvp 86: 9007 < 9005 < 9006 < 9011 < 9010 Both B3LYP: 9005 < 9007 < 9006 < 9011 < 9010 HF: 9010 < 9005 = 9006 < 9011 < 9007 Chemical potential (μ) and electronegativity (χ) are intrinsically related (χ = −μ) thus, one is often reported instead of the other.The overall tendency in ethyl carbamates regard Electronegativity is: 9011 < 9005 < 9007, 9006 < 9010 By comparing the ChelpG charges between optimized and crystallographic geometry there seems to be an overall tendency of the atoms charge whether is positive or negative, except in the atoms that forms the aromatic ring structure.None of the methods estimate with precision the charge of those C sp2 atoms, thus to simplify the results the total charge sum of the aromatic ring was taken into account for the analysis (Table 16).The results showed that both DFT/B3LYP are the best methods to describe atomic charges despite their geometric optimization are not the best.
The Molecular Volume is a measure of the molecular contact surface (Table 17).The results display an overall tendency as expected: the more atoms the molecule has, higher the molecular volume.Mv: 9005 < 9006, 9007 < 9011 < 9010.Comparing the data between Gas phase and Polarizable continuum solvation model we observe that the molecular volume given by PCM is on average higher except with the B3LYP/6-311+G(d,p) method (Table 18).
Comparing the general results given by each of the basis functions, we could say that HF is the most accurate method to predict the geometry structure of a LQM9000 molecule but is unable to predict electronic data with high accuracy.
Usually describing a chemical structure using theoretical calculations, take into consideration the quantum methods that offer a high level of theory, however the objective of carrying out theoretical calculations is to obtain optimized geometric structures so that the physicochemical parameters that are obtained are the most attached to a chemical reality (Table 19).
Although HF methods have ceased to be in force for its level of theory, to some chemical systems, carbamates in question, these methods showed better describe its geometry and validating them with critical parameters such as bond lengths, bond angles and frequencies of molecular vibration.For these chemical systems, carbamates, HF method showed a better description.Although the differences with the DFT methods are not large, we always take into account a method which describes the aforementioned parameters with better accuracy.

LD50 Prediction Calculations
LD50 prediction showed that LQM9005, LQM9006, LQM9007 and LQM9011 compounds are relatively safe.This is because of the high levels needed to induce lethality.However, the LQM9010 showed to be more unsafe than the other compounds.These results are supported with the statistical parameters needed to validate a QSAR model, such as an R 2 higher than 0.7 and the difference between R 2 and Q 2 shouldn't be higher than 0.3.Also, the number of compounds used to generate the model are 78 points.And the model showed to be robust enough according with the validation parameter Q 2 .The test set corroborates how robust is the model.In this case there is an R 2 = 0.9647, showing how does the predictive level of the equation correlates with a high level of confidence.

Synthesis
The Aliphatic carbamates synthesis has been made during the last years  19.
The methodology described in this work had not previously been used and has advantages over those previously reported, specifically focused on the synthesis of ethylbenzyl carbamate and very similar compounds, which consisted in using a stronger basis (K 2 CO 3 ), with no need to increase temperature, making the process under ambient conditions and the yields obtained are excellent from 86% to 96%.The reason why a stronger basis than NaHCO 3 is required, is because of the difference in constants of acidity of aromatic and aliphatic amines: CO − = 10.2. 39Even though the values are referred to the water environment, it is useful to compare the relative basicity of each kind of compounds, being, in case of the amines shown with a deeply great difference, as methylamine is 105.97 times more basic than aniline.It can also be observed because of the HCO 3− ion is an efficient base by using it in the aromatic carbamate synthesis, being more basic than the aromatic amines.Besides it is being considerably less basic than aliphatic amine, is not useful to function as a base (on the contrary, it acts as an acid) with the mentioned substrates, requiring the IR radiation appliance to be used.On the contrary, CO 3 2− ion has a relative basicity very like aliphatic amines, therefore they can be neutralized, nevertheless, to benefit its efficiency it is necessary to use a higher quantity according to the stoichiometric balance (Equation ( 2)).
The previous equation corresponds to the neutralization of an ammoniac salt by the potassium carbonate.Even though the basicity that the amine and the carbonate present is very similar, according to Le Châtelier principle, it benefited the reaction of neutralization in our favor, by increasing the stoichiometric quantity of potassium carbonate.Another fact that benefited the reaction, was the gradual consumption of amine as the synthesis of carbamate passed.Five compounds were chosen for its synthesis, modifying the structure with the increase of one carbon atom such as benzylamine and its derivatives, with electron donors and electron acceptors, which could allow the effect evaluation of this substituent upon its properties, considering that the synthetized structures are hydrophobic and none of them represent a meaningful steric effect.

Conclusion
The vibrational frequencies of synthesized compounds were calculated, using the PM3 semi empirical method, the method ab initio HF with 6-31+G(d) function, the DFT method, using the hybrid functional B3LYP with 6-31+G(d) and 6-311+G(d,p) functions, the basis pbepbe/6-31+G(d) and bvp86/6-31+G(d).The calculated frequencies were scaled with the appropriate factors and were compared with the ones observed experimentally.The level of theory HF/6-31+G(d) was calculated with the intention of providing closer values to the experimental ones (%E Abs = 1.88% ± 0.26%, RM Sov = 32.88 ± 0.23 cm −1 ) and requiring less time for the computational procedure, which is the best prediction of vibrational frequencies of the synthesized compounds.Additionally, lengths and theoretical and experimental angles of the bonds were compared from the LQM 9006 model [(4-fluorobenzyl) ethyl carbamate]; once again, the ab initio method was the most effective for the calculation of such properties, even it did not provide the results of minor error (for bond lengths), the values caused by the other methods were not very different, so therefore by considering the computational time HF is the most suitable method.Synthesis and relative safety of new 5 different ethyl benzyl carbamates was carried out, and it was clear that the synthetic methodology applied was better; in addition the reaction is performed at room temperature with excellent reaction yields.Also the computed endpoint (LD50) shows us that the proposed compounds are safe enough due to the high doses needed to trigger an acute toxicity response.

Figure 2 .
Figure 2. ORTEP representation of the LQM 9006 compound obtained through X ray diffraction.The labels shown are the same used in the result charts.

[ 37 ]
-[63].In these examples different changes are made, changing the nature of the base, solvent or temperature; giving a different efficiency and compares the different variations upon the analyzed reaction, as it is shown in some examples are described in Table

Table 3 .
Average error and Standard Deviation of theory levels used.

Table 6 .
Descriptors calculated for the LQM studied.

Table 8 .
Statistic of the prediction model.

Table 9 .
Test set used to validate the QSAR model.

Table 16 .
ChelpG charged computed for the LQM series.