An evolutionary algorithm for the surface structure problem

Many macroscopic properties: hardness, corrosion, catalytic activity, etc. are directly related to the surface structure, that is, to the position and chemical identity of the outermost atoms of the material. Current experimental techniques for its determination produce a “signature” from which the structure must be inferred by solving an inverse problem: a solution is proposed, its corresponding signature computed and then compared to the experiment. This is a challenging optimization problem where the search space and the number of local minima grows exponentially with the number of atoms, hence its solution cannot be achieved for arbitrarily large structures. Nowadays, it is solved by using a mixture of human knowledge and local search techniques: an expert proposes a solution that is refined using a local minimizer. If the outcome does not fit the experiment, a new solution must be proposed again. Solving a small surface can take from days to weeks of this trial and error method. Here we describe our ongoing work in its solution. We use an hybrid algorithm that mixes evolutionary techniques with trusted region methods and reuses knowledge gained during the execution to avoid repeated search of structures. Its parallelization produces good results even when not requiring the gathering of the full population, hence it can be used in loosely coupled environments such as grids. With this algorithm, the solution of test cases that previously took weeks of expert time can be automatically solved in a day or two of uniprocessor time.

SVis: A computational steering visualization environment for surface structure determination

The arrangement of atoms at the surface of a solid accounts for many of its properties: Hardness, chemical activity, corrosion, etc. are dictated by the precise surface structure. Hence, finding it, has a broad range of technical and industrial applications. The ability to solve this problem opens the possibility of designing by computer materials with properties tailored to specific applications. Since the search space grows exponentially with the number of atoms, its solution cannot be achieved for arbitrarily large structures. Presently, a trial and error procedure is used: an expert proposes an structure as a candidate solution and tries a local optimization procedure on it. The solution relaxes to the local minimum in the attractor basin corresponding to the initial point, that might be the one corresponding to the global minimum or not. This procedure is very time consuming and, for reasonably sized surfaces, can take many iterations and much effort from the expert. Here we report on a visualization environment designed to steer this process in an attempt to solve bigger structures and reduce the time needed. The idea is to use an immersive environment to interact with the computation. It has immediate feedback to assess the quality of the proposed structure in order to let the expert explore the space of candidate solutions. The visualization environment is also able to communicate with the de facto local solver used for this problem. The user is then able to send trial structures to the local minimizer and track its progress as they approach the minimum. This allows for simultaneous testing of candidate structures. The system has also proved very useful as an educational tool for the field.

The solution structure of the Zα domain of the human RNA editing enzyme ADAR1 reveals a prepositioned binding surface for Z-DNA

Double-stranded RNA deaminase I (ADAR1) contains the Z-DNA binding domain Zα. Here we report the solution structure of free Zα and map the interaction surface with Z-DNA, confirming roles previously assigned to residues by mutagenesis. Comparison with the crystal structure of the (Zα)2/Z-DNA complex shows that most Z-DNA contacting residues in free Zα are prepositioned to bind Z-DNA, thus minimizing the entropic cost of binding. Comparison with homologous (α+β)helix–turn–helix/B-DNA complexes suggests that binding of Zα to B-DNA is disfavored by steric hindrance, but does not eliminate the possibility that related domains may bind to both B- and Z-DNA.

Linking topography of its potential surface with the dynamics of folding of a protein model

The “3-color, 46-bead” model of a folding polypeptide is the vehicle for adapting to proteins a mode of analysis used heretofore for atomic clusters, to relate the topography of the potential surface to the dynamics that lead to formation of selected structures. The analysis is based on sequences of stationary points—successive minima, joined by saddles—that rise monotonically in energy from basin bottoms. Like structure-seeking clusters, the potential surface of the model studied here is staircase-like, rather than sawtooth-like, with highly collective motions required for passage from one minimum to the next. The surface has several deep basins whose minima correspond to very similar structures, but which are separated by high energy barriers.

Controlled unzipping of a bacterial surface layer with atomic force microscopy

