Structure and Mechanical Response of Polybutylcarbosilane Dendrimers Confined in a Flat Slit: Effect of Molecular Architecture and Generation Number.

Due to the absence of specific interactions, carbosilane dendrimers are ideal models to study the effect of a hyperbranched regular structure on the molecular response to external influences. In this work, we have studied the conformational behavior of single polybutylcarbosilane dendrimers under confinement between impermeable flat surfaces using atomistic molecular dynamics simulations. Dendrimers of different generations belonging to two homologous series with a tetra-functional core and three- and four-functional branches were simulated. The analysis of the dependence of the internal energy of the dendrimers on the wall distance allowed us to determine the critical degree of compression at which the dendrimers are able to change their shape without energy loss. The effects of generation number and branching functionality on the number of wall contacts, density distribution and shape changes were elucidated. It was found that for high generation dendrimers, the inner layers are not accessible for external interaction. It was shown that the excess stresses occurring at high compressions are concentrated in the structural center of the dendrimer. The nature of the elastic response, which is strongly nonlinear, was analyzed at different compressions depending on the dendrimer architecture and generation. We believe that our results are useful for further studies of dendrimer films under compression and can also serve as a basis for developing model concepts to describe the dynamics of dendrimer melts.

Is It Possible to Find an Antimicrobial Peptide That Passes the Membrane Bilayer with Minimal Force Resistance? An Attempt at a Predictive Approach by Molecular Dynamics Simulation.

There is still no answer to the mechanism of penetration of AMP peptides through the membrane bilayer. Several mechanisms for such a process have been proposed. It is necessary to understand whether it is possible, using the molecular dynamics method, to determine the ability of peptides of different compositions and lengths to pass through a membrane bilayer. To explain the passage of a peptide through a membrane bilayer, a method for preparing a membrane phospholipid bilayer was proposed, and 656 steered molecular dynamics calculations were carried out for pulling 7 amyloidogenic peptides with antimicrobial potential, and monopeptides (homo-repeats consisting of 10 residues of the same amino acid: Poly (Ala), Poly (Leu), Poly (Met), Poly (Arg), and Poly (Glu)) with various sequences through the membrane. Among the 15 studied peptides, the peptides exhibiting the least force resistance when passing through the bilayer were found, and the maximum reaction occurred at the boundary of the membrane bilayer entry. We found that the best correlation between the maximum membrane reaction force and the calculated parameters corresponds to the instability index (the correlation coefficient is above 0.9). One of the interesting results of this study is that the 10 residue amyloidogenic peptides and their extended peptides, with nine added residue cell-penetrating peptides and four residue linkers, both with established antimicrobial activity, have the same bilayer resistance force. All calculated data are summarized and posted on the server.

Signal transfer in human protein tyrosine phosphatase PTP1B from allosteric inhibitor P00058.

Protein tyrosine phosphatases constitute a family of cytosolic and receptor-like signal transducing enzymes that catalyze the hydrolysis of phospho-tyrosine residues of phosphorylated proteins. PTP1B, encoded by PTPN1, is a key negative regulator of insulin and leptin receptor signaling, linking it to two widespread diseases: type 2 diabetes mellitus and obesity. Here, we present crystal structures of the PTP1B apo-enzyme and a complex with a newly identified allosteric inhibitor, 2-(2,5-dimethyl-pyrrol-1-yl)-5-hydroxy-benzoic acid, designated as P00058. The inhibitor binding site is located about 18 Å away from the active center. However, the inhibitor causes significant re-arrangements in the active center of enzyme: residues 45-50 of catalytic Tyr-loop are shifted at their Cα-atom positions by 2.6 to 5.8 Å. We have identified an event of allosteric signal transfer from the inhibitor to the catalytic area using molecular dynamic simulation. Analyzing change of complex structure along the fluctuation trajectory we have found the large Cα-atom shifts in external strand, residues 25-40, which occur at the same time with the shifts in adjacent catalytic p-Tyr-loop. Coming of the signal to this loop arises due to dynamic fluctuation of protein structure at about 4.0 nanoseconds after the inhibitor takes up its space. Communicated by Ramaswamy H. Sarma.

Adsorption of Silicon-Containing Dendrimers: Effects of Chemical Composition, Structure, and Generation Number.

