Vibrational Properties of the Phosphate Group Investigated by Molecular Dynamics and Density Functional Theory

What are the minimal folding seeds in proteins? Experimental and theoretical assessment of secondary structure propensities of small peptide fragments

Certain peptide sequences, some of them as short as amino acid triplets, are significantly overpopulated in specific secondary structure motifs in folded protein structures. For example, 74% of the EAM triplet is found in α-helices, and only 3% occurs in the extended parts of proteins (typically β-sheets). In contrast, other triplets (such as VIV and IYI) appear almost exclusively in extended parts (79% and 69%, respectively). In order to determine whether such preferences are structurally encoded in a particular peptide fragment or appear only at the level of a complex protein structure, NMR, VCD, and ECD experiments were carried out on selected tripeptides: EAM (denoted as pro-'α-helical' in proteins), KAM(α), ALA(α), DIC(α), EKF(α), IYI(pro-β-sheet or more generally, pro-extended), and VIV(β), and the reference α-helical CATWEAMEKCK undecapeptide. The experimental data were in very good agreement with extensive quantum mechanical conformational sampling. Altogether, we clearly showed that the pro-helical vs. pro-extended propensities start to emerge already at the level of tripeptides and can be fully developed at longer sequences. We postulate that certain short peptide sequences can be considered minimal "folding seeds". Admittedly, the inherent secondary structure propensity can be overruled by the large intramolecular interaction energies within the folded and compact protein structures. Still, the correlation of experimental and computational data presented herein suggests that the secondary structure propensity should be considered as one of the key factors that may lead to understanding the underlying physico-chemical principles of protein structure and folding from the first principles.

Experimental detection of conformational transitions between forms of DNA: problems and prospects

Under different conditions, the DNA double helix can take different geometric forms. Of the large number of its conformations, in addition to the "canonical" B form, the A, C, and Z forms are widely known, and the D, Hoogsteen, and X forms are less known. DNA locally takes the A, C, and Z forms in the cell, in complexes with proteins. We compare different methods for detecting non-canonical DNA conformations: X-ray, IR, and Raman spectroscopy, linear and circular dichroism in both the infrared and ultraviolet regions, as well as NMR (measurement of chemical shifts and their anisotropy, scalar and residual dipolar couplings and inter-proton distances from NOESY (nuclear Overhauser effect spectroscopy) data). We discuss the difficulties in applying these methods, the problems of theoretical interpretation of the experimental results, and the prospects for reliable identification of non-canonical DNA conformations.

Methods to functionalize gold nanoparticles with tandem-phosphorothioate DNA: role of physicochemical properties of the phosphorothioate backbone in DNA adsorption to gold nanoparticles

Perception of the differences in the physicochemical properties of phosphorothioate DNA (PS-DNA) and phosphodiester DNA (PO-DNA) greatly aids in understanding the AuNP-DNA binding process. Replacing one non-bridging oxygen atom of the anionic phosphodiester backbone with a sulfur atom leads to a major change in the DNA adsorption mechanism of AuNPs. In this work, we investigated and compared salt-aging, low pH-assisted, and freeze-thaw methods for conjugating phosphorothioate-modified oligonucleotides to AuNPs. The results obtained clearly demonstrate that only the pH-assisted method can successfully bind tandem phosphorothioate DNA to gold nanoparticles and sufficiently maintain the colloidal stability of AuNPs. When a phosphate group is converted to a phosphorothioate group, the negative charge of the phosphate group is located on the sulfur atom. Due to the soft nature of sulfur (a very weak H-bond acceptor), the negative charge on the sulfur atom cannot be shielded even with the gradual addition of salt to increase the ionic strength, so, the pH-assisted based method is the best for the functionalization of AuNPs with tandem-PS DNA.

Improving IDP theoretical chemical shift accuracy and efficiency through a combined MD/ADMA/DFT and machine learning approach

