FROG: Exploiting all-atom molecular dynamics trajectories to calculate linear and non-linear optical responses of molecular liquids within Dalton's QM/MM polarizable embedding scheme.

Quantum mechanical/molecular mechanics (QM/MM) methods are interesting to model the impact of a complex environment on the spectroscopic properties of a molecule. In this context, a FROm molecular dynamics to second harmonic Generation (FROG) code is a tool to exploit molecular dynamics trajectories to perform QM/MM calculations of molecular optical properties. FROG stands for "FROm molecular dynamics to second harmonic Generation" since it was developed for the calculations of hyperpolarizabilities. These are relevant to model non-linear optical intensities and compare them with those obtained from second harmonic scattering or second harmonic generation experiments. FROG's specificity is that it is designed to study simple molecular liquids, including solvents or mixtures, from the bulk to the surface. For the QM/MM calculations, FROG relies on the Dalton package: its electronic-structure models, response theory, and polarizable embedding schemes. FROG helps with the global workflow needed to deal with numerous QM/MM calculations: it permits the user to separate the system into QM and MM fragments, to write Dalton's inputs, to manage the submission of QM/MM calculations, to check whether Dalton's calculation finished successfully, and finally to perform averages on relevant QM observables. All molecules within the simulation box and several time steps are tackled within the same workflow. The platform is written in Python and installed as a package. Intermediate data such as local electric fields or individual molecular properties are accessible to the users in the form of Python object arrays. The resulting data are easily extracted, analyzed, and visualized using Python scripts that are provided in tutorials.

Inviting C5-Trifluoromethylated Pseudoprolines into Collagen Mimetic Peptides.

Numerous collagen mimetic peptides (CMPs) have been engineered using proline derivatives substituted at their C(3) and/or C(4) position in order to stabilize or functionalize collagen triple-helix mimics. However, no example has been reported so far with C(5) substitutions. Here, we introduce a fluorinated CMP incorporating trifluoromethyl groups at the C(5) position of pseudoproline residues. In tripeptide models, our CD, NMR, and molecular dynamics (MD) studies have shown that, when properly arranged, these residues meet the structural requirements for a triple-helix assembly. Two host-guest CMPs were synthesized and analyzed by CD spectroscopy. The NMR analysis in solution of the most stable confirmed the presence of structured homotrimers that we interpret as triple helices. MD calculations showed that the triple-helix model remained stable throughout the simulation with all six trifluoromethyl groups pointing outward from the triple helix. Pseudoprolines substituted at the C(5) positions appeared as valuable tools for the design of new fluorinated collagen mimetic peptides.

First hyperpolarizability of water in bulk liquid phase: long-range electrostatic effects included via the second hyperpolarizability.

The molecular first hyperpolarizability ß contributes to second-order optical non-linear signals collected from molecular liquids. For the Second Harmonic Generation (SHG) response, the first hyperpolarizability ß(2ω, ω, ω) often depends on the molecular electrostatic environment. This is especially true for water, due to its large second hyperpolarizability γ(2ω, ω, ω,0). In this study we compute the electronic γ(2ω, ω, ω,0) and ß(2ω, ω, ω) for water molecules in liquid water using QM/MM calculations. The average value of γ(2ω, ω, ω,0) is smaller than the one for the gaz phase, and its standard deviation among the molecules is relatively small. In addition, we demonstrate that the average bulk second hyperpolarizability 〈γ(2ω, ω, ω,0)〉 can be used to describe the electrostatic effects of the distant neighborhood on the first hyperpolarizability ß(2ω, ω, ω). In comparison with more complex schemes to take into account long-range effects, the approximation is simple, and does not require any modifications of the QM/MM implementation. The long-range correction can be added explicitly, using an average value of γ for water in the condensed phase. It can also be easily added implicitly in QM/MM calculations through an additional embedding electric field, without the explicit calculation of γ.

First hyperpolarizability of water at the air-vapor interface: a QM/MM study questions standard experimental approximations.

