Interfacial interactions and structure of organic-inorganic nanohybrids

Understanding the interfacial interactions and structure is important to better design and application of organic-inorganic nanohybrids. This paper presents our recent molecular dynamic studies on organoclays and polymer nanocomposites, including the layering behavior of organoclays, structural and dynamic properties of dioctadecyldimethyl ammoniums in organoclays, and interfacial interactions and structure of polyurethane nanocomposites. The results demonstrate that the layering behaviors of organoclays are closely related to the chain length of quaternary alkyl ammoniums and cation exchangeable capacity of clays. In addition to typical layered structures such as monolayer, bilayer and pseudo-trilayer, a pseudo-quadrilayer structure was also observed in organoclays modified with dioctadecyldimethyl ammoniums (DODDMA). In such a structure, alkyl chains do not lie flat within a single layer but interlace, and also jump to the next layer or even the next nearest layer. Moreover, the diffusion constants of nitrogen and methylene atoms increase with the temperature and methelene towards the tail groups. For polyurethane nanocomposite, the van der Waals interaction between apolar alkyl chains and soft segments of polyurethane predominates the interactions between organoclay and polyurethane. Different from most bulk polyurethane systems, there is no distinct phase-separated structure for the polyurethane.

Mapping of suramin binding sites on the group IIA human secreted phospholipase A(2)

Suramin is a polysulphonated napthylurea used as an antiprotozoal/anthelminitic drug, which also inhibits a broad range of enzymes. Suramin binding to recombinant human secreted group IIA phospholipase A(2) (hsPLA(2)GIIA) was investigated by molecular dynamics simulations (MD) and isothermal titration calorimetry (ITC). MD indicated two possible bound suramin conformations mediated by hydrophobic and electrostatic interactions with amino-acids in three regions of the protein. namely the active-site and residues located in the N- and C-termini, respectively. All three binding sites are located on the phospholipid membrane recognition surface, suggesting that suramin may inhibit the enzyme, and indeed a 90% reduction in hydrolytic activity was observed in the presence of 100 nM suramin. These results correlated with ITC data, which demonstrated 2.7 suramin binding sites on the hsPLA(2)GIIA, and indicates that suramin represents a novel class of phosphohpase A(2) inhibitor. (C) 2009 Elsevier Inc. All rights reserved.

Macromolecular assembly of polycystin-2 intracytosolic C-terminal domain

Mutations in PKD2 are responsible for approximately 15% of the autosomal dominant polycystic kidney disease cases. This gene encodes polycystin-2, a calcium-permeable cation channel whose C-terminal intracytosolic tail (PC2t) plays an important role in its interaction with a number of different proteins. In the present study, we have comprehensively evaluated the macromolecular assembly of PC2t homooligomer using a series of biophysical and biochemical analyses. Our studies, based on a new delimitation of PC2t, have revealed that it is capable of assembling as a homotetramer independently of any other portion of the molecule. Our data support this tetrameric arrangement in the presence and absence of calcium. Molecular dynamics simulations performed with a modified all-atoms structure-based model supported the PC2t tetrameric assembly, as well as how different populations are disposed in solution. The simulations demonstrated, indeed, that the best-scored structures are the ones compatible with a fourfold oligomeric state. These findings clarify the structural properties of PC2t domain and strongly support a homotetramer assembly of PC2.

Stability improvement of the fatty acid binding protein Sm14 from S. mansoni by Cys replacement: Structural and functional characterization of a vaccine candidate

