The Plasmodium falciparum Malaria M1 Alanyl Aminopeptidase (PfA-M1): Insights of Catalytic Mechanism and Function from MD Simulations

Malaria caused by several species of Plasmodium is major parasitic disease of humans, causing 1-3 million deaths worldwide annually. The widespread resistance of the human parasite to current drug therapies is of major concern making the identification of new drug targets urgent. While the parasite grows and multiplies inside the host erythrocyte it degrades the host cell hemoglobin and utilizes the released amino acids to synthesize its own proteins. The P. falciparum malarial M1 alanyl-aminopeptidase (PfA-M1) is an enzyme involved in the terminal stages of hemoglobin digestion and the generation of an amino acid pool within the parasite. The enzyme has been validated as a potential drug target since inhibitors of the enzyme block parasite growth in vitro and in vivo. In order to gain further understanding of this enzyme, molecular dynamics simulations using data from a recent crystal structure of PfA-M1 were performed. The results elucidate the pentahedral coordination of the catalytic Zn in these metallo-proteases and provide new insights into the roles of this cation and important active site residues in ligand binding and in the hydrolysis of the peptide bond. Based on the data, we propose a two-step catalytic mechanism, in which the conformation of the active site is altered between the Michaelis complex and the transition state. In addition, the simulations identify global changes in the protein in which conformational transitions in the catalytic domain are transmitted at the opening of the N-terminal 8 angstrom-long channel and at the opening of the 30 angstrom-long C-terminal internal chamber that facilitates entry of peptides to the active site and exit of released amino acids. The possible implications of these global changes with regard to enzyme function are discussed.

Free energy and molecular dynamics calculations for the cubic-tetragonal phase transition in zirconia

The high-temperature cubic-tetragonal phase transition of pure stoichiometric zirconia is studied by molecular dynamics (MD) simulations and within the framework of the Landau theory of phase transformations. The interatomic forces are calculated using an empirical, self-consistent, orthogonal tight-binding model, which includes atomic polarizabilities up to the quadrupolar level. A first set of standard MD calculations shows that, on increasing temperature, one particular vibrational frequency softens. The temperature evolution of the free-energy surfaces around the phase transition is then studied with a second set of calculations. These combine the thermodynamic integration technique with constrained MD simulations. The results seem to support the thesis of a second-order phase transition but with unusual, very anharmonic behavior above the transition temperature.

Kinetics of phase transformations in a model with metastable fluid-fluid separation: A molecular dynamics study

By molecular dynamics (MD) simulations we study the crystallization process in a model system whose particles interact by a spherical pair potential with a narrow and deep attractive well adjacent to a hard repulsive core. The phase diagram of the model displays a solid-fluid equilibrium, with a metastable fluid-fluid separation. Our computations are restricted to fairly small systems (from 2592 to 10368 particles) and cover long simulation times, with constant energy trajectories extending up to 76x10(6) MD steps. By progressively reducing the system temperature below the solid-fluid line, we first observe the metastable fluid-fluid separation, occurring readily and almost reversibly upon crossing the corresponding line in the phase diagram. The nucleation of the crystal phase takes place when the system is in the two-fluid metastable region. Analysis of the temperature dependence of the nucleation time allows us to estimate directly the nucleation free energy barrier. The results are compared with the predictions of classical nucleation theory. The critical nucleus is identified, and its structure is found to be predominantly fcc. Following nucleation, the solid phase grows steadily across the system, incorporating a large number of localized and extended defects. We discuss the relaxation processes taking place both during and after the crystallization stage. The relevance of our simulation for the kinetics of protein crystallization under normal experimental conditions is discussed. (C) 2002 American Institute of Physics.

Small angle neutron scattering from 1-alkyl-3-methylimidazolium hexafluorophosphate ionic liquids ([C(n)mim][PF6], n=4, 6, and 8)

