Efficient exploration of reactive potential energy surfaces using Car-Parrinello molecular dynamics

Refined metadynamics through canonical sampling using time-invariant bias potential: A study of polyalcohol dehydration in hot acidic solutions

We propose a canonical sampling method to refine metadynamics simulations a posteriori, where the hills obtained from metadynamics are used as a time-invariant bias potential. In this way, the statistical error in the computed reaction barriers is reduced by an efficient sampling of the collective variable space at the free energy level of interest. This simple approach could be useful particularly when two or more free energy barriers are to be compared among chemical reactions in different or competing conditions. The method was then applied to study the acid dependence of polyalcohol dehydration reactions in high-temperature aqueous solutions. It was found that the reaction proceeds consistently via an S_N 2 mechanism, whereby the free energy of protonation of the hydroxyl group created as an intermediate is affected significantly by the acidic species. Although demonstration is shown for a specific problem, the computational method suggested herein could be generally used for simulations of complex reactions in the condensed phase.

Accessing slow diffusion in solids by employing metadynamics simulation

Molecular dynamics (MD) is a powerful tool to investigate microscopic transport of atoms and molecules in condensed matter. However, there lies a large class of systems wherein atomic diffusion is too slow a process relative to the feasible time scales of typical atomistic simulations. Here, we demonstrate that with judicial implementation of a metadynamics (MTD) technique, the microscopic mechanism of atomic transport in solids can be accessed within a reasonable computational time. The calculations are carried out on the two end members of the true NASICON solid solutions, namely NaZr2(PO4)3 and Na4Zr2(SiO4)3, wherein Na+ diffusion is too slow to be accessed through standard MD simulations. The study also provides fresh insights on correlated ion hops and their implications on the effective diffusion barrier. The results are compared with climbing image nudged elastic band (CI-NEB) calculation, and available experimental results.

Precision Molecular Threading/Dethreading

The general principles guiding the design of molecular machines based on interlocked structures are well known. Nonetheless, the identification of suitable molecular components for a precise tuning of the energetic parameters that determine the mechanical link is still challenging. Indeed, what are the reasons of the "all-or-nothing" effect, which turns a molecular "speed-bump" into a stopper in pseudorotaxane-based architectures? Here we investigate the threading and dethreading processes for a representative class of molecular components, based on symmetric dibenzylammonium axles and dibenzo[24]crown-8 ether, with a joint experimental-computational strategy. From the analysis of quantitative data and an atomistic insight, we derive simple rules correlating the kinetic behaviour with the substitution pattern, and provide rational guidelines for the design of modules to be integrated in molecular switches and motors with sophisticated dynamic features.

Identification of potential drugs against SARS-CoV-2 non-structural protein 1 (nsp1)

Non-structural protein 1 (nsp1) is found in all Betacoronavirus genus, an important viral group that causes severe respiratory human diseases. This protein has significant role in pathogenesis and it is considered a probably major virulence factor. As it is absent in humans, it becomes an interesting target of study, especially when it comes to the rational search for drugs, since it increases the specificity of the target and reduces possible adverse effects that may be caused to the patient. Using approaches in silico we seek to study the behavior of nsp1 in solution to obtain its most stable conformation and find possible drugs with affinity to all of them. For this purpose, complete model of nsp1 of SARS-CoV-2 were predicted and its stability analyzed by molecular dynamics simulations in five different replicas. After main pocket validation using two control drugs and the main conformations of nsp1, molecular docking based on virtual screening were performed to identify novel potential inhibitors from DrugBank database. It has been found 16 molecules in common to all five nsp1 replica conformations. Three of them was ranked as the best compounds among them and showed better energy score than control molecules that have in vitro activity against nsp1 from SARS-CoV-2. The results pointed out here suggest new potential drugs for therapy to aid the rational drug search against COVID-19.

Communicated by Ramaswamy H. Sarma

Mechanism of Nucleotidyltransfer Reaction and Role of Mg2+ Ion in Sugar Nucleotidyltransferases

Sugar nucleotidyl transferases (SNTs) catalyze nucleotidyltransfer reactions to form sugar-nucleotides and pyrophosphate in the presence of two Mg²⁺ ions (Mg²⁺_A and Mg²⁺_B). We unveil the mechanism and free energetics of nucleotidyl transfer reaction in an SNT called GlmU through hybrid quantum mechanics-molecular mechanics molecular dynamics simulations and free energy calculations. The study identifies the roles of the active site residues and the Mg²⁺ ions in catalyzing the reaction. Of great significance, we are able to compare the free energy barrier for the reaction with that for the Mg²⁺-assisted release of the product (i.e., pyrophosphate) into the solution, shedding light on the general mechanistic and kinetic aspects of catalysis by SNTs.

