In this work, the solvation free energies of 20 organic molecules from the 4th Statistical Assessment of the Modeling of Proteins and Ligands (SAMPL4) have been calculated. The sampling of phase space is carried out at a molecular mechanical level, and the associated free energy changes are estimated using the Bennett Acceptance Ratio (BAR). Then the quantum mechanical (QM) corrections are computed through the indirect Non-Boltzmann Bennett's acceptance ratio (NBB) or the thermodynamics perturbation (TP) method. We show that BAR+TP gives a minimum analytic variance for the calculated solvation free energy at the Gaussian limit and performs slightly better than NBB in practice. Furthermore, the expense of the QM calculations in TP is only half of that in NBB. We also show that defining the biasing potential as the difference of the solute-solvent interaction energy, instead of the total energy, can converge the calculated solvation free energies much faster but possibly to different values. Based on the experimental solvation free energies which have been published before, it is discovered in this study that BLYP yields better results than MP2 and some other later functionals such as B3LYP, M06-2X, and ωB97X-D.

The partitioning of solute molecules between immiscible solvents with significantly different polarities is of great importance. The polarization between the solute and solvent molecules plays an essential role in determining the solubility of the solute, which makes computational studies utilizing molecular mechanics (MM) rather difficult. In contrast, quantum mechanics (QM) can provide more reliable predictions. In this work, the partition coefficients of the side chain analogs of some amino acids between water and chloroform were computed. The QM solvation free energies were calculated indirectly via a series of MM states using the multistate Bennett acceptance ratio (MBAR) and the MM-to-QM corrections were applied at the two endpoints using thermodynamic perturbation (TP). Previously, it has been shown (Jia et al. J. Chem. Theory Comput. 2016, 12, 499-511) that this method provides the minimal variance in the results without running QM simulations. However, if there is insufficient overlap in phase space between the MM and QM Hamiltonians, this method fails. In this work, we propose, for the first time, a quantity termed the reweighting entropy that serves as a metric for the reliability of the TP calculations. If the reweighting entropy is below a certain threshold (0.65 for the solvation free energy calculations in this work), this MM-to-QM correction should be avoided and two alternative methods can be employed by either introducing a semiempirical state or conducting nonequilibrium simulations. However, the results show that the QM methods are not guaranteed to yield better results than the MM methods. Further improvement of the QM methods are imperative, especially the treatment of the van der Waals and the electrostatic interactions between the QM region and the MM region in the first shell. We also propose a scheme for the calculation of the van der Waals parameters for the solute molecules in nonaqueous solvent, which improves the quality of the computed thermodynamic properties. Furthermore, the force field parameters for the sulfur-containing molecules are also optimized.

Previous cross-sectional studies have suggested an association between migraine and rheumatoid arthritis (RA), but no longitudinal study has been performed to evaluate the temporal relationship between the two conditions. The purpose of the present population-based, propensity score-matched cohort study was to investigate whether migraineurs are at a higher risk of developing RA. A total of 58,749 subjects aged between 20 and 90 years with at least two ambulatory visits with a diagnosis of migraine were recruited in the migraine group. We fit a logistic regression model that included age, sex, comorbid conditions, and socioeconomic status as covariates to compute the propensity score. The non-migraine group consisted of 58,749 propensity score-matched, randomly sampled subjects without migraine. The RA-free survival curves were generated using the Kaplan-Meier method. Stratified Cox proportional hazard regression was used to estimate the effect of migraine on the risk of RA. During follow-up, 461 subjects in the migraine group and 220 in the non-migraine group developed RA. The incidence rate of RA was 3.18 (95% confidence interval [CI] 2.90-3.49) per 1000 person-years in the migraine group and 1.54 (95% CI 1.34-1.76) per 1000 person-years in the non-migraine group. Compared to the non-migraine group, the crude hazard ratio of RA for the migraine group was 2.15 (95% CI 1.82-2.56, P < 0.0001), and the multivariable-adjusted hazard ratio was 1.91 (95% CI 1.58-2.31, P < 0.0001). This study showed that patients with migraine had an increased risk of developing RA.

Calculating binding free energies with quantum mechanical (QM) methods is notoriously time consuming. In this work, we study whether such calculations can be accelerated by using non-equilibrium (NE) molecular dynamics simulations employing Jarzynski's equality. We study the binding of nine cyclic carboxylate ligands to the octa-acid deep-cavity host from the SAMPL4 challenge with the reference-potential approach. The binding free energies were first calculated at the molecular-mechanics (MM) level with free-energy perturbation, using the generalised Amber force field with restrained electrostatic-potential charges for the host and the ligands. Then, the freeenergy corrections for going from the MM Hamiltonian to a hybrid QM/MM Hamiltonian were estimated by averaging over many short NE molecular dynamics simulations. In the QM/MM calculations, the ligand is described at the semiempirical PM6-DH+ level. We show that this approach yields MM→QM/MM free energy corrections that agree with those from other approaches within statistical uncertainties. The desired precision can be obtained by running a proper number of independent NE simulations. For the systems studied in this work, a total simulation length of 20 ps was appropriate for most ligands and 36-324 simulations were necessary in order to reach a precision of 0.3 kJ/mol.

