Introducing KICK-MEP: exploring potential energy surfaces in systems with significant non-covalent interactions.

CONTEXT: Exploring potential energy surfaces (PES) is fundamental in computational chemistry, as it provides insights into the relationship between molecular energy, geometry, and chemical reactivity. We introduce Kick-MEP, a hybrid method for exploring the PES of atomic and molecular clusters, particularly those dominated by non-covalent interactions. Kick-MEP computes the Coulomb integral between the maximum and minimum electrostatic potential values on a 0.001 a.u. electron density isosurface for two interacting fragments. This approach efficiently estimates interaction energies and selects low-energy configurations at reduced computational cost. Kick-MEP was evaluated on silicon-lithium clusters, water clusters, and thymol encapsulated within Cucurbit[7]uril, consistently identifying the lowest energy structures, including global minima and relevant local minima. METHODS: Kick-MEP generates an initial population of molecular structures using the stochastic Kick algorithm, which combines two molecular fragments (A and B). The molecular electrostatic potential (MEP) values on a 0.001 a.u. electron density isosurface for each fragment are used to compute the Coulomb integral between them. Structures with the lowest Coulomb integral are selected and refined through gradient-based optimization and DFT calculations at the PBE0-D3/Def2-TZVP level. Molecular docking simulations for the thymol-Cucurbit[7]uril complex using AutoDock Vina were performed for benchmarking. Kick-MEP was validated across different molecular systems, demonstrating its effectiveness in identifying the lowest energy structures, including global minima and relevant local minima, while maintaining a low computational cost.

A mechanistic insight for the biosynthesis of N,N-dimethyltryptamine: An ONIOM theoretical approach.

Psychoactive natural products are potent serotonergic agonists capable of modulating brain functions such as memory and cognition. These substances have shown therapeutic potential for treating various mental disorders. The fact that N,N-dimethyltryptamine (DMT) is produced endogenously in several plants and animals, including humans, makes it particularly attractive. As an amino acid-derived alkaloid, the DMT biosynthetic pathway is part of the L-tryptophan biochemical cascade and can be divided into the decarboxylation by an aromatic L-amino acid decarboxylase (AADC) for tryptamine formation and the subsequent double-methylation by the indolethylamine-N-methyltransferase (INMT) through the cofactor S-adenosyl-L-methionine (SAM), a methyl donor. Unlike the decarboxylation mechanism of L-tryptophan, the molecular details of the double methylation of tryptamine have not been elucidated. Therefore, we propose an in silico model using molecular dynamics (MD), non-covalent interaction index (NCI) and density functional theory (DFT) calculations with the ONIOM QM:MM B3LYP/6-31+G(d,p):MM/UFF level of theory. Based on the obtained energetic data, the potential energy surface (PES) indicates an SN2 mechanism profile, with the second methylation energy barrier being the rate-limiting step with δG‡=60kJ∙mol-1 larger than the previous methylation, following the NCI analysis showing more repulsive interactions for the second transition state. In addition, the hybridization information of each reaction step provides geometric details about the double-methylation.

An approach based on genetic algorithms and machine learning coupled for studying alloy and molecular clusters by optimizing quantum energy surfaces.

A new genetic algorithm has been proposed focusing on direct ab initio potential energy surface (PES) global minima search. Besides the commonly used operators, this new approach uses an operator to: improve the initial cluster generation, classify and compare all generated clusters, and use machine learning to model the quantum PES used in parallel optimization. Part of the validation process for this methodology was done with C u n A u m ( n + m ≤ X for X = 14 , 19 , 38 , 55 ) and A u n A g n ( n = 10 , 20 , 30 , 40 , 50 , 60 , 70 , and 75). The results are in fair agreement with the literature and led to a new global minimum for C u 12 A u 7 . A search has been done for the lowest energies of L i n nanoclusters with 2-8 atoms using the DFT approach and for L i 3 , L i 4 , L i 2 H , L i 3 H using DLPNO-CCSD(T) approach. NQGA successfully performed the MP2 optimizations for ( H 2 O ) 11 cluster. In all cases, the proposed genetic algorithm located the previously reported global minima with very efficient performance. The new proposed methodology makes it possible to optimize cluster geometries directly using high-level ab initio methods relinquishing any bias introduced by a classical approach. Our results show that this proposed method has great potential applications due to its flexibility and efficiency in identifying global minima in the tested atomic systems.

Long-range parameter optimization for a better description of potential energy surfaces using Density Functional Theory.

The advance of computing and the development of modern quantum chemistry models such as Density Functional Theory (DFT) have allowed scientists to perform fast in silico studies with accurate results. It also allowed for the achievement of empirically unattainable quantities such as Potential Energy Surfaces (PES), a fundamental construct in various applications, such as the study of weakly bound systems. One of DFT's current weaknesses is a reliable description of PESs, due to a lack of suitable exchange-correlation functionals. In general, other post-Hartree-Fock methods are employed, such as n th-order Møller-Plesset's Perturbation Theory (MPn) or Coupled Cluster Theory (CCSD(T)) with large basis sets. Despite producing good results, these methods demand much computational power when applied to large systems. This work presents a novel approach of PES description of the H2O2-Kr system using DFT by optimizing a long-range parameter present in some DFT functionals, obtaining results similar to those of the MPn methods with somewhat less computational time necessary. Graphical Abstract By optimizing the omega value of certain functionals, DFT can describe PESs with accuracy comparable to MP4 methods.

Planar Elongated B12 Structure in M3B12 Clusters (M = Cu-Au).

