Explore the vast range of titles available.

Allosteric drugs have the potential to revolutionize biomedicine due to their enhanced selectivity and protection against overdosage. However, we need to better understand allosteric mechanisms in order to fully harness their potential in drug discovery. In this study, molecular dynamics simulations and nuclear magnetic resonance spectroscopy are used to investigate how increases in temperature affect allostery in imidazole glycerol phosphate synthase. Results demonstrate that temperature increase triggers a cascade of local amino acid-to-amino acid dynamics that remarkably resembles the allosteric activation that takes place upon effector binding. The differences in the allosteric response elicited by temperature increase as opposed to effector binding are conditional to the alterations of collective motions induced by either mode of activation. This work provides an atomistic picture of temperature-dependent allostery, which could be harnessed to more precisely control enzyme function. Using a combination of MD simulations and NMR, the authors investigate how temperature affects allostery in imidazole glycerol phosphate synthase (IGPS), revealing that increase of temperature triggers local amino acid dynamics and providing insights into mechanism of allosteric regulation.

Membrane proteins mediate processes that are fundamental for the flourishing of biological cells. Membrane-embedded transporters move ions and larger solutes across membranes; receptors mediate communication between the cell and its environment and membrane-embedded enzymes catalyze chemical reactions. Understanding these mechanisms of action requires knowledge of how the proteins couple to their fluid, hydrated lipid membrane environment. We present here current studies in computational and experimental membrane protein biophysics, and show how they address outstanding challenges in understanding the complex environmental effects on the structure, function, and dynamics of membrane proteins.

We calculate the absolute binding free energies of tetra-methylated octa-acids host–guest systems as a part of the SAMPL6 blind challenge (receipt ID vq30p). We employed two different free energy simulation methods, i.e., the umbrella sampling (US) and double decoupling method (DDM). The US method was used with the weighted histogram analysis method (WHAM) (US-WHAM scheme). In the DDM scheme, Hamiltonian replica-exchange method (HREM) was combined with the Bennett acceptance ratio (BAR) (HREM-BAR scheme). We obtained initial binding poses via molecular docking using GalaxyDock-HG program, which is developed for the SAMPL challenge. The root mean square deviation (RMSD) and the mean absolute deviations (MAD) using US-WHAM scheme were 1.33 and 1.02 kcal/mol, respectively. The MAD was the top among all submissions, however the correlation with respect to experiment was unexceptional. While the RMSD and MAD via HREM-BAR scheme were greater than US-WHAM scheme, (i.e., 2.09 and 1.76 kcal/mol), their correlations were slightly better than US-WHAM. The correlation between the two methods was high. Further discussion on the DDM method can be found in a companion paper by Han et al. (receipt ID 3z83m) in the same issue.

SAMHD1 is a deoxynucleoside triphosphate triphosphohydrolase (dNTPase) that depletes cellular dNTPs in noncycling cells to promote genome stability and to inhibit retroviral and herpes viral replication. In addition to being substrates, cellular nucleotides also allosterically regulate SAMHD1 activity. Recently, it was shown that high expression levels of SAMHD1 are also correlated with significantly worse patient responses to nucleotide analog drugs important for treating a variety of cancers, including acute myeloid leukemia (AML). In this study, we used biochemical, structural, and cellular methods to examine the interactions of various cancer drugs with SAMHD1. We found that both the catalytic and the allosteric sites of SAMHD1 are sensitive to sugar modifications of the nucleotide analogs, with the allosteric site being significantly more restrictive. We crystallized cladribine-TP, clofarabine-TP, fludarabine-TP, vidarabine-TP, cytarabine-TP, and gemcitabine-TP in the catalytic pocket of SAMHD1. We found that all of these drugs are substrates of SAMHD1 and that the efficacy of most of these drugs is affected by SAMHD1 activity. Of the nucleotide analogs tested, only cladribine-TP with a deoxyribose sugar efficiently induced the catalytically active SAMHD1 tetramer. Together, these results establish a detailed framework for understanding the substrate specificity and allosteric activation of SAMHD1 with regard to nucleotide analogs, which can be used to improve current cancer and antiviral therapies.

This study reports the results of binding free energy calculations for CB[8] host–guest systems in the SAMPL6 blind challenge (receipt ID 3z83m). Force-field parameters were developed specific for each of host and guest molecules to improve configurational sampling. We used quantum mechanical (QM) implicit solvent calculations and QM force matching to determine non-bonded (partial atomic charges) and bonded terms, respectively. Free energy calculations were carried out using the double-decoupling method (DDM) combined with Hamiltonian replica exchange method (HREM) and Bennett acceptance ratio (BAR) method. The root mean square error (RMSE) of the predicted values using DDM with respect to the experimental results was 4.32 kcal/mol. The coefficient of determination (R2) and Kendall rank coefficient (τ) were 0.49 and 0.52, respectively, highest of all submissions. In addition, these were compared to the results obtained by umbrella sampling (US) and weighted histogram analysis method (WHAM). Overall, DDM achieved a higher prediction accuracy than the US method. Results are discussed in terms of parameterization and free energy simulations.