A highly unsaturated six-vertex amido-substituted silicon cluster

Thermal treatment of the bicyclo[1.1.0]tetrasilatetraamide [Si4{N(SiMe3)Dipp}4] 1 resulted in the formation of a highly unsaturated six-vertex silicon cluster [Si6{N(SiMe3)Dipp}4] 2 with only four amine-substituents and two ligand-free silicon atoms. In solution, a major and a minor conformer of this cluster are in equilibrium according to multinuclear NMR spectroscopy, lineshape analysis, DFT calculations and molecular dynamics simulations. The bonding situation in the highly unsaturated cluster features lone pair type character at the ligand-free silicon atoms and partial single and double bond character in the upper butterfly-shaped ring of 2. This allows to consider 2 as the silicon analogue of a butalene isomer.

What Coordinate Best Describes the Affinity of the Hydrated Excess Proton for the Air-Water Interface?

Molecular dynamics simulations and free energy sampling are employed in this work to investigate the surface affinity of the hydrated excess proton with two definitions of the interface: The Gibbs dividing interface (GDI) and the Willard-Chandler interface (WCI). Both the multistate empirical valence bond (MS-EVB) reactive molecular dynamics method and the density functional theory-based ab initio molecular dynamics (AIMD) were used to describe the hydrated excess proton species, including "vehicular" (standard diffusion) transport and (Grotthuss) proton hopping transport and associated structures of the hydrated excess proton net positive charge defect. The excess proton is found to exhibit a similar trend and quantitative free energy behavior in terms of its surface affinity as a function of the GDI or WCI. Importantly, the definitions of the two interfaces in terms of the excess proton charge defect are highly correlated and far from independent of one another, thus undermining the argument that one interface is superior to the other when describing the proton interface affinity. Moreover, the hydrated excess proton and its solvation shell significantly influence the location and local curvature of the WCI, making it difficult to disentangle the interfacial thermodynamics of the excess proton from the influence of that species on the instantaneous surface curvature.

Deciphering the Self-Cleavage Reaction Mechanism of Hairpin Ribozyme

Hairpin ribozyme catalyzes the reversible self-cleavage of phosphodiester bonds which plays prominent roles in key biological processes involving RNAs. Despite impressive advances on ribozymatic self-cleavage, critical aspects of its molecular reaction mechanism remain controversially debated. Here, we generate and analyze the multidimensional free energy landscape that underlies the reaction using extensive QM/MM metadynamics simulations to investigate in detail the full self-cleavage mechanism. This allows us to answer several pertinent yet controversial questions concerning activation of the 2'-OH group, the mechanistic role of water molecules present in the active site, and the full reaction pathway including the structures of transition states and intermediates. Importantly, we find that a sufficiently unrestricted reaction subspace must be mapped using accelerated sampling methods in order to compute the underlying free energy landscape. It is shown that lower-dimensional sampling where the bond formation and cleavage steps are coupled does not allow the system to sufficiently explore the landscape. On the basis of a three-dimensional free energy surface spanned by flexible generalized coordinates, we find that 2'-OH is indirectly activated by adjacent G8 nucleobase in conjunction with stabilizing H-bonding involving water. This allows the proton of the 2'-OH group to directly migrate toward the 5'-leaving group via a nonbridging oxygen of the phosphodiester link. At variance with similar enzymatic processes where water wires connected to protonable side chains of the protein matrix act as transient proton shuttles, no such de/reprotonation events of water molecules are found to be involved in this ribozymatic transesterification. Overall, our results support an acid-catalyzed reaction mechanism where A38 nucleobase directly acts as an acid whereas G8, in stark contrast, participates only indirectly via stabilizing the nascent nucleophile for subsequent attack. Moreover, we conclude that self-cleavage of hairpin ribozyme follows an A_N + D_N two-step associative pathway where the rate-determining step is the cleavage of the phosphodiester bond. These results provide a major advancement in our understanding of the unique catalytic mechanism of hairpin ribozyme which will fruitfully impact on the design of synthetic ribozymes.

