The role of glycerol-water mixtures in the stability of FKBP12-rapalog-FRB complexes

The thermodynamic and biophysical implications of the introduction of a co-solvent during protein-ligand binding remain elusive. Using ternary complexes of 12-kDa FK506 binding protein (FKBP12), FKBP-rapamycin binding (FRB) domain of the mammalian/mechanistic target of rapamycin (mTOR) kinase, and rapamycin analogs (rapalogs) in glycerol-water mixtures, the influence of solvent composition on ligand binding dynamics was explored. The pharmaceutical potential of rapalogs and the utility of glycerol as a co-solvent in drug delivery applications were critical in deciding the system to be studied. Consolidation of existing studies on rapamycin modification was first performed to strategically design a new rapalog called T1. The results from 100-ns dual-boost Gaussian accelerated molecular dynamics simulations showed that protein stability was induced in the presence of glycerol. Reweighting of the trajectories revealed that the glycerol-rich solvent system lowers the energy barrier in the conformational space of the protein while also preserving native contacts between the ligand and the residues in the binding site. Calculated binding free energies using MM/GBSA also showed that electrostatic energy and polar contribution of solvation energy are heavily influenced by the changes in solvation. Glycerol molecules are preferentially excluded through electrostatic interactions from the solvation shell which induce complex stability as seen in existing experiments. Hence, using glycerol as a co-solvent in rapamycin delivery has a significant role in maintaining stability. In addition, compound T1 is a potential mTORC1-selective inhibitor with strong affinity for the FKBP12-FRB complex. This study aims to provide insights on the design of new rapalogs, and the applicability of glycerol as co-solvent for FKBP12-rapalog-FRB complexes.

Reparameterization of Non-Bonded Parameters for Copper Ions in Plastocyanin: An Adaptive Force Matching Study

Molecular mechanics rely on existing experimental and theoretical inputs to confidently calculate the trajectories of molecular systems. These calculations, however, are often hindered by missing force field parameters. A notable subject of this problem is metal centers of proteins. This study parameterized, through an adaptive force matching (AFM) workflow, the copper cofactor of plastocyanin in its two oxidation states. New 12-6 Lennard-Jones (LJ) parameters and atomic partial charges were generated to complete the non-bonded description of the copper site. Our models show uniform distorted tetrahedral structures for reduced plastocyanin, Cu(I), and oxidized plastocyanin, Cu(II). These structures align with the QM/MM MD results and existing crystallography studies. TD-DFT calculations, meanwhile, showed that conformations with elongated axial Cu-SMet and shortened equatorial Cu-SCys bonds retain the experimental UV-Vis profile of blue copper (BC) proteins, thus signifying the importance of Cu-S interactions on BC proteins' unique spectroscopic properties.