Normal Mode Analysis of KRas4B Reveals Partner Specific Dynamics

Pioneer in Molecular Biology: Conformational Ensembles in Molecular Recognition, Allostery, and Cell Function

In 1978, for my PhD, I developed the efficient O(n3) dynamic programming algorithm for the-then open problem of RNA secondary structure prediction. This algorithm, now dubbed the "Nussinov algorithm", "Nussinov plots", and "Nussinov diagrams", is still taught across Europe and the U.S. As sequences started coming out in the 1980s, I started seeking genome-encoded functional signals, later becoming a bioinformatics trend. In the early 1990s I transited to proteins, co-developing a powerful computer vision-based docking algorithm. In the late 1990s, I proposed the foundational role of conformational ensembles in molecular recognition and allostery. At the time, conformational ensembles and free energy landscapes were viewed as physical properties of proteins but were not associated with function. The classical view of molecular recognition and binding was based on only two conformations captured by crystallography: open and closed. I proposed that all conformational states preexist. Proteins always have not one folded form-nor two-but many folded forms. Thus, rather than inducing fit, binding can work by shifting the ensembles between states, and this shifting, or redistributing the ensembles to maintain equilibrium, is the origin of the allosteric effect and protein, thus cell, function. This transformative paradigm impacted community views in allosteric drug design, catalysis, and regulation. Dynamic conformational ensemble shifts are now acknowledged as the origin of recognition, allostery, and signaling, underscoring that conformational ensembles-not proteins-are the workhorses of the cell, pioneering the fundamental idea that dynamic ensembles are the driving force behind cellular processes. Nussinov was recognized as pioneer in molecular biology by JMB.

X-ray crystallographic and hydrogen deuterium exchange studies confirm alternate kinetic models for homolog insulin monomers

Despite the crucial role of various insulin analogs in achieving satisfactory glycemic control, a comprehensive understanding of their in-solution dynamic mechanisms still holds the potential to further optimize rapid insulin analogs, thus significantly improving the well-being of individuals with Type 1 Diabetes. Here, we employed hydrogen-deuterium exchange mass spectrometry to decipher the molecular dynamics of newly modified and functional insulin analog. A comparative analysis of H/D dynamics demonstrated that the modified insulin exchanges deuterium atoms faster and more extensively than the intact insulin aspart. Additionally, we present new insights derived from our 2.5 Å resolution X-ray crystal structure of modified hexamer insulin analog at ambient temperature. Furthermore, we obtained a distinctive side-chain conformation of the Asn3 residue on the B chain (AsnB3) by operating a comparative analysis with a previously available cryogenic rapid-acting insulin structure (PDB_ID: 4GBN). The experimental conclusions have demonstrated compatibility with modified insulin's distinct cellular activity, comparably to aspart. Additionally, the hybrid structural approach combined with computational analysis employed in this study provides novel insight into the structural dynamics of newly modified and functional insulin vs insulin aspart monomeric entities. It allows further molecular understanding of intermolecular interrelations driving dissociation kinetics and, therefore, a fast action mechanism.

In Silico Identification of Putative Allosteric Pockets and Inhibitors for the KRASG13D-SOS1 Complex in Cancer Therapy

