Conformations of Three Types of Ultra-Long-Chain Fatty Acids in Multicomponent Lipid Bilayers.

Ultra-long-chain fatty acids (ULCFAs) are biosynthesized in certain types of tissues, but their biological roles remain unknown. Here, we report how the conformation of ULCFAs depends on the length and unsaturated-bond ratio of the ultra-long chains and the composition of the host bilayer membrane using molecular dynamics simulations. The ultra-long chain of ULCFAs flips between the two leaflets and fluctuates among three conformations: elongated, L-shaped, and turned. Furthermore, we found that the saturated ultra-long chain exhibited an elongated conformation more frequently than the unsaturated chain. In addition, the truncation of the ultra-long chain at C26 had little effect on the remaining ULCFAs. ULCFAs respond to lipid-density differences in the two leaflets, and the ratio of the elongated and turned conformations changed to reduce this difference. However, in cholesterol-containing membranes, ULCFAs exhibit no density difference after the flip-flop of cholesterol removes the difference.

Conformation of ultra-long-chain fatty acid in lipid bilayer: Molecular dynamics study.

Ultra-long-chain fatty acids (ULCFAs) are biosynthesized in the restricted tissues such as retina, testis, and skin. The conformation of a single ULCFA, in which the sn-1 unsaturated chain has 32 carbons, in three types of phospholipid bilayers is studied by molecular dynamics simulations. It is found that the ultra-long tail of the ULCFA flips between two leaflets and fluctuates among an elongation into the opposite leaflet, lies between two leaflets, and turns back. As the number ratio of lipids in the opposite leaflet increases, the ratio of the elongated shape linearly decreases in all three cases. Thus, ULCFAs can sense the density differences between the two leaflets and respond to these changes.

Effects of Hydrophilic Residues and Hydrophobic Length on Flip-Flop Promotion by Transmembrane Peptides.

Peptide-induced phospholipid flip-flop (scrambling) was evaluated using transmembrane model peptides in which the central residue was substituted with various amino acid residues (sequence: Ac-GKK(LA) nXW(LA) nLKKA-CONH2). Peptides with a strongly hydrophilic residue (X = Q, N, or H) had higher scramblase activity than that of other peptides, and the activity was also dependent on the length of the peptides. Peptides with a hydrophobic stretch of 17 residues showed high flip-promotion propensity, whereas those of 21 and 25 residues did not, suggesting that membrane thinning under negative mismatch conditions promotes the flipping. Interestingly, a hydrophobic stretch of 19 residues intensively promoted phospholipid scrambling and membrane leakage. The distinctive characteristics of the peptide were ascribed by long-term molecular dynamics simulation to the arrangement of central glutamine and terminal four lysine residues on the same side of the helix. The combination of simulated and experimental data enables understanding of the mechanisms by which transmembrane helices, and ultimately unidentified scramblases in biomembranes, cause lipid scrambling.

MODYLAS: A Highly Parallelized General-Purpose Molecular Dynamics Simulation Program for Large-Scale Systems with Long-Range Forces Calculated by Fast Multipole Method (FMM) and Highly Scalable Fine-Grained New Parallel Processing Algorithms.

Our new molecular dynamics (MD) simulation program, MODYLAS, is a general-purpose program appropriate for very large physical, chemical, and biological systems. It is equipped with most standard MD techniques. Long-range forces are evaluated rigorously by the fast multipole method (FMM) without using the fast Fourier transform (FFT). Several new methods have also been developed for extremely fine-grained parallelism of the MD calculation. The virtually buffering-free methods for communications and arithmetic operations, the minimal communication latency algorithm, and the parallel bucket-relay communication algorithm for the upper-level multipole moments in the FMM realize excellent scalability. The methods for blockwise arithmetic operations avoid data reload, attaining very small cache miss rates. Benchmark tests for MODYLAS using 65 536 nodes of the K-computer showed that the overall calculation time per MD step including communications is as short as about 5 ms for a 10 million-atom system; that is, 35 ns of simulation time can be computed per day. The program enables investigations of large-scale real systems such as viruses, liposomes, assemblies of proteins and micelles, and polymers.