We studied the conformational behavior of silicon-containing dendrimers during their adsorption onto a flat impenetrable surface by molecular dynamics (MD) simulations. Four homologous series of dendrimers from the 4th up to the 7th generations were modeled, namely, two types of carbosilane dendrimers differing by the functionality of the core Si atom and two types of siloxane dendrimers with different lengths of the spacers. Comparative analysis of the fractions of adsorbed atoms belonging to various structural layers within dendrimers as well as density profiles allowed us to elucidate not only some general trends but also the effects determined by dendrimer specificity. In particular, it was found that in contrast to the carbosilane dendrimers interacting with the adsorbing surface mainly by their peripheral layers, the siloxane dendrimers with the longer -O-Si(CH3)2-O spacers expose atoms from their interior to the surface spreading out on it. These findings are important for the design of functional materials on the basis of silicon-containing dendrimers.

Hybrid Polycarbosilane-Siloxane Dendrimers: Synthesis and Properties.

A series of carbosilane dendrimers of the 4th, 6th, and 7th generations with a terminal trimethylsilylsiloxane layer was synthesized. Theoretical models of these dendrimers were developed, and equilibrium dendrimer conformations obtained via molecular dynamics simulations were in a good agreement with experimental small-angle X-ray scattering (SAXS) data demonstrating molecule monodispersity and an almost spherical shape. It was confirmed that the glass transition temperature is independent of the dendrimer generation, but is greatly affected by the chemical nature of the dendrimer terminal groups. A sharp increase in the zero-shear viscosity of dendrimer melts was found between the 5th and the 7th dendrimer generations, which was qualitatively identical to that previously reported for polycarbosilane dendrimers with butyl terminal groups. The viscoelastic properties of high-generation dendrimers seem to follow some general trends with an increase in the generation number, which are determined by the regular branching structure of dendrimers.

Effects of generation number, spacer length and temperature on the structure and intramolecular dynamics of siloxane dendrimer melts: molecular dynamics simulations.

The structure and properties of polysiloxane dendrimer melts are studied by extensive atomistic molecular dynamics simulations. Two homologous series differing in the spacer length are considered. In the first series the dendrimer spacers are the shortest ones, comprising only one oxygen atom, while in the second series the spacers consist of two oxygen atoms with the silicon atom in between. Melts of the dendrimers from the 3rd up to the 6th generation number are modelled in a wide temperature range from 273 to 600 K. A comparative study of the macroscopic melt characteristics such as the melt density and thermal expansion coefficients is performed for the two series. Analysis of the dendrimer structure in melts and in the isolated state shows that intermolecular interactions and interpenetration of dendrimer molecules in melts hardly affect the dendrimer interior organization. However, the presence of neighboring molecules significantly slows down their intramolecular dynamics in melts in comparison with that of isolated dendrimers. An increasing generation number causes an increase of the radius of the dendrimer interior region unavailable for neighboring molecules, which starts to exceed the length of the peripheral interpenetration layer for high-generation dendrimers; this fact could lead to different mechanisms of melt dynamics for lower and higher generation dendrimers.

A Comparative Study of Intramolecular Mobility of Single Siloxane and Carbosilane Dendrimers via Molecular Dynamics Simulations.

A comparative analysis of intramolecular dynamics of four types of isolated dendrimers from the fourth to the seventh generations belonging to the siloxane and carbosilane families, differing in spacer length, core functionality, and the type of chemical bonds, has been performed via atomic molecular dynamics simulations. The average radial and angular positions of all Si branching atoms of various topological layers within the dendrimer interior, as well as their variations, have been calculated, and the distributions of the relaxation times of their radial and angular motions have been found. It has been shown that the dendrons of all the dendrimers elongate from the center and decrease in a solid angle with an increasing generation number. The characteristic relaxation times of both angular and radial motions of Si atoms are of the order of a few nanoseconds, and they increase with an increasing generation number and decrease with temperature, with the angular relaxation times being larger than the radial ones. The relaxation times in the carbosilanes are larger than those in the siloxanes. The rotational angle dynamics of the carbosilane dendrimers show that the chain bending is mainly realized via trans-gauche transitions in the Si branching bonds.

Comparative mechanical unfolding studies of spectrin domains R15, R16 and R17.

