Abil E. Aliev

Motional timescale predictions by molecular dynamics simulations: case study using proline and hydroxyproline sidechain dynamics.

We propose a new approach for force field optimizations which aims at reproducing dynamics characteristics using biomolecular MD simulations, in addition to improved prediction of motionally averaged structural properties available from experiment. As the source of experimental data for dynamics fittings, we use 13C NMR spin‐lattice relaxation times T1 of backbone and sidechain carbons, which allow to determine correlation times of both overall molecular and intramolecular motions. For structural fittings, we use motionally averaged experimental values of NMR J couplings. The proline residue and its derivative 4‐hydroxyproline with relatively simple cyclic structure and sidechain dynamics were chosen for the assessment of the new approach in this work. Initially, grid search and simplexed MD simulations identified large number of parameter sets which fit equally well experimental J couplings. Using the Arrhenius‐type relationship between the force constant and the correlation time, the available MD data for a series of parameter sets were analyzed to predict the value of the force constant that best reproduces experimental timescale of the sidechain dynamics. Verification of the new force‐field (termed as AMBER99SB‐ILDNP) against NMR J couplings and correlation times showed consistent and significant improvements compared to the original force field in reproducing both structural and dynamics properties. The results suggest that matching experimental timescales of motions together with motionally averaged characteristics is the valid approach for force field parameter optimization. Such a comprehensive approach is not restricted to cyclic residues and can be extended to other amino acid residues, as well as to the backbone. Proteins 2014; 82:195–215.

Experimental verification of force fields for molecular dynamics simulations using Gly-Pro-Gly-Gly.

Experimental NMR verification of MD simulations using 12 different force fields (AMBER, CHARMM, GROMOS, and OPLS-AA) and 5 different water models has been undertaken to identify reliable MD protocols for structure and dynamics elucidations of small open chain peptides containing Gly and Pro. A conformationally flexible tetrapeptide Gly-Pro-Gly-Gly was selected for NMR (3)J-coupling, chemical shift, and internuclear distance measurements, followed by their calculations using 2 μs long MD simulations in water. In addition, Ramachandran population maps for Pro-2 and Gly-3 residues of GPGG obtained from MD simulations were used for detailed comparisons with similar maps from the protein data bank (PDB) for large number of Gly and Pro residues in proteins. The MD simulations revealed strong dependence of the populations and geometries of preferred backbone and side chain conformations, as well as the time scales of the peptide torsional transitions on the force field used. On the basis of the analysis of the measured and calculated data, AMBER99SB is identified as the most reliable force field for reproducing NMR measured parameters, which are dependent on the peptide backbone and the Pro side chain geometries and dynamics. Ramachandran maps showing the dependence of conformational populations as a function of backbone ϕ/ψ angles for Pro-2 and Gly-3 residues of GPGG from MD simulations using AMBER99SB, AMBER03, and CHARMM were found to resemble similar maps for Gly and Pro residues from the PDB survey. Three force fields (AMBER99, AMBER99ϕ, and AMBER94) showed the least satisfactory agreement with both the solution NMR and the PDB survey data. The poor performance of these force fields is attributed to their propensity to overstabilize helical peptide backbone conformations at the Pro-2 and Gly-3 residues. On the basis of the similarity of the MD and PDB Ramachandran plots, the following sequence of transitions is suggested for the Gly backbone conformation: α(L) ⇆ β(PR) ⇆ β(S) ⇆ β(P) ⇆ α, where backbone secondary structures α(L) and α are associated with helices and turns, β(P) and β(PR) correspond to the left- and right-handed polyproline II structures and β(S) denotes the fully stretched backbone conformation. Compared to the force field dependence, less significant, but noteworthy, variations in the populations of the peptide backbone conformations were observed. For different solvent models considered, a correlation was noted between the number of torsional transitions in GPGG and the water self-diffusion coefficient on using TIP3P, TIP4P, and TIP5P models. In addition to MD results, we also report DFT derived Karplus relationships for Gly and Pro residues using B972 and B3LYP functionals.

Photodetachment Spectra of Deprotonated Fluorescent Protein Chromophore Anions

