Victor Volkov

Electrostatic interactions in phospholipid membranes revealed by coherent 2D IR spectroscopy

The inter- and intramolecular interactions of the carbonyl moieties at the polar interface of a phospholipid membrane are probed by using nonlinear femtosecond infrared spectroscopy. Two-dimensional IR correlation spectra separate homogeneous and inhomogeneous broadenings and show a distinct cross-peak pattern controlled by electrostatic interactions. The inter- and intramolecular electrostatic interactions determine the inhomogeneous character of the optical response. Using molecular dynamics simulation and the nonlinear exciton equations approach, we extract from the spectra short-range structural correlations between carbonyls at the interface.

What are the Sites Water Occupies at the Interface of a Phospholipid Membrane

We explore the two-dimensional infrared response of D(2)O and of OD impurity at the interface of phospholipid membrane fragments. The spectra of the two aqueous states are inhomogeneously broadened due to the absorption of water molecule associated with the membrane interface in different structural sites. The nonlinear spectra of the two species allow disentangling of the spectral contributions of the aqueous states of two types: where the stretching modes are under effective mixing and where the stretching modes are uncoupled. By reviewing the results of the experimental studies with the support of molecular dynamic simulation we identify the spectral signatures of the main structural motives responsible for the inhomogeneous distribution of resonances in the infrared OD stretching region. Our analysis provides a quantitative estimate of the statistical population of the different aqueous species at the polar interface of the bilayer.

Structural Properties of a Membrane Associated Anchor Dipeptide

The association of peptides to phospholipid membranes through the insertion of an anchoring hydrocarbon tail is common to some viruses and to several anticancer drugs. We investigate the association of an anchor dipeptide, N-myristoylated methyl glycine (MrG), to phospholipid membrane fragments made of 1-palmitoyl-2-linoleyl phosphatidylcholine (PLPC). Here we report on the experimental findings of two-dimensional infrared spectroscopy of an MrG backbone in the 6 μm wavelength region. The experimental outcomes are supported by ab initio calculations and by a molecular dynamics simulation accomplished with the replica exchange method. We find that the guest molecule has a preferential unfolded conformation, with dihedral angles Φ = -90 ± 20° and Ψ = -180 ± 20°, while the average orientational distribution of the amide I transition dipole moments with respect to the neighbor PLPC carbonyls is peaked at angles in the range 21-33°. The depth of penetration of MrG inside the membrane corresponds rather well to the one estimated in our previous paper [J. Phys. Chem. B, 2009, 113, 16246], where we found that the backbone moieties of MrG are localized slightly above the carbonyl groups of PLPC. According to the simulation results, the anchor tail is completely inserted in the hydrophobic region of the bilayer, with a largely prevalent extended conformation and a preferential alignment along the average direction of the PLPC hydrocarbon tails.

Two-dimensional infrared spectroscopy of a structured liquid: Neat formamide

Vibrational dynamics of liquid formamide is studied in the spectral region of the amide I mode by means of linear and two-dimensional infrared spectroscopies. The two-dimensional spectrum has a complex structure to be connected to the partially excitonic nature of the vibrational states. The measurements performed on a 1:10 (12)C:(13)C formamide isotopic mixture allow separating the broadening contribution due to the inhomogeneous frequency distribution of the local oscillators from that of excitonic origin. A model based on the Kubo picture of the line broadening is used, together with the dynamic information obtained from a molecular dynamics simulation, to fit the spectra of the (12)C formamide impurity in the isotopic mixture. The relevant dynamical information, such as the amplitude of the frequency fluctuations, lifetime of the second vibrational excited state, and anharmonicity, is thus recovered. By appropriately combining the outcomes of experiments and molecular dynamics simulation, we demonstrate that motional narrowing determines the line shape of the amide I resonance to a large extent. The same analysis provides an estimate of the transition dipole moment of formamide, which results in good agreement with an ab initio calculation. The calculated frequency fluctuation correlation time is found to be comparable to the hydrogen-bond lifetime, which defines the basic structural relaxation rate of the networked liquid.

Retrieval of spectral and dynamic properties from two-dimensional infrared pump-probe experiments

We have developed a fitting algorithm able to extract spectral and dynamic properties of a three level oscillator from a two‐dimensional infrared spectrum (2D‐IR) detected in time resolved nonlinear experiments. Such properties go from the frequencies of the ground‐to‐first and first‐to‐second vibrational transitions (and hence anharmonicity) to the frequency‐fluctuation correlation function. This last is represented through a general expression that allows one to approach the various strategies of modeling proposed in the literature. The model is based on the Kubo picture of stochastic fluctuations of the transition frequency as a result of perturbations by a fluctuating surrounding. To account for the line‐shape broadening due to pump pulse spectral width in double‐resonance measurements, we supply the fitting algorithm with the option to perform the convolution of the spectral signal with a Lorentzian function in the pump‐frequency dimension. The algorithm is tested here on 2D‐IR pump‐probe spectra of a Gly‐Ala dipeptide recorded at various pump‐probe delay times. Speedup benchmarks have been performed on a small Beowulf cluster. The program is written in FORTRAN language for both serial and parallel architectures and is available free of charge to the interested reader.

Binding Free Energies of Host-Guest Systems by Nonequilibrium Alchemical Simulations with Constrained Dynamics: Theoretical Framework

