Abstract:The substitution of methyl (Me or −CH3) by trifluoromethyl (TFM or −CF3) is frequently used in medicinal chemistry. However, the exact effect of −CH3/–CF3 substitution on bioactivity is still controversial. We compiled a data set containing 28 003 pairs of compounds with the only difference that −CH3 is substituted by −CF3, and the statistical results showed that the replacement of −CH3 with −CF3 does not improve bioactivity on average. Yet, 9.19% substitution of −CH3 by −CF3 could increase the biological activity by at least an order. A PDB survey revealed that −CF3 prefers Phe, Met, Leu, and Tyr, while −CH3 prefers Leu, Met, Cys, and Ile. If we substitute the −CH3 by −CF3 near Phe, His, and Arg, the bioactivity is most probably improved. We performed QM/MM calculations for 39 −CH3/–CF3 pairs of protein–ligand complexes and found that the −CH3/–CF3 substitution does achieve a large energy gain in some systems, although the mean energy difference is subtle, which is consistent with the statistical survey. The −CF3 substitution on the benzene ring could be particularly effective at gaining binding energy. The maximum improvements in energy achieved −4.36 kcal/mol by QM/MM calculation. Moreover, energy decompositions from MM/GBSA calculations showed that the large energy gains for the −CH3/–CF3 substitution are largely driven by the electrostatic energy or the solvation free energy. These findings may shed some light on the biological activity profile for −CH3/–CF3 substitution, which should be useful for further drug discovery and drug design.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acs.jcim.0c00898?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acs.jcim.0c00898</a>.Additional figures and tables containing −CH3/–CF3 group ligands structures and related information from statistical and energy calculations (<a class="ext-link" href="/doi/suppl/10.1021/acs.jcim.0c00898/suppl_file/ci0c00898_si_001.pdf">PDF</a>)28 003 pairs of small molecule compounds with bioactivity information (<a class="ext-link" href="/doi/suppl/10.1021/acs.jcim.0c00898/suppl_file/ci0c00898_si_002.zip">ZIP</a>)This article has not yet been cited by other publications.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: What are the specific effects of trifluoromethyl (-CF₃) replacing methyl (-CH₃) on the biological activity in medicinal chemistry? Although trifluoromethyl is often used to replace methyl in medicinal chemistry to optimize drug performance, the specific effect of this substitution on biological activity is still controversial. Specifically, the research aims to answer the following questions: 1. **Does substitution generally improve biological activity?** Can replacing methyl with trifluoromethyl generally improve the biological activity of compounds? 2. **The effect of substitution in specific cases**: In which cases can trifluoromethyl replacing methyl significantly improve biological activity? 3. **The effect of substitution on binding affinity**: How does trifluoromethyl replacing methyl affect the binding affinity between ligands and proteins? ### Overview of research methods To answer these questions, the author conducted the following studies: 1. **Statistical analysis**: A data set of 28,003 pairs of compounds was collected. The only difference between these compounds was that -CH₃ at one position was replaced by -CF₃. Through statistical analysis, the effect of this substitution on biological activity was evaluated. 2. **PDB investigation**: The interactions between -CF₃ and -CH₃ with protein residues were investigated from the Protein Data Bank (PDB) to understand their preferences in protein - binding pockets. 3. **Quantum mechanics/molecular mechanics (QM/MM) calculation**: 39 pairs of protein - ligand complexes were selected for QM/MM calculations to evaluate the energy change of -CH₃ to -CF₃ substitution. 4. **Molecular mechanics/generalized Born surface area (MM/GBSA) calculation**: Further through MM/GBSA calculations, each component of the binding free energy was decomposed to understand the main source of energy gain. ### Main findings 1. **Statistical results**: The statistical results of 28,003 pairs of compounds show that the substitution from -CH₃ to -CF₃ does not guarantee an increase in biological activity in the average sense. However, in some cases, 9.19% of the substitutions can increase the biological activity by at least one order of magnitude. 2. **PDB investigation results**: -CF₃ is more likely to interact with phenylalanine (Phe), methionine (Met), leucine (Leu) and tyrosine (Tyr), while -CH₃ is more likely to interact with leucine, methionine, cysteine (Cys) and isoleucine (Ile). In particular, substitution near phenylalanine (Phe), histidine (His) and arginine (Arg) may improve the binding affinity. 3. **Energy calculation results**: - QM/MM calculations show that although the average energy difference is small (-0.35 kcal/mol), in some systems, substitution can bring significant energy gain (up to -4.36 kcal/mol). - MM/GBSA calculations indicate that electrostatic energy and solvation free energy are the main driving forces for energy gain. ### Conclusion The effect of trifluoromethyl replacing methyl on biological activity is complex and does not always bring a general improvement. However, under specific conditions (such as near certain amino acid residues), this substitution may significantly improve biological activity. These findings provide valuable references for drug design and development, especially for understanding the effect of -CH₃ to -CF₃ substitution on biological activity. ### Formula display Some of the formulas involved in the above process are as follows: - Binding energy calculation formula: \[ \Delta E_{\text{bind}} = E_{\text{complex}} - (E_{\text{ligand}} + E_{\text{protein}} + BSSE) \] where: - \( E_{\text{complex}} \) is the energy of the complex structure - \( E_{\text{ligand}} \) is the energy of the ligand - \( E_{\text{protein}} \) is the energy of protein atoms - \( BSSE \) is the basis set superposition error - MM/GBSA binding free energy decomposition: \[ \Delta G_{\text{bind}} = \Delta E_{\text{vdw}} + \Delta E_{\text{elec}}+ \Delta G_{\text{solv}}+ \Delta G_{\text{non - pol}} \]