Isolated model anion chromophores of the green and cyan fluorescent proteins were generated in an electrospray ion source, and their photodetachment spectra were recorded using photoelectron imaging. Vertical photodetachment energies of 2.85(10) and 4.08(10) eV have been measured for the model green fluorescent protein chromophore anion, corresponding to photodetachment from the ground electronic state of the anion to the ground and first excited electronic states of the radical, respectively. For the model cyan fluorescent protein chromophore anion, vertical photodetachment energies of 2.88(10) and 3.96(10) eV have been measured, corresponding to detachment from the ground electronic state of the anion to the ground and first excited electronic states of the neutral radical, respectively. We also find evidence suggesting that autoionization of electronically excited states of the chromophore anions competes with direct photodetachment. For comparison and to benchmark our measurements, the vertical photodetachment energies of deprotonated phenol and indole anions have also been recorded and presented. Quantum chemistry calculations support our assignments. We discuss our results in the context of the isolated protein chromophore anions acting as electron donors, one of their potential biological functions.

Dynamic mechanical analysis (DMA), 13c solid state nmr and micro-thermomechanical studies of historical parchment

DMA and solid state 13C NMR techniques were used to measure historical parchment samples within the framework of the project (MAP) Micro Analysis of Parchment (EC contract No. SMT4-96-2101) in collaboration with the School of Conservation in Copenhagen. DMA was used in both thermal scan and creep modes. Thermal scans provided information on the transitions associated with the collagen polymer. Microthermal analysis was also used to obtain information on the topography and thermal conductivity of sample areas of 100 μm. Localised heating enabled measurements of softening transitions in the sample. This behaviour is influenced by the chemical composition of parchment. 13C NMR provided information on the carbon atoms associated with the polypeptide chains of the collagen in parchment. The behaviour of samples immersed in water and measured in DMA creep mode was used to measure the shrinkage behaviour of the parchment samples. The different but complementary techniques provided a means for characterising the physicochemical state of parchment samples.

An NMR Method for the Quantitative Assessment of Intramolecular Hydrogen Bonding; Application to Physicochemical, Environmental, and Biochemical Properties

(1)H NMR chemical shifts have been obtained in the solvents deuterochloroform and dimethyl sulfoxide. The difference in the chemical shifts of an OH or NH group in these two solvents, Δδ = δ(DMSO) - δ(CDCl3), can be converted into the hydrogen bond acidity, A, of the group using the equation A = 0.0065 + 0.133Δδ. The NMR A value, ANMR, can be used as a quantitative assessment of intramolecular hydrogen bonding. We list values of Δδ and ANMR for 55 compounds containing an OH group and 60 compounds with an NH group. For the hydroxy compounds, if A > 0.5 then the OH group is not part of an intramolecular hydrogen bond, but if A < 0.1 then the OH group forms part of an intramolecular hydrogen bond. For NH compounds, if A > 0.16 the NH group is not part of an intramolecular hydrogen bond, and if A < 0.05 the NH group is part of an intramolecular hydrogen bond. No comparison compounds are needed, and the method is extremely simple. We further show how it is possible to relate intramolecular hydrogen bonding to the actual effect on values of a number of physicochemical, environmental, and biochemical properties.

Regioselective Dihalohydration Reactions of Propargylic Alcohols: Gold‐Catalyzed and Noncatalyzed Reactions

The regioselective conversion of propargylic alcohols into previously unreported α,α-diiodo-β-hydroxyketones was achieved by treatment with N-iodosuccinimide in the presence of a gold catalyst. The corresponding α,α-dichloro-β-hydroxyketones were obtained by treatment with trichloroisocyanuric acid in the absence of a catalyst. The latter reaction can be extended to other alkynols. These transformations can be used to prepare potentially useful halogenated building blocks. Preliminary mechanistic studies suggest that the reaction involves participation of the acetonitrile solvent in the formation of a 5-halo-1,3-oxazine intermediate.

High-resolution solid-state 2H NMR spectroscopy of polymorphs of glycine

