l-Alanine in a Droplet of Water: A Density-Functional Molecular Dynamics Study

Relaxation dynamics of l-alanine in water medium investigated by dielectric relaxation spectroscopy

The complex dielectric permittivity of L-alanine in aqueous medium at different concentrations and different temperatures were measured in the microwave (0.02 < ν/GHz < 20) frequency region by using open-ended coaxial probe technique. From the reflection coefficient and impedance data, the real and imaginary part of the dielectric permittivity values is determined. It is observed that there is a decrease in the real part of the dielectric permittivity up to certain frequency and an increase in the imaginary part of the dielectric permittivity with increase in the molar concentration of L-alanine in water medium. Based on the experimental data the average relaxation time values are calculated and its behavior is analyzed in terms of bound water and free water molecules. The theoretical dipole moment of L-alanine is calculated at gaseous state as well as in aqueous medium by using PCM and IEFPCM model at HF, DFT/B3LYP and MP2 calculations using 6-311G* basis set. Analysis between experimentally determined parameters and computed dipole moments were discussed. The mean molecular polarizability is calculated from the Lippincott δ function potential model and compared with the Le Fèvre method of polarizability values.

Molecular Dynamics Study of the Effect of l-Alanine Chiral Dopants on Diluted Chromonic Solutions

Atomistic molecular dynamics simulations have been performed for disodium cromoglycate (DSCG) chromonic solutions mixed with l-alanine chiral dopants. We study the fundamental molecular mechanisms induced by low concentrations of l-alanine on diluted DSCG solutions, including their effect on the chromonic aggregates, the solvent, and sodium counterions. Simulations reveal that l-alanine molecules primarily interact with DSCG stacks establishing salt bridges between their respective ammonium and carboxylate groups. Our results demonstrate that l-alanine and sodium counterions jointly establish an intricate network of noncovalent interactions around DSCG aggregates that decreases the global electrostatic repulsion of the chromonic system. Two possible structural effects in DSCG aggregates arise from this electronic stabilization: the increment of the total number of consecutively stacked aromatic planes per DSCG aggregate (intracolumnar effect) or the partial separation reduction between neighboring DSCG columnar sections due to the simultaneous bridging of intercolumnar DSCG carboxylate sites by sodium counterions, forming sodium bridges (intercolumnar effect). Sodium bridges may be responsible for the formation of stacking faults in DSCG aggregates in the form of lateral overlap junctions. This mechanism would explain the difference between lower X-ray correlation lengths with the expected persistence length in chromonics.

Fluorescence Spectral Studies on the Interaction of Alanine and Valine with Resorcinol-Based Acridinedione Dyes in Aqueous Solution: A Comparative Study with Glycine

Photophysical studies were carried out for simple amino acids like alanine and valine with resorcinol-based aqueous acridinedione (ADDR) dyes. ADDR dyes exhibit interesting excited-state characteristics on altering the substituents at the 9th and 10th sites (Scheme 1). The longest-wavelength absorption maxima remain the same on adding the amino acids to the fluorophore, whereas the excited-state behavior varies significantly mostly based on the nature of the substituent at the 9th position. The absence of fluorescence enhancement was observed with addition of β-alanine, l-alanine, and l-valine to ADDR1 dye (photoinduced electron transfer, PET), whereas addition of glycine exhibits enhancement accompanied with a shift toward a longer-wavelength region. Interestingly, the addition of amino acids to non-PET dyes results in a fluorescence quenching accompanied with a larger shift toward the shorter-wavelength region. The properties of fluorophore and nonfluorophore dyes in the presence of alanine or valine are found to be entirely different from those of glycine. The interaction of alanine with ADDR dyes is predominantly through H-bonding, but the structural aspects of H-bonding interactions of alanine and water are completely different from those of glycine and water. The time-correlated single-photon counting method portrays the existence of fluorophore in two distinguishable microenvironments in the presence of amino acids. The fluorescence spectral technique used as a tool in elucidating the mode of interaction of dye with neutral amino acids in aqueous solution is illustrated in the present study.

