Structure optimization via free energy gradient method: Application to glycine zwitterion in aqueous solution

Calculation of the geometry, absorption spectrum, and first hyperpolarizability of 4,5-dicyanoimidazole derivatives in solution. A multiscale ASEC-FEG study

We investigate the effects of solvents on the geometry, absorption spectrum, and first hyperpolarizability of six push-pull molecules, each containing a 4,5-dicyanoimidazole group as an electron acceptor and a N,N-dimethylamino group as an electron donor, with systematically extended π-conjugated systems. Geometry optimizations in dichloromethane, methanol, water, and formamide under normal thermodynamic conditions were performed using the average solvent electrostatic configuration-free energy gradient method, which employs a discrete solvent model. The conformational structure of molecules is moderately affected by the environment, with the π-conjugated system becoming more planar in protic solvents. Solvent effects on the first hyperpolarizability result in marked increases that are in line with the red shifts of the absorption spectrum. The hyperpolarizability of smaller molecules within the set may be significantly influenced by the effects of geometry relaxation in highly polar protic solvents. The results illustrate the role of hydrogen bonds in the structure and electronic properties of push-pull molecules in protic environments. For smaller molecules, hydrogen bonds significantly contribute to enhancing the hyperpolarizability, but the effect of these specific interactions becomes less significant with the length of the π-conjugated system.

Cartesian constraints in QM/MM optimizations

With the rise of quantum mechanical/molecular mechanical (QM/MM) methods, the interest in the calculation of molecular assemblies has increased considerably. The structures and dynamics of such assemblies are usually governed to a large extend by intermolecular interactions. As a result, the corresponding potential energy surfaces are topological rich and possess many shallow minima. Therefore, local structure optimizations of QM/MM molecular assemblies can be challenging, in particular if optimization constraints are imposed. To overcome this problem, structure optimization in normal coordinate space is advocated. To do so, the external degrees of freedom of a molecule are separated from the internal ones by a projector matrix in the space of the Cartesian coordinates. Here we extend this approach to Cartesian constraints. To this end, we devise an algorithm that adds the Cartesian constraints directly to the projector matrix and in this way eliminates them from the reduced coordinate space in which the molecule is optimized. To analyze the performance and stability of the constrained optimization algorithm in normal coordinate space, we present constrained minimizations of small molecular systems and amino acids in gas phase as well as water employing QM/MM constrained optimizations. All calculations are performed in the framework of auxiliary density functional theory as implemented in the program deMon2k.

Alternative Strategy for Spectral Tuning of Flavin-Binding Fluorescent Proteins

iLOV is an engineered flavin-binding fluorescent protein (FbFP) with applications for in vivo cellular imaging. To expand the range of applications of FbFPs for multicolor imaging and FRET-based biosensing, it is desirable to understand how to modify their absorption and emission wavelengths (i.e., through spectral tuning). There is particular interest in developing FbFPs that absorb and emit light at longer wavelengths, which has proven challenging thus far. Existing spectral tuning strategies that do not involve chemical modification of the flavin cofactor have focused on placing positively charged amino acids near flavin's C4a and N5 atoms. Guided by previously reported electrostatic spectral tunning maps (ESTMs) of the flavin cofactor and by quantum mechanical/molecular mechanical (QM/MM) calculations reported in this work, we suggest an alternative strategy: placing a negatively charged amino acid near flavin's N1 atom. We predict that a single-point mutant, iLOV-Q430E, has a slightly red-shifted absorption and fluorescence maximum wavelength relative to iLOV. To validate our theoretical prediction, we experimentally expressed and purified iLOV-Q430E and measured its spectral properties. We found that the Q430E mutation results in a slight change in absorption and a 4-8 nm red shift in the fluorescence relative to iLOV, in good agreement with the computational predictions. Molecular dynamics simulations showed that the carboxylate side chain of the glutamate in iLOV-Q430E points away from the flavin cofactor, which leads to a future expectation that further red shifting may be achieved by bringing the side chain closer to the cofactor.

Ionic Atmosphere Effect on the Absorption Spectrum of a Flavoprotein: A Reminder to Consider Solution Ions