High-resolution solid-state (2)H MAS NMR studies of the α and γ polymorphs of fully deuterated glycine (glycine-d(5)) are reported. Analysis of spinning sideband patterns is used to determine the (2)H quadrupole interaction parameters, and is shown to yield good agreement with the corresponding parameters determined from single-crystal (2)H NMR measurements (the maximum deviation in quadrupole coupling constants determined from these two approaches is only 1%). From analysis of simulated (2)H MAS NMR sideband patterns as a function of reorientational jump frequency (κ) for the -N(+)D(3) group in glycine-d(5), the experimentally observed differences in the (2)H MAS NMR spectrum for the -N(+)D(3) deutrons in the α and γ polymorphs is attributed to differences in the rate of reorientation of the -N(+)D(3) group. These simulations show severe broadening of the (2)H MAS NMR signal in the intermediate motion regime, suggesting that deuterons undergoing reorientational motions at rates in the range κ ≈ 10(4)-10(6) s(-1) are likely to be undetectable in (2)H MAS NMR measurements for materials with natural isotopic abundances. The (1)H NMR chemical shifts for the α and γ polymorphs of glycine have been determined from the (2)H MAS NMR results, taking into account the known second-order shift. Further quantum mechanical calculations of (2)H quadrupole interaction parameters and (1)H chemical shifts reveal the structural dependence of these parameters in the two polymorphs and suggest that the existence of two short intermolecular C-H···O contacts for one of the H atoms of the >CH(2) group in the α polymorph have a significant influence on the (2)H quadrupole coupling and (1)H chemical shift for this site.

Dynamic properties of dioctanoyl peroxide guest molecules constrained within the urea tunnel structure: A combined incoherent quasielastic neutron scattering and solid state 2H nuclear magnetic resonance investigation

The dynamic properties of dioctanoyl peroxide guest molecules within the urea host tunnel structure in the dioctanoyl peroxide/urea inclusion compound have been investigated by incoherent quasielastic neutron scattering (IQNS) and solid state 2H nuclear magnetic resonance (NMR) techniques. The IQNS investigations were carried out on samples of urea inclusion compounds containing perdeuterated urea to ensure that the incoherent scattering is dominated by the dioctanoyl peroxide guest molecules. Using semioriented polycrystalline samples, translational motions of the guest molecules along the tunnel were investigated separately from reorientational motions of the guest molecules about the tunnel axis. The 2H NMR experiments used dioctanoyl peroxide deuterated selectively in both the α CD2 groups and urea with natural isotopic abundance. The dynamic models that have been found to describe the translational and reorientational motions of the guest molecules from the IQNS and 2H NMR data are discussed in detai...

Quantum Mechanical and NMR Studies of Ring Puckering and cis/trans-Rotameric Interconversion in Prolines and Hydroxyprolines

Nuclear magnetic resonance (NMR) and quantum mechanical (QM) studies have been carried out for proline (Pro) containing peptides: N-acetyl-l-proline (AcProOH) and N-acetyl-4-hydroxy-l-proline (AcHypOH). Preliminary results of variable temperature NMR measurements for Gly-Pro-Gly-Gly (GPGG), Val-Ala-Pro-Gly (VAPG), and Ala-Pro-Gly-Trp amide acetate salt (APGW) are also reported. The effect of solvent (D(2)O, DMSO-d(6) and CD(3)CN) on the pyrrolidine ring conformation and cis/trans-rotamerisation along the amide bond preceding Pro was investigated by temperature dependent NMR followed by detailed transition state (TS) searches for both conformational equilibria using QM methods. The results revealed the energetic characteristics of the TS, which were in satisfactory agreement with NMR, and the corresponding TS geometries, which are not available from experiment. The most remarkable feature of the cis/trans-rotamerisation is that the amide nitrogen in AcProOH and AcHypOH adopts a tetrahedral geometry in the TS. Various HF, DFT, and MP2 calculations together with implicit solvation modeling were employed in order to identify the most suitable QM protocols for reliable predictions of the geometry and the relative energies of the conformations of Pro and Hyp containing peptides in aqueous solution. Solution NMR results were used for the verification of the reliability of the QM predictions. The results indicate that the MP2 calculations combined with implicit solvation models are reasonably accurate in reproducing NMR measured populations of four different conformations of either AcProOH or AcHypOH in different solvents, whereas HF and DFT B3LYP calculations were significantly less accurate.

Cytosine modules in quadruple hydrogen bonded arrays

Cytosine modules have been investigated for applications in supramolecular quadruple hydrogen bonded arrays. Notably, the importance of the C-5–H in the formation of unfolded and folded arrays by substitution to C-5–F was established. In addition, the incorporation of different alkyl chain lengths at N-1 and N-9 indicated that longer alkyl chains give rise to more of the unfolded rotamer, with the chain length and degree of unsaturation at N-1 having the major effect. Methyl cytosine modules were also able to readily form hetero-associated Upy–UCyt dimers as efficiently as the hexyl cytosine modules and a polyadipate telechelic polymer was used to prepare cytosine polymers.