Toward theoretical terahertz spectroscopy of glassy aqueous solutions: partially frozen solute-solvent couplings of glycine in water

The molecular-level understanding of THz spectra of aqueous solutions under ambient conditions has been greatly advanced in recent years. Here, we go beyond previous analyses by performing ab initio molecular dynamics simulations of glycine in water with artificially frozen solute or solvent molecules, respectively, while computing the total THz response as well as its decomposition into mode-specific resonances based on the "supermolecular solvation complex" technique. Clamping the water molecules and keeping glycine moving breaks the coupling of glycine to the structural dynamics of the solvent, however, the polarization and dielectric solvation effects in the static solvation cage are still at work since the full electronic structure of the quenched solvent is taken into account. The complementary approach of fixing glycine reveals both the dynamical and electronic response of the solvation cage at the level of its THz response. Moreover, to quantitatively account for the electronic contribution solely due to solvent embedding, the solute species is "vertically desolvated", thus preserving the fully coupled solute-solvent motion in terms of the solute's structural dynamics in solution, while its electronic structure is no longer subject to solute-solvent polarization and charge transfer effects. When referenced to the free simulation of Gly(aq), this three-fold approach allows us to decompose the THz spectral contributions due to the correlated solute-solvent dynamics into entirely structural and purely electronic effects. Beyond providing hitherto unknown insights, the observed systematic changes of THz spectra in terms of peak shifts and lineshape modulations due to conformational freezing and frozen solvation cages might be useful to investigate the solvation of molecules in highly viscous H-bonding solvents such as ionic liquids and even in cryogenic ices as relevant to polar stratospheric and dark interstellar clouds.

Influence of Solvent Representation on Nuclear Shielding Calculations of Protonation States of Small Biological Molecules

In this study, we assess the influence of solvation on the accuracy and reliability of isotropic nuclear magnetic shielding calculations for amino acids in comparison to experimental data. We focus particularly on the performance of solvation methods for different protonation states, as biological molecules occur almost exclusively in aqueous solution and are subject to protonation with pH. We identify significant shortcomings of current implicit solvent models and present a hybrid solvation approach that improves agreement with experimental data by taking into account the presence of direct interactions between amino acid protonation state and water molecules.

Surface Tension and Adsorption Studies by Drop Profile Analysis Tensiometry

Surface tension and dilational viscoelasticity of solutions of various surfactants measured with bubble and drop profile analysis tensiometry are discussed. The study also includes experiments on the co-adsorption of surfactant molecules from a solution drop and alkane molecules from saturated alkane vapor phase. Using experimental data for 12 surfactants with different surface activities, it is shown that depletion due to adsorption of surfactant from the drop bulk can be significant. An algorithm is proposed quantitatively to take into consideration the depletion effect which is required for a correct description of the co-adsorption of alkanes on the solution drop surface and the correct analysis of experimental dynamic surface tension data to determine the adsorption mechanism. Bubble and drop profile analysis tensiometry is also the method of choice for measuring the dilational viscoelasticity of the adsorbed interfacial layer. The same elasticity moduli are obtained with the bubble and drop method only when the equilibrium surface pressures are sufficiently small (Π < 15 mN m-1). When the surface pressure for a surfactant solution is larger than this value, the viscoelasticity moduli determined from drop profile experiments become significantly larger than those obtained from bubble profile measurements.

Dynamic interfacial tension of surfactant solutions

