Probing Liquid-Ordered and Disordered Phases in Lipid Model Membranes: A Combined Theoretical and Spectroscopic Study of a Fluorescent Molecular Rotor.

An integrated theoretical/experimental strategy has been applied to the study of environmental effects on the spectroscopic parameters of 4-(diphenylamino)phtalonitrile (DPAP), a fluorescent molecular rotor. The computational part starts from the development of an effective force field for the first excited electronic state of DPAP and proceeds through molecular dynamics simulations in solvents of different polarities toward the evaluation of Stokes shifts by quantum mechanics/molecular mechanics (QM/MM) approaches. The trends of the computed results closely parallel the available experimental results thus giving confidence to the interpretation of new experimental studies of the photophysics of DPAP in lipid bilayers. In this context, results show unambiguously that both flexible dihedral angles and global rotations are significantly retarded in a cholesterol/DPPC lipid matrix with respect to the DOPC matrix, thus confirming the sensitivity of DPAP to probe different environments and, therefore, its applicability as a probe for detecting different structures and levels of plasma membrane organization.

Development of Nonbonded Models for Metal Cations Using the Electronic Continuum Correction.

The parametrization of classical nonbonded models of metal ions has been widely addressed in the recent years. Despite the continuous development of novel and more physically inspired functional forms, the 12-6 Lennard-Jones plus Coulomb potential is still the most adopted force field in molecular dynamics (MD) codes, owing to its simple form and easy implementation. However, due to the integer formal charge, unpolarizable force fields of ions may suffer from overestimated interatomic electrostatic interactions, leading to nonphysical clustering or repulsion between such full charges. The electronic continuum correction (ECC) can fix this problem through a simple inclusion of solvent polarization effects via ionic charge rescaling. In this work, the development of novel nonbonded models for mono, divalent, and highly charged metal ions is presented. For each metal species, the ionic charge has been scaled, according to the ECC. Lennard-Jones parameters have been optimized using experimental structural and thermodynamic properties as target quantities. Performances of the proposed models are discussed and compared with the literature data, while transferability attitudes among different and well-known water models are evaluated. © 2019 Wiley Periodicals, Inc.

A combined experimental and theoretical study of optical rotatory dispersion for (R)-glycidyl methyl ether in aqueous solution.

The dispersive optical activity for aqueous solutions of non-rigid (R)-glycidyl methyl ether (R-GME) has been explored synergistically from experimental and theoretical perspectives. Density functional theory analyses performed with the polarizable continuum model for implicit solvation identified nine low-lying stable conformers that are interconverted by rotation about two large-amplitude torsional coordinates. The antagonistic chiroptical signatures predicted for these structural isomers were averaged under a Boltzmann-weighting ansatz to estimate the behavior expected for a thermally equilibrated ensemble. This led to optical rotatory dispersion profiles that reproduced the overall shape of observations but failed to achieve uniform agreement with measured specific-rotation values even when anharmonic vibrational corrections were applied. A mixed QM/FQ paradigm, whereby quantum-mechanical (QM) calculations of optical activity were combined with classical molecular dynamics simulations of explicit solvation that included mutual-polarization effects by means of fluctuating charges (FQ), was enlisted to elucidate the microsolvation environment and gauge its impact upon conformer distributions and response properties. Although quantitative accord with experiments remained elusive, this approach revealed strong variations in the magnitude and sign of rotatory powers for R-GME as the configuration of surrounding water molecules evolved, thereby highlighting the inherently dynamical nature of the solvated chiroptical response, calling into question the validity of "static" descriptions based on the presumption of distinct energy minima, and giving insight into the inherent complexity posed by the modeling of such properties for solvated systems.

Effective computational route towards vibrational optical activity spectra of chiral molecules in aqueous solution.

We present a computational methodology based on a polarizable Quantum Mechanical (QM)/Molecular Mechanics (MM) approach to accurately compute the Vibrational Optical Activity (VOA) spectra of chiral systems. This approach is applied for the calculation of Infrared (IR), Vibrational Circular Dichroism (VCD), Raman and Raman Optical Activity (ROA) spectra of aqueous solutions of (l)-methyl lactate and (S)-glycidol. Remarkable agreement between calculations and experiments is reported, showing the reliability and accuracy of the methodology, especially with respect to standard continuum solvation approaches.

Force Field Parametrization of Metal Ions from Statistical Learning Techniques.

