Characterization of the structure of RAMP1 by mutagenesis and molecular modeling

Receptor activity modifying proteins (RAMPs) are a family of single-pass transmembrane proteins that dimerize with G-protein-coupled receptors. They may alter the ligand recognition properties of the receptors (particularly for the calcitonin receptor-like receptor, CLR). Very little structural information is available about RAMPs. Here, an ab initio model has been generated for the extracellular domain of RAMP1. The disulfide bond arrangement (Cys 27-Cys82, Cys40-Cys72, and Cys 57-Cys104) was determined by site-directed mutagenesis. The secondary structure (a-helices from residues 29-51, 60-80, and 87-100) was established from a consensus of predictive routines. Using these constraints, an assemblage of 25,000 structures was constructed and these were ranked using an all-atom statistical potential. The best 1000 conformations were energy minimized. The lowest scoring model was refined by molecular dynamics simulation. To validate our strategy, the same methods were applied to three proteins of known structure; PDB:1HP8, PDB:1V54 chain H (residues 21-85), and PDB:1T0P. When compared to the crystal structures, the models had root mean-square deviations of 3.8 Å, 4.1 Å, and 4.0 Å, respectively. The model of RAMP1 suggested that Phe93, Tyr 100, and Phe101 form a binding interface for CLR, whereas Trp74 and Phe92 may interact with ligands that bind to the CLR/RAMP1 heterodimer. © 2006 by the Biophysical Society.

Molecular phase space transport in water:non-stationary random walk model

Molecular transport in phase space is crucial for chemical reactions because it defines how pre-reactive molecular configurations are found during the time evolution of the system. Using Molecular Dynamics (MD) simulated atomistic trajectories we test the assumption of the normal diffusion in the phase space for bulk water at ambient conditions by checking the equivalence of the transport to the random walk model. Contrary to common expectations we have found that some statistical features of the transport in the phase space differ from those of the normal diffusion models. This implies a non-random character of the path search process by the reacting complexes in water solutions. Our further numerical experiments show that a significant long period of non-stationarity in the transition probabilities of the segments of molecular trajectories can account for the observed non-uniform filling of the phase space. Surprisingly, the characteristic periods in the model non-stationarity constitute hundreds of nanoseconds, that is much longer time scales compared to typical lifetime of known liquid water molecular structures (several picoseconds).

21-residue peptide's dynamics at and between elementary structural transitions

Elementary conformational changes of the backbone of a 21-residue peptide A5(A3RA)3A are studied using molecular dynamics simulations in explicit water. The processes of the conformational transitions and the regimes of stationary fluctuations between them are investigated using minimal perturbations of the system. The perturbations consist of a few degrees rotation of the velocity of one of the systems' atoms and keep the system on the same energy surface. It is found that (i) the system dynamics is insignificantly changed by the perturbations in the regimes between the transitions; (ii) it is very sensitive to the perturbations just before the transitions that prevents the peptide from making the transitions; and (iii) the perturbation of any atom of the system, including distant water molecules is equally effective in preventing the transition. The latter implies strongly collective dynamics of the peptide and water during the transitions.

Identifying and correcting non-Markov states in peptide conformational dynamics

Conformational transitions in proteins define their biological activity and can be investigated in detail using the Markov state model. The fundamental assumption on the transitions between the states, their Markov property, is critical in this framework. We test this assumption by analyzing the transitions obtained directly from the dynamics of a molecular dynamics simulated peptide valine-proline-alanine-leucine and states defined phenomenologically using clustering in dihedral space. We find that the transitions are Markovian at the time scale of ˜ 50 ps and longer. However, at the time scale of 30–40 ps the dynamics loses its Markov property. Our methodology reveals the mechanism that leads to non-Markov behavior. It also provides a way of regrouping the conformations into new states that now possess the required Markov property of their dynamics.

Ion interactions with the carbon nanotube surface in aqueous solutions:understanding the molecular mechanisms

We study the molecular mechanisms of alkali halide ion interactions with the single-wall carbon nanotube surface in water by means of fully atomistic molecular dynamics simulations. We focus on the basic physical-chemical principles of ion–nanotube interactions in aqueous solutions and discuss them in light of recent experimental findings on selective ion effects on carbon nanotubes.

Water-peptide dynamics during conformational transitions

Transitions between metastable conformations of a dipeptide are investigated using classical molecular dynamics simulation with explicit water molecules. The distribution of the surrounding water at different moments before the transitions and the dynamical correlations of water with the peptide's configurational motions indicate that the water molecules represent an integral part of the molecular system during the conformational changes, in contrast to the metastable periods when water and peptide dynamics are essentially decoupled.

Dynamic complexity of chaotic transitions in high-dimensional classical dynamics:Leu-Enkephalin folding