Surface Second-Harmonic Generation (S-SHG) experiments provide a unique approach to probe interfaces. One important issue for S-SHG is how to interpret the S-SHG intensities at the molecular level. Established frameworks commonly assume that each molecule emits light according to an average molecular hyperpolarizability tensor ß(-2ω,ω,ω). However, for water molecules, this first hyperpolarizability is known to be extremely sensitive to their environment. We have investigated the molecular first hyperpolarizability of water molecules within the liquid-vapor interface, using a quantum description with explicit, inhomogeneous electrostatic embedding. The resulting average molecular first hyperpolarizability tensor depends on the distance relative to the interface, and it practically respects the Kleinman symmetry everywhere in the liquid. Within this numerical approach, based on the dipolar approximation, the water layer contributing to the Surface Second Harmonic Generation (S-SHG) intensity is less than a nanometer. The results reported here question standard interpretations based on a single, averaged hyperpolarizability for all molecules at the interface. Not only the molecular first hyperpolarizability tensor significantly depends on the distance relative to the interface, but it is also correlated to the molecular orientation. Such hyperpolarizability fluctuations may impact the S-SHG intensity emitted by an aqueous interface.

Molecular electrometer and binding of cations to phospholipid bilayers.

Despite the vast amount of experimental and theoretical studies on the binding affinity of cations - especially the biologically relevant Na+ and Ca2+ - for phospholipid bilayers, there is no consensus in the literature. Here we show that by interpreting changes in the choline headgroup order parameters according to the 'molecular electrometer' concept [Seelig et al., Biochemistry, 1987, 26, 7535], one can directly compare the ion binding affinities between simulations and experiments. Our findings strongly support the view that in contrast to Ca2+ and other multivalent ions, Na+ and other monovalent ions (except Li+) do not specifically bind to phosphatidylcholine lipid bilayers at sub-molar concentrations. However, the Na+ binding affinity was overestimated by several molecular dynamics simulation models, resulting in artificially positively charged bilayers and exaggerated structural effects in the lipid headgroups. While qualitatively correct headgroup order parameter response was observed with Ca2+ binding in all the tested models, no model had sufficient quantitative accuracy to interpret the Ca2+:lipid stoichiometry or the induced atomistic resolution structural changes. All scientific contributions to this open collaboration work were made publicly, using nmrlipids.blogspot.fi as the main communication platform.

Mixed Mechanism of Lubrication by Lipid Bilayer Stacks.

Although the key role of lipid bilayer stacks in biological lubrication is generally accepted, the mechanisms underlying their extreme efficiency remain elusive. In this article, we report molecular dynamics simulations of lipid bilayer stacks undergoing load and shear. When the hydration level is reduced, the velocity accommodation mechanism changes from viscous shear in hydration water to interlayer sliding in the bilayers. This enables stacks of hydrated lipid bilayers to act as efficient boundary lubricants for various hydration conditions, structures, and mechanical loads. We also propose an estimation for the friction coefficient; thanks to the strong hydration forces between lipid bilayers, the high local viscosity is not in contradiction with low friction coefficients.

Optical properties of prodigiosin and obatoclax: action spectroscopy and theoretical calculations.

Prodiginine molecules (prodigiosin and obatoclax) are well-known pH-chromic dyes with promising anti-tumor properties. They present multiple tautomeric and rotameric forms. The protonation state and the structure of such flexible ligands in interaction with a protein are crucial to understand and to model the protein's biological activities. The determination of the protonation state via UV/vis absorption is possible if the ligand spectra of the neutral and protonated states are sufficiently different, and also if we can eliminate other factors potentially impacting the spectrum. Upon measuring the absorption spectra of the ligand in solution, varying solvents and pH values, we have determined that the optical properties of prodigiosin and obatoclax depend on the protonation state and not on the solvent permittivity constant. In parallel, action spectroscopy (using tunable lasers coupled to ion traps) in the gas phase of protonated and sodiated prodigiosin and obatoclax molecules has been performed to evaluate the sensitivity of the charge and the conformational state to their optical properties free of solvent. The spectra are interpreted using computational simulations of molecular structures and electronic excitations. The excitation energies are only slightly sensitive to various isomerizations, and may be used to distinguish between protonated and deprotonated states, even in the presence of a sodium counter-ion.

Multi-scale modeling of mycosubtilin lipopeptides at the air/water interface: structure and optical second harmonic generation.

