Krystel El Hage

A Toolkit to Fit Nonbonded Parameters from and for Condensed Phase Simulations

The quality of atomistic simulations depends decisively on the accuracy of the underlying energy function (force field). Of particular importance for condensed-phase properties are nonbonded interactions, including the electrostatic and Lennard-Jones terms. Permanent atomic multipoles (MTPs) are an extension to common point-charge (PC) representations in atomistic simulations. MTPs are commonly determined from and fitted to an ab initio Electrostatic Potential (ESP), and Lennard-Jones (LJ) parameters are obtained from comparison of experimental and computed observables using molecular dynamics (MD) simulations. For this a set of thermodynamic observables such as density, heat of vaporization, and hydration free energy is chosen, to which the parametrization is fitted. The current work introduces a comprehensive computing environment (Fitting Wizard (FW)) for optimizing nonbonded interactions for atomistic force fields of different qualities. The FW supports fitting of standard PC-based force fields and more physically motivated multipolar (MTP) force fields. A broader study including 20 molecules ranging from N-methyl-acetamide and benzene to halogenated benzenes, phenols, anilines, and pyridines yields a root mean squared deviation for hydration free energies of 0.36 kcal/mol over a range of 8 kcal/mol. It is furthermore shown that PC-based force fields are not necessarily inferior compared to MTP parametrizations depending on the molecule considered.

Extending Halogen-based Medicinal Chemistry to Proteins: IODO-INSULIN AS A CASE STUDY

Insulin, a protein critical for metabolic homeostasis, provides a classical model for protein design with application to human health. Recent efforts to improve its pharmaceutical formulation demonstrated that iodination of a conserved tyrosine (TyrB26) enhances key properties of a rapid-acting clinical analog. Moreover, the broad utility of halogens in medicinal chemistry has motivated the use of hybrid quantum- and molecular-mechanical methods to study proteins. Here, we (i) undertook quantitative atomistic simulations of 3-[iodo-TyrB26]insulin to predict its structural features, and (ii) tested these predictions by X-ray crystallography. Using an electrostatic model of the modified aromatic ring based on quantum chemistry, the calculations suggested that the analog, as a dimer and hexamer, exhibits subtle differences in aromatic-aromatic interactions at the dimer interface. Aromatic rings (TyrB16, PheB24, PheB25, 3-I-TyrB26, and their symmetry-related mates) at this interface adjust to enable packing of the hydrophobic iodine atoms within the core of each monomer. Strikingly, these features were observed in the crystal structure of a 3-[iodo-TyrB26]insulin analog (determined as an R6 zinc hexamer). Given that residues B24–B30 detach from the core on receptor binding, the environment of 3-I-TyrB26 in a receptor complex must differ from that in the free hormone. Based on the recent structure of a “micro-receptor” complex, we predict that 3-I-TyrB26 engages the receptor via directional halogen bonding and halogen-directed hydrogen bonding as follows: favorable electrostatic interactions exploiting, respectively, the halogens electron-deficient σ-hole and electronegative equatorial band. Inspired by quantum chemistry and molecular dynamics, such “halogen engineering” promises to extend principles of medicinal chemistry to proteins.

Correction to A Toolkit to Fit Nonbonded Parameters from and for Condensed Phase Simulations

J. Chem. Inf. Model. 2016, 56 (8), 1479−1489. DOI: 10.1021/acs.jcim.6b00280 I a recent publication we presented a fitting environment for parametrizing point charge (PC) and multipolar (MTP) force fields for condensed-phase simulations. After publication of this work it came to our attention that one of the scripts contained an error which caused an energy component in the free energy simulations to return incorrect values. This affects the optimization of the parameter when scaling the LennardJones (LJ) parameters according to ε* = ε and Rmin * /2 = Rmin/2 but not the MTP terms. Hence, all compounds considered were reparametrized according to the procedure described in ref 1. The corresponding correlations between experiment and the optimized parametrizations are reported in Figures 1 and 2.

Sequential Proton Coupled Electron Transfer (PCET): Dynamics Observed over 8 Orders of Magnitude in Time

Charge transfer mechanisms lay at the heart of chemistry and biochemistry. Proton coupled electron transfers (PCET) are central in biological processes such as photosynthesis and in the respiratory chain, where they mediate long-range charge transfers. These mechanisms are normally difficult to harness experimentally due to the intrinsic complexity of the associated biological systems. Metal-peptide cations experience both electron and proton transfers upon photoexcitation, proving an amenable model system to study PCET. We report on a time-resolved experiment designed to follow this dual charge transfer kinetics in [HG3W+Ag](+) (H = histidine, G = glycine, W = tryptophan) on time scales ranging from femtoseconds to milliseconds. While electron transfer completes in less than 4 ps, it triggers a proton transfer lasting over hundreds of microseconds. Molecular dynamics simulations show that conformational dynamic plays an important role in slowing down this reaction. This combined experimental and computational approach provides a view of PCET as a single phenomenon despite its very wide time-domain span.

Implications of short time scale dynamics on long time processes