RAS mutations occur in about 30% of human cancers, leading to enhanced RAS signaling and tumor growth. KRAS is the most commonly mutated oncogene in human tumors, especially lung, pancreatic, and colorectal cancers. Direct targeting of KRAS is difficult due to its highly conserved sequence; but, its complex with the guanine nucleotide exchange factor Son of Sevenless (SOS) 1 promises an attractive target for inhibiting RAS-mediated signaling. Here, we first revealed putative allosteric binding sites of the SOS1, KRASG12C-SOS1 complex, and the ternary KRASG13D-SOS1 complex structures using two network-based models, the essential site scanning analysis and the residue interaction network model. The results enabled us to identify two new putative allosteric pockets for the ternary KRASG13D-SOS1 complex. These were then screened together with the known ligand binding site against the natural compounds in the InterBioScreen (IBS) database using the Glide software package developed by Schrödinger, Inc. The docking poses of seven hit compounds were assessed using 400 ns long molecular dynamics (MD) simulations with two independent replicas using Desmond, coupled with thermal MM-GBSA calculations for the estimation of the binding free energy values. The structural skeleton of the seven proposed compounds consists of different functional groups and heterocyclic rings that possess anti-cancer activity and exhibit persistent interactions with key residues in binding pockets throughout the MD simulations. STOCK1N-09823 was determined as the most promising hit that promoted the disruption of the interactions R73 (chain A)/N879 and R73 (chain A)/Y884, which are key for SOS1-mediated KRAS activation.

Computational insights into the allosteric behavior of mini proinsulin driven by C peptide mobility

The production of recombinant insulin remains challenging, particularly in enhancing refolding efficiency and bioactivity. Mini-proinsulin analogs, which involve reducing the length of the C-peptide, offer potential improvements in insulin production. This study aims to evaluate mini-proinsulin analogs' design and receptor binding dynamics to optimize recombinant insulin production in E. coli. Mini-proinsulin analogs were engineered by replacing the 33-residue C-peptide with a pentapeptide sequence to improve refolding. The three-dimensional structure of mini-proinsulin was predicted using AlphaFold and performed docking analysis of mini-proinsulin analogs to the insulin receptor using AutoDock Tools, with comparisons made to previously available NMR-determined analog and the native insulin-insulin receptor complex. Normal Mode Analyses (GNM and ANM) were performed in detail to assess binding dynamics. In silico analyses revealed that mini-proinsulin analogs closely replicate the structural features of native insulin and display receptor binding dynamics similar to native insulin, though they follow distinct receptor interaction paths. All analysis suggests that C-peptide mobility may contribute to the allosteric behavior observed in mini-proinsulin analogs during receptor interaction.

mTOR Variants Activation Discovers PI3K-like Cryptic Pocket, Expanding Allosteric, Mutant-Selective Inhibitor Designs

mTOR plays a crucial role in PI3K/AKT/mTOR signaling. We hypothesized that mTOR activation mechanisms driving oncogenesis can advise effective therapeutic designs. To test this, we combined cancer genomic analysis with extensive molecular dynamics simulations of mTOR oncogenic variants. We observed that conformational changes within mTOR kinase domain are associated with multiple mutational activation events. The mutations disturb the α-packing formed by the kαAL, kα3, kα9, kα9b, and kα10 helices in the kinase domain, creating cryptic pocket. Its opening correlates with opening of the catalytic cleft, including active site residues realignment, favoring catalysis. The cryptic pocket created by disrupted α-packing coincides with the allosteric pocket in PI3Kα can be harmoniously fitted by the PI3Kα allosteric inhibitor RLY-2608, suggesting that analogous drugs designed based on RLY-2608 can restore the packed α-structure, resulting in mTOR inactive conformation. Our results exemplify that knowledge of detailed kinase activation mechanisms can inform innovative allosteric inhibitor development.

Conformations of KRAS4B Affected by Its Partner Binding and G12C Mutation: Insights from GaMD Trajectory-Image Transformation-Based Deep Learning

Binding of partners and mutations highly affects the conformational dynamics of KRAS4B, which is of significance for deeply understanding its function. Gaussian accelerated molecular dynamics (GaMD) simulations followed by deep learning (DL) and principal component analysis (PCA) were carried out to probe the effect of G12C and binding of three partners NF1, RAF1, and SOS1 on the conformation alterations of KRAS4B. DL reveals that G12C and binding of partners result in alterations in the contacts of key structure domains, such as the switch domains SW1 and SW2 together with the loops L4, L5, and P-loop. Binding of NF1, RAF1, and SOS1 constrains the structural fluctuation of SW1, SW2, L4, and L5; on the contrary, G12C leads to the instability of these four structure domains. The analyses of free energy landscapes (FELs) and PCA also show that binding of partners maintains the stability of the conformational states of KRAS4B while G12C induces greater mobility of the switch domains SW1 and SW2, which produces significant impacts on the interactions of GTP with SW1, L4, and L5. Our findings suggest that partner binding and G12C play important roles in the activity and allosteric regulation of KRAS4B, which may theoretically aid in further understanding the function of KRAS4B.