Spectrins belong to repetitive three-helix bundle proteins that have vital functions in multicellular organisms and are of potential value in nanotechnology. To reveal the unique physical features of repeat proteins we have studied the structural and mechanical properties of three repeats of chicken brain α-spectrin (R15, R16 and R17) at the atomic level under stretching at constant velocities (0.01, 0.05 and 0.1 ŷps-1) and constant forces (700 and 900 pN) using molecular dynamics (MD) simulations at T = 300 K. 114 independent MD simulations were performed and their analysis has been done. Despite structural similarity of these domains we have found that R15 is less mechanically stable than R16, which is less stable than R17. This result is in agreement with the thermal unfolding rates. Moreover, we have observed the relationship between mechanical stability, flexibility of the domains and the number of aromatic residues involved in aromatic clusters.

Interfacial and topological effects on the glass transition in free-standing polystyrene films.

United-atom molecular-dynamics computer simulations of atactic polystyrene (PS) were performed for the bulk and free-standing films of 2 nm-20 nm thickness, for both linear and cyclic polymers comprised of 80 monomers. Simulated volumetric glass-transition temperatures (Tg) show a strong dependence on the film thickness below 10 nm. The glass-transition temperature of linear PS is 13% lower than that of the bulk for 2.5 nm-thick films, as compared to less than 1% lower for 20 nm films. Our studies reveal that the fraction of the chain-end groups is larger in the interfacial layer with its outermost region approximately 1 nm below the surface than it is in the bulk. The enhanced population of the end groups is expected to result in a more mobile interfacial layer and the consequent dependence of Tg on the film thickness. In addition, the simulations show an enrichment of backbone aliphatic carbons and concomitant deficit of phenyl aromatic carbons in the interfacial film layer. This deficit would weaken the strong phenyl-phenyl aromatic (π-π) interactions and, hence, lead to a lower film-averaged Tg in thin films, as compared to the bulk sample. To investigate the relative importance of the two possible mechanisms (increased chain ends at the surface or weakened π-π interactions in the interfacial region), the data for linear PS are compared with those for cyclic PS. For the cyclic PS, the reduction of the glass-transition temperature is also significant in thin films, albeit not as much as for linear PS. Moreover, the deficit of phenyl carbons in the film interface is comparable to that observed for linear PS. Therefore, chain-end effects alone cannot explain the observed pronounced Tg dependence on the thickness of thin PS films; the weakened phenyl-phenyl interactions in the interfacial region seems to be an important cause as well.

Mechanical stability analysis of the protein L immunoglobulin-binding domain by full alanine screening using molecular dynamics simulations.

This article is the first to study the mechanical properties of the immunoglobulin-binding domain of protein L (referred to as protein L) and its mutants at the atomic level. In the structure of protein L, each amino acid residue (except for alanines and glycines) was replaced sequentially by alanine. Thus, 49 mutants of protein L were obtained. The proteins were stretched at their termini at constant velocity using molecular dynamics simulations in water, i.e. by forced unfolding. 19 out of 49 mutations resulted in a large decrease of mechanical protein stability. These amino acids were affecting either the secondary structure (11 mutations) or loop structures (8 mutations) of protein L. Analysis of mechanical unfolding of the generated protein that has the same topology as protein L but consists of only alanines and glycines allows us to suggest that the mechanical stability of proteins, and specifically protein L, is determined by interactions between certain amino acid residues, although the unfolding pathway depends on the protein topology. This insight can now be used to modulate the mechanical properties of proteins and their unfolding pathways in the desired direction for using them in various biochips, biosensors and biomaterials for medicine, industry, and household purposes.

Right- and left-handed three-helix proteins. II. Similarity and differences in mechanical unfolding of proteins.

Here, we study mechanical properties of eight 3-helix proteins (four right-handed and four left-handed ones), which are similar in size under stretching at a constant speed and at a constant force on the atomic level using molecular dynamics simulations. The analysis of 256 trajectories from molecular dynamics simulations with explicit water showed that the right-handed three-helix domains are more mechanically resistant than the left-handed domains. Such results are observed at different extension velocities studied (192 trajectories obtained at the following conditions: v = 0.1, 0.05, and 0.01 Å ps(-1) , T = 300 K) and under constant stretching force (64 trajectories, F = 800 pN, T = 300 K). We can explain this by the fact, at least in part, that the right-handed domains have a larger number of contacts per residue and the radius of cross section than the left-handed domains.

Computer simulation of gas-phase neutralization of electrospray-generated protein macroions.