The dynamics of surfactant interfacial layers was first discussed more than a century ago. In 1946 the most important work by Ward and Tordai was published which is still the theoretical basis of all new models to describe the time dependence of interfacial properties. In addition to the diffusion controlled adsorption mechanism, many other models have been postulated in literature, however, well performed experiments with well defined surfactant systems have shown that the diffusional transport is the main process governing the entire formation of surfactant adsorption layers. The main prerequisite, in addition to the diffusional transport, is the consideration of the right boundary condition at the interface, given by a respective equation of state. In addition to the classical models of Langmuir and Frumkin, also the so-called reorientation or interfacial aggregation models are to be assumed to reach a quantitative description of respective experimental data. Moreover, the adsorption of surfactants at the interface between water and a gas phase different from air can be strongly influenced by the type of molecules within the gas phase, such as alkane vapours. These oil molecules co-adsorb from the gas phase and change the adsorption kinetics strongly. Besides the discussion of how to apply theoretical adsorption kinetics models correctly, a large number of experimental data are presented and the way of a quantitative analysis of the adsorption mechanism and the main characteristic parameters is presented. This includes micellar solutions as well as mixtures of surfactants of ionic and non-ionic nature.

Conformational analysis of short polar side-chain amino-acids through umbrella sampling and DFT calculations

Molecular and quantum mechanics calculations were carried out in a series of tripeptides (GXG, where X = D, N and C) as models of the unfolded states of proteins. The selected central amino acids, especially aspartic acid (D) and asparagine (N) are known to present significant average conformations in partially allowed areas of the Ramachandran plot, which have been suggested to be important in unfolded protein regions. In this report, we present the calculation of the propensity values through an umbrella sampling procedure in combination with the calculation of the NMR J-coupling constants obtained by a DFT model. The experimental NMR observations can be reasonably explained in terms of a conformational distribution where PP_II and β basins sum up propensities above 0.9. The conformational analysis of the side chain dihedral angle (χ₁), along with the computation of ³J(H^αH^β), revealed a preference for the g _- and g ₊ rotamers. These may be connected with the presence of intermolecular H-bonding and carbonyl-carbonyl interactions sampled in the PP_II and β basins. Taking into account all those results, it can be established that these residues show a similar behavior to other amino acids in short peptides regarding backbone φ,ψ dihedral angle distribution, in agreement with some experimental analysis of capped dipeptides.

Electronic structures of intermolecular hydrogen bond contacts with solute in aqueous solution: glycine as a working prototype

The intermolecular hydrogen bond (H-bond) interactions play vital roles in many biological systems. Despite continued interest, the nature of their electronic structures has remained elusive. Based on the unique features of aqueous solution, a simple model depicting the H-bond electronic states by orbital hybridizations is developed. The model is demonstrated by reproducing the experimental IR data and yielding favorable solute-solvent interactions for the prototype glycine. The H-bond state for solute H, O and N atoms is found to be characterized by sp(1), sp(2), and sp(3) hybridizations, respectively. The model provides a new way for probing the intricate solute-solvent contacts.

Different proton transfer channels for the transformation of zwitterionic alanine-(H₂O)(n=2-4) to nonzwitterionic alanine-(H₂O)(n=2-4): a density functional theory study

We report here the various possibilities of proton transfer between the zwitterionic and the non-zwetterionic form of alanine (Ala) via (H₂O)(n=2-4) clusters by calculating the transition state structures of zwitterionic alanine (ZAla)-(H₂O)(n=2-4) and non-zwitterionic alanine (Ala)-(H₂O)(n=2-4) complexes at B3LYP/6-311++G(d,p) and CAM-B3LYP/6-311++G(d,p) level of theory. In order to determine the most feasible channel for proton transfer, the barrier energy corresponding to each channel was calculated. For the transformation of ZAla-(H₂O)(n=2) to Ala-(H₂O)(n=2), we identified eight channels for proton transfer. The lowest barrier energy (2.57 kcal mol⁻¹) channel, where ZAla-(H₂O)(n=2) transforms to Ala-(H₂O)(n=2) via triple proton transfer, is said to be the energetically most feasible channel. The values of barrier energy corresponding to the least energy pathway for proton transfer were calculated to be 1.14 and 9.82 kcal mol⁻¹ for n = 3 and n = 4 complexes, respectively, at B3LYP/6-311++G(d,p) level of theory. For complex n = 3, the structure where proton transfer takes place directly from -NH₃⁺ to -COO⁻ has the lowest energy pathway. However, the complexes for n = 2 and 3--the channels where proton transferred from -NH₃⁺ to -COO⁻ via two water molecules have the lowest barrier energy. For each n, the values of barrier energy calculated at CAM-B3LYP/6-311++G(d,p) level of theory were always less compared those calculated at B3LYP/6-311++G(d,p) level of theory. The value of rate constants corresponding to each proton transfer channel was also calculated.