Molecular Mechanism of Phosphorylation-Mediated Impacts on the Conformation Dynamics of GTP-Bound KRAS Probed by GaMD Trajectory-Based Deep Learning

The phosphorylation of different sites produces a significant effect on the conformational dynamics of KRAS. Gaussian accelerated molecular dynamics (GaMD) simulations were combined with deep learning (DL) to explore the molecular mechanism of the phosphorylation-mediated effect on conformational dynamics of the GTP-bound KRAS. The DL finds that the switch domains are involved in obvious differences in conformation contacts and suggests that the switch domains play a key role in the function of KRAS. The analyses of free energy landscapes (FELs) reveal that the phosphorylation of pY32, pY64, and pY137 leads to more disordered states of the switch domains than the wild-type (WT) KRAS and induces conformational transformations between the closed and open states. The results from principal component analysis (PCA) indicate that principal motions PC1 and PC2 are responsible for the closed and open states of the phosphorylated KRAS. Interaction networks were analyzed and the results verify that the phosphorylation alters interactions of GTP and magnesium ion Mg2+ with the switch domains. It is concluded that the phosphorylation pY32, pY64, and pY137 tune the activity of KRAS through changing conformational dynamics and interactions of the switch domains. We anticipated that this work could provide theoretical aids for deeply understanding the function of KRAS.

Mutations influence the conformational dynamics of the GDP/KRAS complex

Mutations near allosteric sites can have a significant impact on the function of KRAS. Three specific mutations, K104Q, G12D/K104Q, and G12D/G75A, which are located near allosteric positions, were selected to investigate the molecular mechanisms behind mutation-induced influences on the activity of KRAS. Gaussian accelerated molecular dynamics (GaMD) simulations followed by the principal component analysis (PCA) were performed to improve the sampling of conformational states. The results revealed that these mutations significantly alter the structural flexibility, correlated motions, and dynamic behavior of the switch regions that are essential for KRAS binding to effectors or regulators. Furthermore, the mutations have a significant impact on the hydrogen bonding interactions between GDP and the switch regions, as well as on the electrostatic interactions of magnesium ions (Mg2+) with these regions. Our results verified that these mutations strongly influence the binding of KRAS to its effectors or regulators and allosterically regulate the activity. We believe that this work can provide valuable theoretical insights into a deeper understanding of KRAS function.

Impacts of Mutations in the P-Loop on Conformational Alterations of KRAS Investigated with Gaussian Accelerated Molecular Dynamics Simulations

G12 mutations heavily affect conformational transformation and activity of KRAS. In this study, Gaussian accelerated molecular dynamics (GaMD) simulations were performed on the GDP-bound wild-type (WT), G12A, G12D, and G12R KRAS to probe mutation-mediated impacts on conformational alterations of KRAS. The results indicate that three G12 mutations obviously affect the structural flexibility and internal dynamics of the switch domains. The analyses of the free energy landscapes (FELs) suggest that three G12 mutations induce more conformational states of KRAS and lead to more disordered switch domains. The principal component analysis shows that three G12 mutations change concerted motions and dynamics behavior of the switch domains. The switch domains mostly overlap with the binding region of KRAS to its effectors. Thus, the high disorder states and concerted motion changes of the switch domains induced by G12 mutations affect the activity of KRAS. The analysis of interaction network of GDP with KRAS signifies that the instability in the interactions of GDP and magnesium ion with the switch domain SW1 drives the high disordered state of the switch domains. This work is expected to provide theoretical aids for understanding the function of KRAS.