Supramolecular Catalysis of m-Xylene Isomerization by Cucurbiturils: Transition State Stabilization, Vibrational Coupling, and Dynamic Binding Equilibrium

The ability of cucurbit[6]uril (CB6) and cucurbit[7]uril (CB7) to catalyze the thermally activated 1,2-methyl shift isomerization pathway of m-xylene in vacuum is investigated using infrequent metadynamics. CB6 is predicted to effectively and selectively catalyze the meta-to-para isomerization through stabilization of the transition state (TS) by van der Waals push (packing coefficient ≈74%), while inhibiting the meta-to-ortho pathway by molding effects of the cavity. Interestingly, despite the snug binding, a very low rate of host-guest vibrational energy transfer is revealed using a novel approach of host-guest partition of the mode-specific anharmonic relaxation rates and ab initio molecular dynamics. The weak vibrational coupling suggests that CB can act as a thermal buffer, possibly shielding encapsulated guests from outside vibrational perturbations such as solvent effects. This dynamic effect could provide an additional boost to the reaction rate by blocking the occurrence of reaction barrier recrossing caused by the friction with surrounding molecules. Finally, mean residence times of xylene into the hosts' cavity were estimated for a range of host-guest complexes, revealing a highly dynamic equilibrium allowing very high guest turnover rates that could minimize catalyst inhibition effects commonly suffered by other supramolecular catalysts.

Density-Functional Tight-Binding Study of Carbonaceous Species Diffusion on the (100)-γ-Al2O3 Surface

Carbonaceous or oxy-carbon species are intermediates formed during C _x H _y combustion on a Pt _n /Al₂O₃ catalyst, which contain carbon, hydrogen, and oxygen atoms. The accumulation of the carbonaceous species, arguably, leads to catalytic deactivation; therefore, their removal is of importance. As the diffusion process is occasionally the rate-determining step in the growth of carbonaceous species, the present study aims to reveal the diffusion mechanisms. The free energy barriers of acetate, formate, and methoxy diffusion on the (100)-γ-Al₂O₃ surface were evaluated through extensive metadynamics simulations at the density-functional tight-binding level. The present work deduces that each adopted carbonaceous species exhibits different diffusion mechanisms and supports experimental evidence that the acetate species exhibits the slowest diffusivity among the adopted carbonaceous species.

Hydration, Solvation, and Isomerization of Methylglyoxal at the Air/Water Interface: New Mechanistic Pathways

Aqueous-phase processing of methylglyoxal (MG) has been suggested to play a key role in the formation of secondary organic aerosols and catalyze particle growth in the atmosphere. However, the details of these processes remain speculative owing to the lack of a complete description of the physicochemical behavior of MG on atmospheric aerosols. Here, the solvation and hydrolysis of MG at the air/liquid water interface is studied via classical and first-principles molecular dynamics simulations combined with free-energy methods. Our results reveal that the polarity of the water solvent catalyzed the trans-to-cis isomerization of MG at the air/liquid water interface relative to the gas phase. Despite the presence of a hydrophobic group, MG often solvates with both the ketone and methyl groups parallel to the water interface. Analysis of the instantaneous water surface reveals that when MG is in the trans state, the methyl group repels interfacial water to maintain the planarity of the molecule, indicating that lateral and temporal inhomogeneities of interfacial environments are important for fully characterizing the solvation of MG. The counterintuitive behavior of the hydrophobic group is ascribed to a tendency to maximize the number of hydrogen bonds between MG and interfacial water while minimizing the torsional free energy. This drives MG hydration, and our simulations indicate that the formation of MG diol is catalyzed at the air/liquid water interface compared to the gas phase and occurs through nucleophilic attack of water on the carbonyl carbon.

Adoption of a Turn Conformation Drives the Binding Affinity of p53 C-Terminal Domain Peptides to 14-3-3σ

The interaction between the adapter protein 14-3-3σ and transcription factor p53 is important for preserving the tumor-suppressor functions of p53 in the cell. A phosphorylated motif within the C-terminal domain (CTD) of p53 is key for binding to the amphipathic groove of 14-3-3. This motif is unique among 14-3-3 binding partners, and the precise dynamics of the interaction is not yet fully understood. Here, we investigate this interaction at the molecular level by analyzing the binding of different length p53 CTD peptides to 14-3-3σ using ITC, SPR, NMR, and MD simulations. We observed that the propensity of the p53 peptide to adopt turn-like conformation plays an important role in the binding to the 14-3-3σ protein. Our study contributes to elucidate the molecular mechanism of the 14-3-3-p53 binding and provides useful insight into how conformation properties of a ligand influence protein binding.

