The creation of a database of odorous compounds focused on molecular rigidity and analysis of the molecular features of the compounds in the database

It is important to select odorous molecules for experiments on olfaction and for the development of an electronic nose. Odorous molecules having a small number of conformers, namely structurally rigid molecules, are assumed to interact with a small number of types of olfactory receptor proteins or to interact with the proteins in a simpler manner than that of fairly flexible molecules. Focusing on the rigidity of molecular structures, we collected 287 odorous molecules from data sources, which included 1205 chemicals in total and a database of the 287 odorous molecules (DB_odMOL) was created using CS ChemFinder Pro (version 5.0). The logarithmic value of the octanol/water partition coefficient (log P) and melting point, boiling point and vapour pressure of the molecules were estimated using CS Chem3D Pro. The database DB_odMOL accumulates these estimated data in addition to literature values for odour quality, odour detection thresholds and the safety of molecules. The rigidity of the 287 molecules was further analysed by conformational analysis performed by molecular mechanics using Conformer in CS Chem3D Pro (version 5) and 72 rigid odorous molecules were selected. The 287 molecules were also analysed based on atomic composition, substructure and molecular size. Sixty-two odorous molecules among the 72 rigid odorous molecules were further selected based on their atomic composition. The 62 rigid molecules with simple atomic composition that were finally selected should be useful for researchers in choosing odorous molecules for the study of olfaction, including the fields of molecular biology, physiology, structure-odour relationships and other fields of the study concerning odour.     

Intermolecular interactions between protein and other molecules including hydration effects

The structural aspects of protein functions, e.g., molecular recognition such as enzyme-substrate and antibody-antigen interactions, are elucidated in terms of dehydration and atomic interactions. When a protein interacts with some target molecule, water molecules at the interacting regions of both molecules are removed, with loss of the hydration free energy, but gaining atomic interactions between atoms of the contact sites in both molecules. The free energies of association originating from the dehydration and interactions between the atoms can be computed from changes in the accessible surface areas of the atoms involved. The free energy due to interactions between atomic groups at the contact sites is estimated as the sum of those estimated from the changes in the accessible surface area of 7 atomic groups, assuming that the interactions are proportional to the change of the area. The chain enthalpies and entropies evaluated from experimental thermodynamic properties and hydration quantities at the standard temperature for 10 proteins were available to determine the proportional constants for the atomic groups. This method was applied to the evaluation of association constants for the dimerization of proteins and the formation of proteolytic enzyme-inhibitor complexes, and the computed constants were in agreement with the experimental ones. However, the method is not accurate enough to account quantitatively for the change in the thermal stability of mutants of T4 lysozyme. Nevertheless, this method provides a way to elucidate the interactions between molecules in solution.     

Probing protein-protein interactions in real time

We have used a prototype small cantilever atomic force microscope to observe, in real time, the interactions between individual protein molecules. In particular, we have observed individual molecules of the chaperonin protein GroES binding to and then dissociating from individual GroEL proteins, which were immobilized on a mica support. This work suggests that the small cantilever atomic force microscope is a useful tool for studying protein dynamics at the single molecule level.     

Fingerprinting polysaccharides with single-molecule atomic force microscopy

We report the use of an atomic force microscopy (AFM)-based force spectroscopy technique to identify, at the single-molecule level, the components of mixtures of polysaccharides. Previously, we showed that the elasticity of certain types of polysaccharides is governed by force-induced conformational transitions of the pyranose ring. These transitions produce atomic fingerprints in the force-extension spectrum that are characteristic of the ground-energy conformation of the pyranose ring and the type of glycosidic linkages. Using this approach we find that commercially available agarose and lambda-carrageenan contain molecules that, when stretched in an atomic force microscope, produce a force spectrum characteristic of alpha-(1-->4) d-glucans. We have identified these molecules as amylopectin or floridean starch, a storage polysaccharide in algae. Our methodology can identify individual polysaccharide molecules in solution, which is not possible by any other spectroscopic technique, and therefore is an important addition to the arsenal of analytical techniques used in carbohydrate research.     

Incorporation of lone pairs into the electrostatic term in the empirical intra- and intermolecular energy calculations

The method hitherto used for estimating the electrostatic term in empirical intramolecular calculations of stable conformations of biologically important molecules and macromolecules and intermolecular calculations of molecular associations or packing energy in molecular crystals had been analyzed. It has been shown that the contribution of atomic hybridization moments is omitted in the calculation of electrostatic interactions from net atomic charges localized on nuclei which have been determined by standard quantum-chemical methods. This contribution plays an important part in determining electrostatic interactions, mainly in molecules containing atoms with lone pairs. Simultaneously, a modified method for calculating the electrostatic term comprising the interaction of the lone pairs, which are represented by atomic hybridization moments, has been proposed. The relationship between the atomic hybridization moment and the bond angle has been expressed for some typical configurations occurring in biologically important molecules. Finally, this new approach is illustrated by results of the conformational analysis of some model compounds for biomolecules and compared with the approach used so far for the estimation of the electrostatic interaction in empirical methods of calculation of the intra- and intermolecular energy.     