Classification of GTP-dependent K-Ras4B active and inactive conformational states

Classifying reliably active and inactive molecular conformations of wildtype (WT) and mutated oncogenic proteins is a key, ongoing challenge in molecular cancer studies. Here, we probe the GTP-bound K-Ras4B conformational dynamics using long-time atomistic molecular dynamics (MD) simulations. We extract and analyze the detailed underlying free energy landscape of WT K-Ras4B. We use two key reaction coordinates, labeled d₁ and d₂ (i.e., distances coordinating the P_β atom of the GTP ligand with two key residues, T35 and G60), shown to correlate closely with activities of WT and mutated K-Ras4B. However, our new K-Ras4B conformational kinetics study reveals a more complex network of equilibrium Markovian states. We show that a new reaction coordinate is required to account for the orientation of acidic K-Ras4B sidechains such as D38 with respect to the interface with binding effector RAF1 and rationalize the activation/inactivation propensities and the corresponding molecular binding mechanisms. We use this understanding to unveil how a relatively conservative mutation (i.e., D33E, in the switch I region) can lead to significantly different activation propensities compared with WT K-Ras4B. Our study sheds new light on the ability of residues near the K-Ras4B-RAF1 interface to modulate the network of salt bridges at the binding interface with the RAF1 downstream effector and, thus, to influence the underlying GTP-dependent activation/inactivation mechanism. Altogether, our hybrid MD-docking modeling approach enables the development of new in silico methods for quantitative assessment of activation propensity changes (e.g., due to mutations or local binding environment). It also unveils the underlying molecular mechanisms and facilitates the rational design of new cancer drugs.

Binding modes of GDP, GTP and GNP to NRAS deciphered by using Gaussian accelerated molecular dynamics simulations

Probing binding modes of GDP, GTP and GNP to NRAS are of significance for understanding the regulation mechanism on the activity of RAS proteins. Four separate Gaussian accelerated molecular dynamics (GaMD) simulations were performed on the apo, GDP-, GTP- and GNP-bound NRAS. Dynamics analyses suggest that binding of three ligands highly affects conformational states of the switch domains from NRAS, which disturbs binding of NRAS to its effectors. The analyses of free energy landscapes (FELs) indicate that binding of GDP, GTP and GNP induces more energetic states of NRAS compared to the apo NRAS but the presence of GNP makes the switch domains more ordered than binding of GDP and GNP. The information of interaction networks of ligands with NRAS reveals that the π-π interaction of residue F28 and the salt bridge interactions of K16 and D119 with ligands stabilize binding of GDP, GTP and GNP to NRAS. Meanwhile magnesium ion plays a bridge role in interactions of ligands with NRAS, which is favourable for associations of GDP, GTP and GNP with NRAS. This work is expected to provide useful information for deeply understanding the function and activity of NRAS.

Deciphering Conformational Changes of the GDP-Bound NRAS Induced by Mutations G13D, Q61R, and C118S through Gaussian Accelerated Molecular Dynamic Simulations

The conformational changes in switch domains significantly affect the activity of NRAS. Gaussian-accelerated molecular dynamics (GaMD) simulations of three separate replicas were performed to decipher the effects of G13D, Q16R, and C118S on the conformational transformation of the GDP-bound NRAS. The analyses of root-mean-square fluctuations and dynamics cross-correlation maps indicated that the structural flexibility and motion modes of the switch domains involved in the binding of NRAS to effectors are highly altered by the G13D, Q61R, and C118Smutations. The free energy landscapes (FELs) suggested that mutations induce more energetic states in NRAS than the GDP-bound WT NRAS and lead to high disorder in the switch domains. The FELs also indicated that the different numbers of sodium ions entering the GDP binding regions compensate for the changes in electrostatic environments caused by mutations, especially for G13D. The GDP-residue interactions revealed that the disorder in the switch domains was attributable to the unstable hydrogen bonds between GDP and two residues, V29 and D30. This work is expected to provide information on the energetic basis and dynamics of conformational changes in switch domains that can aid in deeply understanding the target roles of NRAS in anticancer treatment.