This work extends the multi-scale computational scheme for the quantum mechanics (QM) calculations of Nuclear Magnetic Resonance (NMR) chemical shifts (CSs) in proteins that lack a well-defined 3D structure. The scheme couples the sampling of an intrinsically disordered protein (IDP) by classical molecular dynamics (MD) with protein fragmentation using the adjustable density matrix assembler (ADMA) and density functional theory (DFT) calculations. In contrast to our early investigation on IDPs (Pavlíková Přecechtělová et al., J. Chem. Theory Comput., 2019, 15, 5642-5658) and the state-of-the art NMR calculations for structured proteins, a partial re-optimization was implemented on the raw MD geometries in vibrational normal mode coordinates to enhance the accuracy of the MD/ADMA/DFT computational scheme. In addition, machine-learning based cluster analysis was performed on the scheme to explore its potential in producing protein structure ensembles (CLUSTER ensembles) that yield accurate CSs at a reduced computational cost. The performance of the cluster-based calculations is validated against results obtained with conventional structural ensembles consisting of MD snapshots extracted from the MD trajectory at regular time intervals (REGULAR ensembles). CS calculations performed with the refined MD/ADMA/DFT framework employed the 6-311++G(d,p) basis set that outperformed IGLO-III calculations with the same density functional approximation (B3LYP) and both explicit and implicit solvation. The partial geometry optimization did not universally improve the agreement of computed CSs with the experiment but substantially decreased errors associated with the ensemble averaging. A CLUSTER ensemble with 50 structures yielded ensemble averages close to those obtained with a REGULAR ensemble consisting of 500 MD frames. The cluster based calculations thus required only a fraction of the computational time.

The mutual interactions of RNA, counterions and water - quantifying the electrostatics at the phosphate-water interface

The structure and dynamics of polyanionic biomolecules, like RNA, are decisively determined by their electric interactions with the water molecules and the counterions in the environment. The solvation dynamics of the biomolecules involves a subtle balance of non-covalent and many-body interactions with structural fluctuations due to thermal motion occurring in a femto- to subnanosecond time range. This complex fluctuating many particle scenario is crucial in defining the properties of biological interfaces with far reaching significance for the folding of RNA structures and for facilitating RNA-protein interactions. Given the inherent complexity, suited model systems, carefully calibrated and benchmarked by experiments, are required to quantify the relevant interactions of RNA with the aqueous environment. In this feature article we summarize our recent progress in the understanding of the electrostatics at the biological interface of double stranded RNA (dsRNA) and transfer RNA (tRNA). Dimethyl phosphate (DMP) is introduced as a viable and rigorously accessible model system allowing the interaction strength with water molecules and counterions, their relevant fluctuation timescales and the spatial reach of interactions to be established. We find strong (up to ≈90 MV cm-1) interfacial electric fields with fluctuations extending up to ≈20 THz and demonstrate how the asymmetric stretching vibration νAS(PO2)- of the polarizable phosphate group can serve as the most sensitive probe for interfacial interactions, establishing a rigorous link between simulations and experiment. The approach allows for the direct interfacial observation of interactions of biologically relevant Mg2+ counterions with phosphate groups in contact pair geometries via the rise of a new absorption band imposed by exchange repulsion interactions at short interatomic distances. The systematic extension to RNA provides microscopic insights into the changes of the hydration structure that accompany the temperature induced melting of the dsRNA double helix and quantify the ionic interactions in the folded tRNA. The results show that pairs of negatively charged phosphate groups and Mg2+ ions represent a key structural feature of RNA embedded in water. They highlight the importance of binding motifs made of contact pairs in the electrostatic stabilization of RNA structures that have a strong impact on the surface potential and enable the fine tuning of the local electrostatic properties which are expected to be relevant for mediating the interactions between biomolecules.

How important are the intermolecular hydrogen bonding interactions in methanol solvent for interpreting the chiroptical properties?

Two crispine A analogs and tetrahydrofuro[2,3-b]furan-3,3a(6aH)-diol, endowed with hydroxyl groups that can participate in intramolecular hydrogen bonding, have been synthesized and experimental vibrational circular dichroism (VCD) spectra and optical rotatory dispersion (ORD) data have been measured in CD₃OD/CH₃OH solvents. The absolute configurations (ACs) of these compounds have been determined using their synthetic schemes, supplemented wherever possible with X-ray diffraction data. The ACs are also analyzed with quantum chemical (QC) calculations of VCD and ORD utilizing implicit solvation as well as explicit solvation models, with the later employing classical molecular dynamics (MD) simulations. It is found that VCD calculations with implicit solvation model are adequate for determining the ACs, despite propensity of studied compounds for intermolecular hydrogen bonding between solute and solvent molecules. This observation is important because time-consuming MD simulations may not be necessary in the type of situations studied here. Additionally, it is found that the QC predicted VCD spectra provided enough diastereomer discrimination for determining the correct AC of studied compounds independently. The same observation did not apply to ORD.

