Accurate force fields are one of the major pillars on which successful molecular dynamics simulations of complex biomolecular processes rest. They have been optimized for ambient conditions, whereas high-pressure simulations become increasingly important in pressure perturbation studies, using pressure as an independent thermodynamic variable. Here, we explore the design of non-polarizable force fields tailored to work well in the realm of kilobar pressures--while avoiding complete reparameterization. Our key is to first compute the pressure-induced electronic and structural response of a solute by combining an integral equation approach to include pressure effects on solvent structure with a quantum-chemical treatment of the solute within the embedded cluster reference interaction site model (EC-RISM) framework. Next, the solute's response to compression is taken into account by introducing pressure-dependence into selected parameters of a well-established force field. In our proof-of-principle study, the full machinery is applied to N,N,N-trimethylamine-N-oxide (TMAO) in water being a potent osmolyte that counteracts pressure denaturation. EC-RISM theory is shown to describe well the charge redistribution upon compression of TMAO(aq) to 10 kbar, which is then embodied in force field molecular dynamics by pressure-dependent partial charges. The performance of the high pressure force field is assessed by comparing to experimental and ab initio molecular dynamics data. Beyond its broad usefulness for designing non-polarizable force fields for extreme thermodynamic conditions, a good description of the pressure-response of solutions is highly recommended when constructing and validating polarizable force fields.

Receptor tyrosine kinases represent one of the prime targets in cancer therapy, as the dysregulation of these elementary transducers of extracellular signals, like the epidermal growth factor receptor (EGFR), contributes to the onset of cancer, such as non-small cell lung cancer (NSCLC). Strong efforts were directed to the development of irreversible inhibitors and led to compound CO-1686, which takes advantage of increased residence time at EGFR by alkylating Cys797 and thereby preventing toxic effects. Here, we present a structure-based approach, rationalized by subsequent computational analysis of conformational ligand ensembles in solution, to design novel and irreversible EGFR inhibitors based on a screening hit that was identified in a phenotype screen of 80 NSCLC cell lines against approximately 1500 compounds. Using protein X-ray crystallography, we deciphered the binding mode in engineered cSrc (T338M/S345C), a validated model system for EGFR-T790M, which constituted the basis for further rational design approaches. Chemical synthesis led to further compound collections that revealed increased biochemical potency and, in part, selectivity toward mutated (L858R and L858R/T790M) vs nonmutated EGFR. Further cell-based and kinetic studies were performed to substantiate our initial findings. Utilizing proteolytic digestion and nano-LC-MS/MS analysis, we confirmed the alkylation of Cys797.

Biophysics under extreme conditions is the fundamental platform for scrutinizing life in unusual habitats, such as those in the deep sea or continental subsurfaces, but also for putative extraterrestrial organisms. Therefore, an important thermodynamic variable to explore is pressure. It is shown that the combination of infrared spectroscopy with simulation is an exquisite approach for unraveling the intricate pressure response of the solvation pattern of TMAO in water, which is expected to be transferable to biomolecules in their native solvent. Pressure-enhanced hydrogen bonding was found for TMAO in water. TMAO is a molecule known to stabilize proteins against pressure-induced denaturation in deep-sea organisms.

No abstract

High-pressure (HP) NMR spectroscopy is an important method for detecting rare functional states of proteins by analyzing the pressure response of chemical shifts. However, for the analysis of the shifts it is mandatory to understand the origin of the observed pressure dependence. Here we present experimental HP NMR data on the (15) N-enriched peptide bond model, N-methylacetamide (NMA), in water, combined with quantum-chemical computations of the magnetic parameters using a pressure-sensitive solvation model. Theoretical analysis of NMA and the experimentally used internal reference standard 4,4-dimethyl-4-silapentane-1-sulfonic (DSS) reveal that a substantial part of observed shifts can be attributed to purely solvent-induced electronic polarization of the backbone. DSS is only marginally responsive to pressure changes and is therefore a reliable sensor for variations in the local magnetic field caused by pressure-induced changes of the magnetic susceptibility of the solvent.