The fast-switching decoupling method is a powerful nonequilibrium technique to compute absolute binding free energies of ligand-receptor complexes (Sandberg et al., J. Chem. Theory Comput. 2014, 11, 423-435). Inspired by the theory of noncovalent binding association of Gilson and co-workers (Biophys. J. 1997, 72, 1047-1069), we develop two approaches, termed binded-domain and single-point alchemical-path schemes (BiD-AP and SiP-AP), based on the possibility of performing alchemical trajectories during which the ligand is constrained to fixed positions relative to the receptor. The BiD-AP scheme exploits a recent generalization of nonequilibrium work theorems to estimate the free energy difference between the coupled and uncoupled states of the ligand-receptor complex. With respect to the fast-switching decoupling method without constraints, BiD-AP prevents the ligand from leaving the binding site, but still requires an estimate of the positional binding-site volume, which may not be a simple task. On the other side, the SiP-AP scheme allows avoidance of the calculation of the binding-site volume by introducing an additional equilibrium simulation of ligand and receptor in the bound state. In the companion article (DOI: 10.1021/acs.jctc.7b00595), we show that the extra computational effort required by SiP-AP leads to a significant improvement of accuracy in the free energy estimates.

Binding Free Energies of Host–Guest Systems by Nonequilibrium Alchemical Simulations with Constrained Dynamics: Illustrative Calculations and Numerical Validation

In the companion article (Giovannelli et al., 10.1021/acs.jctc.7b00594), we presented an alchemical approach, based on nonequilibrium molecular dynamics simulations, to compute absolute binding free energies of a generic host-guest system. Two alternative computational routes, called binded-domain and single-point alchemical-path schemes, have been proposed. This study is addressed to furnish numerical validation and illustrative examples of the above-mentioned alchemical schemes. Validation is provided by comparing binding free-energy data relative to two poses of a Zn(II)·anion complex with those recovered from an alternative approach, based on steered molecular dynamics simulations. We illustrate important technical and theoretical aspects for a good practice in applying both alchemical schemes, not only through the calculations on the Zn(II)·anion complex, but also estimating absolute binding free energies of 1:1 complexes of β-cyclodextrin with aromatic compounds (benzene and naphthalene). Comparison with experimental data and previous molecular dynamics simulation studies further confirms the validity of the present nonequilibrium-alchemical methodology.

Sum frequency generation image reconstruction: Aliphatic membrane under spherical cap geometry

The article explores an opportunity to approach structural properties of phospholipid membranes using Sum Frequency Generation microscopy. To establish the principles of sum frequency generation image reconstruction in such systems, at first approach, we may adopt an idealistic spherical cap uniform assembly of hydrocarbon molecules. Quantum mechanical studies for decanoic acid (used here as a representative molecular system) provide necessary information on transition dipole moments and Raman tensors of the normal modes specific to methyl terminal - a typical moiety in aliphatic (and phospholipid) membranes. Relative degree of localization and frequencies of the normal modes of methyl terminals make nonlinearities of this moiety to be promising in structural analysis using Sum Frequency Generation imaging. Accordingly, the article describes derivations of relevant macroscopic nonlinearities and suggests a mapping procedure to translate amplitudes of the nonlinearities onto microscopy image plane according to geometry of spherical assembly, local molecular orientation, and optical geometry. Reconstructed images indicate a possibility to extract local curvature of bilayer envelopes of spherical character. This may have practical implications for structural extractions in membrane systems of practical relevance.

Sum-frequency generation echo and grating from interface

The work addresses spectroscopy of fourth-order Sum Frequency Generation Echo and Grating responses as an experimental tool to study structure and dynamics at interfaces. First, it addresses experimental geometry to extract background-free fourth-order Echo and Grating responses. Further, the article provides the analytical expressions of the response functions for these nonlinearities. The derived expressions are used to model the χ((4)) two-dimensional spectral responses of a hydrated methyl acetate, which resembles a hydrated carbonyl moiety at the polar outer side of a phospholipid membrane. Orientation, transition dipole moments, and Raman tensors are obtained from the results of classical and quantum calculations, respectively. The numerical studies for the nonlinear responses under different polarization schemes and timings suggest the possibility of securely factoring of spectral contributions of χ(YYYZX) and χ(YYYZY) macroscopic susceptibilities. As such, the nonlinearities provide an experimental perspective on orientation of a generic (low-symmetry) molecular system at interfaces. Besides, the spectral properties of the tensors may reflect correlations of the in-plane and out-of-plane field components specific to the interface. For the case of a phospholipid membrane, the experiment would address in-plane and out-of-plane anisotropy of hydrogen bonding and related dynamics.

Partitioning of an anchor dipeptide in a phospholipid membrane.

We explore the localization of a guest N-myristoylated methyl glycine anchor dipeptide in a phospholipid environment. The dipeptide is part of a conservative sequence, which ensures proper association of a wide group of proteins in living organisms with a cellular phospholipid membrane. Using linear and two-color anharmonic infrared spectroscopy, we measure relative degrees of hydration of the amide I modes of the dipeptide and of phospholipid carbonyls. The atomic density of water in dependence of the distance from the hydrophobic center of the bilayer (a result of an independent Neutron scattering experiment) allows us to determine the relative altitudes of the peptide carbonyls with respect to those of the phospholipid ones. Considering this, and the dimensions of the dipeptide molecular frame, we anticipate the average angle between the backbone of the dipeptide and the normal to the membrane surface. The results provide a descriptive picture of the depth and geometry of partitioning of a guest N-myristoylated methyl glycine anchor dipeptide into a phospholipid membrane.