Microscopic Mechanism and Kinetics of Ice Formation at Complex Interfaces: Zooming in on Kaolinite

Most ice in nature forms because of impurities which boost the exceedingly low nucleation rate of pure supercooled water. However, the microscopic details of ice nucleation on these substances remain largely unknown. Here, we have unraveled the molecular mechanism and the kinetics of ice formation on kaolinite, a clay mineral playing a key role in climate science. We find that the formation of ice at strong supercooling in the presence of this clay is about 20 orders of magnitude faster than homogeneous freezing. The critical nucleus is substantially smaller than that found for homogeneous nucleation and, in contrast to the predictions of classical nucleation theory (CNT), it has a strong two-dimensional character. Nonetheless, we show that CNT describes correctly the formation of ice at this complex interface. Kaolinite also promotes the exclusive nucleation of hexagonal ice, as opposed to homogeneous freezing where a mixture of cubic and hexagonal polytypes is observed.

A disappearing act performed by magnesium: the nucleotide exchange mechanism of Ran GTPase by Quantum mechanics/Molecular mechanics studies

The Ran GTPase protein is a guanine nucleotide-binding protein (GNBP) with an acknowledged profile in cancer onset, progression and metastases. The complex mechanism adopted by GNBPs in exchanging GDP for GTP is an intriguing process and crucial for Ran viability. The successful completion of the process is a fundamental aspect of propagating downstream signalling events. QM/MM molecular dynamics simulations were employed in this study to provide a deeper mechanistic understanding of the initiation of nucleotide exchange in Ran. Results indicate significant disruption of the metal-binding site upon interaction with RCC1 (the Ran guanine nucleotide exchange factor), overall culminating in the prominent shift of the divalent magnesium ion. The observed ion drifting is reasoned to occur as a consequence of the complex formation between Ran and RCC1 and is postulated to be a critical factor in the exchange process adopted by Ran. This is the first report to observe and detail such intricate dynamics for a protein in Ras superfamily.

Self-healing of cracks during ductile regime machining of silicon: Insights from molecular dynamics simulation

During nanoindentation and ductile-regime machining of silicon, a phenomenon known as “self-healing” takes place in that the microcracks, microfractures, and small spallings generated during the machining are filled by the plastically flowing ductile phase of silicon. However, this phenomenon has not been observed in simulation studies. In this work, using a long-range potential function, molecular dynamics simulation was used to provide an improved explanation of this mechanism. A unique phenomenon of brittle cracking was discovered, typically inclined at an angle of 45° to 55° to the cut surface, leading to the formation of periodic arrays of nanogrooves being filled by plastically flowing silicon during cutting. This observation is supported by the direct imaging. The simulated X-ray diffraction analysis proves that in contrast to experiments, Si-I to Si-II (beta tin) transformation during ductile-regime cutting is highly unlikely and solid-state amorphisation of silicon caused solely by the machining stress rather than the cutting temperature is the key to its brittle-ductile transition observed during the MD simulations

Intraclonal heterogeneity and distinct molecular mechanisms characterize the development of t(4;14) and t(11;14) myeloma.

We have used whole exome sequencing to compare a group of presentation t(4;14) with t(11;14) cases of myeloma to define the mutational landscape. Each case was characterized by a median of 24.5 exonic nonsynonymous single-nucleotide variations, and there was a consistently higher number of mutations in the t(4;14) group, but this number did not reach statistical significance. We show that the transition and transversion rates in the 2 subgroups are similar, suggesting that there was no specific mechanism leading to mutation differentiating the 2 groups. Only 3% of mutations were seen in both groups, and recurrently mutated genes include NRAS, KRAS, BRAF, and DIS3 as well as DNAH5, a member of the axonemal dynein family. The pattern of mutation in each group was distinct, with the t(4;14) group being characterized by deregulation of chromatin organization, actin filament, and microfilament movement. Recurrent RAS pathway mutations identified subclonal heterogeneity at a mutational level in both groups, with mutations being present as either dominant or minor subclones. The presence of subclonal diversity was confirmed at a single-cell level using other tumor-acquired mutations. These results are consistent with a distinct molecular pathogenesis underlying each subgroup and have important impacts on targeted treatment strategies. The Medical Research Council Myeloma IX trial is registered under ISRCTN68454111.

Phosphorylation Mechanism of Phosphomevalonate Kinase: Implications for Rational Engineering of Isoprenoid Biosynthetic Pathway Enzymes

Mevalonate pathway is of important clinical, pharmaceutical and biotechnological relevance. However, lack of the understanding of the phosphorylation mechanism of the kinases in this pathway has limited rationally engineering the kinases in industry. Here the phosphorylation reaction mechanism of a representative kinase in the mevalonate pathway, phosphomevalonate kinase, was studied by using molecular dynamics and hybrid QM/MM methods. We find that a conserved residue (Ser106) is reorientated to anchor ATP via a stable H-bond interaction. In addition, Ser213 located on the α-helix at the catalytic site is repositioned to further approach the substrate, facilitating the proton transfer during the phosphorylation. Furthermore, we elucidate that Lys101 functions to neutralize the negative charge developed at the β-, γ-bridging oxygen atom of ATP during phosphoryl transfer. We demonstrate that the dissociative catalytic reaction occurs via a direct phosphorylation pathway. This is the first study on the phosphorylation mechanism of a mevalonate pathway kinase. The elucidation of the catalytic mechanism not only sheds light on the common catalytic mechanism of GHMP kinase superfamily, but also provides the structural basis for engineering the mevalonate pathway kinases to further exploit their applications in the production of a wide range of fine chemicals such as biofuels or pharmaceuticals. 

Insight into the mechanism of galactokinase: role of a critical glutamate residue and helix/coil transitions

Galactokinase, the enzyme which catalyses the first committed step in the Leloir pathway, has attracted interest due to its potential as a biocatalyst and as a possible drug target in the treatment of type I galactosemia. The mechanism of the enzyme is not fully elucidated. Molecular dynamics (MD) simulations of galactokinase with the active site residues Arg-37 and Asp-186 altered predicted that two regions (residues 174-179 and 231-240) had different dynamics as a consequence. Interestingly, the same two regions were also affected by alterations in Arg-105, Glu-174 and Arg- 228. These three residues were identified as important in catalysis in previous computational studies on human galactokinase. Alteration of Arg-105 to methionine resulted in a modest reduction in activity with little change in stability. When Arg-228 was changed to methionine, the enzyme’s interaction with both ATP and galactose was affected. This variant was significantly less stable than the wild-type protein. Changing Glu-174 to glutamine (but not to aspartate) resulted in no detectable activity and a less stable enzyme. Overall, these combined in silico and in vitro studies demonstrate the importance of a negative charge at position 174 and highlight the critical role of the dynamics in to key regions of the protein. We postulate that these regions may be critical for mediating the enzyme’s structure and function.