Leu-Enkephalin in explicit water is simulated using classical molecular dynamics. A ß-turn transition is investigated by calculating the topological complexity (in the "computational mechanics" framework [J. P. Crutchfield and K. Young, Phys. Rev. Lett., 63, 105 (1989)]) of the dynamics of both the peptide and the neighbouring water molecules. The complexity of the atomic trajectories of the (relatively short) simulations used in this study reflect the degree of phase space mixing in the system. It is demonstrated that the dynamic complexity of the hydrogen atoms of the peptide and almost all of the hydrogens of the neighbouring waters exhibit a minimum precisely at the moment of the ß-turn transition. This indicates the appearance of simplified periodic patterns in the atomic motion, which could correspond to high-dimensional tori in the phase space. It is hypothesized that this behaviour is the manifestation of the effect described in the approach to molecular transitions by Komatsuzaki and Berry [T. Komatsuzaki and R.S. Berry, Adv. Chem. Phys., 123, 79 (2002)], where a "quasi-regular" dynamics at the transition is suggested. Therefore, for the first time, the less chaotic character of the folding transition in a realistic molecular system is demonstrated. © Springer-Verlag Berlin Heidelberg 2006.

Non-Markov state model of peptide dynamics

A hidden Markov state model has been applied to classical molecular dynamics simulated small peptide in explicit water. The methodology allows increasing the time resolution of the model and describe the dynamics with the precision of 0.3 ps (comparing to 6 ps for the standard methodology). It also permits the investigation of the mechanisms of transitions between the conformational states of the peptide. The detailed description of one of such transitions for the studied molecule is presented. © 2012 Elsevier B.V. All rights reserved.

Molecular mechanisms of salt effects on carbon nanotube dispersions in an organic solvent (N-methyl-2-pyrrolidone)

We consider the effects of salt (sodium iodide) on pristine carbon nanotube (CNT) dispersions in an organic solvent, N-methyl-2-pyrrolidone (NMP). We investigate the molecular-scale mechanisms of ion interactions with the nanotube surface and we show how the microscopic ion-surface interactions affect the stability of CNT dispersions in NMP. In our study we use a combination of fully atomistic Molecular Dynamics simulations of sodium and iodide ions at the CNT-NMP interface with direct experiments on the CNT dispersions. In the experiments we analyze the effects of salt on the stability of the dispersions by photoluminescence (PL) and optical absorption spectroscopy of the samples as well as by visual inspection. By fully atomistic Molecular Dynamics simulations we investigate the molecular-scale mechanisms of sodium and iodide ion interactions with the nanotube surface. Our simulations reveal that both ions are depleted from the CNT surface in the CNT-NMP dispersions mainly due to the two reasons: (1) there is a high energy penalty for the ion partial desolvation at the CNT surface; (2) NMP molecules form a dense solvation layer at the CNT surface that prevents ions to come close to the CNT surface. As a result, an increase of the salt concentration increases the "osmotic" stress in the CNT-NMP system and, thus, decreases the stability of the CNT dispersions in NMP. Direct experiments confirm the simulation results: addition of NaI salt into the NMP dispersions of pristine CNTs leads to precipitation of CNTs (bundle formation) even at very small salt concentration (∼10 -3 mol L -1). In line with the simulation predictions, the effect increases with the increase of the salt concentration. Overall, our results show that dissolved salt ions have strong effects on the stability of CNT dispersions. Therefore, it is possible to stimulate the bundle formation in the CNT-NMP dispersions and regulate the overall concentration of nanotubes in the dispersions by changing the NaI concentration in the solvent. © 2012 The Royal Society of Chemistry.

Molecular alignment relaxation in polymer optical fibers for sensing applications

A systematic study of annealing behavior of drawn PMMA fibers was performed. Annealing dynamics were investigated under different environmental conditions by fiber longitudinal shrinkage monitoring. The shrinkage process was found to follow a stretched exponential decay function revealing the heterogeneous nature of the underlying molecular dynamics. The complex dependence of the fiber shrinkage on initial degree of molecular alignment in the fiber, annealing time and temperature was investigated and interpreted. Moreover, humidity was shown to have a profound effect on the annealing process, which was not recognized previously. Annealing was also shown to have considerable effect on the fiber mechanical properties associated with the relaxation of molecular alignment in the fiber. The consequences of fiber annealing for the climatic stability of certain polymer optical fiber-based sensors are discussed, emphasizing the importance of fiber controlled pre-annealing with respect to the foreseeable operating conditions.

Dynamical frustration of protein's environment at the nanoseconds time scale

A 21-residue peptide in explicit water has been simulated using classical molecular dynamics. The system's trajectory has been analysed with a novel approach that quantifies the process of how atom's environment trajectories are explored. The approach is based on the measure of Statistical Complexity that extracts complete dynamical information from the signal. The introduced characteristic quantifies the system's dynamics at the nanoseconds time scale. It has been found that the peptide exhibits nanoseconds long periods that significantly differ in the rates of the exploration of the dynamically allowed configurations of the environment. During these periods the rates remain the same but different from other periods and from the rate for water. Periods of dynamical frustration are detected when only limited routes in the space of possible trajectories of the surrounding atoms are realised.