The process of neutralizing hydrated multicharged gas-phase protein ions with small counterions was simulated using a molecular dynamics (MD) technique. Hen egg white lysozyme (HEWL) molecules with different numbers of positive charges, both dry and solvated by up to 1500 water molecules, were first equilibrated. Simulations revealed that the hydration layer over a highly charged protein surface adapted a spiny structure with water protrusions composed of oriented water dipoles. MD simulations of the neutralization process showed that the impact of a small dehydrated single-charged counterion with a dehydrated HEWL ion bearing eight uncompensated charges resulted in a short local increase in temperature by 600-1000 K, which quickly (in 3-5 ps) dissipated over the whole protein molecule, increasing its average temperature by 20-25 K. When the protein ion was solvated, no drastic local increase in the temperature of the protein atoms was observed, because the impact energy was dissipated among the water molecules near the collision site.

Two-, three-, and four-state events occur in the mechanical unfolding of small protein L using molecular dynamics simulation.

Mechanical properties of (protein L)(5) have been recently investigated by single-molecule force spectroscopy. It has been demonstrated that the unfolding of individual domains proceeds through a two-state mechanism. Here, we study mechanical properties of protein L at the atomic level under stretching at constant velocity using molecular dynamics simulations. We have found that the unfolding process of protein L can occur either in a single step or through short living and quite native like intermediate states, which was not observed in previous studies. Analysis of the 24 trajectories from molecular dynamics simulations with explicit water showed that the mechanical unfolding of protein L occurs through at least two pathways. These pathways coincide in two- and multi-state events and at different extension velocities studied (0.125, 0.0625 and 0.005 A x ps(-1)).

Molecular dynamics simulation of thermomechanical properties of montmorillonite crystal. 3. montmorillonite crystals with PEO oligomer intercalates.

We present the results of molecular dynamics (MD) simulation of the structure and thermomechanical behavior of Wyoming-type Na+-montmorillonite (MMT) with poly(ethylene oxide) (PEO) oligomer intercalates. Periodic boundary conditions in all three directions and simulation cells containing two MMT lamellae [Si248Al8][Al112Mg16]O640[OH]128 oriented parallel to the XY-plane were used. The interlamellar space, or gallery, between neighboring MMT lamellae was populated by 24 Na+ counterions and PEO macromolecules of different lengths, ranging from 2 up to 240 repeat units. We considered three different loadings of PEO within the gallery: 80, 160, and 240 repeat units, corresponding to 13, 23, and 31 wt % PEO based on total mass of the nanocomposite, respectively. In the cases of 13 and 23 wt %, the polymer chains formed one or two well-defined amorphous layers with interlayer distances of 1.35 and 1.8 nm, respectively. We have observed also formation of a wider monolayer gallery with interlayer distances of 1.6 nm. Three-layer PEO films formed in the case of 31 wt % loading. The thermal properties were analyzed over the range 300-400 K, and the isothermal linear compressibility, transversal moduli, and shear moduli were calculated at 300 K. These properties are compared with the results of our simulation of thermal and mechanical properties of MMT crystal with galleries filled by one or two water layers as well as with those of an isolated clay nanoplate.

Molecular dynamics simulation of thermomechanical properties of montmorillonite crystal. 1. Isolated clay nanoplate.

The structure and mechanical properties of clay nanoparticles is a subject of growing interest because of their numerous applications in engineering. We present the results of molecular dynamics simulation for a single nanoplate of pyrophyllite - a 2:1 clay mineral consisting of two tetrahedral sheets of SiO4 and an intervening octahedral AlO6 sheet. Simulations were performed in the temperature interval from 5 to 750 K using the ionic-type potentials of Cygan et al. On this basis the temperature dependences of structural parameters, characterizing both tetrahedral and octahedral sheets as well as single lamella, have been studied. Two slightly different structures were observed in this wide temperature interval. The mechanical properties of the nanoplate were calculated from stress-strain diagrams, which have been obtained at relatively slow rates of deformation (for molecular simulations). Using different types of loading, we calculated the full elasticity tensor and estimated the influence of temperature on its components. We estimated also the bending and torsion stiffnesses of the nanoplate as specific characteristics of this type of particle. Because the nanoplate is atomically thin, a reasonable determination of the thickness is a nontrivial problem, both in the modeling of mechanical properties and in physical interpretation of the obtained data. We propose a procedure for its calculation.