Monolayers of the lipopeptide mycosubtilin are studied at the air/water interface. Their structure is investigated using molecular dynamics simulations. All-atom models suggest that the lipopeptide is flexible and aggregates at the interface. To achieve simulation times of several microseconds, a coarse-grained (CG) model based on the MARTINI force field was also used. These CG simulations describe the formation of half-micelles at the interface for surface densities up to 1 lipopeptide per nm(2). In these aggregates, the tyrosine side chain orientation is found to be constrained: on average, its main axis, as defined along the C-OH bond, aligns along the interface normal and points towards the air side. The origin of the optical second harmonic generation (SHG) from mycosubtilin monolayers at the air/water interface is also investigated. The molecular hyperpolarizability of the lipopeptide is obtained from quantum chemistry calculations. The tyrosine side chain contribution to the hyperpolarizability is found to be dominant. The orientation distribution of tyrosine, associated with a dominant hyperpolarizability component along the C-OH bond of the tyrosine, yields a ratio of the susceptibility elements χ((2))(ZZZ)/χ((2))(ZXX) consistent with the experimental measurements recently reported by M. N. Nasir et al. [Phys. Chem. Chem. Phys., 2013, 15, 19919].

Interleaflet sliding in lipidic bilayers under shear flow: comparison of the gel and fluid phases using reversed non-equilibrium molecular dynamics simulations.

The friction between two rubbing surfaces lubricated by water can be diminished if they are coated with phospholipidic bilayers or brushes of polyelectrolytes. In the case of a coating by lipid membranes, the friction is lower when the lipids are in the gel phase rather than in the liquid phase. We investigated the response of fluid or gel bilayers to a mechanical load or under shear using non-equilibrium molecular dynamics simulations (NEMD) to understand whether this difference could come from intermonolayer sliding. The system is composed of a single fully hydrated bilayer of coarse grained phospholipids under a parallel shear with vorticity parallel to the bilayer. In both the liquid and the gel phases, an intermonolayer slip was measured in the velocity profile. In the liquid phase this slip is proportional to the shear stress. In the tilted gel phase of our model the stress is not systematically linear and relaxes differently when the shear is in the direction of the tilt or perpendicular to it. The impact of surface tension (or load) on the friction is different for the liquid and gel phases, but grossly the slip remains of the same order of magnitude.

Influence of the tyrosine environment on the second harmonic generation of iturinic antimicrobial lipopeptides at the air-water interface.

The second harmonic generation (SHG) response at the air-water interface from the tyrosine-containing natural iturinic cyclo-lipopeptides mycosubtilin, iturin A and bacillomycin D is reported. It is shown that this response is dominated by the single tyrosine residue present in these molecules owing to the large first hyperpolarizability arising from the non-centrosymmetric aromatic ring structure of this amino acid. The SHG response of these iturinic antibiotics is also compared to the response of surfactin, a cyclo-lipopeptide with a similar l,d-amino acid sequence but lacking a tyrosine residue, and PalmATA, a synthetic linear lipopeptide possessing a single tyrosine residue but lacking the amino acid sequence structuring the cycle of the iturinic antibiotics. From the light polarization analysis of the SHG response, it is shown that the tyrosine local environment is critical in defining the SHG response of these peptides at the air-water interface. Our results demonstrate that tyrosine, similar to tryptophan, can be used as an endogenous molecular probe of peptides and proteins for SHG at the air-water interface, paving the way for SHG studies of other tyrosine-containing bioactive molecules.

Liquid Water: When Hyperpolarizability Fluctuations Boost and Reshape the Second Harmonic Scattering Intensities.

Second harmonic scattering (SHS) is a method of choice to investigate the molecular structure of liquids. While a clear interpretation of SHS intensity exists for diluted solutions of dyes, the scattering due to solvents remains difficult to interpret quantitatively. Here, we report a quantum mechanics/molecular mechanics (QM/MM) approach to model the polarization-resolved SHS intensity of liquid water, quantifying different contributions to the signal. We point out that the molecular hyperpolarizability fluctuations and correlations cannot be neglected. The intermolecular orientational and hyperpolarizability correlations up to the third solvation layer strongly increase the scattering intensities and modulate the polarization-resolved oscillation that is predicted here by QM/MM without fitting parameters. Our approach can be generalized to other pure liquids to provide a quantitative interpretation of SHS intensities in terms of short-range molecular ordering.