Inverted micelle formation of cell-penetrating peptide studied by coarse-grained simulation: importance of attractive force between cell-penetrating peptides and lipid head group.

Arginine-rich peptide and Antennapedia are cell-penetrating peptides (CPPs) which have the ability to permeate plasma membrane. Deformation of the plasma membrane with CPPs is the key to understand permeation mechanism. We investigate the dynamics of CPP and the lipid bilayer membrane by coarse-grained simulation. We found that the peptide makes inverted micelle in the lipid bilayer membrane, when the attractive potential between the peptide and lipid heads is strong. The inverted micelle is formed to minimize potential energy of the peptide. For vesicle membrane, the peptide moves from the outer vesicle to the inner vesicle through the membrane. The translocation of the peptide suggests inverted micelle model as a possible mechanism of CPPs.

In silico and in vitro approaches to elucidate the thermal stability of human UDP-glucuronosyltransferase (UGT) 1A9.

UDP-Glucuronosyltransferases (UGTs) are predominant drug metabolizing enzymes in the liver and extrahepatic tissues. Human UGT1A9 is uniquely stable against heat treatment. To understand the unique properties of UGT1A9, the three-dimensional structure was constructed by homology modeling using a crystal structure of TDP-epi-vancosaminyltransferase as template. Sequence alignment analysis revealed that 13 amino acid residues (Arg42, Lys91, Ala92, Tyr106, Gly111, Tyr113, Asp115, Asn152, Leu173, Leu219, His221, Arg222, and Glu241) are unique to UGT1A9 as compared with UGT1A7, UGT1A8 and UGT1A10. To examine the roles of these residues in the conformational stability of UGT1A9, molecular dynamics simulation of the structures was carried out at 310 K and 360 K in aqueous solution for 3.0 nanoseconds. Root mean square deviation analyses revealed that Arg42, Leu173, Leu219, His221 and Arg222 were responsible for the thermal stability. Root mean square fluctuation analyses and a dynamical cross correlation map revealed that Lys91, Ala92, Tyr106, Gly111, Tyr113, Asp115, Leu219, His221, Arg222 and Glu241 were responsible for the thermal stability. In vitro study using mutants of these residues demonstrated that all these amino acids may be collectively involved in the thermal stability of UGT1A9. The results presented here provide a molecular basis for the thermal stability of human UGT1A9.

Electronic and geometric structures of the blue copper site of azurin investigated by QM/MM hybrid calculations.

The electronic and geometric structures of the copper-binding site in a fully solvated azurin were investigated using quantum mechanics (QM) and molecular mechanics (MM) hybrid calculations. Two types of computational models were applied to evaluate the effects of the environment surrounding the active site. In model I, long-distance electrostatic interactions between QM region atoms and partial point charges of the surrounding protein moieties and solvent water were calculated in a QM Hamiltonian, for which the spin-unrestricted Hartree-Fock (UHF)/density functional theory (DFT) hybrid all-electron calculation with the B3LYP functional was adopted. In model II, the QM Hamiltonian was not allowed to be polarized by those partial point charges. Models I and II provided different descriptions of the copper coordination structure, particularly for the coordinative bonds including a large dipole. In fact, the Cu-O(Gly45) and Cu-S(Cys112) bonds are sensitive to the treatment of long-distance electrostatic interactions in the QM Hamiltonian. This suggests that biological processes occurring in the active site are regulated by the surrounding structures of protein and solvent, and therefore the effects of long-range electrostatic interactions involved in the QM Hamiltonian are crucial for accurate descriptions of electronic structures of the copper active site of metalloenzymes.

Theoretical study of above-threshold dissociation on diatomic molecules by using nonresonant intense laser pulses.

The above-threshold dissociation of the ground state of a OH molecule under intense nonresonant laser pulses has been studied using the time-dependent Schrödinger equation with discrete variable representation. The applied field is assumed as a two-color mixed nonresonant laser pulses which has the nonresonant frequency omega and the overtone 2omega. After modulating the relative phase factor between the omega and 2omega pulse, we extracted a three-photon absorption peak or a five-photon absorption peak in the ATD spectrum.