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.

Cation complexation by chemically modified calixarenes. 11. Complexation and extraction of alkali cations by calix[5]- and -[6]arene ketones. Crystal and molecular structures of calix[5]arene ketones and Na+ and Rb+ complexes

A series of four calix[5]arenes and three calix[6]arenes (R-calixarene-OCH2COR1) (R = H or Bu-t) with alkyl ketone residues (R-1 = Me or Bu-t) on the lower rim have been synthesized, and their affinity for complexation of alkali cations has been assessed through phase-transfer experiments and stability constant measurements. The conformations of these ketones have been probed by H-1 NMR and X-ray diffraction analysis, and by molecular mechanics calculations. Pentamer 3 (R R-1 = Bu-t) possesses a symmetrical cone conformation in solution and a very distorted cone conformation in the solid state. Pentamer 5 (R = H, R-1 = Bu-t) exists in a distorted 1,2-alternate conformation in the solid state, but in solution two slowly interconverting conformations, one a cone and the other presumed to be 1,2-alternate, can be detected. X-ray structure analysis of the sodium and rubidium perchlorate complexes of 3 reveal the cations deeply encapsulated by the ethereal and carbonyl oxygen atoms in distorted cone conformations which can be accurately reproduced by molecular mechanics calculations. The phase-transfer and stability constant data reveal that the extent of complexation depends on calixarene size and the nature of the alkyl residues adjacent to the ketonic carbonyls with tert-butyl much more efficacious than methyl.

Molecular modeling of calixarenes with group I metal ions

Molecular mechanics calculations have been used to model the geometries of the complexes of Group I metal ions with calix[n]arenes (n = 4,5). A simple procedure in which the calixarene atoms are assigned partial charges on the basis of AM1 calculations and the metal ions are allowed to bind electrostatically to the calixarenes produces surprising good results when the resulting structures are compared to known crystallographic data on the complexes. Encapsulated solvent molecules and/or counterions can be included in the calculations and, indeed, are necessary to reproduce the X-ray data.

The Effect of Tacticity on the Conformational Properties of Poly(1-olefin sulfone)s

The effect of tacticity on the conformational properties of poly(olefin sulfone)s was studied. Tactic polymers, prepared from racemic thiirane monomers using chiral inititators were compared with atactic polymers prepared by free radical co-polymerisation of the 1-olefin with sulfur dioxide. Analysis of the XRD patterns showed that the tactic polymers formed more ordered structures in the bulk with longer layer spacings, consistent with a model in which their side chains meet at the tips in contrast with the atactic polymers whose side chains interdigitate. 13C MAS nmr experiments suggest that as tacticity increases so too does the proportion of C-S bonds in the gauche conformation, however the proportion of S-C bonds in the trans conformation falls, in contrast to a reported molecular mechanics study. Finally, DSC measurements on the polymers with longer side chains showed the presence of two endotherms on heating, illustrating definite liquid crystalline behaviour.

Suggested Binding Mechanism of the HIV-gp120 to its CD4 Receptor

The molecular recognition and attachment of the CD4 molecule and the HIV envelope glycoprotein (gp120) might be described as a consecutive three-step molecular recognition process. 1. (a) Long range interaction: electrostatic pre-orientation, 2. (b) short range interaction: electronic attachment followed by a ‘Locking-in’ (via aromatic ring orientation) and 3. (c) internal interaction (induced fit): conformational readjustment of the protein molecules. On the basis of the preliminary investigations (X-ray structures of CD4 and biological studies of CD4 and gp120 point mutants) we described a computational model. This approach consists of empirical calculations as well as ab initio level of quantum chemistry. The conformational analysis of the wild type and mutant CD4 molecules was supported by molecular mechanics and dynamics (Amber force field). The latter analysis involves the application of a novel method, the Amino Acid Conformation Assignment of Proteins (ACAP) software, developed for the notation of secondary protein structures. According to the cardinal role of the electrostatic factors during this interaction, several ab initio investigations were performed for better understanding of the recognition process on submolecular level. Using the above mentioned computational model, we could interpret the basic behaviours and predict some additional features of CD4-gp120 interaction, in spite of the missing gp120 X-ray structure.

Modelling of the sodium complex of a calixarene tetraester in the 1,3-alternate conformation