This study utilizes the FMN-dependent NADH:quinone oxidoreductase from Pseudomonas aeruginosa PAO1 to investigate the effect of introducing an active site negative charge on the flavin absorption spectrum both in the absence and presence of a long-range electrostatic potential coming from solution ions. There were no observed changes in the flavin UV-visible spectrum when an active site tyrosine (Y277) becomes deprotonated in vitro. These results could only be reproduced computationally using average solvent electrostatic configuration (ASEC) QM/MM simulations that include both positive and negative solution ions. The same calculations performed with minimal ions to neutralize the total protein charge predicted that deprotonating Y277 would significantly alter the flavin absorption spectrum. Analyzing the distribution of solution ions indicated that the ions reorganize around the protein surface upon Y277 deprotonation to cancel the effect of the tyrosinate on the flavin absorption spectrum. Additional biochemical experiments were performed to test this hypothesis.

Free Energy Computation for an Isomerizing Chromophore in a Molecular Cavity via the Average Solvent Electrostatic Configuration Model: Applications in Rhodopsin and Rhodopsin-Mimicking Systems

We present a novel technique for computing the free energy differences between two chromophore "isomers" hosted in a molecular environment (a generalized solvent). Such an environment may range from a relatively rigid protein cavity to a flexible solvent environment. The technique is characterized by the application of the previously reported "average electrostatic solvent configuration" method, and it is based on the idea of using the free energy perturbation theory along with a chromophore annihilation procedure in thermodynamic cycle calculations. The method is benchmarked by computing the ground-state room-temperature relative stabilities between (i) the cis and trans isomers of prototypal animal and microbial rhodopsins and (ii) the analogue isomers of a rhodopsin-like light-driven molecular switch in methanol. Furthermore, we show that the same technology can be used to estimate the activation free energy for the thermal isomerization of systems i-ii by replacing one isomer with a transition state. The results show that the computed relative stability and isomerization barrier magnitudes for the selected systems are in line with the available experimental observation in spite of their widely diverse complexity.

Applicability of DFT functionals for evaluating the first hyperpolarizability of phenol blue in solution

The first electronic hyperpolarizability (β) of phenol blue (PB) in several solvents in a wide range of dielectric constants is investigated using the density functional theory (DFT). The reliability of various exchange-correlation functionals is assessed by a comparison to reference Møller-Plesset second-order perturbation theory (MP2) calculations. The equilibrium geometry of PB in each solvent is obtained by using the average solvent electrostatic configuration/free energy gradient method, which performs optimizations on the free energy hyper-surface by employing iteratively the sequential quantum mechanics/molecular mechanics methodology. The dependence of β on the bond length alternation (BLA) coordinate is rationalized by means of the two-level model. Within the employed exchange-correlation functionals, the LC-BLYP functional shows the best performance for describing the static and dynamic MP2 results of β, which increases as the BLA diminishes, reaching a maximum in an intermediate value of BLA. The results also illustrate the role played by the difference between the ground- and excited-state dipole moments (Δμ) in determining the hyperpolarizability behavior in solution. Particularly, in the aqueous solution case, Δμ goes to around zero when BLA is near zero, leading to an abrupt decline in the β value. The DFT results of this study, therefore, indicate a clear relationship between the first hyperpolarizability and the BLA coordinate for the PB in solution, in agreement with experiment.

Ab Initio Evaluation of the Redox Potential of Cytochrome c

Various biochemical activities of metabolism and biosynthesis are fulfilled by redox processes with explicit electron exchange, which furnish redox enzymes with high chemical reactivity. However, theoretical investigation of a redox process, which simultaneously involves a complex electronic change at a redox metal center and conformational reorganization of the surrounding protein environment coupled to the electronic change, requires computationally conflicting approaches, highly accurate quantum chemical calculations, and long-time molecular dynamics (MD) simulations, limiting the physicochemical understanding of biological redox processes. Here, we theoretically examined a redox process of cytochrome c by means of a hybrid molecular simulation technique, which enables one to consistently treat the redox center at the ab initio quantum chemistry level of theory and the protein reorganization with long-time MD simulations on the microsecond timescale. The calculations successfully evaluated a large absolute redox potential, 4.34 eV, with errors of only 0.03 to 0.34 eV to the experimental ones without any problem-specific empirical parameters. Through the long-time MD sampling, large and nonlinear reorganization of the protein environment was unveiled and the molecular determinants for the redox potential were identified. The present ab initio approach significantly expands the applicability of theoretical investigation to biological redox systems with more electronically complicated redox centers such as polynuclear transition metal complexes.