We have combined high-resolution atomic force microscopy (AFM) imaging and force spectroscopy to gain insight into the interaction forces between the individual protomers of the hexagonally packed intermediate (HPI) layer of Deinococcus radiodurans. After imaging the HPI layer, the AFM stylus was attached to individual protomers by enforced stylus-sample contact to allow force spectroscopy experiments. Imaging of the HPI layer after recording force-extension curves allowed adhesion forces to be correlated with structural alterations. By using this approach, individual protomers of the HPI layer were found to be removed at pulling forces of ≈300 pN. Furthermore, it was possible to sequentially unzip entire bacterial pores formed by six HPI protomers. The combination of high-resolution AFM imaging of individual proteins with the determination of their intramolecular forces is a method of studying the mechanical stability of supramolecular structures at the level of single molecules.

Structural activity of a cloned potassium channel (ROMK1) monitored with the atomic force microscope: The “molecular-sandwich” technique

The atomic force microscope (AFM) was used to continuously follow height changes of individual protein molecules exposed to physiological stimuli. A AFM tip was coated with ROMK1 (a cloned renal epithelial potassium channel known to be highly pH sensitive) and lowered onto atomically flat mica surface until the protein was sandwiched between AFM tip and mica. Because the AFM tip was an integral part of a highly flexible cantilever, any structural alterations of the sandwiched molecule were transmitted to the cantilever. This resulted in a distortion of the cantilever that was monitored by means of a laser beam. With this system it was possible to resolve vertical height changes in the ROMK1 protein of ≥0.2 nm (approximately 5% of the molecule’s height) with a time resolution of ≥1 msec. When bathed in electrolyte solution that contained the catalytic subunit of protein kinase A and 0.1 mM ATP (conditions that activate the native ion channel), we found stochastically occurring height fluctuations in the ROMK1 molecule. These changes in height were pH-dependent, being greatest at pH 7.6, and lowering the pH (either by titration or by the application of CO2) reduced their magnitude. The data show that overall changes in shape of proteins occur stochastically and increase in size and frequency when the proteins are active. This AFM “molecular-sandwich” technique, called MOST, measures structural activity of proteins in real time and could prove useful for studies on the relationship between structure and function of proteins at the molecular level.

The crystal structure of the immunity protein of colicin E7 suggests a possible colicin-interacting surface.

The immunity protein of colicin E7 (ImmE7) can bind specifically to the DNase-type colicin E7 and inhibit its bactericidal activity. Here we report the 1.8-angstrom crystal structure of the ImmE7 protein. This is the first x-ray structure determined in the superfamily of colicin immunity proteins. The ImmE7 protein consists of four antiparallel alpha-helices, folded in a topology similar to the architecture of a four-helix bundle structure. A region rich in acidic residues is identified. This negatively charged area has the greatest variability within the family of DNase-type immunity proteins; thus, it seems likely that this area is involved in specific binding to colicin. Based on structural, genetic, and kinetic data, we suggest that all the DNase-type immunity proteins, as well as colicins, share a "homologous-structural framework" and that specific interaction between a colicin and its cognate immunity protein relies upon how well these two proteins' charged residues match on the interaction surface, thus leading to specific immunity of the colicin.

Mapping leaf surface landscapes.

Leaf surfaces provide the ecologically relevant landscapes to those organisms that encounter or colonize the leaf surface. Leaf surface topography directly affects microhabitat availability for colonizing microbes, microhabitat quality and acceptability for insects, and the efficacy of agricultural spray applications. Prior detailed mechanistic studies that examined particular fungi-plant and pollinator-plant interactions have demonstrated the importance of plant surface topography or roughness in determining the outcome of the interactions. Until now, however, it has not been possible to measure accurately the topography--i.e., the three-dimensional structure--of such leaf surfaces or to record precise changes in patterns of leaf surface elevation over time. Using contact mode atomic force microscopy, we measured three-dimensional coordinates of upper leaf surfaces of Vaccinium macrocarpon (cranberry), a perennial plant, on leaves of two age classes. We then produced topographic maps of these leaf surfaces, which revealed striking differences between age classes of leaves: old leaves have much rougher surfaces than those of young leaves. Atomic force microscope measurements were analyzed by lag (1) autocorrelation estimates of leaf surfaces by age class. We suggest that the changes in topography result from removal of epicuticular lipids and that the changes in leaf surface topography influence phylloplane ecology. Visualizing and mapping leaf surfaces permit detailed investigations into leaf surface-mediated phenomena, improving our understanding of phylloplane interactions.