Origins of Optical Activity in an Oxo-Helicene: Experimental and Computational Studies

Helicenes are known to provide extremely strong optical activity. Prediction of the properties of helicenes may facilitate their design and synthesis for analytical or materials sciences. On a model 7,12,17-trioxa[11]helicene molecule, experimental results from multiple spectroscopic techniques are analyzed on the basis of density functional theory (DFT) simulations to test computational methodology and analyze the origins of chirality. Infrared (IR), vibrational circular dichroism (VCD), electronic circular dichroism (ECD), magnetic circular dichroism (MCD), and Raman optical activity (ROA, computations only) spectra are compared. Large dissymmetry factors are predicted both for vibrational (ROA/Raman ∼ VCD/IR ∼ 10^-3) and electronic (ECD/Abs ∼10^-2) optical activity, which could be verified experimentally except for ROA. Largest VCD signals come from a strong vibrational coupling of the C-H in-plane and out-of-plane bending modes in stacked helicene rings. The sum-over-states (SOS) approach appeared convenient for simulation of MCD spectra. Our results demonstrated that selected computational methods can be successfully used for reliable modeling of spectral and chiroptical properties of large helicenes. In particular, they can be used for guiding rational design of strongly chiral chromophores.

Deconvolution of conformational exchange from Raman spectra of aqueous RNA nucleosides

Ribonucleic acids (RNAs) are key to the central dogma of molecular biology. While Raman spectroscopy holds great potential for studying RNA conformational dynamics, current computational Raman prediction and assignment methods are limited in terms of system size and inclusion of conformational exchange. Here, a framework is presented that predicts Raman spectra using mixtures of sub-spectra corresponding to major conformers calculated using classical and ab initio molecular dynamics. Experimental optimization allowed purines and pyrimidines to be characterized as predominantly syn and anti, respectively, and ribose into exchange between equivalent south and north populations. These measurements are in excellent agreement with Raman spectroscopy of ribonucleosides, and previous experimental and computational results. This framework provides a measure of ribonucleoside solution populations and conformational exchange in RNA subunits. It complements other experimental techniques and could be extended to other molecules, such as proteins and carbohydrates, enabling biological insights and providing a new analytical tool.

Effects of counterions and solvents on the geometrical and vibrational features of dinucleoside-monophosphate (dNMP): case of 3',5'-dideoxycytidine-monophosphate (dDCMP)

The effects of the interaction of the monovalent (Li⁺, Na⁺, K⁺) and divalent (Mg²⁺) counterions hexahydrated (6H₂O), with the PO₂^- group, on the geometrical and vibrational characteristics of 3', 5'-dDCMP, were studied using the DFT/B3LYP/6-31++G(d) method. These calculations were performed using the explicit (6H₂O) and hybrid (6H₂O/Continuum) solvation models. The optimizations reveal that in the conformation g^-g^- and in the explicit model of solvation, the small ions (Li⁺, Na⁺) deviate from the bisector plane of the angle O1-P-O2 and the large ions (K⁺ and Mg²⁺) remain in this plane, whereas in the hybrid model of solvation, the counterions deviate from this plane. However, when the conformer is g⁺g⁺, the monovalent counterions deviate and divide the remainder of the plane regardless of the type of solvation model. In addition, the g^-g^- conformer is the most stable in the presence of the explicit solvent, while the g⁺g⁺ conformer is the most stable in the presence of the hybrid solvent. Finally, the normal modes of the conformers g^-g^- and g⁺g⁺ in the presence of the counterions in the hybrid model show a better agreement with the available experimental data of the DNA forms A, B (g^-g^-), and Z (g⁺g⁺) relatively to the explicit model. This very good agreement is illustrated by the very small deviations ≤ 0.08% (g^-g^-) and ≤ 0.41% (g⁺g⁺) observed between the calculated and experimental data for the PO₂^- (asymmetric) stretching mode in the presence of the counterion K⁺ in the hybrid model. Graphical abstract.

Luminescence properties of tricalcium phosphate doped with dysprosium