This paper describes the use of molecular mechanics to model the geometry of the sodium complex of a calix[4] arene tetraester, in the 1,3-alternate conformation 1. Partial charges were assigned to the calixarene on the basis of semi-empirical (AM1, PM3, MNDO, INDO, CNDO and ZINDO) calculations and the binding of the sodium ion to the calixarene was modelled using molecular mechanics. Agreement between the optimised and X-ray structures of the complex was very good. The effect of placing the cation in different starting positions on the energy-minimised geometry of the complex is described.

Structural and kinetic studies of spin crossover in an Iron(II) complex with a novel tripodal ligand

Configurational and ligand conformational influences on the kinetics of (1)A(1) reversible arrow T-5(2) spin crossover in the Fe(II) complex with the novel tripodal ligand, 1,1,1-tris((N-(2-pyridylmethyl)-N-methylamino)methyl)ethane (tptMetame), have been explored. Despite having six chelate rings and three chiral nitrogen atoms, only one enantiomeric pair of isomers, Delta, SSS, and Lambda, RRR, of the complex ion is observed. The conformation of the three rings forming the upper ''cap'' of the complex structure can be assigned delta or lambda with respect to the 3-fold molecular axis. X-ray data at 300 and 153 K, above and below the critical temperature for the spin transition, show that the conformation of the ligand ''cap'' is the same as the absolute configuration of the complex, with the same Lambda lambda(CAP)(or Delta delta(CAP)) combination prevailing for both the LS ((1)A(1)) and HS (T-5(2)) isomers. Molecular mechanics calculations further show that the ligand energy remains lowest for this Lambda lambda(CAP) (or Delta delta(CAP)) combination at all Fe-N distances over the range spanning the LS and HS isomers. Measurements of the spin crossover relaxation time have been carried out in solution over the temperature range 293-170 K. The observed monophasic relaxation traces are also consistent with the absolute configuration of the complex remaining unaltered during the spin crossover.

The role of the active site residues in human galactokinase: implications for the mechanisms of GHMP kinases.

Galactokinase catalyses the phosphorylation of galactose at the expense of ATP. Like other members of the GHMP family of kinases it is postulated to function through an active site base mechanism in which Asp-186 abstracts a proton from galactose. This asparate residue was altered to alanine and to asparagine by site-directed mutagenesis of the corresponding gene. This resulted in variant enzyme with no detectable galactokinase activity. Alteration of Arg-37, which lies adjacent to Asp-186 and is postulated to assist the catalytic base, to lysine resulted in an active enzyme. However, alteration of this residue to glutamate abolished activity. All the variant enzymes, except the arginine to lysine substitution, were structurally unstable (as judged by native gel electrophoresis in the presence of urea) compared to the wild type. This suggests that the lack of activity results from this structural instability, in addition to any direct effects on the catalytic mechanism. Computational estimations of the pK(a) values of the arginine and aspartate residues, suggest that Arg-37 remains protonated throughout the catalytic cycle whereas Asp-186 has an abnormally high pK(a) value (7.18). Quantum mechanics/molecular mechanics (QM/MM) calculations suggest that Asp-186 moves closer to the galactose molecule during catalysis. The experimental and theoretical studies presented here argue for a mechanism in which the C-1-OH bond in the sugar is weakened by the presence of Asp-186 thus facilitating nucleophilic attack by the oxygen atom on the gamma-phosphorus of ATP.

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.

Molecular dynamics simulation of nanoindentation of Fe3C and tetrahedral Fe4C

Study of nanomechanical response of iron carbides is important because presence of iron carbides greatly influences the performance and longevity of steel components. This work contributes to the literature by exploring nanoindentation of nanocrystalline Fe3C and tetrahedral-Fe4C using molecular dynamics simulation. The chemical interactions of iron and carbon were described through an analytical bond order inter-atomic potential (ABOP) energy function. The indentations were performed at an indentation speed of 50 m/sec and a repeat trial was performed at 5 m/sec. Load-displacement (P-h) curve for both these carbides showed residual indentation depth and maximum indentation depth (hf/hmax) ratio to be higher than 0.7 i.e. a circumstance where Oliver and Pharr method was not appropriate to be applied to evaluate the material properties. Alternate evaluation revealed Fe3C to be much harder than Fe4C. Gibbs free energy of formation and radial distribution function, coupled with state of the average local temperature and von Mises stresses indicate the formation of a new phase of iron-carbide. Formation of this newer phase was found to be due to deviatoric strain rather than the high temperature induced in the substrate during nanoindentation