Physical properties and phase separation structure of Antheraea pernyi/Bombyx mori silk fibroin blend films

The physical structure and compatibility of solution-cast Antheraea pernyi/Bombyx mori silk fibroin blend films were stuided by differential scanning calorimetry (DSC), thermomechanical (TMA) and thermogravimetric (TGA) analysis, dynamic viscoelastic measurement, infrared spectroscopy, and x-ray diffractometry. The DSC curves of the blend films showed independent endotherms at 280 and 358°C, corresponding to the thermal decomposition of B. mori and A. pernyi silk fibroins with random coil conformation. The intensity was roughly proportionate to the amount of each component in the blend. The thermal behavior corresponding to the conformational transitions induced by heating on A. pernyi and B. mori silk fibroins overlapped in the temperature range 190–230°C. Thermal expansion and contraction properties, as well as weight retention behavior of the blend films were intermediate between the pure components, as shown by the TMA and TGA curves. The onset temperature of the storage modulus curve decreased markedly, approaching that of B. mori silk fibroin film when the amount of this component in the blend increased. The loss modulus curve of the blend films showed two peaks at ca. 190 and 210°C, the former corresponding to B. mori, and the latter to A. pernyi silk fibroin. Infrared spectra of the blends exhibited absorption bands characteristic of the pure components overlapping in the spectral region 2000–400 cm^?1. The x-ray diffraction peaks at 23 and 21.5°, attributed to the crystalline spacings of A. pernyi and B. mori fibroins, respectively, overlapped in the diffraction curves of the blends, while the peak at 11.4°, of A. pernyi, increased as the content of this fibroin in the blend increased. The degree of crystallinity, calculated from the x-ray diffraction curves, diminished as the amount of B. mori silk fibroin decreased. A low degree of compatibility exists between the two fibroins when they are cast from aqueous solution in the experimental conditions adopted in this work. © 1994 John Wiley & Sons, Inc.     

Origin of attraction and directionality of the pi/pi interaction: model chemistry calculations of benzene dimer interaction.

A model chemistry for the evaluation of intermolecular interaction between aromatic molecules (AIMI Model) has been developed. The CCSD(T) interaction energy at the basis set limit has been estimated from the MP2 interaction energy near the basis set limit and the CCSD(T) correction term obtained by using a medium size basis set. The calculated interaction energies of the parallel, T-shaped,and slipped-parallel benzene dimers are -1.48, -2.46, and -2.48 kcal/mol, respectively. The substantial attractive interaction in benzene dimer, even where the molecules are well separated, shows that the major source of attraction is not short-range interactions such as charge-transfer but long-range interactions such as electrostatic and dispersion. The inclusion of electron correlation increases attraction significantly. The dispersion interaction is found to be the major source of attraction in the benzene dimer. The orientation dependence of the dimer interaction is mainly controlled by long-range interactions. Although electrostatic interaction is considerably weaker than dispersion interaction, it is highly orientation dependent. Dispersion and electrostatic interactions are both important for the directionality of the benzene dimer interaction.     

Molecular dynamics study of the methanol effect on the membrane morphology of perfluorosulfonic ionomers

The membranes of a perfluorosulfonic acid polymer swollen in 10-80 wt % methanol solution were investigated to elucidate the methanol effect on their morphologies, such as size of the solvent cluster, solvent location, and polymer structure, by using isothermal-isobaric molecular dynamics simulations. In higher methanol concentrations, we found less-spherical solvent aggregation and a more spread polymer structure because of the ampholytic nature of methanol. The partial radial distribution functions between solvent oxygen and fluorocarbons, which are composed of the main chain, clearly show that methanol is located closer to the polymer matrix than water. On the other hand, water is preferentially located in the vicinity of an acidic headgroup, SO(3)(-), compared with methanol, although both have similar attractive interaction energies to the acidic group. Furthermore, we discussed solvent dynamics and hydrogen bonding between sulfonic oxygen and solvent O-H groups.     

Rigid-body dynamics in the isothermal-isobaric ensemble: a test on the accuracy and computational efficiency

We have developed a time-reversible rigid-body (rRB) molecular dynamics algorithm in the isothermal-isobaric (NPT) ensemble. The algorithm is an extension of rigid-body dynamics [Matubayasi and Nakahara, J Chem Phys 1999, 110, 3291] to the NPT ensemble on the basis of non-Hamiltonian statistical mechanics [Martyna, G. J. et al., J Chem Phys 1994, 101, 4177]. A series of MD simulations of water as well as fully hydrated lipid bilayer systems have been undertaken to investigate the accuracy and efficiency of the algorithm. The rRB algorithm was shown to be superior to the state-of-the-art constraint-dynamics algorithm SHAKE/RATTLE/ROLL, with respect to computational efficiency. However, it was revealed that both algorithms produced accurate trajectories of molecules in the NPT as well as NVT ensembles, as long as a reasonably short time step was used. A couple of multiple time-step (MTS) integration schemes were also examined. The advantage of the rRB algorithm for computational efficiency increased when the MD simulation was carried out using MTS on parallel processing computer systems; total computer time for MTS-MD of a lipid bilayer using 64 processors was reduced by about 40% using rRB instead of SHAKE/RATTLE/ROLL.     