Here, it is shown that the M3B12 (M = Cu-Au) clusters' global minima consist of an elongated planar B12 fragment connected by an in-plane linear M3 fragment. This result is striking since this B12 planar structure is not favored in the bare cluster, nor when one or two metals are added. The minimum energy structures were revealed by screening the potential energy surface using genetic algorithms and density functional theory calculations. Chemical bonding analysis shows that the strong electrostatic interactions with the metal compensate for the high energy spent in the M3 and B12 fragment distortion. Furthermore, metals participate in the delocalized π-bonds, which infers an aromatic character to these species.

Spherical harmonics representation of the potential energy surface for the H2⋯H2 van der Waals complex.

We perform a study of the molecular anisotropy for the H2⋯H2 van der Waals system using a spherical harmonics expansion. We use six leading stable configurations to construct our analytical potential energy surface (PES) from ab initio calculations guided qualitatively by the symmetry-adapted perturbation theory (SAPT) analyses. We extrapolate the energies of the PES performed at the CCSD(T)/aug-cc-pVnZ (n = 2 and 3) levels to the complete basis set (CBS) limit. To best fit the shallow potential energy surface of each leading configuration with the intermolecular distance, it was employed an extended version of the Rydberg potential. To assess the quality of our extrapolated analytical PES, we calculate the second virial coefficients, which are in relatively good agreement with the experimental data. As a result, the spherical harmonics coefficients obtained might be of considerable relevance in spectroscopy and dynamics applications.

Conformational and electronic study of dopamine interacting with the D2 dopamine receptor.

We report an exhaustive conformational and electronic study on dopamine (DA) interacting with the D2 dopamine receptor (D2 DR). For the first time, the complete surface of the conformational potential energy of the complex DA/D2 DR is reported. Such a surface was obtained through the use of QM/MM calculations. A detailed study of the molecular interactions that stabilize and destabilize the different molecular complexes was carried out using two techniques: Quantum Theory of Atoms in Molecules computations and nuclear magnetic shielding constants calculations. A comparative study of the behavior of DA in the gas phase, aqueous solution, and in the active site of D2 DR has allowed us to evaluate the degree of deformation suffered by the ligand and, therefore, analyze how rustic are the lock-key model and the induced fit theory in this case. Our results allow us to propose one of the conformations obtained as the "biologically relevant" conformation of DA when it is interacting with the D2 DR.

A Fukui function-guided genetic algorithm. Assessment on structural prediction of Sin (n = 12-20) clusters.

Theoretical studies are essential for the structural characterization of clusters, when it comes to rationalize their unique size-dependent properties and composition. However, the rapid growth of local minima on the potential energy surface (PES), with respect to cluster size, makes the candidate identification a challenging undertaking. In this article, we introduce a hybrid strategy to explore the PES of clusters. This proposal involves the use of a biased initial population of a genetic algorithm procedure. Each individual in this population is built by assembling small fragments, according to the best matching of the Fukui function. The performance of a genetic algorithm procedure. The performance of the method is assessed on the PES exploration of medium-sized Sin clusters (n = 12-20). The most relevant results are: (a) the method converges at almost half of the time used by the canonical version of the GA and, (b) in all the studied cases, with the exception of Si13 and Si16 , the method allowed to identify the global minimum (GM) and other important low-lying structures. Additionally, the apparent deficiency of the proposal to identify the GM was corrected when a Si atom, or other low-lying isomers, were considered to build the clusters. © 2017 Wiley Periodicals, Inc.

Inversion of the stereochemical configuration (3S, 5S)-clavaminic acid into (3R, 5R)-clavulanic acid: A computationally-assisted approach based on experimental evidence.

Clavulanic acid (CA), a potent inhibitor of ß-lactamase enzymes, is produced by Streptomyces clavuligerus (Sc) cultivation processes, for which low yields are commonly obtained. Improved knowledge of the clavam biosynthetic pathway, especially the steps involved in the inversion of 3S-5S into 3R-5R stereochemical configuration, would help to eventually identify bottlenecks in the pathway. In this work, we studied the role of acetate in CA biosynthesis by a combined continuous culture and computational simulation approach. From this we derived a new model for the synthesis of N-acetyl-glycyl-clavaminic acid (NAG-clavam) by Sc. Acetylated compounds, such as NAG-clavam and N-acetyl-clavaminic acid, have been reported in the clavam pathway. Although the acetyl group is present in the ß-lactam intermediate NAG-clavam, it is unknown how this group is incorporated. Hence, under the consideration of the experimentally proven accumulation of acetate during CA biosynthesis, and the fact that an acetyl group is present in the NAG-clavam structure, a computational evaluation of the tentative formation of NAG-clavam was performed for the purpose of providing further understanding. The proposed reaction mechanism consists of two steps: first, acetate reacts with ATP to produce a reactive acylphosphate intermediate; second, a direct nucleophilic attack of the terminal amino group of N-glycyl-clavaminic on the carbonyl carbon of the acylphosphate intermediate leads to a tetrahydral intermediate, which collapses and produces ADP and N-acetyl-glycyl-clavaminic acid. The calculations suggest that for the proposed reaction mechanism, the reaction proceeds until completion of the first step, without the direct action of an enzyme, where acetate and ATP are involved. For this step, the computed activation energy was ≅2.82kcal/mol while the reaction energy was ≅2.38kcal/mol. As this is an endothermic chemical process with a relatively small activation energy, the reaction rate should be considerably high. The calculations offered in this work should not be considered as a definite characterization of the potential energy surface for the reaction between acetate and ATP, but rather as a first approximation that provides valuable insight about the reaction mechanism. Finally, a complete route for the inversion of the stereochemical configuration from (3S, 5S)-clavaminic acid into (3R, 5R)-clavulanic acid is proposed, including a novel alternative for the double epimerization using proline racemase and NAG-clavam formation.