Comparative effects of oncogenic mutations G12C, G12V, G13D, and Q61H on local conformations and dynamics of K-Ras

K-Ras is the most frequently mutated protein in human cancers. However, until very recently, its oncogenic mutants were viewed as undruggable. To develop inhibitors that directly target oncogenic K-Ras mutants, we need to understand both their mutant-specific and pan-mutant dynamics and conformations. Recently, we have investigated how the most frequently observed K-Ras mutation in cancer patients, G12D, changes its local dynamics and conformations (Vatansever et al., 2019). Here, we extend our analysis to study and compare the local effects of other frequently observed oncogenic mutations, G12C, G12V, G13D and Q61H. For this purpose, we have performed Molecular Dynamics (MD) simulations of each mutant when active (GTP-bound) and inactive (GDP-bound), analyzed their trajectories, and compared how each mutant changes local residue conformations, inter-protein distance distributions, local flexibility and residue pair correlated motions. Our results reveal that in the four active oncogenic mutants we have studied, the α2 helix moves closer to the C-terminal of the α3 helix. However, P-loop mutations cause α3 helix to move away from Loop7, and only G12 mutations change the local conformational state populations of the protein. Furthermore, the motions of coupled residues are mutant-specific: G12 mutations lead to new negative correlations between residue motions, while Q61H destroys them. Overall, our findings on the local conformational states and protein dynamics of oncogenic K-Ras mutants can provide insights for both mutant-selective and pan-mutant targeted inhibition efforts.

The dynamical Matryoshka model: 1. Incoherent neutron scattering functions for lipid dynamics in bilayers

Fluid lipid bilayers are the building blocks of biological membranes. Although there is a large amount of experimental data using incoherent quasi-elastic neutron scattering (QENS) techniques to study membranes, very little theoretical works have been developed to study the local dynamics of membranes. The main objective of this work is to build a theoretical framework to study and describe the local dynamics of lipids and derive analytical expressions of intermediate scattering functions (ISF) for QENS. As results, we developed the dynamical Matryoshka model which describes the local dynamics of lipid molecules in membrane layers as a nested hierarchical convolution of three motional processes: (i) individual motions described by the vibrational motions of H-atoms; (ii) internal motions including movements of the lipid backbone, head groups and tails, and (iii) molecule movements of the lipid molecule as a whole. The analytical expressions of the ISF associated with these movements are all derived. For use in analyzing the QENS experimental data, we also derived an analytical expression for the aggregate ISF of the Matryoshka model which involves an elastic term plus three inelastic terms of well-separated time scales and whose amplitudes and rates are functions of the lipid motions. And as an illustrative application, we used the aggregated ISF to analyze the experimental QENS data on a lipid sample of multilamellar bilayers of DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine). It is clear from this analysis that the dynamical Matryoshka model describes very well the experimental data and allow extracting the dynamical parameters of the studied system.

Towards the energy landscape of adenylate kinase in crowded milieu: Activity, conformation, structure and dynamics in sequence

Enzyme function is governed by a complex network of conformational changes and internal dynamics, with the same getting more convoluted in the crowded cellular environment. Here, we have explored an intricate interplay amongst activity, structure, conformation, and dynamics of a multidomain enzyme, AK3L1 (UniProtKB: Q9UIJ7) in the crowded milieu. We have monitored changes in the enzyme landscape in response to the chemical denaturant, urea, under the influence of different concentrations of macromolecular crowders. Extensive experimental analyses using FRET-based domain displacement measurements, sub-nanosecond time scale local dynamics, and global structural changes, along with enzymatic activity studies, has been carried out to get deeper insights into the factors that may modulate the functional landscape of adenylate kinase (AK3L1). It was observed that AK3L1 gets activated at low urea concentrations, whereas higher urea concentrations unfold and thereby deactivate the enzyme. A sequential response of AK3L1 is observed towards external perturbation (urea) occurring through a series of well-defined steps. Incorporation of crowders not only shift the maximum activity of enzyme to a higher urea concentration, but also enhance domain compaction, as revealed by FRET studies. The modulation in enzyme activity and solvation dynamics acting as local response, precede global unfolding of the enzyme, indicating that the structural alterations around the active site are quite decoupled from the large amplitude global transitions.