Visualizing the transiently populated closed-state of human HSP90 ATP binding domain.

HSP90 are abundant molecular chaperones, assisting the folding of several hundred client proteins, including substrates involved in tumor growth or neurodegenerative diseases. A complex set of large ATP-driven structural changes occurs during HSP90 functional cycle. However, the existence of such structural rearrangements in apo HSP90 has remained unclear. Here, we identify a metastable excited state in the isolated human HSP90α ATP binding domain. We use solution NMR and mutagenesis to characterize structures of both ground and excited states. We demonstrate that in solution the HSP90α ATP binding domain transiently samples a functionally relevant ATP-lid closed state, distant by more than 30 Å from the ground state. NMR relaxation enables to derive information on the kinetics and thermodynamics of this interconversion, while molecular dynamics simulations establish that the ATP-lid in closed conformation is a metastable exited state. The precise description of the dynamics and structures sampled by human HSP90α ATP binding domain provides information for the future design of new therapeutic ligands.

Rheology of sliding leaflets in coarse-grained DSPC lipid bilayers.

Amphiphilic lipid bilayers modify the friction properties of the surfaces on top of which they are deposited. In particular, the measured sliding friction coefficient can be significantly reduced compared with the native surface. We investigate in this work the friction properties of a numerical coarse-grained model of DSPC (1,2-distearoyl-sn-glycero-3-phosphocholine) lipid bilayer subject to longitudinal shear. The interleaflet friction coefficient is obtained from out-of-equilibrium pulling or from relaxation simulations. In particular, we gain access to the transient viscoelastic response of a sheared bilayer. The bilayer mechanical response is found to depend significantly on the membrane physical state, with evidence in favor of a linear response regime in the fluid but not in the gel region. The linear response validity domain is established, and the timescales appearing in the membrane response discussed.

Coarse-Grain Simulations of Solid Supported Lipid Bilayers with Varying Hydration Levels.

Supported lipid bilayers (SLBs) are a very popular system for the study of biomimetic membranes. Understanding of the interactions between the solid substrate and the lipid membrane opens pathways to the design of new materials with fine-tunable properties. While it is possible to study SLBs via molecular dynamics (MD) simulations, difficulties still remain for these strategies; in particular, the confined water layer thickness and structure are difficult to reproduce in simulations. We have explored different coarse-grained (CG) models for the membrane/support interaction, and their impact on the substrate hydration level. Our results highlight the relevance of including long-range interactions in CG-MD simulations of fluid SLBs. Modeled neutron reflectivity curves are deduced from the structures obtained by molecular simulations, and substrate parameters are optimized to match the experimental and modeled reflectivity curves. We expect our coarse-grained approach to open new perspectives for the simulations of SLBs of increasing complexity, including lipid layers of complex compositions, or adsorbed lipidic layers on patterned surfaces.

Anomeric memory of the glycosidic bond upon fragmentation and its consequences for carbohydrate sequencing.

Deciphering the carbohydrate alphabet is problematic due to its unique complexity among biomolecules. Strikingly, routine sequencing technologies-which are available for proteins and DNA and have revolutionised biology-do not exist for carbohydrates. This lack of structural tools is identified as a crucial bottleneck, limiting the full development of glycosciences and their considerable potential impact for the society. In this context, establishing generic carbohydrate sequencing methods is both a major scientific challenge and a strategic priority. Here we show that a hybrid analytical approach integrating molecular spectroscopy with mass spectrometry provides an adequate metric to resolve carbohydrate isomerisms, i.e the monosaccharide content, anomeric configuration, regiochemistry and stereochemistry of the glycosidic linkage. On the basis of the spectroscopic discrimination of MS fragments, we report the unexpected demonstration of the anomeric memory of the glycosidic bond upon fragmentation. This remarkable property is applied to de novo sequencing of underivatized oligosaccharides.Establishing generic carbohydrate sequencing methods is both a major scientific challenge and a strategic priority. Here the authors show a hybrid analytical approach integrating molecular spectroscopy and mass spectrometry to resolve carbohydrate isomerism, anomeric configuration, regiochemistry and stereochemistry.