A novel statistical procedure has been developed to optimize the parameters of nonbonded force fields of metal ions in soft matter. The criterion for the optimization is the minimization of the deviations from ab initio forces and energies calculated for model systems. The method exploits the combination of the linear ridge regression and the cross-validation techniques with the differential evolution algorithm. Wide freedom in the choice of the functional form of the force fields is allowed since both linear and nonlinear parameters can be optimized. In order to maximize the information content of the data employed in the fitting procedure, the composition of the training set is entrusted to a combinatorial optimization algorithm which maximizes the dissimilarity of the included instances. The methodology has been validated using the force field parametrization of five metal ions (Zn2+, Ni2+, Mg2+, Ca2+, and Na+) in water as test cases.

Computational study of the DPAP molecular rotor in various environments: from force field development to molecular dynamics simulations and spectroscopic calculations.

Fluorescent molecular rotors (FMRs) belong to an important class of environment-sensitive dyes capable of acting as nanoprobes in the measurement of viscosity and polarity of their micro-environment. FMRs have found widespread applications in various research fields, ranging from analytical to biochemical sciences, for example in intracellular imaging studies or in volatile organic compound detection. Here, a computational investigation of a recently proposed FMR, namely 4-(diphenylamino)phthalonitrile (DPAP), in various chemical environments is presented. A purposely developed molecular mechanics force field is proposed and then applied to simulate the rotor in a high- and low-polar solvent (i.e., acetonitrile, tetrahydrofuran, o-xylene and cyclohexane), a polymer matrix and a lipid membrane. Subtle effects of the molecular interactions with the embedding medium, the structural fluctuations of the rotor and its rotational dynamics are analyzed in some detail. The results correlate with a previous work, thus supporting the reliability of the model, and provide further insights into the environment-specific properties of the dye. In particular, it is shown how molecular diffusion and rotational correlation times of the FMR are affected by the surrounding medium and how the molecular orientation of the dye becomes anisotropic once immersed in the lipid bilayer. Moreover, a qualitative correlation between the FMR rotational dynamics and the fluorescence lifetime is detected, a result in line with the observed viscosity dependence of its emission. Finally, optical absorption spectra are computed and successfully compared with their experimental counterparts.

Immersive virtual reality in computational chemistry: Applications to the analysis of QM and MM data.

The role of Virtual Reality (VR) tools in molecular sciences is analyzed in this contribution through the presentation of the Caffeine software to the quantum chemistry community. Caffeine, developed at Scuola Normale Superiore, is specifically tailored for molecular representation and data visualization with VR systems, such as VR theaters and helmets. Usefulness and advantages that can be gained by exploiting VR are here reported, considering few examples specifically selected to illustrate different level of theory and molecular representation.

Tuning of dye optical properties by environmental effects: a QM/MM and experimental study.

The present work is aimed at a deeper investigation of two recently synthesized heteroaromatic fluorophores by means of a computational multilayer approach, integrating quantum mechanics (QM) and molecular mechanics (MM). In particular, dispersion of the title dyes in a polymer matrix is studied in connection with potential applications as photoactive species in luminescent solar concentrators (LSCs). Molecular dynamics simulations, based on accurate QM-derived force fields, reveal increased stiffness of these organic dyes when going from CHCl3 solution to the polymer matrix. QM/MM computations of UV spectra for snapshots extracted from MD simulations show that this different flexibility permits explaining the different spectral shapes obtained experimentally for the two different environments. Moreover, the general spectroscopic trends are reproduced well by static computations employing a polarizable continuum description of environmental effects.

Development and validation of a docking-based virtual screening platform for the identification of new lactate dehydrogenase inhibitors.

The human muscle isoform of lactate dehydrogenase (hLDH5) is one of the key enzymes of the glycolytic process. It is overexpressed in metastatic cancer cells and is linked to the vitality of tumors in hypoxic conditions. With the aim of identifying new hLDH5 inhibitors, a fully automated docking-based virtual screening platform was developed by considering different protein conformations and the consensus docking strategy. In order to verify the reliability of the reported platform, a small database of about 10,000 compounds was filtered by using this method, and the top-ranked compounds were tested for their hLDH5 inhibition activity. Enzymatic assays revealed that, among the ten selected compounds, two proved to efficiently inhibit enzyme activity with IC50 values in the micromolar range. These results demonstrate the validity of the methodologies we followed, encouraging the application of larger virtual screening studies and further refinements of the platform. Furthermore, the two active compounds herein described may be considered as interesting leads for the development of new and more efficient LDH inhibitors.