Ab initio metadynamics calculations of dimethylamine for probing pKb variations in bulk vs. surface environments

Free energy landscape obtained from ab initio metadynamics calculations for dimethylamine protonation at the air–water interface.

Mechanism of allosteric activation of human mRNA cap methyltransferase (RNMT) by RAM: insights from accelerated molecular dynamics simulations

The RNA guanine-N7 methyltransferase (RNMT) in complex with RNMT-activating miniprotein (RAM) catalyses the formation of a N7-methylated guanosine cap structure on the 5' end of nascent RNA polymerase II transcripts. The mRNA cap protects the primary transcript from exonucleases and recruits cap-binding complexes that mediate RNA processing, export and translation. By using microsecond standard and accelerated molecular dynamics simulations, we provide for the first time a detailed molecular mechanism of allosteric regulation of RNMT by RAM. We show that RAM selects the RNMT active site conformations that are optimal for binding of substrates (AdoMet and the cap), thus enhancing their affinity. Furthermore, our results strongly suggest the likely scenario in which the cap binding promotes the subsequent AdoMet binding, consistent with the previously suggested cooperative binding model. By employing the network community analyses, we revealed the underlying long-range allosteric networks and paths that are crucial for allosteric regulation by RAM. Our findings complement and explain previous experimental data on RNMT activity. Moreover, this study provides the most complete description of the cap and AdoMet binding poses and interactions within the enzyme's active site. This information is critical for the drug discovery efforts that consider RNMT as a promising anti-cancer target.

Radiocesium interaction with clay minerals: Theory and simulation advances Post-Fukushima

Insights at the microscopic level of the process of radiocesium adsorption and interaction with clay mineral particles have improved substantially over the past several years, triggered by pressing social issues such as management of huge amounts of waste soil accumulated after the Fukushima Dai-ichi nuclear power plant accident. In particular, computer-based molecular modeling supported by advanced hardware and algorithms has proven to be a powerful approach. Its application can now generally encompass the full complexity of clay particle adsorption sites from basal surfaces to interlayers with inserted water molecules, to edges including fresh and weathered frayed ones. On the other hand, its methodological schemes are now varied from traditional force-field molecular dynamics on large-scale realizations composed of many thousands of atoms including water molecules to first-principles methods on smaller models in rather exacting fashion. In this article, we overview new understanding enabled by simulations across methodological variations, focusing on recent insights that connect with experimental observations, namely: 1) the energy scale for cesium adsorption on the basal surface, 2) progress in understanding the structure of clay edges, which is difficult to probe experimentally, 3) cesium adsorption properties at hydrated interlayer sites, 4) the importance of the size relationship between the ionic radius of cesium and the interlayer distance at frayed edge sites, 5) the migration of cesium into deep interlayer sites, and 6) the effects of nuclear decay of radiocesium. Key experimental observations that motivate these simulation advances are also summarized. Furthermore, some directions toward future solutions of waste soil management are discussed based on the obtained microscopic insights.

Interactome of Glyceraldehyde-3-Phosphate Dehydrogenase Points to the Existence of Metabolons in Paracoccidioides lutzii

Paracoccidioides is a dimorphic fungus, the causative agent of paracoccidioidomycosis. The disease is endemic within Latin America and prevalent in Brazil. The treatment is based on azoles, sulfonamides and amphotericin B. The seeking for new treatment approaches is a real necessity for neglected infections. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is an essential glycolytic enzyme, well known for its multitude of functions within cells, therefore categorized as a moonlight protein. To our knowledge, this is the first approach performed on the Paracoccidioides genus regarding the description of PPIs having GAPDH as a target. Here, we show an overview of experimental GAPDH interactome in different phases of Paracoccidioides lutzii and an in silico analysis of 18 proteins partners. GAPDH interacted with 207 proteins in P. lutzii. Several proteins bound to GAPDH in mycelium, transition and yeast phases are common to important pathways such as glycolysis and TCA. We performed a co-immunoprecipitation assay to validate the complex formed by GAPDH with triose phosphate isomerase, enolase, isocitrate lyase and 2-methylcitrate synthase. We found GAPDH participating in complexes with proteins of specific pathways, indicating the existence of a glycolytic and a TCA metabolon in P. lutzii. GAPDH interacted with several proteins that undergoes regulation by nitrosylation. In addition, we modeled the GAPDH 3-D structure, performed molecular dynamics and molecular docking in order to identify the interacting interface between GAPDH and the interacting proteins. Despite the large number of interacting proteins, GAPDH has only four main regions of contact with interacting proteins, reflecting its ancestrality and conservation over evolution.