The visualization of amplicons after polymerase chain reaction

Linear DNA molecules amplified by the polymerase chain reaction were visualized by atomic force microscopy. The measured contour length of the PCR product of 1414 bp sequence was 435 +/- 15 nm. Considering that the calculated value of the distance between the nucleotides along the duplex axis is 0.31 nm, it was assumed that linear DNA molecules on the surface of mica, which serve as a support in the atomic force microscopy method, are in the A form. The influence of surface properties of the mica and the sample drying procedure on the conformation of adsorbed DNA molecules is discussed. Possible reasons for the Gaussian distribution of the contour length of the synthesized amplicon are considered.     

Protein dynamics by the combination of high-speed AFM and computational modeling

High-speed atomic force microscopy (HS-AFM) allows direct observation of biological molecules in dynamic action. However, HS-AFM has no atomic resolution. This article reviews recent progress of computational methods to infer high-resolution information, including the construction of 3D atomistic structures, from experimentally acquired resolution-limited HS-AFM images.     

A computational analysis of some diaryl ureas in relation to their observed crystalline hydrogen bonding patterns

We have carried out anab initio self-consistent-field molecular orbital analysis of the structures, atomic charges and electrostatic potentials of 1,3-bisphenylurea, 1,3-bis(m-nitrophenyl)urea and 1,3-bis(p-nitrophenyl)urea. Our results provide insight into the contrasting hydrogen bond patterns and preferences of these molecules in crystalline environments. The tendency for 1,3-bisphenylurea to form homomeric rather than heteromeric (i.e. host-guest) crystals is attributed to a relatively strong and extended homomeric attractive electrostatic interaction between two like molecules. Our calculated internuclear distances and atomic charges indicate that these molecules have some degree of H...O intramolecular hydrogen bonding between the carbonyl oxygen and the nearest aromatic hydrogens when the aromatic and the urea portions of the molecules are coplanar. This interaction is strongest for the meta nitro derivative, consistent with the latter remaining very nearly planar in its cocrystal structures. Our surface electrostatic potentials for the three diaryl ureas are used to interpret their differing crystallization properties and tendencies to form cocrystals with guest molecules.     

Solvating atomic level fine-grained proteins in supra-molecular level coarse-grained water for molecular dynamics simulations

Simulation of the dynamics of a protein in aqueous solution using an atomic model for both the protein and the many water molecules is still computationally extremely demanding considering the time scale of protein motions. The use of supra-atomic or supra-molecular coarse-grained (CG) models may enhance the computational efficiency, but inevitably at the cost of reduced accuracy. Coarse-graining solvent degrees of freedom is likely to yield a favourable balance between reduced accuracy and enhanced computational speed. Here, the use of a supra-molecular coarse-grained water model that largely preserves the thermodynamic and dielectric properties of atomic level fine-grained (FG) water in molecular dynamics simulations of an atomic model for four proteins is investigated. The results of using an FG, a CG, an implicit, or a vacuum solvent environment of the four proteins are compared, and for hen egg-white lysozyme a comparison to NMR data is made. The mixed-grained simulations do not show large differences compared to the FG atomic level simulations, apart from an increased tendency to form hydrogen bonds between long side chains, which is due to the reduced ability of the supra-molecular CG beads that represent five FG water molecules to make solvent-protein hydrogen bonds. But, the mixed-grained simulations are at least an order of magnitude faster than the atomic level ones.     

Composite structural motifs of binding sites for delineating biological functions of proteins

Most biological processes are described as a series of interactions between proteins and other molecules, and interactions are in turn described in terms of atomic structures. To annotate protein functions as sets of interaction states at atomic resolution, and thereby to better understand the relation between protein interactions and biological functions, we conducted exhaustive all-against-all atomic structure comparisons of all known binding sites for ligands including small molecules, proteins and nucleic acids, and identified recurring elementary motifs. By integrating the elementary motifs associated with each subunit, we defined composite motifs that represent context-dependent combinations of elementary motifs. It is demonstrated that function similarity can be better inferred from composite motif similarity compared to the similarity of protein sequences or of individual binding sites. By integrating the composite motifs associated with each protein function, we define meta-composite motifs each of which is regarded as a time-independent diagrammatic representation of a biological process. It is shown that meta-composite motifs provide richer annotations of biological processes than sequence clusters. The present results serve as a basis for bridging atomic structures to higher-order biological phenomena by classification and integration of binding site structures.     

High-speed atomic force microscopy: Structure and dynamics of single proteins