The presence of local anisotropy in the bulk, isotropic, and ionic liquid phases-leading to local mesoscopic inhomogeneity-with nanoscale segregation and expanding nonpolar domains on increasing the length of the cation alkyl-substituents has been proposed on the basis of molecular dynamics (MD) simulations. However, there has been little conclusive experimental evidence for the existence of intermediate mesoscopic structure between the first/second shell correlations shown by neutron scattering on short chain length based materials and the mesophase structure of the long chain length ionic liquid crystals. Herein, small angle neutron scattering measurements have been performed on selectively H/D-isotopically substituted 1-alkyl-3-methylimidazolium hexafluorophosphate ionic liquids with butyl, hexyl, and octyl substituents. The data show the unambiguous existence of a diffraction peak in the low-Q region for all three liquids which moves to longer distances (lower Q), sharpens, and increases in intensity with increasing length of the alkyl substituent. It is notable, however, that this peak occurs at lower values of Q (longer length scale) than predicted in any of the previously published MD simulations of ionic liquids, and that the magnitude of the scattering from this peak is comparable with that from the remainder of the amorphous ionic liquid. This strongly suggests that the peak arises from the second coordination shells of the ions along the vector of alkyl-chain substituents as a consequence of increasing the anisotropy of the cation, and that there is little or no long-range correlated nanostructure in these ionic liquids.

Molecular dynamics studies of the STAT3 homodimer:DNA complex: relationships between STAT3 mutations and protein-DNA recognition.

Signal Transducers and Activators of Transcription (STAT) proteins are a group of latent cytoplasmic transcription factors involved in cytokine signaling. STAT3 is a member of the STAT family and is expressed at elevated levels in a large number of diverse human cancers and is now a validated target for anticancer drug discovery.. Understanding the dynamics of the STAT3 dimer interface, accounting for both protein-DNA and protein-protein interactions, with respect to the dynamics of the latent unphosphorylated STAT3 monomer, is important for designing potential small-molecule inhibitors of the activated dimer. Molecular dynamics (MD) simulations have been used to study the activated STAT3 homodimer:DNA complex and the latent unphosphorylated STAT3 monomer in an explicit water environment. Analysis of the data obtained from MD simulations over a 50 ns time frame has suggested how the transcription factor interacts with DNA, the nature of the conformational changes, and ways in which function may be affected. Examination of the dimer interface, focusing on the protein-DNA interactions, including involvement of water molecules, has revealed the key residues contributing to the recognition events involved in STAT3 protein-DNA interactions. This has shown that the majority of mutations in the DNA-binding domain are found at the protein-DNA interface. These mutations have been mapped in detail and related to specific protein-DNA contacts. Their structural stability is described, together with an analysis of the model as a starting-point for the discovery of novel small-molecule STAT3 inhibitors.

Anisotropy of single crystal 3C-SiC during nanometric cutting

3C–SiC (the only polytype of SiC that resides in a diamond cubic lattice structure) is a relatively new material that exhibits most of the desirable engineering properties required for advanced electronic applications. The anisotropy exhibited by 3C–SiC during its nanometric cutting is significant, and the potential for its exploitation has yet to be fully investigated. This paper aims to understand the influence of crystal anisotropy of 3C–SiC on its cutting behaviour. A molecular dynamics simulation model was developed to simulate the nanometric cutting of single-crystal 3C–SiC in nine (9) distinct combinations of crystal orientations and cutting directions, i.e. (1?1?1) <-1?1?0>, (1?1?1) <-2?1?1>, (1?1?0) <-1?1?0>, (1?1?0) <0?0?1>, (1?1?0) <1?1?-2>, (0?0?1) <-1?1?0>, (0?0?1) <1?0?0>, (1?1?-2) <1?-1?0> and (1?-2?0) <2?1?0>.

In order to ensure the reliability of the simulation results, two separate simulation trials were carried out with different machining parameters. In the first trial, a cutting tool rake angle of -25°, d/r (uncut chip thickness/cutting edge radius) ratio of 0.57 and cutting velocity of 10 m s-1 were used whereas a second trial was done using a cutting tool rake angle of -30°, d/r ratio of 1 and cutting velocity of 4 m s-1. Both the trials showed similar anisotropic variation.

The simulated orthogonal components of thrust force in 3C–SiC showed a variation of up to 45%, while the resultant cutting forces showed a variation of 37%. This suggests that 3C–SiC is highly anisotropic in its ease of deformation. These results corroborate with the experimentally observed anisotropic variation of 43.6% in Young's modulus of 3C–SiC. The recently developed dislocation extraction algorithm (DXA) [1, 2] was employed to detect the nucleation of dislocations in the MD simulations of varying cutting orientations and cutting directions. Based on the overall analysis, it was found that 3C–SiC offers ease of deformation on either (1?1?1) <-1?1?0>, (1?1?0) <0?0?1>, or (1?0?0) <1?0?0> setups.

