Biomedical subjects
Publications and source records attributed to T Ebisuzaki.
We extend the recently proposed order-N algorithms for calculating linear- and nonlinear-response functions in time domain to the systems described by nonorthonormal basis sets.
The Ras protein and its homolog, Rap1A, have an identical "effector region" (residues 32-40) preceded by Asp30-Glu31 and Glu30-Lys31, respectively. In the complex of the "Ras-like" E30D/K31E mutant Rap1A with the Ras-binding domain (RBD), residues 51-131 of Raf-1, Glu31 in Rap1A forms a tight salt bridge with Lys84 in Raf-1. However, we have recently found that Raf-1 RBD binding of Ras is indeed reduced by the E31K mutation, but is not affected by the E31A mutation. Here, the "Rap1A-like" D30E/E31K mutant of Ras was prepared and shown to bind the Raf-1 RBD less strongly than wild-type Ras, but slightly more tightly than the E31K mutant. The backbone 1H, 13C, and 15N magnetic resonances of the Raf-1 RBD were assigned in complexes with the wild-type and D30E/E31K mutant Ras proteins in the guanosine 5'-O-(beta,gamma-imidotriphosphate)-bound form. The Lys84 residue in the Raf-1 RBD exhibited a large change in chemical shift upon binding wild-type Ras, suggesting that Lys84 interacts with wild-type Ras. The D30E/E31K mutant of Ras caused nearly the same perturbations in Raf-1 chemical shifts, including that of Lys84. We hypothesized that Glu31 in Ras may not be the major salt bridge partner of Lys84 in Raf-1. A molecular dynamics simulation of a model structure of the Raf-1 RBD.Ras.GTP complex suggested that Lys84 in Raf-1 might instead form a tight salt bridge with Asp33 in Ras. Consistent with this, the D33A mutation in Ras greatly reduced its Raf-I RBD binding activity. We conclude that the major salt bridge partner of Lys84 in Raf-1 may be Asp33 in Ras.
We propose an efficient time-dependent algorithm for nonlinear response function that requires CPU time proportional to the system size, and study the size effects in two-photon absorption spectra of Si nanocrystallites by using this algorithm.
Molecular dynamics simulations have been extensively used in research of proteins. Since these simulations are quite computer intensive, their acceleration is of main interest of the research. In molecular dynamics simulations, almost all computing time is consumed in calculating the forces between particles, e.g., Coulomb and van der Waals forces. We have designed and built GRAPE-2A (GRAvity PipE 2A), a special-purpose computer for use in simulations of classical many-body systems. GRAPE-2A calculates forces exerted on a particle from the other particles. GRAPE-2A can calculate force of an arbitrary functional form of a central force. The host computer, which is connected to GRAPE-2A through the VME bus, performs other calculations such as time integration. The peak speed of GRAPE-2A is 180 Mflops. We can also stimulate systems with periodic boundary conditions by the Ewald method, using GRAPE-2A and another special-purpose computer, WINE (Wave space INtegrator for the Ewald method).
Explore the source record for details and available documents.