Efficient and Accurate Modeling of Conformational Transitions in Proteins: The Case of c-Src Kinase

Point Mutations at a Key Site Alter the Cytochrome P450 OleP Structural Dynamics

Substrate binding to the cytochrome P450 OleP is coupled to a large open-to-closed transition that remodels the active site, minimizing its exposure to the external solvent. When the aglycone substrate binds, a small empty cavity is formed between the I and G helices, the BC loop, and the substrate itself, where solvent molecules accumulate mediating substrate-enzyme interactions. Herein, we analyzed the role of this cavity in substrate binding to OleP by producing three mutants (E89Y, G92W, and S240Y) to decrease its volume. The crystal structures of the OleP mutants in the closed state bound to the aglycone 6DEB showed that G92W and S240Y occupied the cavity, providing additional contact points with the substrate. Conversely, mutation E89Y induces a flipped-out conformation of this amino acid side chain, that points towards the bulk, increasing the empty volume. Equilibrium titrations and molecular dynamic simulations indicate that the presence of a bulky residue within the cavity impacts the binding properties of the enzyme, perturbing the conformational space explored by the complexes. Our data highlight the relevance of this region in OleP substrate binding and suggest that it represents a key substrate-protein contact site to consider in the perspective of redirecting its activity towards alternative compounds.

Quaternary Structure Transitions of Human Hemoglobin: An Atomic-Level View of the Functional Intermediate States

Human hemoglobin (HbA) is one of the prototypal systems used to investigate structure-function relationships in proteins. Indeed, HbA has been used to develop the basic concepts of protein allostery, although the atomic-level mechanism underlying the HbA functionality is still highly debated. This is due to the fact that most of the three-dimensional structural information collected over the decades refers to the endpoints of HbA functional transition with little data available for the intermediate states. Here, we report molecular dynamics (MD) simulations by focusing on the relevance of the intermediate states of the protein functional transition unraveled by the crystallographic studies carried out on vertebrate Hbs. Fully atomistic simulations of the HbA T-state indicate that the protein undergoes a spontaneous transition toward the R-state. The inspection of the trajectory structures indicates that the protein significantly populates the intermediate HL-(C) state previously unraveled by crystallography. In the structural transition, it also assumes the intermediate states crystallographically detected in Antarctic fish Hbs. This finding suggests that HbA and Antarctic fish Hbs, in addition to the endpoints of the transitions, also share a similar deoxygenation pathway despite a distace of hundreds of millions of years in the evolution scale. Finally, using the essential dynamic sampling methodology, we gained some insights into the reverse R to T transition that is not spontaneously observed in classic MD simulations.

All-atom adaptively biased path optimization of Src kinase conformational inactivation: Switched electrostatic network in the concerted motion of αC helix and the activation loop

A method to optimize a conformational pathway through a space of well-chosen reduced variables is employed to advance our understanding of protein conformational equilibrium. The adaptively biased path optimization strategy utilizes unrestricted, enhanced sampling in the region of a path in the reduced-variable space to identify a broad path between two stable end-states. Application to the inactivation transition of the Src tyrosine kinase catalytic domain reveals new insight into this well studied conformational equilibrium. The mechanistic description gained from identifying the motions and structural features along the path includes details of the switched electrostatic network found to underpin the transition. The free energy barrier along the path results from rotation of a helix, αC, that is tightly correlated with motions in the activation loop (A-loop) as well as distal regions in the C-lobe. Path profiles of the reduced variables clearly demonstrate the strongly correlated motions. The exchange of electrostatic interactions among residues in the network is key to these interdependent motions. In addition, the increased resolution from an all-atom model in defining the path shows multiple components for the A-loop motion and that different parts of the A-loop contribute throughout the length of the path.

Quantitative Characterization of Binding Pockets and Binding Complementarity by Means of Zernike Descriptors

In this work, we describe the application of the Zernike formalism to quantitatively characterize the binding pockets of two sets of biologically relevant systems. Such an approach, when applied to molecular dynamics trajectories, is able to pinpoint the subtle differences between very similar molecular regions and their impact on the local propensity to ligand binding, allowing us to quantify such differences. The statistical robustness of our procedure suggests that it is very suitable to describe protein binding sites and protein-ligand interactions within a rigorous and well-defined framework.