Dataset Information

Modeling the vibrational couplings of nucleobases.

ABSTRACT: Vibrational spectroscopy, in particular infrared spectroscopy, has been widely used to probe the three-dimensional structures and conformational dynamics of nucleic acids. As commonly used chromophores, the C=O and C=C stretch modes in the nucleobases exhibit distinct spectral features for different base pairing and stacking configurations. To elucidate the origin of their structural sensitivity, in this work, we develop transition charge coupling (TCC) models that allow one to efficiently calculate the interactions or couplings between the C=O and C=C chromophores based on the geometric arrangements of the nucleobases. To evaluate their performances, we apply the TCC models to DNA and RNA oligonucleotides with a variety of secondary and tertiary structures and demonstrate that the predicted couplings are in quantitative agreement with the reference values. We further elucidate how the interactions between the paired and stacked bases give rise to characteristic IR absorption peaks and show that the TCC models provide more reliable predictions of the coupling constants as compared to the transition dipole coupling scheme. The TCC models, together with our recently developed through-bond coupling constants and vibrational frequency maps, provide an effective theoretical strategy to model the vibrational Hamiltonian, and hence the vibrational spectra of nucleic acids in the base carbonyl stretch region directly from atomistic molecular simulations.

SUBMITTER: Jiang Y

PROVIDER: S-EPMC7046491 | biostudies-literature | 2020 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Modeling the vibrational couplings of nucleobases.

Jiang Yaoyukun Y Wang Lu L

The Journal of chemical physics 20200201 8

Vibrational spectroscopy, in particular infrared spectroscopy, has been widely used to probe the three-dimensional structures and conformational dynamics of nucleic acids. As commonly used chromophores, the C=O and C=C stretch modes in the nucleobases exhibit distinct spectral features for different base pairing and stacking configurations. To elucidate the origin of their structural sensitivity, in this work, we develop transition charge coupling (TCC) models that allow one to efficiently calcu ...[more]

PMID: 32113367

Similar Datasets

Project description:The solution conformation behavior of a dermatan-derived tetrasaccharide--Delta HexA-(1-->3)-GalNAc4S-beta-(1-->4)-IdoA-alpha-(1-->3)-red-GalNAc4S (S is a sulfate group)--has been explored by means of NMR spectroscopy, especially by NOE-based conformational analysis. The tetrasaccharide was present as four species, two of which are chemically different in the anomeric orientation of the reducing 2-deoxy-2-acetamido-galactose (red-GalNAc) residue, while the other two are the result of different conformations of the iduronic acid (IdoA) unit. The two alpha-beta-interconverting anomers were present in a 0.6:1 ratio. Ring conformations have been defined by analysis of (3)J(H,H) coupling constants and interresidual NOE contacts. Both 2-deoxy-2-acetamido-galactose (GalNAc) residues were found in the (4)C(1) chair conformation, the unsaturated uronic acid (Delta-Hex A) adopts a strongly predominant half-chair (1)H(2) conformation, while the IdoA residue exists either in the (1)C(4) chair or in the (2)S(0) skewed boat geometries, in a 4:1 ratio. There is a moderate flexibility of Phi and Psi torsions as suggested by nuclear Overhauser effects (NOEs), molecular modeling (MM), and molecular dynamics (MD) studies. This was further investigated by residual dipolar couplings (RDCs). One-bond C--H RDCs ((1)D(C,H)) and long-range H-H ((3)D(H,H)) RDCs were measured for the tetrasaccharide in a phage solution and interpreted in combination with restrained MD simulation. The RDC-derived data substantially confirmed the validity of the conformer distribution resulting from the NOE-derived simulations, but allowed an improved definition of the conformational behavior of the oligosaccharides in solution. In summary, the data show a moderate flexibility of the four tetrasaccharide species at the central glycosidic linkage. Differences in the shapes of species with the IdoA in skew and in chair conformations and in the distribution of the sulfate groups have also been highlighted.

Project description:Infrared spectroscopy is playing an important role in the elucidation of amyloid fiber formation, but the coupling models that link spectra to structure are not well tested for parallel ?-sheets. Using a synthetic macrocycle that enforces a two stranded parallel ?-sheet conformation, we measured the lifetimes and frequency for six combinations of doubly (13)C?(18)O labeled amide I modes using 2D IR spectroscopy. The average vibrational lifetime of the isotope labeled residues was 550 fs. The frequencies of the labels ranged from 1585 to 1595 cm(-1), with the largest frequency shift occurring for in-register amino acids. The 2D IR spectra of the coupled isotope labels were calculated from molecular dynamics simulations of a series of macrocycle structures generated from replica exchange dynamics to fully sample the conformational distribution. The models used to simulate the spectra include through-space coupling, through-bond coupling, and local frequency shifts caused by environment electrostatics and hydrogen bonding. The calculated spectra predict the line widths and frequencies nearly quantitatively. Historically, the characteristic features of ?-sheet infrared spectra have been attributed to through-space couplings such as transition dipole coupling. We find that frequency shifts of the local carbonyl groups due to nearest neighbor couplings and environmental factors are more important, while the through-space couplings dictate the spectral intensities. As a result, the characteristic absorption spectra empirically used for decades to assign parallel ?-sheet secondary structure arises because of a redistribution of oscillator strength, but the through-space couplings do not themselves dramatically alter the frequency distribution of eigenstates much more than already exists in random coil structures. Moreover, solvent exposed residues have amide I bands with >20 cm(-1) line width. Narrower line widths indicate that the amide I backbone is solvent protected inside the macrocycle. This work provides calculated and experimentally verified couplings for parallel ?-sheets that can be used in structure-based models to simulate and interpret the infrared spectra of ?-sheet containing proteins and protein assemblies, such as amyloid fibers.

Project description:Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.

Dataset Information

Modeling the vibrational couplings of nucleobases.

Publications

Modeling the vibrational couplings of nucleobases.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets