ABSTRACT
The de novo computational design of proteins with predefined three-dimensional structure is becoming much more routine due to advancements both in force fields and algorithms. However, creating designs with functions beyond folding is more challenging. In that regard, the recent design of small beta barrel proteins that activate the fluorescence of an exogenous small molecule chromophore (DFHBI) is noteworthy. These proteins, termed mini fluorescence activating proteins (mFAPs), have been shown to increase the brightness of the chromophore more than 100-fold upon binding to the designed ligand pocket. The design process created a large library of variants with different brightness levels but gave no rational explanation for why one variant was brighter than another. Here, we use quantum mechanics and molecular dynamics simulations to investigate how molecular flexibility in the ground and excited states influences brightness. We show that the ability of the protein to resist dihedral angle rotation of the chromophore is critical for predicting brightness. Our simulations suggest that the mFAP/DFHBI complex has a rough energy landscape, requiring extensive ground-state sampling to achieve converged predictions of excited-state kinetics. While computationally demanding, this roughness suggests that mFAP protein function can be enhanced by reshaping the energy landscape toward conformations that better resist DFHBI bond rotation.