The Schistosoma mansoni fatty acid binding protein (FABP), SmA, is a vaccine candidate against, S. mansoni and F hepatica. Previously, we demonstrated the importance of a correct fold to achieve protection in immunized animals after cercariae challenge [[10]. C.R.R. Ramos, R.C.R. Figueredo, T.A. Pertinhez, M.M. Vilar, A.L.T.O. Nascimento, M. Tendler, I. Raw, A. Spisni, P.L. Ho, Gene structure and M20T polymorphism of the Schistosoma mansoni Sm14 fatty acid-binding protein: structural, functional and immunoprotection analysis. J. Biol. Chem. 278 (2003) 12745-12751]. Here we show that the reduction of vaccine efficacy over time is due to protein dimerization and subsequent aggregation. We produced the mutants Sm14-M20(C62S) and Sm14M20(C62V) that, as expected, did not dimerize in SDS-PAGE. Molecular dynamics calculations and unfolding experiments highlighted a higher structural stability of these mutants with respect to the wild-type. In addition, we found that the mutated proteins, after thermal denaturation, refolded to their active native molecular architecture as proved by the recovery of the fatty acid binding ability. Sm14-M20(C62V) turned out to be the more stable form over time, providing the basis to determine the first 3D solution structure of a Sm14 protein in its apo-form. Overall, Sm14-M20(C62V) possesses an improved structural stability over time, an essential feature to preserve its immunization capability and, in experimentally immunized animals, it exhibits a protection effect against S. mansoni cercariae infections comparable to the one obtained with the wild-type protein. These facts indicate this protein as a good lead molecule for large-scale production and for developing an effective Sm14 based anti-helminthes vaccine. (C) 2008 Elsevier B.V. All rights reserved.

Solid-state NMR structure determination of melittin in a lipid environment

Solid-state C-13 NMR spectroscopy was used to investigate the three-dimensional structure of melittin as lyophilized powder and in ditetradecylphosphatidylcholine (DTPC) membranes. The distance between specifically labeled carbons in analogs [1-C-13]Gly3-[2-C-13]Ala4, [1-C-13]Gly3-[2-C-13]Leu6, [1-C-13]Leu13-[2-C-13]Ala15, [2-C-13]Leu13-[1-C-13]Ala15, and [1-C-13]Leu13-[2-C-13]Leu16 was measured by rotational resonance. As expected, the internuclear distances measured in [1-C-13]Gly3-[2-C-13]Ala4 and [1-C-13]Gly3-[2-C-13]Leu6 were consistent with alpha -helical structure in the N-terminus irrespective of environment. The Internuclear distances measured in [1-C-13]Leu13-[2-C-13]Ala15, [2-C-13]Leu13-[1-C-13]Ala15, and [1-C-13]Leu13-[2-C-13]Leu16 revealed, via molecular modeling, some dependence upon environment for conformation in the region of the bend in helical structure induced by Pro14. A slightly larger interhelical angle between the N- and C-terminal helices was indicated for peptide in dry or hydrated gel state DTPC (139 degrees -145 degrees) than in lyophilized powder (121 degrees -139 degrees) or crystals (129 degrees). The angle, however, is not as great as deduced for melittin in aligned bilayers of DTPC in the liquid-crystalline state (similar to 160 degrees) (R. Smith, F. Separovic, T. J. Milne, A. Whittaker, F. M. Bennett, B. A. Cornell, and A. Makriyannis, 1994, J. Mol, Biol 241:456-466). The study illustrates the utility of rotational resonance in determining local structure within peptide-lipid complexes.

Ring flexibility within tricyclic antidepressant drugs

The internal flexibility of the central seven-membered ring of a series of tricyclic antidepressant drugs (TCAs), imipramine {l}, amitriptyline {2}, doxepin {3}, and dothiepin {4}, has been investigated by H-1 and C-13 nuclear magnetic (NMR) techniques. Two dynamic processes were examined: ring inversion and bridge flexing. H-1 NMR lineshape analysis was used to obtain ring inversion barriers for 2-4. These studies yielded energy barriers of 14.3, 16.7, and 15.7 +/- 0.6 kcal/mol for the hydrochloride salts of doxepin, dothiepin, and amitriptyline, respectively. The barriers for the corresponding free bases were lower by 0.6 kcal/mol on average. (CT1)-C-13 relaxation measurements were used to determine the degree of bridge flexing associated with the central seven-membered ring for all four compounds. By fitting the T-1 data to a two-state jump model, lifetimes and amplitudes of rapid bridge flexing motions were determined. The results show that imipramine has the fastest rate of bridge flexing, followed by amitriptyline, doxepin, and dothiepin. The pharmacological profiles of the TCAs are complex and they interact with many receptor sites, resulting in numerous side effects and a general lack of understanding of their precise mode of action in different anxiety-related disorders. They all have similar three-dimensional structures, which makes it difficult to rationalize their differing relative potency in different assays/clinical settings. However, the clear finding here that there are significantly different degrees of internal mobility suggests that molecular dynamics should be an additional factor considered when trying to understand the mode of action of this clinically important family of molecules. (C) 2001 Wiley-Liss, Inc. and the American Pharmaceutical Association J Pharm Sci 90:713-721, 2001.