The role of water-based hydrogen atom wires in long-range electron-transfer reactions in aqueous media for the FeII-FeIII self-exchange and related systems

Extended reach: A calculated mechanism for long-range proton-coupled electron transfer (PCET, see picture) through an array of structured water molecules between Fe(III)-Fe(III) complexes accounts for the reaction enthalpy and kinetic isotope effect previously measured for this reaction. This mechanism may be general and occur for other hydroxo-metal complexes.

13C NMR chemical shift calculations of charged surfactants in water — A combined density functional theory (DFT) and molecular dynamics (MD) methodological study

Structural and magnetic properties of one anionic and one cationic amphiphile molecule (sodium octanoate and hexadecyltrimethylammonium chloride, respectively) in water are studied comparing different methods to account for the presence of the solvent. Calculated 13C NMR chemical shifts are used as the probe for accuracy of the theoretical electronic structures obtained with different descriptions of the surfactants in water solution. The best agreement with the experimental data are obtained by averaging 13C NMR isotropic chemical shifts over a large number of conformational structures of sodium octanoate while considering the electronic structure of the solvent molecules. The 13C chemical shift values of the hexadecyltrimethylammonium alkane chain are systematically overestimated by 10–15 ppm even if an extensive conformational sampling and water as the polarized continuum medium have been taken into consideration. The role of the basis set quality has been studied and discussed as well.

QM and QM/MM simulations of proteins

Molecular dynamics simulations of biomolecules have matured into powerful tools of structural biology. In addition to the commonly used empirical force field potentials, quantum mechanical descriptions are gaining popularity for structure optimization and dynamic simulations of peptides and proteins. In this chapter, we introduce methodological developments such as the QM/MM framework and linear-scaling QM that make efficient calculations on large biomolecules possible. We identify the most common scenarios in which quantum descriptions of peptides and proteins are employed, such as structural refinement, force field development, treatment of unusual residues, and predicting spectroscopic and exited state properties. The benefits and shortcomings of QM potentials, in comparison to classical force fields, are discussed, with special emphasis on the sampling problems of protein conformational space. Finally, recent examples of QM/MM calculations in light-sensitive membrane proteins illustrate typical applications of the reviewed methods.

Presence or absence of a novel charge-transfer complex in the base-catalyzed hydrolysis of N-ethylbenzamide or ethyl benzoate

Reaction paths of base-catalyzed hydrolyses of isoelectronic substrates, Ph–C(=O)–X–Et [X = O (ethyl benzoate) and X = NH (N-ethylbenzamide)], were traced by DFT calculations. To simulate bond interchanges accompanied by proton transfers, a cluster model of Ph–C(=O)–X–Et + OH⁻(H₂O)₁₆ was employed. For X = O, three elementary processes and for X = NH four ones were obtained. The rate-determining step of X = O is the first TS (TS1, the OH⁻ addition step), while that of X = NH is TS2. TS2 of X = NH leads to a novel Mulliken charge-transfer complex, Ph–(OH)(O=)C∙∙∙N(H₂)–Et. The superiority or inferiority between the direct nucleophilic process or the general base-catalyzed process for TS1 was examined with the model Ph–C(=O)–X–Et + OH⁻(H₂O)_n, n = 3, 5, 8, 12, 16, 24 and 32. The latter process was calculated to be more favorable regardless of the number (n, except n = 3) of water molecules. The counter ion Na⁺ works unfavorably on the ester hydrolysis, particularly on TS1. A minimal model of TS1 was proposed and was found to be insensitive to n.

