Volume 10, Issue 11-12, Pages 887-984(November 2001)
Original Paper
Imaging of atomic orbitals with the Atomic Force Microscope - experiments and simulations
NASA Astrophysics Data System (ADS)
Giessibl, F. J.; Bielefeldt, H.; Hembacher, S.; Mannhart, J.
2001-11-01
Atomic force microscopy (AFM) is a mechanical profiling technique that allows to image surfaces with atomic resolution. Recent progress in reducing the noise of this technique has led to a resolution level where previously undetectable symmetries of the images of single atoms are observed. These symmetries are related to the nature of the interatomic forces. The Si(111)-(7 × 7) surface is studied by AFM with various tips and AFM images are simulated with chemical and electrostatic model forces. The calculation of images from the tip-sample forces is explained in detail and the implications of the imaging parameters are discussed. Because the structure of the Si(111)-(7 × 7) surface is known very well, the shape of the adatom images is used to determine the tip structure. The observability of atomic orbitals by AFM and scanning tunneling microscopy is discussed.
Spinello, A; Barone, G; Grunenberg, J
2016-01-28
In depth Monte Carlo conformational scans in combination with molecular dynamics (MD) simulations and electronic structure calculations were applied in order to study the molecular recognition process between tetrasubstituted naphthalene diimide (ND) guests and G-quadruplex (G4) DNA receptors. ND guests are a promising class of telomere stabilizers due to which they are used in novel anticancer therapeutics. Though several ND guests have been studied experimentally in the past, the protonation state under physiological conditions is still unclear. Based on chemical intuition, in the case of N-methyl-piperazine substitution, different protonation states are possible and might play a crucial role in the molecular recognition process by G4-DNA. Depending on the proton concentration, different nitrogen atoms of the N-methyl-piperazine might (or might not) be protonated. This fact was considered in our simulation in terms of a case by case analysis, since the process of molecular recognition is determined by possible donor or acceptor positions. The results of our simulations show that the electrostatic interactions between the ND ligands and the G4 receptor are maximized in the case of the protonation of the terminal nitrogen atoms, forming compact ND G4 complexes inside the grooves. The influence of different protonation states in terms of the ability to form hydrogen bonds with the sugar-phosphate backbone, as well as the importance of mediated vs. direct hydrogen bonding, was analyzed in detail by MD and relaxed force constant (compliance constant) simulations.
Beckstein, Oliver; Denning, Elizabeth J.; Perilla, Juan R.; Woolf, Thomas B.
2009-01-01
Adenylate kinase (AdK), a phosphotransferase enzyme, plays an important role in cellular energy homeostasis. It undergoes a large conformational change between an open and a closed state, even in the absence of substrate. We investigate the apo-AdK transition at the atomic level both with free energy calculations and our new dynamic importance sampling (DIMS) molecular dynamics (MD) method. DIMS is shown to sample biologically relevant conformations as verified by comparing an ensemble of hundreds of DIMS transitions to AdK crystal structure intermediates. The simulations reveal in atomic detail how hinge regions partially and intermittently unfold during the transition. Conserved salt bridges are seen to have important structural and dynamic roles; in particular four ionic bonds are identified that open in a sequential, zipper-like fashion and thus dominate the free energy landscape of the transition. Transitions between the closed and open conformations only have to overcome moderate free energy barriers. Unexpectedly, the closed and open state encompass broad free energy basins that contain conformations differing in domain hinge motions by up to 40°. The significance of these extended states is discussed in relation to recent experimental FRET measurements. Taken together, these results demonstrate how a small number of cooperative key interactions can shape the overall dynamics of an enzyme and suggest an “all-or-nothing” mechanism for the opening and closing of AdK. Our efficient DIMS-MD computer simulation approach can provide a detailed picture of a functionally important macromolecular transition and thus help to interpret and suggest experiments to probe the conformational landscape of dynamic proteins such as AdK. PMID:19751742
Rudack, Till; Jenrich, Sarah; Brucker, Sven; Vetter, Ingrid R.; Gerwert, Klaus; Kötting, Carsten
2015-01-01
Small GTPases regulate key processes in cells. Malfunction of their GTPase reaction by mutations is involved in severe diseases. Here, we compare the GTPase reaction of the slower hydrolyzing GTPase Ran with Ras. By combination of time-resolved FTIR difference spectroscopy and QM/MM simulations we elucidate that the Mg2+ coordination by the phosphate groups, which varies largely among the x-ray structures, is the same for Ran and Ras. A new x-ray structure of a Ran·RanBD1 complex with improved resolution confirmed this finding and revealed a general problem with the refinement of Mg2+ in GTPases. The Mg2+ coordination is not responsible for the much slower GTPase reaction of Ran. Instead, the location of the Tyr-39 side chain of Ran between the γ-phosphate and Gln-69 prevents the optimal positioning of the attacking water molecule by the Gln-69 relative to the γ-phosphate. This is confirmed in the RanY39A·RanBD1 crystal structure. The QM/MM simulations provide IR spectra of the catalytic center, which agree very nicely with the experimental ones. The combination of both methods can correlate spectra with structure at atomic detail. For example the FTIR difference spectra of RasA18T and RanT25A mutants show that spectral differences are mainly due to the hydrogen bond of Thr-25 to the α-phosphate in Ran. By integration of x-ray structure analysis, experimental, and theoretical IR spectroscopy the catalytic center of the x-ray structural models are further refined to sub-Å resolution, allowing an improved understanding of catalysis. PMID:26272610
Rudack, Till; Jenrich, Sarah; Brucker, Sven; Vetter, Ingrid R; Gerwert, Klaus; Kötting, Carsten
2015-10-02
Small GTPases regulate key processes in cells. Malfunction of their GTPase reaction by mutations is involved in severe diseases. Here, we compare the GTPase reaction of the slower hydrolyzing GTPase Ran with Ras. By combination of time-resolved FTIR difference spectroscopy and QM/MM simulations we elucidate that the Mg(2+) coordination by the phosphate groups, which varies largely among the x-ray structures, is the same for Ran and Ras. A new x-ray structure of a Ran·RanBD1 complex with improved resolution confirmed this finding and revealed a general problem with the refinement of Mg(2+) in GTPases. The Mg(2+) coordination is not responsible for the much slower GTPase reaction of Ran. Instead, the location of the Tyr-39 side chain of Ran between the γ-phosphate and Gln-69 prevents the optimal positioning of the attacking water molecule by the Gln-69 relative to the γ-phosphate. This is confirmed in the RanY39A·RanBD1 crystal structure. The QM/MM simulations provide IR spectra of the catalytic center, which agree very nicely with the experimental ones. The combination of both methods can correlate spectra with structure at atomic detail. For example the FTIR difference spectra of RasA18T and RanT25A mutants show that spectral differences are mainly due to the hydrogen bond of Thr-25 to the α-phosphate in Ran. By integration of x-ray structure analysis, experimental, and theoretical IR spectroscopy the catalytic center of the x-ray structural models are further refined to sub-Å resolution, allowing an improved understanding of catalysis. © 2015 by The American Society for Biochemistry and Molecular Biology, Inc.
Alberini, Giulio; Benfenati, Fabio
2017-01-01
Tight-junctions between epithelial cells of biological barriers are specialized molecular structures that regulate the flux of solutes across the barrier, parallel to cell walls. The tight-junction backbone is made of strands of transmembrane proteins from the claudin family, but the molecular mechanism of its function is still not completely understood. Recently, the crystal structure of a mammalian claudin-15 was reported, displaying for the first time the detailed features of transmembrane and extracellular domains. Successively, a structural model of claudin-15-based paracellular channels has been proposed, suggesting a putative assembly that illustrates how claudins associate in the same cell (via cis interactions) and across adjacent cells (via trans interactions). Although very promising, the model offers only a static conformation, with residues missing in the most important extracellular regions and potential steric clashes. Here we present detailed atomic models of paracellular single and double pore architectures, obtained from the putative assembly and refined via structural modeling and all-atom molecular dynamics simulations in double membrane bilayer and water environment. Our results show an overall stable configuration of the complex with a fluctuating pore size. Extracellular residue loops in trans interaction are able to form stable contacts and regulate the size of the pore, which displays a stationary radius of 2.5–3.0 Å at the narrowest region. The side-by-side interactions of the cis configuration are preserved via stable hydrogen bonds, already predicted by cysteine crosslinking experiments. Overall, this work introduces an improved version of the claudin-15-based paracellular channel model that strengthens its validity and that can be used in further computational studies to understand the structural features of tight-junctions regulation. PMID:28863193
Investigating the Structure and Dynamics of the PIK3CA Wild-Type and H1047R Oncogenic Mutant
Pavlaki, Maria; Lazani, Vasiliki; Christoforidis, Savvas; Agianian, Bogos; Cournia, Zoe
2014-01-01
The PIK3CA gene is one of the most frequently mutated oncogenes in human cancers. It encodes p110α, the catalytic subunit of phosphatidylinositol 3-kinase alpha (PI3Kα), which activates signaling cascades leading to cell proliferation, survival, and cell growth. The most frequent mutation in PIK3CA is H1047R, which results in enzymatic overactivation. Understanding how the H1047R mutation causes the enhanced activity of the protein in atomic detail is central to developing mutant-specific therapeutics for cancer. To this end, Surface Plasmon Resonance (SPR) experiments and Molecular Dynamics (MD) simulations were carried out for both wild-type (WT) and H1047R mutant proteins. An expanded positive charge distribution on the membrane binding regions of the mutant with respect to the WT protein is observed through MD simulations, which justifies the increased ability of the mutated protein variant to bind to membranes rich in anionic lipids in our SPR experiments. Our results further support an auto-inhibitory role of the C-terminal tail in the WT protein, which is abolished in the mutant protein due to loss of crucial intermolecular interactions. Moreover, Functional Mode Analysis reveals that the H1047R mutation alters the twisting motion of the N-lobe of the kinase domain with respect to the C-lobe and shifts the position of the conserved P-loop residues in the vicinity of the active site. These findings demonstrate the dynamical and structural differences of the two proteins in atomic detail and propose a mechanism of overactivation for the mutant protein. The results may be further utilized for the design of mutant-specific PI3Kα inhibitors that exploit the altered mutant conformation. PMID:25340423
Imaging Cold Gas to 1 kpc scales in high-redshift galaxies with the ngVLA
NASA Astrophysics Data System (ADS)
Casey, Caitlin; Narayanan, Desika; Dave, Romeel; Hung, Chao-Ling; Champagne, Jaclyn; Carilli, Chris Luke; Decarli, Roberto; Murphy, Eric J.; Popping, Gergo; Riechers, Dominik; Somerville, Rachel S.; Walter, Fabian
2017-01-01
The next generation Very Large Array (ngVLA) will revolutionize our understanding of the distant Universe via the detection of cold molecular gas in the first galaxies. Its impact on studies of galaxy characterization via detailed gas dynamics will provide crucial insight on dominant physical drivers for star-formation in high redshift galaxies, including the exchange of gas from scales of the circumgalactic medium down to resolved clouds on mass scales of ~10^5 M_sun. In this study, we employ a series of high-resolution, cosmological, hydrodynamic zoom simulations from the MUFASA simulation suite and a CASA simulator to generate mock ngVLA observations. Based on a direct comparison between the inferred results from our mock observations and the cosmological simulations, we investigate the capabilities of ngVLA to constrain the mode of star formation, dynamical mass, and molecular gas kinematics in individual high-redshift galaxies using cold gas tracers like CO(1-0) and CO(2-1). Using the Despotic radiative transfer code that encompasses simultaneous thermal and statistical equilibrium in calculating the molecular and atomic level populations, we generate parallel mock observations of high-J transitions of CO and C+ from ALMA for comparison. The factor of 100 times improvement in mapping speed for the ngVLA beyond the Jansky VLA and the proposed ALMA Band 1 will make these detailed, high-resolution imaging and kinematic studies routine at z=2 and beyond.
Formoso, Elena; Matxain, Jon M; Lopez, Xabier; York, Darrin M
2010-06-03
The mechanisms of enzymes are intimately connected with their overall structure and dynamics in solution. Experimentally, it is considerably challenging to provide detailed atomic level information about the conformational events that occur at different stages along the chemical reaction path. Here, theoretical tools may offer new potential insights that complement those obtained from experiments that may not yield an unambiguous mechanistic interpretation. In this study, we apply molecular dynamics simulations of bovine pancreatic ribonuclease A, an archetype ribonuclease, to study the conformational dynamics, structural relaxation, and differential solvation that occur at discrete stages of the transesterification and cleavage reaction. Simulations were performed with explicit solvation with rigorous electrostatics and utilize recently developed molecular mechanical force field parameters for transphosphorylation and hydrolysis transition state analogues. Herein, we present results for the enzyme complexed with the dinucleotide substrate cytidilyl-3',5'-adenosine (CpA) in the reactant, and transphosphorylation and hydrolysis transition states. A detailed analysis of active site structures and hydrogen-bond patterns is presented and compared. The integrity of the overall backbone structure is preserved in the simulations and supports a mechanism whereby His12 stabilizes accumulating negative charge at the transition states through hydrogen-bond donation to the nonbridge oxygens. Lys41 is shown to be highly versatile along the reaction coordinate and can aid in the stabilization of the dianionic transition state, while being poised to act as a general acid catalyst in the hydrolysis step.
Formoso, Elena; Matxain, Jon M.; Lopez, Xabier; York, Darrin M.
2010-01-01
The mechanisms of enzymes are intimately connected with their overall structure and dynamics in solution. Experimentally it is considerably challenging to provide detailed atomic level information about the conformational events that occur at different stages along the chemical reaction path. Here, theoretical tools may offer new potential insights that complement those obtained from experiments that may not yield an unambiguous mechanistic interpretation. In this study we apply molecular dynamics simulations of bovine pancreatic ribonuclease A, an archetype ribonuclease, in order to study the conformational dynamics, structural relaxation, and differential solvation that occurs at discreet stages of the transesterification and cleavage reaction. Simulations were performed with explicit solvation with rigorous electrostatics, and utilize recently developed molecular mechanical force field parameters for transphosphorylation and hydrolysis transition state analogs. Herein, we present results for the enzyme complexed with the dinucleotide substrate cytidilyl-3′,5′-adenosine (CpA) in the reactant, and transphosphorylation and hydrolysis transition states. A detailed analysis of active site structures and hydrogen bond patterns are presented and compared. The integrity of the overall backbone structure is preserved in the simulations, and support a mechanism whereby His12 stabilizes accumulating negative charge at the transition states through hydrogen bond donation to the non-bridge oxygens. Lys41 is shown to be highly versatile along the reaction coordinate, and can aid in the stabilization of the dianionic transition state, while being poised to act as a general acid catalyst in the hydrolysis step. PMID:20455590
Santos, Elson C; Neto, Abel F G; Maneschy, Carlos E; Chen, James; Ramalho, Teodorico C; Neto, A M J C
2015-05-01
Here we analyzed several physical behaviors through computational simulation of systems consisting of a zig-zag type carbon nanotube and relaxed cold atoms (Rb, Au, Si and Ar). These atoms were chosen due to their different chemical properties. The atoms individually were relaxed on the outside of the nanotube during the simulations. Each system was found under the influence of a uniform electric field parallel to the carbon nanotube and under the thermal effect of the initial temperature at the simulations. Because of the electric field, the cold atoms orbited the carbon nanotube while increasing the initial temperature allowed the variation of the radius of the orbiting atoms. We calculated the following quantities: kinetic energy, potential energy and total energy and in situ temperature, molar entropy variation and average radius of the orbit of the atoms. Our data suggest that only the action of electric field is enough to generate the attractive potential and this system could be used as a selected atoms sensor.
Atomic Detail Visualization of Photosynthetic Membranes with GPU-Accelerated Ray Tracing
Vandivort, Kirby L.; Barragan, Angela; Singharoy, Abhishek; Teo, Ivan; Ribeiro, João V.; Isralewitz, Barry; Liu, Bo; Goh, Boon Chong; Phillips, James C.; MacGregor-Chatwin, Craig; Johnson, Matthew P.; Kourkoutis, Lena F.; Hunter, C. Neil
2016-01-01
The cellular process responsible for providing energy for most life on Earth, namely photosynthetic light-harvesting, requires the cooperation of hundreds of proteins across an organelle, involving length and time scales spanning several orders of magnitude over quantum and classical regimes. Simulation and visualization of this fundamental energy conversion process pose many unique methodological and computational challenges. We present, in two accompanying movies, light-harvesting in the photosynthetic apparatus found in purple bacteria, the so-called chromatophore. The movies are the culmination of three decades of modeling efforts, featuring the collaboration of theoretical, experimental, and computational scientists. We describe the techniques that were used to build, simulate, analyze, and visualize the structures shown in the movies, and we highlight cases where scientific needs spurred the development of new parallel algorithms that efficiently harness GPU accelerators and petascale computers. PMID:27274603
Ballistic Fracturing of Carbon Nanotubes.
Ozden, Sehmus; Machado, Leonardo D; Tiwary, ChandraSekhar; Autreto, Pedro A S; Vajtai, Robert; Barrera, Enrique V; Galvao, Douglas S; Ajayan, Pulickel M
2016-09-21
Advanced materials with multifunctional capabilities and high resistance to hypervelocity impact are of great interest to the designers of aerospace structures. Carbon nanotubes (CNTs) with their lightweight and high strength properties are alternative to metals and/or metallic alloys conventionally used in aerospace applications. Here we report a detailed study on the ballistic fracturing of CNTs for different velocity ranges. Our results show that the highly energetic impacts cause bond breakage and carbon atom rehybridizations, and sometimes extensive structural reconstructions were also observed. Experimental observations show the formation of nanoribbons, nanodiamonds, and covalently interconnected nanostructures, depending on impact conditions. Fully atomistic reactive molecular dynamics simulations were also carried out in order to gain further insights into the mechanism behind the transformation of CNTs. The simulations show that the velocity and relative orientation of the multiple colliding nanotubes are critical to determine the impact outcome.