For surface analysis of biological molecules, atomic force microscopy (AFM) is an appealing technique combining data acquisition under physiological conditions, for example buffer solution, room temperature and ambient pressure, and high resolution. However, a key feature of life, dynamics, could not be assessed until recently because of the slowness of conventional AFM setups. Thus, for observing bio-molecular processes, the gain of image acquisition speed signifies a key progress. Here, we review the development and recent achievements using high-speed atomic force microscopy (HS-AFM). The HS-AFM is now the only technique to assess structure and dynamics of single molecules, revealing molecular motor action and diffusion dynamics. From this imaging data, watching molecules at work, novel and direct insights could be gained concerning the structure, dynamics and function relationship at the single bio-molecule level.     

Computer simulations of the flexibility of a series of synthetic cyclic peptide analogues

The flexibility of a series of cyclic peptides derived from the epitope of a snake toxin is investigated using computer simulation techniques. Molecular dynamics (MD) simulations and vibrational analyses are performed on chemically constrained peptides modeled on the parent protein loop. In the 50 ps MD simulations, large variations in the atomic fluctuations are seen between the peptides, and can be related to the nature of the chemical constraints present in the molecules. Normal mode analyses are performed on energy-minimized configurations derived from the dynamics trajectories. The atomic fluctuations calculated from the normal modes are about 30% of those of the molecular dynamics for the more flexible peptides and 70% for the more constrained molecules. The calculated differences in flexibility between the molecules are much less significant in the harmonic approximation. © 1993 John Wiley & Sons, Inc.     

MANOSK: A graphics program for analyzing and modeling molecular structure and functions

A program written for the Evans and Sutherland PS300 graphic displays is described in this paper. Called MANOSK, it provides a flexible and powerful environment for displaying, manipulating and analyzing small molecules and macromolecules from atomic coordinates. Translations and rotations of up to four independent molecules are available from the dials, and screen-relative orientations of the molecules are used in all geometrical and energetical calculations. A gradual progression of functions complexity makes the program easy to use for occasional works and efficient for more sophisticated studies.     

Atomic coordinates for tosyl-elastase

The full set of non-hydrogen, atomic coordinates of porcine tosyl elastase is given. This includes the tosyl group and the internal water molecules.     

Direct measurement of DNA by means of AFM.

We have demonstrated that the electrostatic stretch-and-positioning method is useful for the analysis of a long DNA molecule by means of atomic force microscopy (AFM). DNA molecules were stretched parallel to the field line, and immobilized onto the aluminum electrodes patterned on a glass plate. Through AFM observation, we confirmed the immobilization of individual DNA molecules, not aggregate.     

On contribution of known atomic partial charges of protein backbone in electrostatic potential density maps

Partial charges of atoms in a molecule and electrostatic potential (ESP) density for that molecule are known to bear a strong correlation. In order to generate a set of point‐field force field parameters for molecular dynamics, Kollman and coworkers have extracted atomic partial charges for each of all 20 amino acids using restrained partial charge‐fitting procedures from theoretical ESP density obtained from condensed‐state quantum mechanics. The magnitude of atomic partial charges for neutral peptide backbone they have obtained is similar to that of partial atomic charges for ionized carboxylate side chain atoms. In this study, the effect of these known atomic partial charges on ESP is examined using computer simulations and compared with the experimental ESP density recently obtained for proteins using electron microscopy. It is found that the observed ESP density maps are most consistent with the simulations that include atomic partial charges of protein backbone. Therefore, atomic partial charges are integral part of atomic properties in protein molecules and should be included in model refinement.     

Construction and optimisation of a computer model for a bacterial membrane

The main steps in the construction of a computer model for a bacterial membrane are described. The membrane has been built of 72 lipid molecules, 54 of which being 1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphatidylethanolamine (POPE) and 18--1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphatidyl-rac-glycerol (POPG) molecules (thus in the proportion of 3:1). The membrane was hydrated with 1955 water molecules (approximately 27 water molecules per lipid). To neutralise the electronic charge (-e) on each POPG molecule, 18 sodium ions (Na+) were added to the membrane close to the POPG phosphate groups. The atomic charges on the POPE and POPG headgroups were obtained from ab initio quantum mechanical restrained electrostatic potential fitting (RESP) (Bayly et al., 1993, J. Phys. Chem. 97, 10269) using the GAMESS program at the 6-31G* level (Schmidt et al., 1993, J. Comput. Chem. 14, 1347). The model constructed in this way provided an initial structure for subsequent molecular dynamics simulation studies intended to elucidate the atomic level interactions responsible for the structure and dynamics of the bacterial membrane.     

Mean-squared atomic displacements in hydrated lysozyme,native and denatured

We use elastic neutron scattering to demonstrate that a sharp increase in the mean-squared atomic displacements, commonly observed in hydrated proteins above 200 K and often referred to as the dynamical transition, is present in the hydrated state of both native and denatured lysozyme. A direct comparison of the native and denatured protein thus confirms that the presence of the transition in the mean-squared atomic displacements is not specific to biologically functional molecules.