A simplified charge projection scheme for long-range electrostatics in ab initio QM/MM calculations

In a previous work [Pan et al., Molecules 23, 2500 (2018)], a charge projection scheme was reported, where outer molecular mechanical (MM) charges [>10 Å from the quantum mechanical (QM) region] were projected onto the electrostatic potential (ESP) grid of the QM region to accurately and efficiently capture long-range electrostatics in ab initio QM/MM calculations. Here, a further simplification to the model is proposed, where the outer MM charges are projected onto inner MM atom positions (instead of ESP grid positions). This enables a representation of the long-range MM electrostatic potential via augmentary charges (AC) on inner MM atoms. Combined with the long-range electrostatic correction function from Cisneros et al. [J. Chem. Phys. 143, 044103 (2015)] to smoothly switch between inner and outer MM regions, this new QM/MM-AC electrostatic model yields accurate and continuous ab initio QM/MM electrostatic energies with a 10 Å cutoff between inner and outer MM regions. This model enables efficient QM/MM cluster calculations with a large number of MM atoms as well as QM/MM calculations with periodic boundary conditions.

Assessing the structure and first hyperpolarizability of Li@B10H14 in solution: a sequential QM/MM study using the ASEC-FEG method

The structure and electronic properties of the lithium decahydroborate (Li@B₁₀H₁₄) complex in chloroform and water in normal thermodynamic conditions have been investigated using sequential QM/MM calculations by means of the average solvent electrostatic configuration (ASEC) and the Free Energy Gradient (FEG) methods. To obtain the structure of the Li@B₁₀H₁₄ complex in each of the solvents considered, we have performed geometry optimizations in solution using the ASEC-FEG method. The results show, for the first time with a realistic model of the molecular environment, that this alkali-metal-borane cluster is stable in chloroform but unstable in water. We have also explored the role of the electronic polarization of the solute due to solvent in the static first hyperpolarizability. The results show that, despite the reduction due to the effect of electrostatic polarization in chloroform, the Li@B₁₀H₁₄ complex still exhibits a large electronic first hyperpolarizability, with potential for application as a second-order nonlinear optical (NLO) material. In water, in contrast, the contribution of the excess electron for NLO responses is significantly affected by the electrostatic polarization effects. Therefore our results reveal that the influence of the environment must be considered in the design of new stable NLO materials.

Quantum Mechanical and Molecular Mechanics Modeling of Membrane-Embedded Rhodopsins

Computational chemistry provides versatile methods for studying the properties and functioning of biological systems at different levels of precision and at different time scales. The aim of this article is to review the computational methodologies that are applicable to rhodopsins as archetypes for photoactive membrane proteins that are of great importance both in nature and in modern technologies. For each class of computational techniques, from methods that use quantum mechanics for simulating rhodopsin photophysics to less-accurate coarse-grained methodologies used for long-scale protein dynamics, we consider possible applications and the main directions for improvement.

Microscopic Origin of Different Hydration Patterns of para-Nitrophenol and Its Anion: A Study Combining Multiconfigurational Calculations and the Free-Energy Gradient Method

A theoretical study of the solvatochromic shifts of para-nitrophenol ( pNP) and para-nitrophenolate anion ( pNP^-) in aqueous solution is presented using a QM/MM methodology with molecular dynamics simulation. The optimized structures in aqueous solution are obtained using both the polarizable continuum and the free-energy gradient methods. For pNP, the calculated redshifts at the CASPT2 (12,10) level are, respectively, 0.71 and 0.94 eV, in good agreement with the experimental ones (0.80-0.83 eV), whereas for pNP^-, they are small. The difference between the solvatochromic shifts of pNP and pNP^- is calculated as 0.71 eV in good agreement with the experimental one (0.79-0.81 eV). Finally, these shifts are understood in terms of the solvent effect on the solute structure, accurately calculated by the present theoretical treatment.

How Methylation Modifies the Photophysics of the Native All- trans-Retinal Protonated Schiff Base: A CASPT2/MD Study in Gas Phase and in Methanol

A comparison between the free-energy surfaces of the all- trans-retinal protonated Schiff base (RPSB) and its 10-methylated derivative in gas phase and methanol solution is performed at CASSCF//CASSCF and CASPT2//CASSCF levels. Solvent effects were included using the average solvent electrostatic potential from molecular dynamics method. This is a QM/MM (quantum mechanics/molecular mechanics) method that makes use of the mean field approximation. It is found that the methyl group bonded to C10 produces noticeable changes in the solution free-energy profile of the S₁ excited state, mainly in the relative stability of the minimum energy conical intersections (MECIs) with respect to the Franck-Condon (FC) point. The conical intersections yielding the 9- cis and 11- cis isomers are stabilized while that yielding the 13- cis isomer is destabilized; in fact, it becomes inaccessible by excitation to S₁. Furthermore, the planar S₁ minimum is not present in the methylated compound. The solvent notably stabilizes the S₂ excited state at the FC geometry. Therefore, if the S₂ state has an effect on the photoisomerization dynamics, it must be because it permits the RPSB population to branch around the FC point. All these changes combine to speed up the photoisomerization in the 10-methylated compound with respect to the native compound.

Zwitterionization of glycine in water environment: Stabilization mechanism and NMR spectral signatures

At physiological conditions, myriads of biomolecules (e.g., amino acids, peptides, and proteins) exist predominantly in the zwitterionic structural form and their biological functions will result in these conditions. However these geometrical structures are inaccessible energetically in the gas phase, and at this point, stabilization of amino-acids in physiological conditions is still under debate. In this paper, the electronic properties of a glycine molecule in the liquid environment were studied by performing a relaxation of the glycine geometry in liquid water using the free energy gradient method combined with a sequential quantum mechanics/molecular mechanics approach. A series of Monte Carlo Metropolis simulations of the glycine molecule embedded in liquid water, followed by only a quantum mechanical calculation in each of them were carried out. Both the local and global liquid environments were emphasized to obtain nuclear magnetic resonance (NMR) parameters for the glycine molecule in liquid water. The results of the equilibrium structure in solution and the systematic study of the hydrogen bonds were used to discard the direct proton transfer from the carboxyl group to the ammonium group of the glycine molecule in water solution. The calculations of the Density Functional Theory (DFT) were performed to study the polarization of the solvent in the parameters of nuclear magnetic resonance of the glycine molecule in liquid water. DFT calculations predicted isotropic chemical changes on the H, C, N, and O atoms of glycine in liquid water solution which agree with the available experimental data.

QM/MM Study of Substituent and Solvent Effects on the Excited State Dynamics of the Photoactive Yellow Protein Chromophore

Substituent and solvent effects on the excited state dynamics of the Photoactive Yellow Protein chromophore are studied using the average solvent electrostatic potential from molecular dynamics (ASEP/MD) method. Four molecular models were considered: the ester and thioester derivatives of the p-coumaric acid anion and their methylated derivatives. We found that the solvent produces dramatic modifications on the free energy profile of the S1 state: 1) Two twisted structures that are minima in the gas phase could not be located in aqueous solution. 2) Conical intersections (CIs) associated with the rotation of the single bond adjacent to the phenyl group are found for the four derivatives in water solution but only for thio derivatives in the gas phase. 3) The relative stability of minima and CIs is reverted with respect to the gas phase values, affecting the prevalent de-excitation paths. As a consequence of these changes, three competitive de-excitation channels are open in aqueous solution: the fluorescence emission from a planar minimum on S1, the trans-cis photoisomerization through a CI that involves the rotation of the vinyl double bond, and the nonradiative, nonreactive, de-excitation through the CI associated with the rotation of the single bond adjacent to the phenyl group. In the gas phase, the minima are the structures with the lower energy, while in solution these are the conical intersections. In solution, the de-excitation prevalent path seems to be the photoisomerization for oxo compounds, while thio compounds return to the initial trans ground state without emission.

Efficient Computational Research Protocol to Survey Free Energy Surface for Solution Chemical Reaction in the QM/MM Framework: The FEG-ER Methodology and Its Application to Isomerization Reaction of Glycine in Aqueous Solution

In solution chemical reaction, we often need to consider a multidimensional free energy (FE) surface (FES) which is analogous to a Born-Oppenheimer potential energy surface. To survey the FES, an efficient computational research protocol is proposed within the QM/MM framework; (i) we first obtain some stable states (or transition states) involved by optimizing their structures on the FES, in a stepwise fashion, finally using the free energy gradient (FEG) method, and then (ii) we directly obtain the FE differences among any arbitrary states on the FES, efficiently by employing the QM/MM method with energy representation (ER), i.e., the QM/MM-ER method. To validate the calculation accuracy and efficiency, we applied the above FEG-ER methodology to a typical isomerization reaction of glycine in aqueous solution, and reproduced quite satisfactorily the experimental value of the reaction FE. Further, it was found that the structural relaxation of the solute in the QM/MM force field is not negligible to estimate correctly the FES. We believe that the present research protocol should become prevailing as one computational strategy and will play promising and important roles in solution chemistry toward solution reaction ergodography.

Solvent effects on de-excitation channels in the p-coumaric acid methyl ester anion, an analogue of the photoactive yellow protein (PYP) chromophore

Environmental effects on the deactivation channels of the PYP chromophore in the gas phase and water solution are compared at the CASPT2//CASSCF/cc-pVDZ level.

A new QM/MM method oriented to the study of ionic liquids

The interest on room temperature ionic liquids has grown in the last decades because of their use as all-purpose solvent and their low environmental impact. In the present work, a new theoretical procedure is developed to study pure ionic liquids within the framework of the quantum mechanics/molecular mechanics method. Each type of ion (cation or anion) is considered as an independent entity quantum mechanically described that follows a differentiated path in the liquid. The method permits, through an iterative procedure, the full coupling between the polarized charge distribution of the ions and the liquid structure around them. The procedure has been tested with 1-ethyl-3-methylimidazolium tetrafluoroborate. It was found that, similar to non-polar liquids and as a consequence of the low value of the reaction field, the cation and anion charge distributions are hardly polarized by the rest of molecules in the liquid. Their structure is characterized by an alternance between anion and cation shells as evidenced by the coincidence of the first maximum of the anion-anion and cation-cation radial distribution functions with the first minimum of the anion-cation. Some degree of stacking between the cations is also found.

Correcting for the free energy costs of bond or angle constraints in molecular dynamics simulations

BACKGROUND

Free energy simulations are an important tool in the arsenal of computational biophysics, allowing the calculation of thermodynamic properties of binding or enzymatic reactions. This paper introduces methods to increase the accuracy and precision of free energy calculations by calculating the free energy costs of constraints during post-processing. The primary purpose of employing constraints for these free energy methods is to increase the phase space overlap between ensembles, which is required for accuracy and convergence.

METHODS

The free energy costs of applying or removing constraints are calculated as additional explicit steps in the free energy cycle. The new techniques focus on hard degrees of freedom and use both gradients and Hessian estimation. Enthalpy, vibrational entropy, and Jacobian free energy terms are considered.

RESULTS

We demonstrate the utility of this method with simple classical systems involving harmonic and anharmonic oscillators, four-atomic benchmark systems, an alchemical mutation of ethane to methanol, and free energy simulations between alanine and serine. The errors for the analytical test cases are all below 0.0007kcal/mol, and the accuracy of the free energy results of ethane to methanol is improved from 0.15 to 0.04kcal/mol. For the alanine to serine case, the phase space overlaps of the unconstrained simulations range between 0.15 and 0.9%. The introduction of constraints increases the overlap up to 2.05%. On average, the overlap increases by 94% relative to the unconstrained value and precision is doubled.

CONCLUSIONS

The approach reduces errors arising from constraints by about an order of magnitude. Free energy simulations benefit from the use of constraints through enhanced convergence and higher precision.

GENERAL SIGNIFICANCE

The primary utility of this approach is to calculate free energies for systems with disparate energy surfaces and bonded terms, especially in multi-scale molecular mechanics/quantum mechanics simulations. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Computing the Free Energy along a Reaction Coordinate Using Rigid Body Dynamics

The calculations of potential of mean force along complex chemical reactions or rare events pathways are of great interest because of their importance for many areas in chemistry, molecular biology, and material science. The major difficulty for free energy calculations comes from the great computational cost for adequate sampling of the system in high-energy regions, especially close to the reaction transition state. Here, we present a method, called FEG-RBD, in which the free energy gradients were obtained from rigid body dynamics simulations. Then the free energy gradients were integrated along a reference reaction pathway to calculate free energy profiles. In a given system, the reaction coordinates defining a subset of atoms (e.g., a solute, or the quantum mechanics (QM) region of a quantum mechanics/molecular mechanics simulation) are selected to form a rigid body during the simulation. The first-order derivatives (gradients) of the free energy with respect to the reaction coordinates are obtained through the integration of constraint forces within the rigid body. Each structure along the reference reaction path is separately subjected to such a rigid body simulation. The individual free energy gradients are integrated along the reference pathway to obtain the free energy profile. Test cases provided demonstrate both the strengths and weaknesses of the FEG-RBD method. The most significant benefit of this method comes from the fast convergence rate of the free energy gradient using rigid-body constraints instead of restraints. A correction to the free energy due to approximate relaxation of the rigid-body constraint is estimated and discussed. A comparison with umbrella sampling using a simple test case revealed the improved sampling efficiency of FEG-RBD by a factor of 4 on average. The enhanced efficiency makes this method effective for calculating the free energy of complex chemical reactions when the reaction coordinate can be unambiguously defined by a small subset of atoms within the system.

Dual Approach to Vibrational Spectra in Solution: Microscopic Influence of Hydrogen Bonding to the State of Motion of Glycine in Water

We have proposed a new theoretical methodology to clarify the microscopic nature of the vibrational properties in solution, which consists of a combination of the vibrational frequency analyses (VFAs) with two kinds of Hessian matrices, that is, the effective Hessian on the free energy surface (free energy Hessian: "FE-Hessian") and the instantaneous one (instantaneous normal mode Hessian: "INM-Hessian") within QM/MM framework. In these VFAs, the Hessians were obtained by the analytical approach, having the advantages from the aspect of both the computational efficiency and accuracy in comparison to those obtained by the numerical one. In the present study, we have applied them to the glycine aqueous solution. First, by using the VFA with the FE-Hessian (VFA-FEH), we estimated the vibrational frequency shifts induced by solvent water molecules. The calculated values were quantitatively in agreement with experimental ones. It was clearly demonstrated that such vibrational shifts are attributed to not only the structural relaxation but also the explicit solute-solvent interactions (i.e., interatomic interactions). Second, by using the VFA with the INM-Hessian (VFA-INMH), the vibrational spectra in solution were investigated through the vibrational INM densities of states (DOS). By the comparison between the spectroscopic features and the microscopic solvation structure around glycine molecule, it was found that the frequency shifts and bandwidths in IR spectra are closely correlated with the hydrogen bonding (HB) network formations. In particular, the instantaneous changes of vibrational states of the hydroxyl group and carbonyl one, showing apparently inverse tendency on the strength of the HB interaction, can be explained very well on the basis of two different mechanisms, that is, the direct change of electron density in the bonding orbitals and the indirect one due to hyperconjugation between the lone electron pair and the antibonding orbitals, respectively. In conclusion, the present dual VFA approach is a quite useful strategy to interpret the microscopic origin of the experimental vibrational spectra.

A SELF-CONSISTENT-CHARGE DENSITY-FUNCTIONAL TIGHT-BINDING THEORY BASED MOLECULAR DYNAMICS SIMULATION OF A ZWITTERIONIC GLYCINE IN AN EXPLICIT WATER ENVIRONMENT

A zwitterionic glycine (zGLY) is adopted as an example to study the impact of water environment (310 H2O molecules) on the molecular structure and energetics using a self-consistent-charge density-functional tight-binding theory based molecular dynamics (SCC-DFTB/MD) method. It is found that maximal eight hydrogen bonds could be formed simultaneously between eight water molecules and the zGLY. The ability of the COO- terminal to adsorb water molecules is stronger than the [Formula: see text] terminal with respect to hydrogen bonding with more water molecules and exhibits lower adiabatic adsorption energies. The zGLY's intramolecular hydrogen bond appeared unpredictably, without involving any proton transfer and generally helpful for enhancing the system stability. Water molecules play an important role to stabilize the zwitterionic amino acids and restrain the proton migration from the [Formula: see text] to the COO− group. Our results show that the SCC-DFTB/MD method could successfully describe geometry dynamical evolutions and energetics of biomolecules in a nanoscale simulation with the presence of a large number of water molecules. Our study not only verified the feasibility of a QM level methodology for describing the aqueous states of biochemical molecules, but also gave a clear evidence for the impact of water environment on amino acids.