Describing the Conformational Landscape of Small Organic Molecules through Gaussian Mixtures in Dihedral Space

Directed kinetic transition network model

Molecular dynamics simulations contain detailed kinetic information related to the functional states of proteins and macromolecules, but this information is obscured by the high dimensionality of configurational space. Markov state models and transition network models are widely applied to extract kinetic descriptors from equilibrium molecular dynamics simulations. In this study, we developed the Directed Kinetic Transition Network (DKTN)-a graph representation of a master equation which is appropriate for describing nonequilibrium kinetics. DKTN models the transition rate matrix among different states under detailed balance. Adopting the mixing time from the Markov chain, we use the half mixing time as the criterion to identify critical state transition regarding the protein conformational change. The similarity between the master equation and the Kolmogorov equation suggests that the DKTN model can be reformulated into the continuous-time Markov chain model, which is a general case of the Markov chain without a specific lag time. We selected a photo-sensitive protein, vivid, as a model system to illustrate the usage of the DKTN model. Overall, the DKTN model provides a graph representation of the master equation based on chemical kinetics to model the protein conformational change without the underlying assumption of the Markovian property.

Determination of the conformational states of strychnine in solution using NMR residual dipolar couplings in a tensor-free approach

Small molecules with rotatable bonds can occupy different conformational states in solution as a consequence of their thermal fluctuations. The accurate determination of the structures of such states, as well as of their statistical weights, has been challenging because of the technical difficulties in extracting information from experimental measurements, which are normally averaged over the conformational space available. Here, to achieve this objective, we present an approach based on a recently proposed tensor-free method for incorporating NMR residual dipolar couplings as structural restraints in replica-averaged molecular dynamics simulations. This approach enables the information provided by the experimental data to be used in the spirit of the maximum entropy principle to determine the structural ensembles of small molecules. Furthermore, in order to enhance the sampling of the conformational space we incorporated the metadynamics method in the simulations. We illustrate the method in the case of strychnine, determining the three major conformational states of this small molecule and their associated occupation probabilities.

Conformational Preference and Spectroscopical Characteristics of the Active Pharmaceutical Ingredient Levetiracetam

The analysis of the possible conformers and the conformational change between solid and liquid states of a particular drug molecule are mandatory not only for describing reliably its spectroscopical properties but also for understanding the interaction with the receptor and its mechanism of action. Therefore, here we investigated the free-energy conformational landscape of levetiracetam (LEV) in gas phase as well as in water and ethanol, aiming to describe the 3-dimensional structure and energetic stability of its conformers. Twenty-two unique conformers were identified, and their energetic stability was determined at density functional theory B3LYP/6-31+G(2d,2p) level of theory. The 6 most stable monomers in water, within a relative free-energy window of 0.71 kcal mol^-1 and clearly separated in energy from the remaining subset of 16 conformers, as well as the 3 most stable dimers were then used to compute the Boltzmann populations-averaged UV-Vis and NMR spectra of LEV. The conformational landscape in solution is distinctly different from that corresponding to gas phase, particularly due to the relative orientations of the butanamide group. Aiming to clarify the stability of the possible dimers of LEV, we also investigated computationally the structure of a set of 11 nonhydrated and hydrated homochiral hydrogen-bonded LEV dimers.

Assessing the suitability of the multilevel strategy for the conformational analysis of small ligands

Predicting the conformational preferences of flexible compounds is a challenging problem in drug design, where the recognition between ligand and receptor is affected by the ability of the interacting partners to adopt a favorable conformation for the binding. To explore the conformational space of flexible ligands and to obtain the relative free energy of the conformation wells, we have recently reported a multilevel computational strategy that relies on the predominant-state approximation-where the conformational space is partitioned into distinct conformational wells-and combines a low-level method for sampling the conformational minima and high-level ab initio calculations for estimating their relative stability. In this study, we assess the performance of the multilevel strategy for predicting the conformational preferences of a series of structurally related phenylethylamines and streptomycin in aqueous solution. The charged nature of these compounds and the chemical complexity of streptomycin make them a challenging test for the multilevel approach. Furthermore, we explore the suitability of using a molecular mechanics approach as a source of approximate ensembles in the first stage of the multilevel strategy. The results support the reliability of the multilevel approach for obtaining an accurate conformational ensemble of small (bio)organic molecules in aqueous solution.