Tricalcium phosphate having effective atomic number Z_eff = 15.785, equivalent to that of bones was studied for its thermoluminescence (TL) and photoluminescence (PL) properties. Different samples with varied concentrations of the dopant Dy³⁺ (0.1, 0.2, 0.3, 0.4, 0.5, and 0.6 mol %) were synthesized by the chemical co-precipitation technique. Phase and compound confirmations were done using X-Ray Diffraction (XRD) and the phosphors' crystallite size was calculated using Scherrer's formula which was found to range between 27 nm and 49 nm. The surface morphological study was done using Field Emission-Scanning Electron Microscopy (FE-SEM). Other characterization techniques used for compound confirmation included Energy Dispersive X-Ray Spectroscopy (EDX) and Fourier Transform Infrared Spectroscopy (FTIR). Samples were further irradiated by gamma rays (emitted from Co-60) with dose varying from 10 Gy to 5 kGy in order to study their TL properties. Concentration optimization of the given phosphor was done in terms of its TL intensity and was found to be 0.5 mol %. The TL dose response of the phosphor was linear over a wide range of dose (10 Gy-3 kGy). Deconvolution was performed on the glow curve for 10 Gy dose, giving six peaks (at 127^o, 153^o, 185^o, 218^o, 313^o and 335^oC) suggesting the presence of six different types of traps. Other characteristics of the TL material i.e. repeatability and fading were also studied. Overall, the phosphor showed promising results for its utility in TL dosimetry. In addition to the TL, PL further confirmed the presence of dopant in the phosphor. Moreover, the dopant concentration was optimized in terms of the nanophosphor's PL intensity. The Commission International de l'Éclairage (CIE) was used to calculate chromaticity coordinates, colour rendering index and colour temperature in order to investigate the phosphor's application in white LEDs.

Influence of the hydrophobic domain on the self-assembly and hydrogen bonding of hydroxy-amphiphiles

The amphiphiles 1-octadecanol (octadecyl (stearyl) alcohol, ODA) and 1,2-dioleoylglycerol (DOG) were studied by IR spectroscopy and X-ray diffraction combined with multiscale theoretical modeling. The computations allowed us to rationalize the experimental findings and deduce the supramolecular structure of the formed assemblies while providing a fairly detailed insight into their hydrogen-bonding patterns. IR spectra revealed that the amphiphilic assemblies dramatically differ in structural order and hydrogen-bond strength, both being high in ODA and low in DOG. On the other hand, both compounds demonstrated common features, namely a splitting of the IR bands arising from O-H stretching vibrations (νOH) as well as complete hydrophobicity. However, the observed phenomena have different origins in the two amphiphiles. While the νOH split in ODA occurs due to a vibrational coupling along the string of inter-layer O-HO hydrogen bonds, in DOG it arises from different types of hydrogen bonds (intra- and intermolecular). The hydrophobicity of ODA stems from the very tight O-HO hydrogen bonding network connecting the opposite monolayers in a densely packed tilted crystalline phase (L_c'), whereas in DOG it occurs because the polar sites are locked inside reverted micellar-like assemblies. ODA and DOG illustrate that, in the assemblies of amphiphilic hydroxyl compounds, hydrogen bonds can be formed in a wide structural latitude, which is primarily governed by the chemical nature of apolar chains. Such a wide structural variability of OH-involving hydrogen bonds can be essential for the biological functioning of relevant molecules, such as glycolipids, acylglycerols, and, potentially, glycoproteins and carbohydrates.

Insight into vibrational circular dichroism of proteins by density functional modeling

Vibrational circular dichroism (VCD) spectroscopy is an excellent method to determine the secondary structure of proteins in solution. Comparison of experimental spectra with quantum-chemical simulations represents a convenient and objective way to extract information on the structure. This has been difficult for such large molecules where approximate theoretical models have to be used. In the present study we applied the Cartesian-coordinate based tensor transfer (CCT) making it possible to extend the density functional theory (DFT) and model spectral intensities of large globular proteins nearly at quantum-chemical precision. Indeed, comparison with experiment provided a better understanding of the dependence of VCD spectral shapes on the geometry, their sensitivity to fine structural details and interactions with the environment. On a model set of globular proteins the simulated spectra correlated well with experimental data and revealed which structural information can (and cannot) be obtained from this kind of spectroscopy. Although the VCD technique has been regarded as being rather insensitive to side-chain variations, we found that the spectra of human and hen lysozyme differing by a few amino acids only are quite distinct. This has been explained by long-distance coupling of the amide vibrations. Likewise, the modeling reproduced some spectral changes caused by protein deuteration even when the protein structure was conserved.