Probing the Conical Intersection Dynamics of the RNA Base Uracil by UV-Pump Stimulated-Raman-Probe Signals; Ab Initio Simulations.
Ontology highlight
ABSTRACT:
SUBMITTER: Fingerhut BP
PROVIDER: S-EPMC3958139 | biostudies-other | 2014 Mar
REPOSITORIES: biostudies-other
Similar Datasets
Project description:The C, N and O 1s XPS spectra of uracil clusters in the gas phase have been measured. A new bottom-up approach, which relies on computational simulations starting from the crystallographic structure of uracil, has been adopted to interpret the measured spectra. This approach sheds light on the different molecular interactions (H-bond, π-stacking, dispersion interactions) at work in the cluster and provides a good understanding of the observed XPS chemical shifts with respect to the isolated molecule in terms of intramolecular and intermolecular screening occurring after the core-hole ionization. The proposed bottom-up approach, reasonably expensive in terms of computational resources, has been validated by finite-temperature molecular dynamics simulations of clusters composed of up to fifty molecules.
| S-EPMC7403593 | biostudies-literature
Project description:Ab initio molecular dynamics simulations have been performed in order to investigate the relaxation dynamics of uracil after UV excitation in gas phase. Intersystem crossing (ISC) has been included for the first time into time-dependent simulations of uracil, allowing the system to relax in the singlet as well as in the triplet states. The results show a qualitatively different picture than similar simulations that include singlet states only. The inclusion of ISC effectively quenches the relaxation to the singlet ground state and instead privileges transitions from the low-lying n?* state (S1) to a ??* triplet state (T2) followed by rapid internal conversion to the lowest triplet state.
| S-EPMC4391640 | biostudies-literature
Project description:Accurate ab initio simulation of protein folding is a critical step toward elucidation of protein-folding mechanisms. Here, we demonstrate highly accurate folding of the 35 residue villin headpiece subdomain (HP35) by all-atom molecular dynamics simulations using AMBER FF03 and the generalized-Born solvation model. In a set of 20 micros long simulations, the protein folded to the native state in multiple trajectories, with the lowest C(alpha) RMSD being 0.39 A for residues 2-34 (excluding residues 1 and 35). The native state had the highest population among all sampled conformations, and the center of most populated cluster had a C(alpha) RMSD of 1.63 A. Folding of this protein can be described as a two-stage process that followed a well-defined pathway. In the first stage, formation of helices II and III as a folding intermediate constituted the rate-limiting step and was initiated at a folding nucleus around residues Phe17 and Pro21. The folding intermediate further acted as a template that facilitated the folding and docking of helix I in the second stage. Detailed descriptions of the folding kinetics and the roles of key residues are presented.
| S-EPMC2701201 | biostudies-literature
Project description:BackgroundPredicting 3-dimensional protein structures from amino-acid sequences is an important unsolved problem in computational structural biology. The problem becomes relatively easier if close homologous proteins have been solved, as high-resolution models can be built by aligning target sequences to the solved homologous structures. However, for sequences without similar folds in the Protein Data Bank (PDB) library, the models have to be predicted from scratch. Progress in the ab initio structure modeling is slow. The aim of this study was to extend the TASSER (threading/assembly/refinement) method for the ab initio modeling and examine systemically its ability to fold small single-domain proteins.ResultsWe developed I-TASSER by iteratively implementing the TASSER method, which is used in the folding test of three benchmarks of small proteins. First, data on 16 small proteins (< 90 residues) were used to generate I-TASSER models, which had an average Calpha-root mean square deviation (RMSD) of 3.8A, with 6 of them having a Calpha-RMSD < 2.5A. The overall result was comparable with the all-atomic ROSETTA simulation, but the central processing unit (CPU) time by I-TASSER was much shorter (150 CPU days vs. 5 CPU hours). Second, data on 20 small proteins (< 120 residues) were used. I-TASSER folded four of them with a Calpha-RMSD < 2.5A. The average Calpha-RMSD of the I-TASSER models was 3.9A, whereas it was 5.9A using TOUCHSTONE-II software. Finally, 20 non-homologous small proteins (< 120 residues) were taken from the PDB library. An average Calpha-RMSD of 3.9A was obtained for the third benchmark, with seven cases having a Calpha-RMSD < 2.5A.ConclusionOur simulation results show that I-TASSER can consistently predict the correct folds and sometimes high-resolution models for small single-domain proteins. Compared with other ab initio modeling methods such as ROSETTA and TOUCHSTONE II, the average performance of I-TASSER is either much better or is similar within a lower computational time. These data, together with the significant performance of automated I-TASSER server (the Zhang-Server) in the 'free modeling' section of the recent Critical Assessment of Structure Prediction (CASP)7 experiment, demonstrate new progresses in automated ab initio model generation. The I-TASSER server is freely available for academic users http://zhang.bioinformatics.ku.edu/I-TASSER.
| S-EPMC1878469 | biostudies-literature
Project description:Full multiple spawning (FMS) offers a strategy to simulate the nonadiabatic dynamics of molecular systems by describing their nuclear wavefunctions by a linear combination of coupled trajectory basis functions (TBFs). Applying a series of controlled approximations to the full multiple spawning (FMS) equations leads to the ab initio multiple spawning (AIMS), which is compatible with an on-the-fly propagation of the TBFs and an accurate description of nonadiabatic processes. The AIMS strategy and its numerical implementations, however, rely on a series of user-defined parameters. Herein, we investigate the influence of these parameters on the electronic-state population of two molecular systems- trans-azomethane and a two-dimensional model of the butatriene cation. This work highlights the stability of AIMS with respect to most of its parameters, underlines the specific parameters that require particular attention from the user of the method, and offers prescriptions for an informed selection of their value.Supplementary informationThe online version contains supplementary material available at 10.1007/s00214-023-03004-w.
| S-EPMC10382418 | biostudies-literature
Project description:We investigated excess electron solvation dynamics in (NH3)n- ammonia clusters in the n = 8-32 size range by performing finite temperature molecular dynamics simulations. In particular, we focused on three possible scenarios. The first case is designed to model electron attachment to small neutral ammonia clusters (n ≤ ∼10) that form hydrogen-bonded chains. The excess electron is bound to the clusters via dipole bound states, and persists with a VDE of ∼160 meV at 100 K for the n = 8 cluster. The coupled nuclear and electronic relaxation is fast (within ∼100 fs) and takes place predominantly by intermolecular librational motions and the intramolecular umbrella mode. The second scenario illustrates the mechanism of excess electron attachment to cold compact neutral clusters of medium size (8 ≤ n ≤ 32). The neutral clusters show increasing tendency with size to bind the excess electron on the surface of the clusters in weakly bound, diffuse, and highly delocalized states. Anionic relaxation trajectories launched from these initial states provide no indication for excess electron stabilization for sizes n < 24. Excess electrons are likely to autodetach from these clusters. The two largest investigated clusters (n = 24 and 32) can accommodate the excess electron in electronic states that are mainly localized on the surface, but they may be partly embedded in the cluster. In the last 500 fs of the simulated trajectories, the VDE of the solvated electron fluctuates around ∼200 meV for n = 24 and ∼500 meV for n = 32, consistent with the values extrapolated from the experimentally observed linear VDE-n-1/3 trend. In the third case, we prepared neutral ammonia cluster configurations, including an n = 48 cluster, that contain possible electron localization sites within the interior of the cluster. Excess electrons added to these clusters localize in cavities with high VDEs up to 1.9 eV for n = 48. The computed VDE values for larger clusters are considerably higher than the experimentally observed photoelectric threshold energy for the solvated electron in bulk ammonia (∼1.4 eV). Molecular dynamics simulations launched from these initial cavity states strongly indicate, however, that these cavity structures exist only for ∼200 fs. During the relaxation the electron leaves the cavity and becomes delocalized, while the cluster loses its initial compactness.
| S-EPMC7458421 | biostudies-literature
Project description:Lead (Pb) free non-toxic perovskite solar cells have become more important in the commercialization of the photovoltaic devices. In this study the structural, electronic, optical and mechanical properties of Pb-free inorganic metal halide cubic perovskites CsBX3 (B = Sn, Ge; X = I, Br, Cl) for perovskite solar cells are simulated using first-principles Density Functional Theory (DFT). These compounds are semiconductors with direct band gap energy and mechanically stable. Results suggest that the materials have high absorption coefficient, low reflectivity and high optical conductivity with potential application in solar cells and other optoelectronic energy devices. On the basis of the optical properties, one can expect that the Germanium (Ge) would be a better replacement of Pb as Ge containing compounds have higher optical absorption and optical conductivity than that of Pb containing compounds. A combinational analysis of the electronic, optical and mechanical properties of the compounds suggests that CsGeI3 based perovskite is the best Pb-free inorganic metal halide semiconductor for the solar cell application. However, the compound with solid solution of CsGe(I0.7Br0.3)3 is found to be mechanically more ductile than CsGeI3. This study will also guide to obtain Pb-free organic perovskites for optoelectronic devices.
| S-EPMC5656601 | biostudies-literature
Project description:We present a mixed quantum-classical simulation of polariton dynamics for molecule-cavity hybrid systems. In particular, we treat the coupled electronic-photonic degrees of freedom (DOFs) as the quantum subsystem and the nuclear DOFs as the classical subsystem and use the trajectory surface hopping approach to simulate non-adiabatic dynamics among the polariton states due to the coupled motion of nuclei. We use the accurate nuclear gradient expression derived from the Pauli-Fierz quantum electrodynamics Hamiltonian without making further approximations. The energies, gradients, and derivative couplings of the molecular systems are obtained from the on-the-fly simulations at the level of complete active space self-consistent field (CASSCF), which are used to compute the polariton energies and nuclear gradients. The derivatives of dipoles are also necessary ingredients in the polariton nuclear gradient expression but are often not readily available in electronic structure methods. To address this challenge, we use a machine learning model with the Kernel ridge regression method to construct the dipoles and further obtain their derivatives, at the same level as the CASSCF theory. The cavity loss process is modeled with the Lindblad jump superoperator on the reduced density of the electronic-photonic quantum subsystem. We investigate the azomethane molecule and its photoinduced isomerization dynamics inside the cavity. Our results show the accuracy of the machine-learned dipoles and their usage in simulating polariton dynamics. Our polariton dynamics results also demonstrate the isomerization reaction of azomethane can be effectively tuned by coupling to an optical cavity and by changing the light-matter coupling strength and the cavity loss rate.
| S-EPMC10134431 | biostudies-literature
Project description:Characterization of enzyme active site structure and interactions at high resolution is important for the understanding of the enzyme catalysis. Vibrational frequency and NMR chemical shift measurements of enzyme-bound ligands are often used for such purpose when X-ray structures are not available or when higher resolution active site structures are desired. This review is focused on how ab initio calculations may be integrated with vibrational and NMR chemical shift measurements to quantitatively determine high-resolution ligand structures (up to 0.001 Å for bond length and 0.01 Å for hydrogen bonding distance) and how interaction energies between bound ligand and its surroundings at the active site may be determined. Quantitative characterization of substrate ionic states, bond polarizations, tautomeric forms, conformational changes and its interactions with surroundings in enzyme complexes that mimic ground state or transition state can provide snapshots for visualizing the substrate structural evolution along enzyme-catalyzed reaction pathway. Our results have shown that the integration of spectroscopic studies with theoretical computation greatly enhances our ability to interpret experimental data and significantly increases the reliability of the theoretical analysis.
| S-EPMC5484042 | biostudies-literature
Project description:Raman spectroscopy is frequently used to identify composition, structure and layer thickness of 2D materials. Here, we describe an efficient first-principles workflow for calculating resonant first-order Raman spectra of solids within third-order perturbation theory employing a localized atomic orbital basis set. The method is used to obtain the Raman spectra of 733 different monolayers selected from the Computational 2D Materials Database (C2DB). We benchmark the computational scheme against available experimental data for 15 known monolayers. Furthermore, we propose an automatic procedure for identifying a material based on an input experimental Raman spectrum and apply it to the cases of MoS2 (H-phase) and WTe2 (T[Formula: see text]-phase). The Raman spectra of all materials at different excitation frequencies and polarization configurations are freely available from the C2DB. Our comprehensive and easily accessible library of ab initio Raman spectra should be valuable for both theoreticians and experimentalists in the field of 2D materials.
| S-EPMC7296020 | biostudies-literature