Combining Free Energy Simulations and NMR Chemical-Shift Perturbation To Identify Transient Cation-π Contacts in Proteins

Mechanism of rotenone binding to respiratory complex I depends on ligand flexibility

Respiratory complex I is a major cellular energy transducer located in the inner mitochondrial membrane. Its inhibition by rotenone, a natural isoflavonoid, has been used for centuries by indigenous peoples to aid in fishing and, more recently, as a broad-spectrum pesticide or even a possible anticancer therapeutic. Unraveling the molecular mechanism of rotenone action will help to design tuned derivatives and to understand the still mysterious catalytic mechanism of complex I. Although composed of five fused rings, rotenone is a flexible molecule and populates two conformers, bent and straight. Here, a rotenone derivative locked in the straight form was synthesized and found to inhibit complex I with 600-fold less potency than natural rotenone. Large-scale molecular dynamics and free energy simulations of the pathway for ligand binding to complex I show that rotenone is more stable in the bent conformer, either free in the membrane or bound to the redox active site in the substrate-binding Q-channel. However, the straight conformer is necessary for passage from the membrane through the narrow entrance of the channel. The less potent inhibition of the synthesized derivative is therefore due to its lack of internal flexibility, and interconversion between bent and straight forms is required to enable efficient kinetics and high stability for rotenone binding. The ligand also induces reconfiguration of protein loops and side-chains inside the Q-channel similar to structural changes that occur in the open to closed conformational transition of complex I. Detailed understanding of ligand flexibility and interactions that determine rotenone binding may now be exploited to tune the properties of synthetic derivatives for specific applications.

Accuracy of Alternate Nonpolarizable Force Fields for the Determination of Protein-Ligand Binding Affinities Dominated by Cation-π Interactions

Modifying pair-specific Lennard-Jones parameters through the nonbonded FIX (NBFIX) feature of the CHARMM36 force field has proven cost-effective for improving the description of cation-π interactions in biological objects by means of pairwise additive potential energy functions. Here, two sets of newly optimized CHARMM36 force-field parameters including NBFIX corrections, coined CHARMM36m-NBF and CHARMM36-WYF, and the original force fields, namely CHARMM36m and Amber ff14SB, are used to determine the standard binding free energies of seven protein-ligand complexes containing cation-π interactions. Compared with precise experimental measurements, our results indicate that the uncorrected, original force fields significantly underestimate the binding free energies, with a mean error of 5.3 kcal/mol, while the mean errors of CHARMM36m-NBF and CHARMM36-WYF amount to 0.8 and 2.1 kcal/mol, respectively. The present study cogently demonstrates that the use of modified parameters jointly with NBFIX corrections dramatically increases the accuracy of the standard binding free energy of protein-ligand complexes dominated by cation-π interactions, most notably with CHARMM36m-NBF.