The multiple dissociation constants of glutathione disulfide: interpreting experimental pH-titration curves with ab initio MD simulations

The hexapeptide glutathione disulfide (GSSG) has six ionizable groups with six associated dissociation constants. The experimentally measured pH-titration curve, however, does not exhibit the six corresponding equivalence points and bears little resemblance to standard textbook examples of acid-base pH-titration curves. The curve highlights the difficulties in determining multiple pKa values of polyprotic acids - typically proteins and peptides - from experiment. The six pKa values of GSSG can, however, be estimated using Car-Parrinello molecular dynamics (CPMD) simulations in conjunction with metadynamics sampling of the underlying free energy landscape of the dissociation reactions. Ab initio MD simulations were performed on a GSSG molecule solvated by 200 water molecules. Using the estimated pKa values the theoretical titration curve was calculated and found to be in good agreement with experiment. The results clearly highlight how dissociation constants estimated from ab initio MD simulations can facilitate the interpretation of the pH-titration curves of complex chemical and biological systems.

Advances in Sustainable Catalysis: A Computational Perspective

The enormous challenge of moving our societies to a more sustainable future offers several exciting opportunities for computational chemists. The first principles approach to "catalysis by design" will enable new and much greener chemical routes to produce vital fuels and fine chemicals. This prospective outlines a wide variety of case studies to underscore how the use of theoretical techniques, from QM/MM to unrestricted DFT and periodic boundary conditions, can be applied to biocatalysis and to both homogeneous and heterogenous catalysts of all sizes and morphologies to provide invaluable insights into the reaction mechanisms they catalyze.

Structural properties of sub-nanometer metallic clusters

At the nanoscale, the investigation of structural features becomes fundamental as we can establish relationships between cluster geometries and their physicochemical properties. The peculiarity lies in the variety of shapes often unusual and far from any geometrical and crystallographic intuition clusters can assume. In this respect, we should treat and consider nanoparticles as a new form of matter. Nanoparticle structures depend on their size, chemical composition, ordering, as well as external conditions e.g. synthesis method, pressure, temperature, support. On top of that, at finite temperatures nanoparticles can fluctuate among different structures, opening new and exciting horizons for the design of optimal nanoparticles for advanced applications. This article aims to overview geometrical features of transition metal clusters and of their various rearrangements.

-60 °C solution synthesis of atomically dispersed cobalt electrocatalyst with superior performance

Temperature can govern morphologies, structures and properties of products from synthesis in solution. A reaction in solution at low temperature may result in different materials than at higher temperature due to thermodynamics and kinetics of nuclei formation. Here, we report a low-temperature solution synthesis of atomically dispersed cobalt in a catalyst with superior performance. By using a water/alcohol mixed solvent with low freezing point, liquid-phase reduction of a cobalt precursor with hydrazine hydrate is realized at -60 °C. A higher energy barrier and a sluggish nucleation rate are achieved to suppress nuclei formation; thus atomically dispersed cobalt is successfully obtained in a catalyst for oxygen reduction with electrochemical performance superior to that of a Pt/C catalyst. Furthermore, the atomically dispersed cobalt catalyst is applied in a microbial fuel cell to obtain a high maximum power density (2550 ± 60 mW m^-2) and no current drop upon operation for 820 h.

Dissociation Process of a MDM2/p53 Complex Investigated by Parallel Cascade Selection Molecular Dynamics and the Markov State Model

Recently, we efficiently generated dissociation pathways of a protein-ligand complex without applying force bias with parallel cascade selection molecular dynamics (PaCS-MD) and showed that PaCS-MD in combination with the Markov state model (MSM) yielded a binding free energy comparable to experimental values. In this work, we applied the same procedure to a complex of MDM2 protein and the transactivation domain of p53 protein (TAD-p53), the latter of which is known to be very flexible in the unbound state. Using 30 independent MD simulations in PaCS-MD, we successfully generated 25 dissociation pathways of the complex, which showed complete or partial unfolding of the helical region of TAD-p53 during the dissociation process within an average simulation time of 154.8 ± 46.4 ns. The standard binding free energy obtained in combination with one-dimensional-, three-dimensional (3D)- or C_α-MSM was in good agreement with those determined experimentally. Using 3D-MSM based on the center of mass position of TAD-p53 relative to MDM2, the dissociation rate constant was calculated, which was comparable to those measured experimentally. C_α-MSM based on all C_α coordinates of TAD-p53 reproduced the experimentally measured standard binding free energy, and dissociation and association rate constants. We conclude that the combination of PaCS-MD and MSM offers an efficient computational procedure to calculate binding free energies and kinetic rates.