Gas phase basicities of polyfunctional molecules. Part 3: Amino acids

The present article is the third part of a general overview of the gas-phase protonation thermochemistry of polyfunctional molecules (first part: Mass Spectrom. Rev., 2007, 26:775-835, second part: Mass Spectrom. Rev., 2011, in press). This review is devoted to the 20 proteinogenic amino acids and is divided in two parts. In the first one, the experimental data obtained during the last 30 years using the equilibrium, thermokinetic and kinetic methods are presented. A general re-assignment of the values originating from these various experiments has been done on the basis of the commonly accepted Hunter & Lias 1998 gas-phase basicity scale in order to provide an homogeneous set of data. In the second part, theoretical investigations on gaseous neutral and protonated amino acids are reviewed. Conformational landscapes of both types of species were examined in order to provide theoretical protonation thermochemistry based on the truly identified most stable conformers. Proton affinities computed at the presently highest levels of theory (i.e. composite methods such as Gn procedures) are presented. Estimates of thermochemical parameters calculated using a Boltzmann distribution of conformers at 298K are also included. Finally, comparison between experiment and theory is discussed and a set of evaluated proton affinities, gas-phase basicities and protonation entropies is proposed.

Assessment of the intrinsic conformational preferences of dipeptide amino acids in aqueous solution by combined umbrella sampling/MBAR statistics. A comparison with experimental results

The propensities of 19 amino acid dipeptides have been calculated by a distributed umbrella sampling molecular dynamics simulation procedure using the OPLS-AA force field. The potential of mean force maps was estimated with the multiple Bennett acceptance ratio statistics. The resulting propensities compare satisfactorily well with very recently published experimental data on equivalent systems. In particular, α conformation-probabilities for all of the dipeptides remain much lower than either β or P(II) propensities. This result is in agreement with most experimental data for dipeptides. However, it is also in contrast with most simulation studies performed so far with other force fields, where α conformations result even more probable than P(II) or β ones. We discuss the behavior of the OPLS-AA force field, which can be useful for the improvement of this model in reproducing the recent experimental observations on amino acid dipeptides.

Theoretical study on conformational preferences of ribose in 2-thiouridine--the role of the 2'OH group

Conformational changes in ribose are well-known to play a significant role in biomolecular identification. The mechanism of selectivity towards C3'-endo conformation (conformer b) in ribose of 2-thiouridine has been studied using DFT (B3LYP) and MP2 methodology, together with 6-31+G(d,p) basis set. The polarity of the C2S2 bond is enhanced due to the orientation of H2' towards the S2 atoms, which leads to a difference in the corresponding bond lengths, the atomic charges and the vO2'H2' stretch vibrations in all the conformers. NBO analysis shows that charge transfer mainly occurs in the C2N3 and C2S2 orbitals. The higher stability of conformer b is attributed to its larger orbital interaction energies within the 2-thiouracil base, and total orbital interaction energies of conformer b. Our conclusion is that the distant electrostatic rather than hydrogen bonding effects between 2'OH and the S2 atoms play the dominant role in the orbital interaction, and enhance the selectivity towards the C3'-endo conformation of ribose.

A new model for mapping the peptide backbone: predicting proton chemical shifts in proteins

This paper describes a methodology that correlates experimental chemical shifts (at the alpha proton) of proteins with their geometrical data (both dihedral angles and distances) obtained from 13 representative proteins, which are taken from the Protein Data Bank (PDB) and the BioMagRes Data Bank (BMRB). To this end, the experimentally measured proton chemical shifts of simple amides have been correlated with DFT-based calculated structures, at the B3PW91/6-31G* level. This results in a series of mathematical relationships, which are extrapolated to the above-mentioned proteins giving rise to a modified equation for such skeleta. It is relevant to note that the equation is also supported by a clear comparison with NMR data of a protein beyond the chosen set, such as insulin, even with lower errors. The model also relates the dependence of chemical shifts on hydrophobic and anisotropic effects at the amino acid residues.