Modeling the Dynamics of Protein-Protein Interfaces, How and Why?

Protein-protein assemblies act as a key component in numerous cellular processes. Their accurate modeling at the atomic level remains a challenge for structural biology. To address this challenge, several docking and a handful of deep learning methodologies focus on modeling protein-protein interfaces. Although the outcome of these methods has been assessed using static reference structures, more and more data point to the fact that the interaction stability and specificity is encoded in the dynamics of these interfaces. Therefore, this dynamics information must be taken into account when modeling and assessing protein interactions at the atomistic scale. Expanding on this, our review initially focuses on the recent computational strategies aiming at investigating protein-protein interfaces in a dynamic fashion using enhanced sampling, multi-scale modeling, and experimental data integration. Then, we discuss how interface dynamics report on the function of protein assemblies in globular complexes, in fuzzy complexes containing intrinsically disordered proteins, as well as in active complexes, where chemical reactions take place across the protein-protein interface.

Q61 mutant-mediated dynamics changes of the GTP-KRAS complex probed by Gaussian accelerated molecular dynamics and free energy landscapes

Understanding the molecular mechanism of the GTP-KRAS binding is significant for improving the target roles of KRAS in cancer treatment. In this work, multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations were applied to decode the effect of Q61A, Q61H and Q61L on the activity of KRAS. Dynamics analyses based on MR-GaMD trajectory reveal that motion modes and dynamics behavior of the switch domain in KRAS are heavily affected by the three Q61 mutants. Information of free energy landscapes (FELs) shows that Q61A, Q61H and Q61L induce structural disorder of the switch domain and disturb the activity of KRAS. Analysis of the interaction network uncovers that the decrease in the stability of hydrogen bonding interactions (HBIs) of GTP with residues V29 and D30 induced by Q61A, Q61H and Q61L is responsible for the structural disorder of the switch-I and that in the occupancy of the hydrogen bond between GTP and residue G60 leads to the structural disorder of the switch-II. Thus, the high disorder of the switch domain caused by three current Q61 mutants produces a significant effect on binding of KRAS to its effectors. This work is expected to provide useful information for further understanding function and target roles of KRAS in anti-cancer drug development.

Conformational transformation of switch domains in GDP/K-Ras induced by G13 mutants: An investigation through Gaussian accelerated molecular dynamics simulations and principal component analysis

Mutations in K-Ras are involved in a large number of all human cancers, thus, K-Ras is regarded as a promising target for anticancer drug design. Understanding the target roles of K-Ras is important for providing insights on the molecular mechanism underlying the conformational transformation of the switch domains in K-Ras due to mutations. In this study, multiple replica Gaussian accelerated molecular (MR-GaMD) simulations and principal component analysis (PCA) were applied to probe the effect of G13A, G13D and G13I mutations on conformational transformations of the switch domains in GDP-associated K-Ras. The results suggest that G13A, G13D and G13I enhance the structural flexibility of the switch domains, change the correlated motion modes of the switch domains and strengthen the total motion strength of K-Ras compared with the wild-type (WT) K-Ras. Free energy landscape analyses not only show that the switch domains of the GDP-bound inactive K-Ras mainly exist as a closed state but also indicate that mutations evidently alter the free energy profile of K-Ras and affect the conformational transformation of the switch domains between the closed and open states. Analyses of hydrophobic interaction contacts and hydrogen bonding interactions show that the mutations scarcely change the interaction network of GDP with K-Ras and only disturb the interaction of GDP with the switch (SW1). In summary, two newly introduced mutations, G13A and G13I, play similar adjustment roles in the conformational transformations of two switch domains to G13D and are possibly utilized to tune the activity of K-Ras and the binding of guanine nucleotide exchange factors.