Exposing catalytic versatility of GTPases: taking reaction detours in mutants of hGBP1 enzyme without additional energetic cost

Our study reveals that the replacement of catalytically competent residues by the inert amino acid alanine, S73A and E99A, in hGBP1 opens a plethora of molecularly different reaction pathways featuring very similar energy barriers as the wild type.

A novel ring-shaped reaction pathway with interconvertible intermediates in chitinase A as revealed by QM/MM simulation combined with a one-dimensional projection technique

Efficient sampling achieved by the use of a one-dimensional projection technique reveals the catalytic mechanism of chitinase A from Serratia marcescens.

Analyzing and Biasing Simulations with PLUMED

This chapter discusses how the PLUMED plugin for molecular dynamics can be used to analyze and bias molecular dynamics trajectories. The chapter begins by introducing the notion of a collective variable and by then explaining how the free energy can be computed as a function of one or more collective variables. A number of practical issues mostly around periodic boundary conditions that arise when these types of calculations are performed using PLUMED are then discussed. Later parts of the chapter discuss how PLUMED can be used to perform enhanced sampling simulations that introduce simulation biases or multiple replicas of the system and Monte Carlo exchanges between these replicas. This section is then followed by a discussion on how free-energy surfaces and associated error bars can be extracted from such simulations by using weighted histogram and block averaging techniques.

Carbon dioxide capture in 2,2'-iminodiethanol aqueous solution from ab initio molecular dynamics simulations

The reaction of carbon dioxide (CO₂) with aqueous 2,2'-iminodiethanol (trivial name is diethanolamine: DEA) has been investigated using both blue moon ensemble and metadynamics approaches combined with ab initio molecular dynamics (AIMD) simulations. A spontaneous direct proton transfer from DEA zwitterion (DEAZW) to DEA but not to H₂O has been observed in straightforward AIMD simulation in the time scale of ps. The ab initio free-energy calculations reproduced the overall free-energy difference, predicting the ionic products DEA carbamate ion (DEAC) and the protonated DEA (DEAH). The computed free-energy barrier for the first reaction step, which is the CO₂ binding (48 kJ mol^-1), is found to agree reasonably well with the available experimental data (52-56 kJ mol^-1). By contrast, the barriers for the next step, the deprotonation of zwitterion realized either via reaction with DEA or H₂O, are underestimated by 25-35 kJ mol^-1 compared to the experimental reference. A part of this error is attributed to the neglected reversible work needed to bring two reactants together, which might significantly contribute to the free-energy of activation of bimolecular reactions in a dilute solution. The computed free-energy profile is compared with our results [Y. Kubota et al., J. Chem. Phys. 146, 094303 (2017)] for the same reaction in 2-aminoethanol (trivial name is monoethanolamine: MEA).

Finding the best density functional approximation to describe interaction energies and structures of ionic liquids in molecular dynamics studies

Ionic liquids raise interesting but complicated questions for theoretical investigations due to the fact that a number of different inter-molecular interactions, e.g., hydrogen bonding, long-range Coulomb interactions, and dispersion interactions, need to be described properly. Here, we present a detailed study on the ionic liquids ethylammonium nitrate and 1-ethyl-3-methylimidazolium acetate, in which we compare different dispersion corrected density functional approximations to accurate local coupled cluster data in static calculations on ionic liquid clusters. The efficient new composite method B97-3c is tested and has been implemented in CP2K for future studies. Furthermore, tight-binding based approaches which may be used in large scale simulations are assessed. Subsequently, ab initio as well as classical molecular dynamics simulations are conducted and structural analyses are presented in order to shed light on the different short- and long-range structural patterns depending on the method and the system size considered in the simulation. Our results indicate the presence of strong hydrogen bonds in ionic liquids as well as the aggregation of alkyl side chains due to dispersion interactions.