Molecular dynamics simulation model for the quantitative assessment of tool wear during single point diamond turning of cubic silicon carbide

Silicon carbide (SiC) is a material of great technological interest for engineering applications concerning hostile environments where silicon-based components cannot work (beyond 623 K). Single point diamond turning (SPDT) has remained a superior and viable method to harness process efficiency and freeform shapes on this harder material. However, it is extremely difficult to machine this ceramic consistently in the ductile regime due to sudden and rapid tool wear. It thus becomes non trivial to develop an accurate understanding of tool wear mechanism during SPDT of SiC in order to identify measures to suppress wear to minimize operational cost.

In this paper, molecular dynamics (MD) simulation has been deployed with a realistic analytical bond order potential (ABOP) formalism based potential energy function to understand tool wear mechanism during single point diamond turning of SiC. The most significant result was obtained using the radial distribution function which suggests graphitization of diamond tool during the machining process. This phenomenon occurs due to the abrasive processes between these two ultra hard materials. The abrasive action results in locally high temperature which compounds with the massive cutting forces leading to sp3–sp2 order–disorder transition of diamond tool. This represents the root cause of tool wear during SPDT operation of cubic SiC. Further testing led to the development of a novel method for quantitative assessment of the progression of diamond tool wear from MD simulations.

The development of a surface defect machining method for hard turning processes

Hard turning (HT) is a material removal process employing a combination of a single point cutting tool and high speeds to machine hard ferrous alloys which exhibit hardness values over 45 HRC. In this paper, a surface defect machining (SDM) method for HT is proposed which harnesses the combined advantages of porosity machining and pulsed laser pre-treatment processing. From previous experimental work, this was shown to provide better controllability of the process and improved quality of the machined surface. While the experiments showed promising results, a comprehensive understanding of this new technique could only be achieved through a rigorous, in depth theoretical analysis. Therefore, an assessment of the SDM technique was carried out using both finite element method (FEM) and molecular dynamics (MD) simulations.
FEM modelling was used to compare the conventional HT of AISI 4340 steel (52 HRC) using an Al2O3 insert with the proposed SDM method. The simulations showed very good agreement with the previously published experimental results. Compared to conventional HT, SDM provided favourable machining outcomes, such as reduced shear plane angle, reduced average cutting forces, improved surface roughness, lower residual stresses on the machined surface, reduced tool–chip interface contact length and increased chip flow velocity. Furthermore, a scientific explanation of the improved surface finish was revealed using a state-of-the-art MD simulation model which suggested that during SDM, a combination of both the cutting action and rough polishing action help improve the machined surface finish.

Influence of test methodology and probe geometry on nanoscale fatigue mechanisms of diamond-like carbon thin film

The aim of this paper is to investigate the mechanism of nanoscale fatigue using nano-impact and multiple-loading cycle nanoindentation tests, and compare it to previously reported findings of nanoscale fatigue using integrated stiffness and depth sensing approach. Two different film loading mechanism, loading history and indenter shapes are compared to comprehend the influence of test methodology on the nanoscale fatigue failure mechanisms of DLC film. An amorphous 100 nm thick DLC film was deposited on a 500 μm silicon substrate using sputtering of graphite target in pure argon atmosphere. Nano-impact and multiple-load cycle indentations were performed in the load range of 100 μN to 1000 μN and 0.1 mN to 100 mN, respectively. Both test types were conducted using conical and Berkovich indenters. Results indicate that for the case of conical indenter, the combination of nano-impact and multiple-loading cycle nanoindentation tests provide information on the life and failure mechanism of DLC film, which is comparable to the previously reported findings using the integrated stiffness and depth sensing approach. However, the comparison of results is sensitive to the applied load, loading mechanism, test-type and probe geometry. The loading mechanism and load history is therefore critical which also leads to two different definitions of film failure. The choice of exact test methodology, load and probe geometry should therefore be dictated by the in-service tribological conditions, and where necessary both test methodologies can be used to provide better insights of failure mechanism. Molecular dynamics (MD) simulations of the elastic response of nanoindentation is reported, which indicates that the elastic modulus of the film measured using MD simulation was higher than that experimentally measured. This difference is attributed to the factors related to the presence of material defects, crystal structure, residual stress, indenter geometry and loading/unloading rate differences between the MD and experimental results.

Diamond machining of silicon: A review of advances in molecular dynamics simulation