Design and synthesis of potential antimicrobial agents

Chorismate mutase is one of the essential enzymes in the shikimate pathway and is key to the survival of the organism Mycobacterium tuberculosis. The x-ray crystal structure of this enzyme from Mycobacterium tuberculosis was manipulated to prepare an initial set of in silico protein models of the active site. Known inhibitors of the enzyme were docked into the active site using the flexible ligand / flexible active site side chains approach implemented in CAChe Worksystem (Fujitsu Ltd). The resulting complexes were refined by molecular dynamics studies in explicit water using Amber 9. This yielded a further set of protein models that were used for additional rounds of ligand docking. A binding hypothesis was established for the enzyme and this was used to screen a database of commercially available drug-like compounds. From these results new potential ligands were designed that fitted appropriately into the active site and matched the functional groups and binding motifs founds therein. Some of these compounds and close analogues were then synthesized and submitted for biological evaluation. As a separate part of this thesis, analogues of very active anti-tuberculosis pyridylcarboxamidrazone were also prepared. This was carried out by the addition and the deletion of the substitutions from the lead compound thereby preparing heteroaryl carboxamidrazone derivatives and related compounds. All these compounds were initially evaluated for biological activity against various gram positive organisms and then sent to the TAACF (USA) for screening against Mycobacterium tuberculosis. Some of the new compounds proved to be at least as potent as the original lead compound but less toxic.

Strain distributionand electronic property modifications in Si/Ge axial nanowires hetrostructures

Molecular dynamics simulations were carried out for Si/Ge axial nanowire heterostructures using modified effective atom method (MEAM) potentials. A Si–Ge MEAM interatomic cross potential was developed based on available experimental data and was used for these studies. The atomic distortions and strain distributions near the Si/Ge interfaces are predicted for nanowires with their axes oriented along the [111] direction. The cases of 10 and 25 nm diameter Si/Ge biwires and of 25 nm diameter Si/Ge/Si axial heterostructures with the Ge disk 1 nm thick were studied. Substantial distortions in the height of the atoms adjacent to the interface were found for the biwires but not for the Ge disks. Strains as high as 3.5% were found for the Ge disk and values of 2%–2.5% were found at the Si and Ge interfacial layers in the biwires. Deformation potential theory was used to estimate the influence of the strains on the band gap, and reductions in band gap to as small as 40% of bulk values are predicted for the Ge disks. The localized regions of increased strain and resulting energy minima were also found within the Si/Ge biwire interfaces with the larger effects on the Ge side of the interface. The regions of strain maxima near and within the interfaces are anticipated to be useful for tailoring band gaps and producing quantum confinement of carriers. These results suggest that nanowire heterostructures provide greater design flexibility in band structure modification than is possible with planar layer growth.

Novel imprinted polymers as artificial enzymes

Derivatives of L-histidine were investigated as suitable models for the Asp-His couple found in the catalytic triad of serine proteases. A combination of molecular dynamics and IH NMR spectroscopy suggested that the most populous conformations of N-acetyl-L-histidine and the N-acetyl-L-histidine anion were predominated by those in which the carboxylate group was gauche to the imidazole ring overcoming steric and electrostatic repulsion, suggesting there is an interaction between the carboxylate group and the imidazole ring. Kinetic studies, using imidazole, N-acetyl-L-histidine and the N-acetyl-L-histidine anion showed that in a DMSO/H20 9: 1 v/v solution, the N-acetyl-L-histidine anion catalysed the hydrolysis of p-nitrophenyl acetate at a greater rate than using either imidazole or N-acetyl-L-histidine as catalyst. This indicates that the carboxylate group affects the nucleophilicity of the unprotonated imidazole ring. 31P MAS NMR spectroscopy was investigated as a new technique for the study of the template molecule environment within the polymer networks. It was found that it was possible to distinguish between template associated with the polymer and that which was precipitated onto the surface, though it was not possible to distinguish between polymer within imprinted cavities and that which was not. Attempts to study the effect of the carboxylate group/imidazole ring interaction in the imprinted cavity of a molecularly imprinted polymer network were hindered by the method used to follow the reaction. It was found though that in a pH 8.0 buffered solution the presence of imprinted cavities increased the rate of reaction for those polymers derived from L-histidine. Some preliminary investigations into the design and synthesis of an MIP which would catalyse the oxy-Cope rearrangement were carried out but the results were inconclusive.

Why are MD simulated protein folding times wrong?

The question of significant deviations of protein folding times simulated using molecular dynamics from experimental values is investigated. It is shown that in the framework of Markov State Model (MSM) describing the conformational dynamics of peptides and proteins, the folding time is very sensitive to the simulation model parameters, such as forcefield and temperature. Using two peptides as examples, we show that the deviations in the folding times can reach an order of magnitude for modest variations of the molecular model. We, therefore, conclude that the folding rate values obtained in molecular dynamics simulations have to be treated with care.