Search

A microscopic study of functional disorder-order folding transitions coupled to binding is performed for the p27 protein, which derives a kinetic advantage from the intrinsically disordered unbound form on binding with the phosphorylated cyclin A-cyclindependent kinase 2 (Cdk2) complex. Hierarchy of structural loss during p27 coupled unfolding and unbinding is simulated by using high-temperature Monte Carlo simulations initiated from the crystal structure of the tertiary complex. Subsequent determination of the transition-state ensemble and the proposed atomic picture of the folding mechanism coupled to binding provide a microscopic rationale that reconciles the initiation recruitment of p27 at the cyclin A docking site with the kinetic benefit for a disordered [alpha]-helix in the unbound form of p27. The emerging structural polarization in the ensemble of unfolding/unbinding trajectories and in the computationally determined transition-state ensemble is not determined by the intrinsic folding preferences of p27 but rather is attributed to the topological requirements of the native intermolecular interface to order [beta]-hairpin and [beta]-strand of p27 that could be critical for nucleating rapid folding transition coupled to binding. In agreement with the experimental data, the disorder-order folding transition for p27 is largely determined by the functional requirement to form a specific intermolecular interface that ultimately dictates the folding mechanism and overwhelms any local folding preferences for creating a stable [alpha]-helix in the p27 structure before overcoming the major free energy barrier.

The search for novel leads is a critical step in the drug discovery process. Computational approaches to identify new lead molecules have focused on discovering complete ligands by evaluating the binding affinity of a large number of candidates, a task of considerable complexity. A new computational method is introduced in this work based on the premise that the primary molecular recognition event in the protein binding site may be accomplished by small core fragments that serve as molecular anchors, providing a structurally stable platform that can be subsequently tailored into complete ligands. To fulfill its role, we show that an effective molecular anchor must meet both the thermodynamic requirement of relative energetic stability of a single binding mode and its consistent kinetic accessibility, which may be measured by the structural consensus of multiple docking simulations. From a large number of candidates, this technique is able to identify known core fragments responsible for primary recognition by the FK506 binding protein (FKBP-12), along with a diverse repertoire of novel molecular cores. By contrast, absolute energetic criteria for selecting molecular anchors are found to be promiscuous. A relationship between a minimum frustration principle of binding energy landscapes and receptor-specific molecular anchors in their role as 'recognition nuclei' is established, thereby unraveling a mechanism of lead discovery and providing a practical route to receptor-biased computational combinatorial chemistry.

A general mean field model of ligand-protein interactions has been developed to investigate the structure of ligand-protein complexes at finite temperature. The model incorporates the gross effects of protein flexibility by employing mean field theory with a finite number of protein conformations. Protein conformers in the model's database were superimposed by least square superposition of the Calpha atoms to derive the energy spectrum.

A new technique to analyze the energy flow in liquid-phase dynamics is described. Some of the energy in the transition is transferred from the reaction coordinate to other internal degrees of the solute. The rate at which the remaining energy is transferred to the solvent is a function of the reaction of coordinate phase space point and its time derivative as well as other velocities characterizing a trajectory of the system.