Analytical exclusion chromatography

This review summarizes the development of exclusion chromatography, also termed gel filtration, molecular-sieve chromatography and gel permeation chromatography, for the quantitative characterization of solutes and solute interactions. As well as affording a means of determining molecular mass and molecular mass distribution, the technique offers a convenient way of characterizing solute selfassociation and solute-ligand interactions in terms of reaction stoichiometry and equilibrium constant. The availability of molecular-sieve media with different selective porosities ensures that very little restriction is imposed on the size of solute amenable to study. Furthermore, access to a diverse array of assay procedures for monitoring the column eluate endows analytical exclusion chromatography with far greater flexibility than other techniques from the viewpoint of solute concentration range that can be examined. In addition to its widely recognized prowess as a means of solute separation and purification, exclusion chromatography thus also possesses considerable potential for investigating the functional roles of the purified solutes. (C) 2003 Elsevier Science B.V. All rights reserved.

Epitope map for a growth hormone receptor agonist monoclonal antibody, MAb 263

Monoclonal antibody (MAb) 263 is a widely used monoclonal antibody that recognizes the extracellular domain (ECD) of the GH receptor. It has been shown to act as a GH agonist both in vitro and in vivo, and we report here that it must be divalent to exert its effect on the full-length receptor. To understand the mechanism of its agonist action, we have determined the precise epitope for this antibody using a novel random PCR mutagenesis approach together with expression screening in yeast. A library of 5200 clones of rabbit GH receptor ECD mutants were screened both with MAb 263 and with an anticarboxy-tag antibody to verify complete ECD expression. Sequencing for clones that expressed complete ECD but were not MAb 263 positive identified 20 epitope residues distributed in a discontinuous manner throughout the ECD. The major part of the epitope, as revealed after mapping onto the crystal structure model of the ECD molecule, was located on the side and upper portion of domain 1, particularly within the D - E strand disulfide loop 79 - 96. Molecular dynamics docking of an antibody of the same isotype as MAb 263 was used to dock the bivalent antibody to the 1528-Angstrom(2) epitope and to visualize the likely consequences of MAb binding. The minimized model enables the antibody to grasp two receptors in a pincer-like movement from opposite sides, facilitating alignment of the receptor dimerization domains in a manner similar to, but not identical with, GH.

Method for determining the shear stress in cylindrical systems

We develop a method for determining the elements of the pressure tensor at a radius r in a cylindrically symmetric system, analogous to the so-called method of planes used in planar systems [B. D. Todd, Denis J. Evans, and Peter J. Daivis, Phys. Rev. E 52, 1627 (1995)]. We demonstrate its application in determining the radial shear stress dependence during molecular dynamics simulations of the forced flow of methane in cylindrical silica mesopores. Such expressions are useful for the examination of constitutive relations in the context of transport in confined systems.

Wall mediated transport in confined spaces: Exact theory for low density

We present a theory for the transport of molecules adsorbed in slit and cylindrical nanopores at low density, considering the axial momentum gain of molecules oscillating between diffuse wall reflections. Good agreement with molecular dynamics simulations is obtained over a wide range of pore sizes, including the regime of single-file diffusion where fluid-fluid interactions are shown to have a negligible effect on the collective transport coefficient. We show that dispersive fluid-wall interactions considerably attenuate transport compared to classical hard sphere theory.

Hydrodynamic origin of diffusion in nanopores

We study the transport of a subcritical Lennard-Jones fluid in a cylindrical nanopore, using a combination of equilibrium and nonequilibrium as well as dual control volume grand canonical molecular dynamics methods. We show that all three techniques yield the same value of the transport coefficient for diffusely reflecting pore walls, even in the presence of viscous transport. We also demonstrate that the classical Knudsen mechanism is not manifested, and that a combination of viscous flow and momentum exchange at the pore wall governs the transport over a wide range of densities.

Multi-scale modeling of the mechanical behavior of polymer-based materials