Toward Atomistic Resolution Structure of Phosphatidylcholine Headgroup and Glycerol Backbone at Different Ambient Conditions.

Phospholipids are essential building blocks of biological membranes. Despite a vast amount of very accurate experimental data, the atomistic resolution structures sampled by the glycerol backbone and choline headgroup in phoshatidylcholine bilayers are not known. Atomistic resolution molecular dynamics simulations have the potential to resolve the structures, and to give an arrestingly intuitive interpretation of the experimental data, but only if the simulations reproduce the data within experimental accuracy. In the present work, we simulated phosphatidylcholine (PC) lipid bilayers with 13 different atomistic models, and compared simulations with NMR experiments in terms of the highly structurally sensitive C-H bond vector order parameters. Focusing on the glycerol backbone and choline headgroups, we showed that the order parameter comparison can be used to judge the atomistic resolution structural accuracy of the models. Accurate models, in turn, allow molecular dynamics simulations to be used as an interpretation tool that translates these NMR data into a dynamic three-dimensional representation of biomolecules in biologically relevant conditions. In addition to lipid bilayers in fully hydrated conditions, we reviewed previous experimental data for dehydrated bilayers and cholesterol-containing bilayers, and interpreted them with simulations. Although none of the existing models reached experimental accuracy, by critically comparing them we were able to distill relevant chemical information: (1) increase of choline order parameters indicates the P-N vector tilting more parallel to the membrane, and (2) cholesterol induces only minor changes to the PC (glycerol backbone) structure. This work has been done as a fully open collaboration, using nmrlipids.blogspot.fi as a communication platform; all the scientific contributions were made publicly on this blog. During the open research process, the repository holding our simulation trajectories and files ( https://zenodo.org/collection/user-nmrlipids ) has become the most extensive publicly available collection of molecular dynamics simulation trajectories of lipid bilayers.

Optical properties of a visible push-pull chromophore covalently bound to carbohydrates: solution and gas-phase spectroscopy combined to theoretical investigations.

The use of visible absorbing and fluorescent tags for sensing and structural analysis of carbohydrates is a promising route in a variety of medical, diagnostic, and therapeutic contexts. Here we report an easy method for covalent attachment of nonfluorescent push-pull chromophores based on the 4-cyano-5-dicyanomethylene-2-oxo-3-pyrroline ring to carbohydrate moieties. The impact of sugar grafting on the optical properties of the push-pull chromophore in the gas phase and in solution was investigated by absorption and action spectroscopy and theoretical methods. The labeled sugars efficiently absorb photons in the visible range, as demonstrated by their intense photodissociation in a quadrupole ion trap. A strong blue shift (-70 nm) of the gas-phase photodissociation intensity maximum is observed upon sugar grafting, whereas no such effect is visible on the solution absorption spectra. Molecular dynamics simulations of labeled maltose in the gas phase describe strong interactions between the sulfonated chromophore and the carbohydrate, which lead to cyclic conformations. These are not observed in the simulations with explicit solvation. Time-dependent density functional theory (TD-DFT) calculations on model molecules permit us to attribute the observed shift to the formation of such cyclic conformations and to the displacement of the negative charge relative to the aromatic moiety of the chromophore.

Microsolvation effects on the optical properties of crystal violet.

We present a joint experimental and theoretical study of the photoabsorption and photodissociation behavior of crystal violet, that is, the tris[p-(dimethylamino)phenyl]methyl cation. The photodissociation spectra of isolated and microsolvated crystal violet have been measured. A single band is observed for the bare cation. This is in good agreement with the calculated vibronic absorption spectrum based on time-dependent density functional theory calculations. The interaction of crystal violet with a single water molecule shifts and broadens the photodissociation spectrum, so that it approaches the spectrum obtained in solution. Theoretical calculations of the structure of the complex suggest that the shift in the absorption spectrum originates from a water molecule bonding with the central carbon atom of crystal violet.