This review provides a comprehensive overview of the structural dynamics in topical gas- and condensed-phase systems on multiple length and time scales. Starting from vibrationally induced dissociation of small molecules in the gas phase, the question of vibrational and internal energy redistribution through conformational dynamics is further developed by considering coupled electron/proton transfer in a model peptide over many orders of magnitude. The influence of the surrounding solvent is probed for electron transfer to the solvent in hydrated I−. Next, the dynamics of a modified PDZ domain over many time scales is analyzed following activation of a photoswitch. The hydration dynamics around halogenated amino acid side chains and their structural dynamics in proteins are relevant for iodinated TyrB26 insulin. Binding of nitric oxide to myoglobin is a process for which experimental and computational analyses have converged to a common view which connects rebinding time scales and the underlying dynamics. Finally, rhodopsin is a paradigmatic system for multiple length- and time-scale processes for which experimental and computational methods provide valuable insights into the functional dynamics. The systems discussed here highlight that for a comprehensive understanding of how structure, flexibility, energetics, and dynamics contribute to functional dynamics, experimental studies in multiple wavelength regions and computational studies including quantum, classical, and more coarse grained levels are required.

Molecular Mechanisms Underlying Solute Retention at Heterogeneous Interfaces

Despite considerable effort, a molecular-level understanding of the mechanisms governing adsorption/desorption in reversed-phase liquid chromatography is still lacking. This impedes rational design of columns and the development of reliable, computationally more efficient approaches to predict the selectivity of a particular column design. Using state-of-the art, validated force fields and free-energy simulations, the adsorption thermodynamics of benzene derivatives is investigated in atomistic detail and provides a quantitative microscopic understanding of retention when compared with experimental data. It is found that pure partitioning or pure adsorption is rather the exception than the rule. Typically, a pronounced ∼1 kcal/mol stabilization on the surface is accompanied by a broad trough indicative of partitioning before the probe molecule incorporates into the mobile phase. The present findings provide a quantitative and rational basis to develop improved effective, coarse-grained computational models and to design columns for specific applications.

Free energy simulations for protein ligand binding and stability

Abstract We summarize several computational techniques to determine relative free energies for condensed-phase systems. The focus is on practical considerations which are capable of making direct contact with experiments. Particular applications include the thermodynamic stability of apo- and holo-myoglobin, insulin dimerization free energy, ligand binding in lysozyme, and ligand diffusion in globular proteins. In addition to provide differential free energies between neighboring states, converged umbrella sampling simulations provide insight into migration barriers and ligand dissociation barriers and analysis of the trajectories yield additional insight into the structural dynamics of fundamental processes. Also, such simulations are useful tools to quantify relative stability changes for situations where experiments are difficult. This is illustrated for NO-bound myoglobin. For the dissociation of benzonitrile from lysozyme it is found that long umbrella sampling simulations are required to approximately converge the free energy profile. Then, however, the resulting differential free energy between the bound and unbound state is in good agreement with estimates from molecular mechanics with generalized Born surface area simulations. Furthermore, comparing the barrier height for ligand escape suggests that ligand dissociation contains a non-equilibrium component.

Valid molecular dynamics simulations of human hemoglobin require a surprisingly large box size

Recent molecular dynamics (MD) simulations of human hemoglobin (Hb) give results in disagreement with experiment. Although it is known that the unliganded (T0) and liganded (R4) tetramers are stable in solution, the published MD simulations of T0 undergo a rapid quaternary transition to an R-like structure. We show that T0 is stable only when the periodic solvent box contains ten times more water molecules than the standard size for such simulations. The results suggest that such a large box is required for the hydrophobic effect, which stabilizes the T0 tetramer, to be manifested. Even in the largest box, T0 is not stable unless His146 is protonated, providing an atomistic validation of the Perutz model. The possibility that extra large boxes are required to obtain meaningful results will have to be considered in evaluating existing and future simulations of a wide range of systems.

Publisher's Note: “Implications of short time scale dynamics on long time processes” (Struct. Dyn. 4, 061507 (2017)]

[This corrects the article DOI: 10.1063/1.4996448.].

The Role of Water in the Stability of Wild-type and Mutant Insulin Dimers

Insulin dimerization and aggregation play important roles in the endogenous delivery of the hormone. One of the important residues at the insulin dimer interface is PheB24, which is an invariant aromatic anchor that packs toward its own monomer inside a hydrophobic cavity formed by ValB12, LeuB15, TyrB16, CysB19, and TyrB26. Using molecular dynamics and free-energy simulations within explicit solvent, the structural and dynamical consequences of mutations of Phe at position B24 to glycine (Gly), alanine (Ala), and d-Ala and the des-PheB25 variant are quantified. Consistent with experiments, it is found that the Gly and Ala modifications lead to insulin dimers with reduced stability by 4 and 5 kcal/mol from thermodynamic integration and 4 and 8 kcal/mol from results using molecular mechanics-generalized Born surface area, respectively. Given the experimental difficulties to quantify the thermodynamic stability of modified insulin dimers, such computations provide a valuable complement. Interestingly, the Gly mutant exists as a strongly and a weakly interacting dimer. Analysis of the molecular dynamics simulations shows that this can be explained by water molecules that replace direct monomer-monomer H-bonding contacts at the dimerization interface involving residues B24 to B26. It is concluded that such solvent molecules play an essential role and must be included in future insulin dimerization studies.