Polymers have become the reference material for high reliability and performance applications. In this work, a multi-scale approach is proposed to investigate the mechanical properties of polymeric based material under strain. To achieve a better understanding of phenomena occurring at the smaller scales, a coupling of a Finite Element Method (FEM) and Molecular Dynamics (MD) modeling in an iterative procedure was employed, enabling the prediction of the macroscopic constitutive response. As the mechanical response can be related to the local microstructure, which in turn depends on the nano-scale structure, the previous described multi-scale method computes the stress-strain relationship at every analysis point of the macro-structure by detailed modeling of the underlying micro- and meso-scale deformation phenomena. The proposed multi-scale approach can enable prediction of properties at the macroscale while taking into consideration phenomena that occur at the mesoscale, thus offering an increased potential accuracy compared to traditional methods.

Mechanical Behavior of the Lamellar Structure in Semi-Crystalline Polymers

We have employed molecular dynamics simulations to study the behavior of virtual polymeric materials under an applied uniaxial tensile load. Through computer simulations, one can obtain experimentally inaccessible information about phenomena taking place at the molecular and microscopic levels. Not only can the global material response be monitored and characterized along time, but the response of macromolecular chains can be followed independently if desired. The computer-generated materials were created by emulating the step-wise polymerization, resulting in self-avoiding chains in 3D with controlled degree of orientation along a certain axis. These materials represent a simplified model of the lamellar structure of semi-crystalline polymers,being comprised of an amorphous region surrounded by two crystalline lamellar regions. For the simulations, a series of materials were created, varying i) the lamella thickness, ii) the amorphous region thickness, iii) the preferential chain orientation, and iv) the degree of packing of the amorphous region. Simulation results indicate that the lamella thickness has the strongest influence on the mechanical properties of the lamella-amorphous structure, which is in agreement with experimental data. The other morphological parameters also affect the mechanical response, but to a smaller degree. This research follows previous simulation work on the crack formation and propagation phenomena, deformation mechanisms at the nanoscale, and the influence of the loading conditions on the material response. Computer simulations can improve the fundamental understanding about the phenomena responsible for the behavior of polymeric materials, and will eventually lead to the design of knowledge-based materials with improved properties.

Modelling the Mechanical Behavior of Polymer-Based Nanocomposites

Molecular dynamics simulations were employed to analyze the mechanical properties of polymer-based nanocomposites with varying nanofiber network parameters. The study was focused on nanofiber aspect ratio, concentration and initial orientation. The reinforcing phase affects the behavior of the polymeric nanocomposite. Simulations have shown that the fiber concentration has a significant effect on the properties, with higher loadings resulting in higher stress levels and higher stiffness, matching the general behavior from experimental knowledge in this field. The results also indicate that, within the studied range, the observed effect of the aspect ratio and initial orientation is smaller than that of the concentration, and that these two parameters are interrelated.

Biochemical and structural characterization of β-lactamases tem-180 and tem-201

With the constant development of new antibiotics, selective pressure is a force to reckon when investigating antibiotic resistance. Although advantageous for medical treatments, it leads to increasing resistance. It is essential to use more potent and toxic antibiotics. Enzymes capable of hydrolyzing antibiotics are among the most common ways of resistance and TEM variants have been detected in several resistant isolates. Due to the rapid evolution of these variants, complex phenotypes have emerged and the need to understand their biological activity becomes crucial. To investigate the biochemical properties of TEM-180 and TEM-201 several computational methodologies have been used, allowing the comprehension of their structure and catalytic activity, which translates into their biological phenotype. In this work we intent to characterize the interface between these proteins and the several antibiotics used as ligands. We performed explicit solvent molecular dynamics (MD) simulations of these complexes and studied a variety of structural and energetic features. The interfacial residues show a distinct behavior when in complex with different antibiotics. Nevertheless, it was possible to identify some common Hot Spots among several complexes – Lys73, Tyr105 and Glu166. The structural changes that occur during the Molecular Dynamic (MD) simulation lead to the conclusion that these variants have an inherent capacity of adapting to the various antibiotics. This capability might be the reason why they can hydrolyze antibiotics that have not been described until now to be degraded by TEM variants. The results obtained with computational and experimental methodologies for the complex with Imipenem have shown that in order to this type of enzymes be able to acylate the antibiotics, they need to be capable to protect the ligand from water molecules.