Structural changes in poly(L-alanine) induced by heat treatment

Structural changes in poly(L -alanine)(PLA), a model compound related to tussah silk fibroin, induced by heat treatment have been studied by differential thermal analysis, x-ray diffraction, and infrared spectroscopy. PLA heated below 300°C shows x-ray patterns very similar to those of the α-helix crystalline phase, in addition to the diffraction patterns due to the β structure. Samples heated at 368°C exhibit predominantly the diffraction patterns due to the β crystalline phase. From infrared spectra, PLA samples heated below 280°C are found to be composed of all three molecular conformations: β sheet, random coil, and α helix. The intensity of the infrared absorption at 650 cm^?1 (amide V) assigned to the random coil conformation is decreased by heat treatment above 325°C. The content of the β-sheet structure remains almost constant when the specimen is heated below 325°C, and increases abruptly on heating to around 335°C, while the random coil content is decreased abruptly by heat treatment above 335°C. The α-helix content does not change, regardless of heat treatment. It is suggested therefore that the random coil conformation of PLA transforms directly into the β sheet on heat treatment above 335°C.     

Unusual stereoselectivity in the synthesis of 4-oxa-δ-valerothiolactones: structure proof of β-hydroxythioamide

Stereoselectivity in the formation of 4-oxa-δ-valerothiolactone 2 from β-hydroxythioamide 1 is discussed. The structure of 1 is deduced from the structure of 2.     

Physical and chemical properties of tussah silk fibroin films

Physical and chemical structure, as well as thermal behavior of solution-cast regenerated films, prepared from tussah (Antheraea pernyi) silk fibroin, were compared with those of solution-cast native films, in order to ascertain whether treatment (degumming, dissolution) used for preparation affected their properties. Regenerated fibroin films exhibited a higher thermal stability than native ones, as shown by differential scanning calorimetry, thermomechanical analysis, and dynamic mechanical behavior. Glass transition temperature and other relevant thermal transitions of the regenerated silk specimen shifted to higher temperatures compared with those of native specimen. Molecular conformation and crystalline structure did not show significant differences between the two kinds of silk films. Amino acid composition and molecular weight, however, distribution changed markedly after dissolving tussah silk fibroin fiber in concentrated LiSCN in polypeptide size was the main features for the regenerated silk fibroin. © 1994 John Wiley & Sons, Inc.     

High field magnetization of Y2Co7 and YCo3 hydrides

The magnetization curves for the Y₂Co₇H_x and the YCo₃H_x systems have been measured at 4.2 K in the pulsed high magnetic fields up to 280 kOe. The metamagnetic transition is observed in the β and the γ hydrides except for Y₂Co₇H₃ and YCo₃H₁, which is interpreted in terms of the itinerant electron metamagnetism.     

Efficient free energy calculation of water across lipid membranes

An efficient free energy (FE) calculation of a water molecule to go across lipid membranes is presented. Both overlapping distribution and cavity insertion Widom methods are complementarily used. The former is useful for a dense region where water molecules are found, i.e., from the interfacial to bulk water region, while the latter works well in the low density region, i.e., the hydrocarbon region. Since both methods evaluate the excess chemical potential of water, the obtained FE profile is free from the fitting problem usually arisen when a FE difference method is used. A diphytanyl phosphatidylcholine bilayer is used for our test calculations. An excellent and fast convergence of the chemical potential is obtained when each method is applied for the appropriate region. The estimated FE barrier using the Ewald method for the electrostatic interaction is approximately 7.2 kcal/mol, which is higher than that using the interaction cutoff of 20 A by about 0.9 kcal/mol.     

Computationally efficient canonical molecular dynamics simulations by using a multiple time‐step integrator algorithm combined with the particle mesh Ewald method and with the fast multipole method

An efficient implementation of the canonical molecular dynamics simulation using the reversible reference system propagator algorithm (r‐RESPA) combined with the particle mesh Ewald method (PMEM) and with the macroscopic expansion of the fast multipole method (MEFMM) was examined. The performance of the calculations was evaluated for systems with 3000, 9999, 30,000, 60,000, and 99,840 particles. For a given accuracy, the optimal conditions for minimizing the CPU time for the implementation of the Ewald method, the PMEM, and the MEFMM were first analyzed. Using the optimal conditions, we evaluated the performance and the reliability of the integrated methods. For all the systems examined, the r‐RESPA with the PMEM was about twice as fast as the r‐RESPA with the MEFMM. The difference arose from the difference in the numerical complexities of the fast Fourier transform in the PMEM and from the transformation of the multipole moments into the coefficients of the local field expansion in the MEFMM. Compared with conventional methods, such as the velocity‐verlet algorithm with the Ewald method, significant speedups were obtained by the integrated methods; the speedup of the calculation was a function of system size, and was a factor of 100 for a system with 3000 particles and increased to a factor of 700 for a system with 99,840 particles. These integrated calculations are, therefore, promising for realizing large‐scale molecular dynamics simulations for complex systems. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 201–217, 2000