Triple hydrogen bonding in a circular arrangement: ab initio, DFT and first-principles MD studies of tris-hydroxyaryl enamines.
Ontology highlight
ABSTRACT:
SUBMITTER: Martyniak A
PROVIDER: S-EPMC3474916 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Recently, AIMD (ab initio molecular dynamics) has been extensively employed to explore the dynamical information of electronic systems. However, it remains extremely challenging to reliably predict the properties of nanosystems with a radical nature using conventional electronic structure methods (e.g., Kohn-Sham density functional theory) due to the presence of static correlation. To address this challenge, we combine the recently formulated TAO-DFT (thermally-assisted-occupation density functional theory) with AIMD. The resulting TAO-AIMD method is employed to investigate the instantaneous/average radical nature and infrared spectra of n-acenes containing n linearly fused benzene rings (n = 2-8) at 300 K. According to the TAO-AIMD simulations, on average, the smaller n-acenes (up to n = 5) possess a nonradical nature, and the larger n-acenes (n = 6-8) possess an increasing radical nature, showing remarkable similarities to the ground-state counterparts at 0 K. Besides, the infrared spectra of n-acenes obtained with the TAO-AIMD simulations are in qualitative agreement with the existing experimental data.
| S-EPMC7674663 | biostudies-literature
Project description:MotivationCysteine-rich proteins cover many important families in nature but there are currently no methods specifically designed for modeling the structure of these proteins. The accuracy of disulfide connectivity pattern prediction, particularly for the proteins of higher-order connections, e.g., >3 bonds, is too low to effectively assist structure assembly simulations.ResultsWe propose a new hierarchical order reduction protocol called Cyscon for disulfide-bonding prediction. The most confident disulfide bonds are first identified and bonding prediction is then focused on the remaining cysteine residues based on SVR training. Compared with purely machine learning-based approaches, Cyscon improved the average accuracy of connectivity pattern prediction by 21.9%. For proteins with more than 5 disulfide bonds, Cyscon improved the accuracy by 585% on the benchmark set of PDBCYS. When applied to 158 non-redundant cysteine-rich proteins, Cyscon predictions helped increase (or decrease) the TM-score (or RMSD) of the ab initio QUARK modeling by 12.1% (or 14.4%). This result demonstrates a new avenue to improve the ab initio structure modeling for cysteine-rich proteins.Availability and implementationhttp://www.csbio.sjtu.edu.cn/bioinf/Cyscon/Contactzhng@umich.edu or hbshen@sjtu.edu.cn.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC5898604 | biostudies-literature
Project description:A fundamental requirement to predict the native conformation, address questions of sequence design and optimization, and gain insights into the folding mechanisms of proteins lies in the definition of an unbiased reaction coordinate that reports on the folding state without the need to compare it to reference values, which might be unavailable for new (designed) sequences. Here, we introduce such a reaction coordinate, which does not depend on previous structural knowledge of the native state but relies solely on the energy partition within the protein: the spectral gap of the pair nonbonded energy matrix (ENergy Gap, ENG). This quantity can be simply calculated along unbiased MD trajectories. We show that upon folding the gap increases significantly, while its fluctuations are reduced to a minimum. This is consistently observed for a diverse set of systems and trajectories. Our approach allows one to promptly identify residues that belong to the folding core as well as residues involved in non-native contacts that need to be disrupted to guide polypeptides to the folded state. The energy gap and fluctuations criteria are then used to develop an automatic detection system which allows us to extract and analyze folding transitions from a generic MD trajectory. We speculate that our method can be used to detect conformational ensembles in dynamic and intrinsically disordered proteins, revealing potential preorganization for binding.
| S-EPMC8009504 | biostudies-literature
Project description:Unlike typical hydrogen-bonded networks such as water, hydrogen bonded ionic liquids display some unusual characteristics due to the complex interplay of electrostatics, polarization, and dispersion forces in the bulk. Protic ionic liquids in particular contain close-to traditional linear hydrogen bonds that define their physicochemical properties. This work investigates whether hydrogen bonded ionic liquids (HBILs) can be differentiated from aprotic ionic liquids with no linear hydrogen bonds using state-of-the-art ab initio calculations. This is achieved through geometry optimizations of a series of single ion pairs of HBILs in the gas phase and an implicit solvent. Using benchmark CCSD(T)/CBS calculations, the electrostatic and dispersion components of the interaction energy of these systems are compared with those of aprotic ionic liquids. The inclusion of the implicit solvent significantly influenced geometries of single ion pairs, with the gas phase shortening the hydrogen bond to reduce electrostatic interactions. HBILs were found to have stronger interactions by at least 10EtMeNH0 kJ mol-1 over aprotic ILs, clearly highlighting the electrostatic nature of hydrogen bonding. Geometric and energetic parameters were found to complement each other in determining the extent of hydrogen bonding present in these ionic liquids.
| S-EPMC6468050 | biostudies-literature
Project description:We introduce a machine learning-based approach called ab initio generalized Langevin equation (AIGLE) to model the dynamics of slow collective variables (CVs) in materials and molecules. In this scheme, the parameters are learned from atomistic simulations based on ab initio quantum mechanical models. Force field, memory kernel, and noise generator are constructed in the context of the Mori-Zwanzig formalism, under the constraint of the fluctuation-dissipation theorem. Combined with deep potential molecular dynamics and electronic density functional theory, this approach opens the way to multiscale modeling in a variety of situations. Here, we demonstrate this capability with a study of two mesoscale processes in crystalline lead titanate, namely the field-driven dynamics of a planar ferroelectric domain wall, and the dynamics of an extensive lattice of coarse-grained electric dipoles. In the first case, AIGLE extends the reach of ab initio simulations to a regime of noise-driven motions not accessible to molecular dynamics. In the second case, AIGLE deals with an extensive set of CVs by adopting a local approximation for the memory kernel and retaining only short-range noise correlations. The scheme is computationally more efficient than molecular dynamics by several orders of magnitude and mimics the microscopic dynamics at low frequencies where it reproduces accurately the dominant far-infrared absorption frequency.
| S-EPMC10998567 | biostudies-literature
Project description:We present an approach for the calculation of spin density distributions for molecules that require very large active spaces for a qualitatively correct description of their electronic structure. Our approach is based on the density-matrix renormalization group (DMRG) algorithm to calculate the spin density matrix elements as a basic quantity for the spatially resolved spin density distribution. The spin density matrix elements are directly determined from the second-quantized elementary operators optimized by the DMRG algorithm. As an analytic convergence criterion for the spin density distribution, we employ our recently developed sampling-reconstruction scheme [J. Chem. Phys.2011, 134, 224101] to build an accurate complete-active-space configuration-interaction (CASCI) wave function from the optimized matrix product states. The spin density matrix elements can then also be determined as an expectation value employing the reconstructed wave function expansion. Furthermore, the explicit reconstruction of a CASCI-type wave function provides insight into chemically interesting features of the molecule under study such as the distribution of α and β electrons in terms of Slater determinants, CI coefficients, and natural orbitals. The methodology is applied to an iron nitrosyl complex which we have identified as a challenging system for standard approaches [J. Chem. Theory Comput.2011, 7, 2740].
| S-EPMC3373140 | biostudies-literature
Project description:Here we present a comprehensive density functional theory (DFT) based ab initio study of copper bismuth oxide CuBi2O4 (CBO) in combination with experimental observations. The CBO samples were prepared following both solid-state reaction (SCBO) and hydrothermal (HCBO) methods. The P4/ncc phase purity of the as-synthesized samples was corroborated by Rietveld refinement of the powdered X-ray diffraction measurements along with Generalized Gradient Approximation of Perdew-Burke-Ernzerhof (GGA-PBE) and the Hubbard interaction U corrected GGA-PBE+U relaxed crystallographic parameters. Scanning and field emission scanning electron micrographs confirmed the particle size of the SCBO and HCBO samples to be ∼250 and ∼60 nm respectively. The GGA-PBE and GGA-PBE+U derived Raman peaks are in better agreement with that of the experimentally observed ones when compared to local density approximation based results. The DFT derived phonon density of states conforms with the absorption bands in Fourier transform infrared spectra. Both structural and dynamic stability criteria of the CBO are confirmed by elastic tensor and density functional perturbation theory-based phonon band structure simulations respectively. The CBO band gap underestimation of GGA-PBE as compared to UV-vis diffuse reflectance derived 1.8 eV was eliminated by tuning the U and the Hartree-Fock exact-exchange mixing parameter αHF in GGA-PBE+U and Heyd-Scuseria-Ernzerhof (HSE06) hybrid functionals respectively. The HSE06 with αHF = 14% yields the optimum linear optical properties of CBO in terms of the dielectric function, absorption, and their derivatives as compared to that of GGA-PBE and GGA-PBE+U functionals. Our as-synthesized HCBO shows ∼70% photocatalytic efficiency in degrading methylene blue dye under 3 h optical illumination. This DFT-guided experimental approach to CBO may help to gain a better understanding of its functional properties.
| S-EPMC10170495 | biostudies-literature
Project description:Reaching the full potential of X-ray nanotomography, in particular for biological samples, is limited by many factors, of which one of the most serious is radiation damage. Although sample deformation caused by radiation damage can be partly mitigated by cryogenic protection, it is still present in these conditions and, as we exemplify here using a specimen extracted from scales of the Cyphochilus beetle, it will pose a limit to the achievable imaging resolution. We demonstrate a generalized tomographic model, which optimally follows the sample morphological changes and attempts to recover the original sample structure close to the ideal, damage-free reconstruction. Whereas our demonstration was performed using ptychographic X-ray tomography, the method can be adopted for any tomographic imaging modality. Our application demonstrates improved reconstruction quality of radiation-sensitive samples, which will be of increasing relevance with the higher brightness of 4th generation synchrotron sources.
| S-EPMC6565693 | biostudies-literature
Project description:AIDA: ab initio domain assembly server, available at http://ffas.burnham.org/AIDA/ is a tool that can identify domains in multi-domain proteins and then predict their 3D structures and relative spatial arrangements. The server is free and open to all users, and there is an option for a user to provide an e-mail to get the link to result page. Domains are evolutionary conserved and often functionally independent units in proteins. Most proteins, especially eukaryotic ones, consist of multiple domains while at the same time, most experimentally determined protein structures contain only one or two domains. As a result, often structures of individual domains in multi-domain proteins can be accurately predicted, but the mutual arrangement of different domains remains unknown. To address this issue we have developed AIDA program, which combines steps of identifying individual domains, predicting (separately) their structures and assembling them into multiple domain complexes using an ab initio folding potential to describe domain-domain interactions. AIDA server not only supports the assembly of a large number of continuous domains, but also allows the assembly of domains inserted into other domains. Users can also provide distance restraints to guide the AIDA energy minimization.
| S-EPMC4086082 | biostudies-literature
Project description:The aim of this study is to describe and compare the supramolecular interactions, in the solid state, of chloro-, bromo-, and iodobenzothiophene diols. The compounds were obtained through organo-catalyzed reactions starting from 3-substituted halobenzothiophene carbaldehydes. Energies of the noncovalent interactions were obtained by density functional theory calculations. Bond distances and angles were found to be in accordance with those determined by X-ray structure analysis. anti-Bromobenzothiophene derivatives showed strong halogen???? interactions between bromine and the heterocyclic phenyl ring, corresponding to an energy of 7.5?kcal?mol(-1). syn-Bromo and syn-iodo derivatives appeared to be isostructural, showing X???O (carbonyl) interactions, ??stacking, and formation of extended hydrogen bonding networks. In contrast, the chloro derivatives displayed no halogen bonding interactions.
| S-EPMC4420588 | biostudies-literature