Tamoxifen is a widely used personalized medicine for estrogen receptor (ER)-positive breast cancer, but approximately 30% of patients receiving the treatment relapse due to tamoxifen resistance (TamR). Recently, several reports have linked lncRNAs to cancer drug resistance. However, the role of lncRNAs in TamR is unclear. To identify TamR-related lncRNAs, we first used a bioinformatic approach to predict whether they have connection with known TamR-associated genes by starBase v2.0 and divided them into two groups. Group A contains lncRNAs that connect with known TamR genes and group B contains lncRNAs that show no predicted interaction. Among the 12 lncRNAs in group A, 58.3% of them are either up- or downregulated in MCF-7/TamR cells compared to the sensitive cells. In contrast, the expression levels of all group B lncRNAs are not changed in MCF-7/TamR cells. LINC00894-002 exhibits the most sophisticated network pattern and is the most downregulated lncRNA in MCF-7/TamR cells. Moreover, we find that LINC00894-002 is directly upregulated by ERα. Knocking down LINC00894-002 downregulates expression of miR-200a-3p and miR-200b-3p, upregulates the expression of TGF-β2 and ZEB1, and finally contributes to TamR. Herein, we report the first case of an inhibitory lncRNA against TamR through the miR-200-TGF-β2-ZEB1 signaling pathway.

No abstract

Relative free energies for the binding of nine cyclic carboxylate ligands to the octaacid deep-cavity host were calculated at the combined density-functional theory and molecular mechanics (DFT/MM) level of theory. The DFT calculations employed the BLYP functional and the 6-31G* basis set for the ligand. We employed free-energy perturbations (FEP) with the reference-potential approach and used molecular dynamics (MD) simulations with the semiempirical quantum mechanical (SQM) PM6-DH+ method for the ligand as an intermediate level between MM and DFT/MM to improve the convergence. Thus, the relative binding free energy of two ligands was first calculated at the MM level by an alchemical transformation from one ligand to another in both the bound and unbound states. Then, for each ligand the free-energy correction for going from the MM to the SQM/MM potentials was calculated using explicit SQM/MM MD simulations. Finally, the free-energy correction for going from the SQM/MM to the DFT/MM potentials was estimated with FEP without running any DFT/MM simulations. Instead, the free energy was calculated by single-step exponential averaging (ssEA) or employing the cumulant approximation to the second order (CA). The results show that CA converges much better than ssEA and with 500-4500 DFT/MM single-point energy calculations, converged free energies with a precision of 0.3 kJ/mol can be obtained. These free energies reproduce the experimental binding free energy differences with a mean absolute deviation of 3.4 kJ/mol, a correlation (R 2) of 0.97 and correct signs for all of the eight free-energy differences. This is appreciably better than the results obtained at the SQM/MM level of theory and also slightly better than those obtained with MM. We show that the convergence of the SQM/MM→DFT/MM perturbations can be monitored by the use of Wu and Kofke's bias metric Π and by the standard deviation of the difference between the SQM/MM and DFT/MM energies. Finally, we show that the use of the intermediate SQM/MM MD simulations improves the convergence of the free energies by a factor of at least two, compared to doing direct MM→DFT/MM perturbations.

We have estimated free energies for the binding of eight carboxylate ligands to two variants of the octa-acid deep-cavity host in the SAMPL6 blind-test challenge (with or without endo methyl groups on the four upper-rim benzoate groups, OAM and OAH, respectively). We employed free-energy perturbation (FEP) for relative binding energies at the molecular mechanics (MM) and the combined quantum mechanical (QM) and MM (QM/MM) levels, the latter obtained with the reference-potential approach with QM/MM sampling for the MM → QM/MM FEP. The semiempirical QM method PM6-DH+ was employed for the ligand in the latter calculations. Moreover, binding free energies were also estimated from QM/MM optimised structures, combined with COSMO-RS estimates of the solvation energy and thermostatistical corrections from MM frequencies. They were performed at the PM6-DH+ level of theory with the full host and guest molecule in the QM system (and also four water molecules in the geometry optimisations) for 10-20 snapshots from molecular dynamics simulations of the complex. Finally, the structure with the lowest free energy was recalculated using the dispersion-corrected density-functional theory method TPSS-D3, for both the structure and the energy. The two FEP approaches gave similar results (PM6-DH+/MM slightly better for OAM), which were among the five submissions with the best performance in the challenge and gave the best results without any fit to data from the SAMPL5 challenge, with mean absolute deviations (MAD) of 2.4-5.2 kJ/mol and a correlation coefficient (R) of 0.77-0.93. This is the first time QM/MM approaches give binding free energies that are competitive to those obtained with MM for the octa-acid host. The QM/MM-optimised structures gave somewhat worse performance (MAD = 3-8 kJ/mol and R = 0.1-0.9), but the results were improved compared to previous studies of this system with similar methods.