Exploring and Engineering the Conformational Landscape of Calmodulin through Specific Interactions

Role of the Residue Q1919 in Increasing Kinase Activity of G2019S LRRK2 Kinase: A Computational Study

Mutations in multi-domain leucine-rich repeat kinase 2 (LRRK2) have been an interest to researchers as these mutations are associated with Parkinson's disease. G2019S mutation in LRRK2 kinase domain leads to the formation of additional hydrogen bonds by S2019 which results in stabilization of the active state of the kinase, thereby increasing kinase activity. Two additional hydrogen bonds of S2019 are reported separately. Here, a mechanistic picture of the formation of additional hydrogen bonds of S2019 with Q1919 (also with E1920) is presented using 'active' Roco4 kinase as a homology model and its relationship with the stabilization of the 'active' G2019S LRRK2 kinase. A conformational flipping of residue Q1919 was found which helped to form stable hydrogen bond with S2019 and made 'active' state more stable in G2019S LRRK2. Two different states were found within the 'active' kinase with respect to the conformational change (flipping) in Q1919. Two doubly-mutated systems, G2019S/Q1919A and G2019S/E1920 K, were studied separately to check the effect of Q1919 and E1920. For both cases, the stable S2 state was not formed, leading to a decrease in kinase activity. These results indicate that both the additional hydrogen bonds of S2019 (with Q1919 and E1920) are necessary to stabilize the active G2019S LRRK2.

Inter-Domain Repulsion of Dumbbell-Shaped Calmodulin during Electrospray Ionization Revealed by Molecular Dynamics Simulations

The mechanisms whereby protein ions are released from nanodroplets at the liquid-gas interface have continued to be controversial since electrospray ionization (ESI) mass spectrometry was widely applied in biomolecular structure analysis in solution. Several viable pathways have been proposed and verified for single-domain proteins. However, the ESI mechanism of multi-domain proteins with more complicated and flexible structures remains unclear. Herein, dumbbell-shaped calmodulin was chosen as a multi-domain protein model to perform molecular dynamics simulations to investigate the structural evolution during the ESI process. For [Ca₄CAM], the protein followed the classical charge residue model. As the inter-domain electrostatic repulsion increased, the droplet was found to split into two sub-droplets, while stronger-repulsive apo-calmodulin unfolded during the early evaporation stage. We designated this novel ESI mechanism as the domain repulsion model, which provides new mechanistic insights into further exploration of proteins containing more domains. Our results suggest that greater attention should be paid to the effect of domain-domain interactions on structure retention during liquid-gas interface transfer when mass spectrometry is used as the developing technique in gas phase structural biology.

Comparative study of the unbinding process of some HTLV-1 protease inhibitors using unbiased molecular dynamics simulations

The HTLV-1 protease is one of the major antiviral targets to overwhelm this virus. Several research groups have developed protease inhibitors, but none has been successful. In this regard, developing new HTLV-1 protease inhibitors to fix the defects in previous inhibitors may overcome the lack of curative treatment for this oncovirus. Thus, we decided to study the unbinding pathways of the most potent (compound 10, PDB ID 4YDF, Ki = 15 nM) and one of the weakest (compound 9, PDB ID 4YDG, Ki = 7900 nM) protease inhibitors, which are very structurally similar. We conducted 12 successful short and long simulations (totaling 14.8 μs) to unbind the compounds from two monoprotonated (mp) forms of protease using the Supervised Molecular Dynamics (SuMD) without applying any biasing force. The results revealed that Asp32 or Asp32′ in the two forms of mp state similarly exert powerful effects on maintaining both potent and weak inhibitors in the binding pocket of HTLV-1 protease. In the potent inhibitor’s unbinding process, His66′ was a great supporter that was absent in the weak inhibitor’s unbinding pathway. In contrast, in the weak inhibitor’s unbinding process, Trp98/Trp98′ by pi-pi stacking interactions were unfavorable for the stability of the inhibitor in the binding site. In our opinion, these results will assist in designing more potent and effective inhibitors for the HTLV-1 protease.