Substitution Effect of the Trifluoromethyl Group on the Bioactivity in Medicinal Chemistry: Statistical Analysis and Energy Calculations

Searching for "Greener" Bioequivalents of CF3 to Lower its Environmental Impact

Automated Fragmentation QM/MM Calculation of NMR Chemical Shifts for Protein-Ligand Complexes.

Computational Study of the Trifluoromethyl Radical Donor Abilities of CF3 Sources

Theoretical Calculations on 1-Substituted-2-aminobenzimidazole Derivatives and Studies on Quantitative Structure-activity Relationship (QSAR)

Exploring the impact of trifluoromethyl (–CF3) functional group on the anti-cancer activity of isoxazole-based molecules: design, synthesis, biological evaluation and molecular docking analysis

Identification of bioactive compounds with popular single-atom modifications: Comprehensive analysis and implications for compound design

Alchemical Binding Free Energy Calculations in AMBER20: Advances and Best Practices for Drug Discovery

Quantum Mechanical-Cluster Approach to Solve the Bioisosteric Replacement Problem in Drug Design

Quantum Descriptors for Predicting and Understanding the Structure–Activity Relationships of Michael Acceptor Warheads

Computational Chemistry: Prediction of Compound Accessibility of Targeted Synthesized Compounds

Computational Exploration of the Impact of Low‐Spin and High‐Spin Ground State on the Chelating Ability of Dimethylglyoxime Ligand on Dihalo Transition Metal: A QTAIM, EDA, and CDA Analysis

QSAR Study of 1-Substitute-2-amido- Benzimidazole

Three-dimensional QSAR of HPPD inhibitors, PSA inhibitors, and anxiolytic agents: effect of tautomerism on the CoMFA models.

Insights into the binding mechanism of 2,5‐substituted 4‐pyrone derivatives as therapeutic agents for fused dimeric interactions: A computational study using QTAIM, dynamics and docking simulations of protein–ligand complexes

Quantum Mechanics/Molecular Mechanics Modeling of Regioselectivity of Drug Metabolism in Cytochrome P450 2C9

Kinetic studies of imidazole and its methyl derivative binding to metmyoglobin: Effects of substitute methyl on the binding affinity

Synthesis and Biological Evaluation of CF 3 Se‐Substituted α‐Amino Acid Derivatives

The rat acute oral toxicity of trifluoromethyl compounds (TFMs): a computational toxicology study combining the 2D-QSTR, read-across and consensus modeling methods

Two halogenated coplanar Schiff base complexes: Insights into the impact of halogen substitution on intermolecular interactions via theoretical calculations

An Energetic Guide for Estimating Trifluoromethyl Cation Donor Abilities of Electrophilic Trifluoromethylating Reagents: Computations of X-Cf3 Bond Heterolytic Dissociation Enthalpies