A Quantitative Model for cAMP Binding to the Binding Domain of MloK1

Spatiotemporal Resolution of Conformational Changes in Biomolecules by Combining Pulsed Electron-Electron Double Resonance Spectroscopy with Microsecond Freeze-Hyperquenching

The function of proteins is linked to their conformations that can be resolved with several high-resolution methods. However, only a few methods can provide the temporal order of intermediates and conformational changes, with each having its limitations. Here, we combine pulsed electron-electron double resonance spectroscopy with a microsecond freeze-hyperquenching setup to achieve spatiotemporal resolution in the angstrom range and lower microsecond time scale. We show that the conformational change of the C_α-helix in the cyclic nucleotide-binding domain of the Mesorhizobium loti potassium channel occurs within about 150 μs and can be resolved with angstrom precision. Thus, this approach holds great promise for obtaining 4D landscapes of conformational changes in biomolecules.

Allosteric Response of DNA Recognition Helices of Catabolite Activator Protein to cAMP and DNA Binding

The homodimeric catabolite activator protein (CAP) regulates the transcription of several bacterial genes based on the cellular concentration of cyclic adenosine monophosphate (cAMP). The binding of cAMP to CAP triggers allosteric communication between the cAMP binding domains (CBD) and DNA binding domains (DBD) of CAP, which entails repositioning of DNA recognition helices (F-helices) in the DBD to dock favorably to the target DNA. Despite considerable progress, much remains to be understood about the mechanistic details of DNA recognition by CAP and about the map of allosteric pathways involved in CAP-mediated gene transcription. The present study uses molecular dynamics and umbrella sampling simulations to investigate the mechanism of cAMP- and DNA-induced changes in the conformation and energetics of F-helices observed during the allosteric regulation of CAP by cAMP and the subsequent binding to the DNA promoter region. Using novel collective variables, the free energy profiles associated with the orientation and dynamics of F-helices in the unliganded, cAMP-bound, and cAMP-DNA-bound states of CAP are calculated and compared. The binding-induced alterations in the resultant free energy profiles reveal important flexibility constraints imposed on DBD upon cAMP and DNA binding. A comprehensive analysis of residue-wise interaction maps reveals potential allosteric pathways between CBD and DBD that facilitate the allosteric transduction of regulatory signals in CAP. The revelation that the predicted allosteric pathways crisscross the intersubunit interface offers important clues on the microscopic origin of the intersubunit cooperativity and dimer stability of CAP.

Heterogeneous and rate-dependent streptavidin-biotin unbinding revealed by high-speed force spectroscopy and atomistic simulations

Receptor-ligand interactions are essential for biological function and their binding strength is commonly explained in terms of static lock-and-key models based on molecular complementarity. However, detailed information on the full unbinding pathway is often lacking due, in part, to the static nature of atomic structures and ensemble averaging inherent to bulk biophysics approaches. Here we combine molecular dynamics and high-speed force spectroscopy on the streptavidin-biotin complex to determine the binding strength and unbinding pathways over the widest dynamic range. Experiment and simulation show excellent agreement at overlapping velocities and provided evidence of the unbinding mechanisms. During unbinding, biotin crosses multiple energy barriers and visits various intermediate states far from the binding pocket, while streptavidin undergoes transient induced fits, all varying with loading rate. This multistate process slows down the transition to the unbound state and favors rebinding, thus explaining the long lifetime of the complex. We provide an atomistic, dynamic picture of the unbinding process, replacing a simple two-state picture with one that involves many routes to the lock and rate-dependent induced-fit motions for intermediates, which might be relevant for other receptor-ligand bonds.

Thermodynamics of camphor migration in cytochrome P450cam by atomistic simulations

Understanding the mechanisms of ligand binding to enzymes is of paramount importance for the design of new drugs. Here, we report on the use of a novel biased molecular dynamics (MD) methodology to study the mechanism of camphor binding to cytochrome P450cam. Microsecond-long MD simulations allowed us to observe reaction coordinates characterizing ligand diffusion from the active site of cytochrome P450cam to solvent via three egress routes. These atomistic simulations were used to estimate thermodynamic quantities along the reaction coordinates and indicate diverse binding configurations. The results suggest that the diffusion of camphor along the pathway near the substrate recognition site (SRS) is thermodynamically preferred. In addition, we show that the diffusion near the SRS is triggered by a transition from a heterogeneous collection of closed ligand-bound conformers to the basin comprising the open conformations of cytochrome P450cam. The conformational change accompanying this switch is characterized by the retraction of the F and G helices and the disorder of the B' helix. These results are corroborated by experimental studies and provide detailed insight into ligand binding and conformational behavior of the cytochrome family. The presented methodology is general and can be applied to other ligand-protein systems.