Molecular dynamics (MD) simulation has enhanced our understanding about ductile-regime machining of brittle materials such as silicon and germanium. In particular, MD simulation has helped understand the occurrence of brittle–ductile transition due to the high-pressure phase transformation (HPPT), which induces Herzfeld–Mott transition. In this paper, relevant MD simulation studies in conjunction with experimental studies are reviewed with a focus on (i) the importance of machining variables: undeformed chip thickness, feed rate, depth of cut, geometry of the cutting tool in influencing the state of the deviatoric stresses to cause HPPT in silicon, (ii) the influence of material properties: role of fracture toughness and hardness, crystal structure and anisotropy of the material, and (iii) phenomenological understanding of the wear of diamond cutting tools, which are all non-trivial for cost-effective manufacturing of silicon. The ongoing developmental work on potential energy functions is reviewed to identify opportunities for overcoming the current limitations of MD simulations. Potential research areas relating to how MD simulation might help improve existing manufacturing technologies are identified which may be of particular interest to early stage researchers.

Molecular dynamics simulation (MDS) to study nanoscale machining processes

Molecular Dynamics Simulations (MDS) are constantly being used to make important contributions to our fundamental understanding of material behaviour, at the atomic scale, for a variety of thermodynamic processes. This chapter shows that molecular dynamics simulation is a robust numerical analysis tool in addressing a range of complex nanofinishing (machining) problems that are otherwise difficult or impossible to understand using other methods. For example the mechanism of nanometric cutting of silicon carbide is influenced by a number of variables such as machine tool performance, machining conditions, material properties, and cutting tool performance (material microstructure and physical geometry of the contact) and all these variables cannot be monitored online through experimental examination. However, these could suitably be studied using an advanced simulation based approach such as MDS. This chapter details how MD simulation can be used as a research and commercial tool to understand key issues of ultra precision manufacturing research problems and a specific case was addressed by studying diamond machining of silicon carbide. While this is appreciable, there are a lot of challenges and opportunities in this fertile area. For example, the world of MD simulations is dependent on present day computers and the accuracy and reliability of potential energy functions [109]. This presents a limitation: Real-world scale simulation models are yet to be developed. The simulated length and timescales are far shorter than the experimental ones which couples further with the fact that contact loading simulations are typically done in the speed range of a few hundreds of m/sec against the experimental speed of typically about 1 m/sec [17]. Consequently, MD simulations suffer from the spurious effects of high cutting speeds and the accuracy of the simulation results has yet to be fully explored. The development of user-friendly software could help facilitate molecular dynamics as an integral part of computer-aided design and manufacturing to tackle a range of machining problems from all perspectives, including materials science (phase of the material formed due to the sub-surface deformation layer), electronics and optics (properties of the finished machined surface due to the metallurgical transformation in comparison to the bulk material), and mechanical engineering (extent of residual stresses in the machined component) [110]. Overall, this chapter provided key information concerning diamond machining of SiC which is classed as hard, brittle material. From the analysis presented in the earlier sections, MD simulation has helped in understanding the effects of crystal anisotropy in nanometric cutting of 3C-SiC by revealing the atomic-level deformation mechanisms for different crystal orientations and cutting directions. In addition to this, the MD simulation revealed that the material removal mechanism on the (111) surface of 3C-SiC (akin to diamond) is dominated by cleavage. These understandings led to the development of a new approach named the “surface defect machining” method which has the potential to be more effective to implement than ductile mode micro laser assisted machining or conventional nanometric cutting.

Ferroelectric Vortices and Related Configurations

This chapter discusses that the theoretical studies, using both atomistic and phenomenological approaches, have made clear predictions about the existence and behaviour of ferroelectric (FE) vortices. Effective Hamiltonians can be implemented within both Monte Carlo (MC) and molecular dynamics (MD) simulations. In contrast to the effective Hamiltonian method, which is atomistic in nature, the phase field method employs a continuum approach, in which the polarization field is the order parameter. Properties of FE nanostructures are largely governed by the existence of a depolarization field, which is much stronger than the demagnetization field in magnetic nanosystems. The topological patterns seen in rare earth manganites are often referred to as vortices and yet this claim never seems to be explicitly justified. By inspection, the form of a vortex structure is such that there is a continuous rotation in the orientation of dipole vectors around the singularity at the centre of the vortex.

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

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.