X-ray structure of clotting factor IXa: active site and module structure related to Xase activity and hemophilia B.

Hereditary deficiency of factor IXa (fIXa), a key enzyme in blood coagulation, causes hemophilia B, a severe X chromosome-linked bleeding disorder afflicting 1 in 30,000 males; clinical studies have identified nearly 500 deleterious variants. The x-ray structure of porcine fIXa described here shows the atomic origins of the disease, while the spatial distribution of mutation sites suggests a structural model for factor X activation by phospholipid-bound fIXa and cofactor VIIIa. The 3.0-A-resolution diffraction data clearly show the structures of the serine proteinase module and the two preceding epidermal growth factor (EGF)-like modules; the N-terminal Gla module is partially disordered. The catalytic module, with covalent inhibitor D-Phe-1I-Pro-2I-Arg-3I chloromethyl ketone, most closely resembles fXa but differs significantly at several positions. Particularly noteworthy is the strained conformation of Glu-388, a residue strictly conserved in known fIXa sequences but conserved as Gly among other trypsin-like serine proteinases. Flexibility apparent in electron density together with modeling studies suggests that this may cause incomplete active site formation, even after zymogen, and hence the low catalytic activity of fIXa. The principal axes of the oblong EGF-like domains define an angle of 110 degrees, stabilized by a strictly conserved and fIX-specific interdomain salt bridge. The disorder of the Gla module, whose hydrophobic helix is apparent in electron density, can be attributed to the absence of calcium in the crystals; we have modeled the Gla module in its calcium form by using prothrombin fragment 1. The arched module arrangement agrees with fluorescence energy transfer experiments. Most hemophilic mutation sites of surface fIX residues occur on the concave surface of the bent molecule and suggest a plausible model for the membrane-bound ternary fIXa-FVIIIa-fX complex structure: fIXa and an equivalently arranged fX arch across an underlying fVIIIa subdomain from opposite sides; the stabilizing fVIIIa interactions force the catalytic modules together, completing fIXa active site formation and catalytic enhancement.

Electronic structure and transport properties of atomic NiO spinvalves

Ab initio quantum transport calculations show that short NiO chains suspended in Ni nanocontacts present a very strong spin-polarization of the conductance.The generalized gradient approximation we use here predicts a similar polarization of the conductance as the one previously computed with non-local exchange, confirming the robustness of the result. Their use as nanoscopic spinvalves is proposed.

Understanding the structure of the first atomic contact in gold

We have studied experimentally jump-to-contact (JC) and jump-out-of-contact (JOC) phenomena in gold electrodes. JC can be observed at first contact when two metals approach each other, while JOC occurs in the last contact before breaking. When the indentation depth between the electrodes is limited to a certain value of conductance, a highly reproducible behaviour in the evolution of the conductance can be obtained for hundreds of cycles of formation and rupture. Molecular dynamics simulations of this process show how the two metallic electrodes are shaped into tips of a well-defined crystallographic structure formed through a mechanical annealing mechanism. We report a detailed analysis of the atomic configurations obtained before contact and rupture of these stable structures and obtained their conductance using first-principles quantum transport calculations. These results help us understand the values of conductance obtained experimentally in the JC and JOC phenomena and improve our understanding of atomic-sized contacts and the evolution of their structural characteristics.

Molecular dynamics simulation of organic-inorganic nanocomposites: Layering behavior and interlayer structure of organoclays