The computation of distribution coefficients between polar and apolar phases requires both an accurate characterization of transfer free energies between phases and proper accounting of ionization and protomerization. We present a protocol for accurately predicting partition coefficients between two immiscible phases, and then apply it to 53 drug-like molecules in the SAMPL5 blind prediction challenge. Our results combine implicit solvent QM calculations with classical MD simulations using the non-Boltzmann Bennett free energy estimator. The OLYP/DZP/SMD method yields predictions that have a small deviation from experiment (RMSD = 2.3 [Formula: see text] D units), relative to other participants in the challenge. Our free energy corrections based on QM protomer and [Formula: see text] calculations increase the correlation between predicted and experimental distribution coefficients, for all methods used. Unfortunately, these corrections are overly hydrophilic, and fail to account for additional effects such as aggregation, water dragging and the presence of polar impurities in the apolar phase. We show that, although expensive, QM-NBB free energy calculations offer an accurate and robust method that is superior to standard MM and QM techniques alone.

One of the central aspects of biomolecular recognition is the hydrophobic effect, which is experimentally evaluated by measuring the distribution coefficients of compounds between polar and apolar phases. We use our predictions of the distribution coefficients between water and cyclohexane from the SAMPL5 challenge to estimate the hydrophobicity of different explicit solvent simulation techniques. Based on molecular dynamics trajectories with the CHARMM General Force Field, we compare pure molecular mechanics (MM) with quantum-mechanical (QM) calculations based on QM/MM schemes that treat the solvent at the MM level. We perform QM/MM with both density functional theory (BLYP) and semi-empirical methods (OM1, OM2, OM3, PM3). The calculations also serve to test the sensitivity of partition coefficients to solute polarizability as well as the interplay of the quantum-mechanical region with the fixed-charge molecular mechanics environment. Our results indicate that QM/MM with both BLYP and OM2 outperforms pure MM. However, this observation is limited to a subset of cases where convergence of the free energy can be achieved.

As part of the SAMPL5 blind prediction challenge, we calculate the absolute binding free energies of six guest molecules to an octa-acid (OAH) and to a methylated octa-acid (OAMe). We use the double decoupling method via thermodynamic integration (TI) or Hamiltonian replica exchange in connection with the Bennett acceptance ratio (HREM-BAR). We produce the binding poses either through manual docking or by using GalaxyDock-HG, a docking software developed specifically for this study. The root mean square deviations for our most accurate predictions are 1.4 kcal mol −1 for OAH with TI and 1.9 kcal mol −1 for OAMe with HREM-BAR. Our best results for OAMe were obtained for systems with ionic concentrations corresponding to the ionic strength of the experimental solution. The most problematic system contains a halogenated guest. Our attempt to model the σ -hole of the bromine using a constrained off-site point charge, does not improve results. We use results from molecular dynamics simulations to argue that the distinct binding affinities of this guest to OAH and OAMe are due to a difference in the flexibility of the host. We believe that the results of this extensive analysis of host-guest complexes will help improve the protocol used in predicting binding affinities for larger systems, such as protein-substrate compounds.

Herein, we report the absolute binding free energy calculations of CBClip complexes in the SAMPL5 blind challenge. Initial conformations of CBClip complexes were obtained using docking and molecular dynamics simulations. Free energy calculations were performed using thermodynamic integration (TI) with soft-core potentials and Bennett’s acceptance ratio (BAR) method based on a serial insertion scheme. We compared the results obtained with TI simulations with soft-core potentials and Hamiltonian replica exchange simulations with the serial insertion method combined with the BAR method. The results show that the difference between the two methods can be mainly attributed to the van der Waals free energies, suggesting that either the simulations used for TI or the simulations used for BAR, or both are not fully converged and the two sets of simulations may have sampled difference phase space regions. The penalty scores of force field parameters of the 10 guest molecules provided by CHARMM Generalized Force Field can be an indicator of the accuracy of binding free energy calculations. Among our submissions, the combination of docking and TI performed best, which yielded the root mean square deviation of 2.94 kcal/mol and an average unsigned error of 3.41 kcal/mol for the ten guest molecules. These values were best overall among all participants. However, our submissions had little correlation with experiments.