Isothermal-isobaric (NPT) molecular dynamics simulation has been performed to investigate the layering behavior and structure of nanoconfined quaternary alkylammoniums in organoclays. This work is focused on systems consisting of two clay layers and a number of alkylammoniums, and involves the use of modified Dreiding force field. The simulated basal spacings of organoclays agree satisfactorily with the experimental results in the literature. The atomic density profiles in the direction normal to the clay surface indicate that the alkyl chains within the interlayer space of montmorillonite exhibit an obvious layering behavior. The headgroups of long alkyl chains are distributed within two layers close to the clay surface, whereas the distributions of methyl and methylene groups are strongly dependent on the alkyl chain length and clay layer charge. Monolayer, bilayer, and pseudo-trilayer structures are found in organoclays modified with single long alkyl chains, which are identical to the structural models based on the measured basal spacings. A pseudo-quadrilayer structure, for the first time to our knowledge, is also identified in organoclays with double long alkyl chains. In the mixture structure of paraffin-type and multilayer, alkyl chains do not lie flat within a single layer but interlace, and also jump to the next layer in pseudo-trilayer as well as next nearest layer in pseudo-quadrilayer.

Evidence for energy relaxation via a radiative cascade in surface-passivated PbS quantum dots

Multiple emission peaks have been observed from surface passivated PbS nanocrystals displaying strong quantum confinement. The emission spectra are shown to be strongly dependent on the excited-state parity. We also find that intraband energy relaxation from initial states excited far above the band-edge is nearly three orders of magnitude slower than that found in other nanocrystal quantum dots, providing evidence of inefficient energy relaxation via phonon emission. The initial-state parity dependence of the photoluminescent emission properties suggests that energy relaxation from the higher excited states occurs via a radiative cascade, analogous to energy relaxation in atomic systems. Such radiative cascade emission is possible from ideal zero-dimensional semiconductors, where electronic transitions can be decoupled from phonon modes.

Structure and thermal stability of Langmuir-Blodgett films of barium arachidate

The structures of multilayer Langmuir-Blodgett films of barium arachidate before and after heat treatment have been investigated using both atomic force microscopy (AFM) and grazing incidence synchrotron X-ray diffraction (GIXD). AFM gave information on surface morphology at molecular resolution while GIXD provided quantitative details of the lattice structures of the films with their crystal symmetries and lattice constants. As-prepared films contained three coexisting structures: two triclinic structures with the molecularchains tilted by about 20degrees from the film normal and with 3 x 1 or 2 x 2 super-lattice features arising from height modulation of the molecules in the films; a rectangular structure with molecules perpendicular to the film surface. Of these, the 3 x 1 structure is dominant with a loose correlation between the bilayers. In the film plane both superstructures are commensurate with the local structures, having different oblique symmetries. The lattice constants for the 3 x 1 structure are a(s) = 3a = 13.86 Angstrom, b(s) = b = 4.31 Angstrom and gamma(s) = gamma = 82.7degrees; for the 2 x 2 structure a(s) = 2a = 16.54 Angstrom, b(s) = 2b = 9.67 Angstrom, gamma(s) = gamma = 88degrees. For the rectangular structure the lattice constants are a = 7.39 Angstrom, b = 4.96 Angstrom and gamma = 90degrees. After annealing, the 2 x 2 and rectangular structures were not observed, while the 3 x 1 structure had developed over the entire film. For the annealed films the correlation length in the film plane is about twice that in the unheated films, and in the out-of-plane direction covers two bilayers. The above lattice parameters, determined by GIXD, differed significantly from the values obtained by AFM, due possibly to distortion of the films by the scanning action of the AFM tip. (C) 2004 Published by Elsevier B.V.

The structure of ultrathin Langmuir-Blodgett films of cadmium behenate

Grazing incidence x-ray-diffraction investigations of the structures of Langmuir-Blodgett films of cadmium behenate with 1, 2, 3, 5, and 21 monolayers are reported. The single monolayer film, deposited on a hydrophilic substrate, showed a hexagonal structure, whereas the bilayer film, deposited on a hydrophobic substrate, had a rectangular structure with herringbone orientation of the acyl chains. With multilayer films formed on a hydrophilic substrate, it was possible to detect that the hexagonal structure of the first layer was retained when additional layers were deposited and that the additional layers had the same rectangular structure as the bilayer. (c) 2005 American Institute of Physics.