Swelling behavior of the cellulose Ibeta crystal models by molecular dynamics

The various crystal models of cellulose Ibeta, each differing in crystal size, have been studied by computer simulation using the amber molecular-dynamics package and the GLYCAM parameters. The four types of crystal model were constructed by a combination of two base-plane sizes, consisting of either 24 or 48 chains and two chain lengths having either 10 or 20 residues. The base planes of the crystal models were composed by the edges of the [1,1,0], [1,-1,0], and [1,0,0] crystal planes, where the [1,1,0] plane was assigned to the longest edge. The crystal models were soaked in water boxes to investigate their swelling behavior. Unexpectedly, the crystal models twisted quickly to form a slightly right-handed shape during the initial approximately 50 ps and that, in a steady, swollen state, the twisted forms remained for the rest of the simulation time. In spite of such overall deformation, the inner part of the swollen model fairly reproduced the important structural features of the original crystal structure, such as the rotational positions of the substituent groups and the hydrogen-bonding scheme. On heating the crystal model up to 550 K, the twisted shape was conserved in most of the temperature range, while the initial conformations of the substituent groups deviated above approximately 430 K, followed by appreciable disordering in chain sheets at higher temperatures. It is suggested that some internal tensions are involved within a chain sheet of the initial structure. In the course of swelling, some of these tensions were released to introduce a twisted shape in the crystal models.     

Refinement of homology-based protein structures by molecular dynamics simulation techniques

The use of classical molecular dynamics simulations, performed in explicit water, for the refinement of structural models of proteins generated ab initio or based on homology has been investigated. The study involved a test set of 15 proteins that were previously used by Baker and coworkers to assess the efficiency of the ROSETTA method for ab initio protein structure prediction. For each protein, four models generated using the ROSETTA procedure were simulated for periods of between 5 and 400 nsec in explicit solvent, under identical conditions. In addition, the experimentally determined structure and the experimentally derived structure in which the side chains of all residues had been deleted and then regenerated using the WHATIF program were simulated and used as controls. A significant improvement in the deviation of the model structures from the experimentally determined structures was observed in several cases. In addition, it was found that in certain cases in which the experimental structure deviated rapidly from the initial structure in the simulations, indicating internal strain, the structures were more stable after regenerating the side-chain positions. Overall, the results indicate that molecular dynamics simulations on a tens to hundreds of nanoseconds time scale are useful for the refinement of homology or ab initio models of small to medium-size proteins.     

The molecular signals that regulate activity-dependent synapse refinement in the brain

The formation of appropriate synaptic connections is critical for the proper functioning of the brain. Early in development, neurons form a surplus of immature synapses. To establish efficient, functional neural networks, neurons selectively stabilize active synapses and eliminate less active ones. This process is known as activity-dependent synapse refinement. Defects in this process have been implicated in neuropsychiatric disorders such as schizophrenia and autism. Here we review the manner and mechanisms by which synapse elimination is regulated through activity-dependent competition. We propose a theoretical framework for the molecular mechanisms of synapse refinement, in which three types of signals regulate the refinement. We then describe the identity of these signals and discuss how multiple molecular signals interact to achieve appropriate synapse refinement in the brain.     

Full-matrix least-squares refinement of lysozymes and analysis of anisotropic thermal motion

Crystal structures of turkey egg lysozyme (TEL) and human lysozyme (HL) were refined by full-matrix least-squares method using anisotropic temperature factors. The refinement converged at the conventional R-values of 0.104 (TEL) and 0.115 (HL) for reflections with F_o > 0 to the resolution of 1.12 Å and 1.15 Å, respectively. The estimated r.m.s. coordinate errors for protein atoms were 0.031 Å (TEL) and 0.034 Å (HL). The introduction of anisotropic temperature factors markedly reduced the R-value but did not significantly affect the main chain coordinates. The degree of anisotropy of atomic thermal motion has strong positive correlation with the square of distance from the molecular centroid. The ratio of the radial component of thermal ellipsoid to the r.m.s. magnitude of three principal components has negative correlation with the distance from the molecular centroid, suggesting the domination of libration rather than breathing motion. The TLS model was applied to elucidate the characteristics of the rigid-body motion. The TLS tensors were determined by the least-squares fit to observed temperature factors. The profile of the magnitude of reproduced temperature factors by the TLS method well fitted to that of observed B_eqv. However, considerable disagreement was observed in the shape and orientation of thermal ellipsoid for atoms with large temperature factors, indicating the large contribution of local motion. The upper estimate of the external motion, 67% (TEL) and 61% (HL) of B_eqv, was deduced from the plot of the magnitude of TLS tensors determined for main chain atoms which were grouped into shells according to the distance from the center of libration. In the external motion, the translational portion is predominant and the contribution of libration and screw motion is relatively small. The internal motion, estimated by subtracting the upper estimate of the external motion from the observed temperature factor, is very similar between TEL and HL in spite of the difference in 54 of 130 amino acid residues and in crystal packing, being suggested to reflect the intrinsic internal motion of chicken-type lysozymes. Proteins 30:232–243, 1998. © 1998 Wiley-Liss, Inc.     

Functional implications of interleukin-1 beta based on the three-dimensional structure.

The molecular structure of interleukin-1 beta, a hormone-like cytokine with roles in several disease processes, has been determined at 2.0 A resolution and refined to a crystallographic R-factor of 0.19. The framework of this molecule consists of 12 antiparallel beta-strands exhibiting pseudo-3-fold symmetry. Six of the strands make up a beta-barrel with polar residues concentrated at either end. Analysis of the three-dimensional structure, together with results from site-directed mutagenesis and biochemical and immunological studies, suggest that the core of the beta-barrel plays an important functional role. A large patch of charged residues on one end of the barrel is proposed as the binding surface with which IL-1 interacts with its receptor.     

The dynamic NMR structure of the TΨC-loop: Implications for the specificity of tRNA methylation

tRNA (m5U54)-methyltransferase (RUMT) catalyzes the S-adenosylmethionine-dependentmethylation of uridine-54 in the TC-loop of all transfer RNAs in E. coli to form the 54-ribosylthymine residue. However, in all tRNA structures, residue 54 is completely buried andthe question arises as to how RUMT gains access to the methylation site. A 17-mer RNAhairpin consisting of nucleotides 49–65 of the T-loop is a substrate for RUMT.Homonuclear NMR methods in conjunction with restrained molecular dynamics (MD)methods were used to determine the solution structure of the 17-mer T-arm fragment. Theloop of the hairpin exhibits enhanced flexibility which renders the conventional NMR averagestructure less useful compared to the more commonly found situation where a molecule existsin predominantly one major conformation. However, when resorting to softer refinementmethods such as MD with time-averaged restraints, the conflicting restraints in the loop canbe satisfied much better. The dynamic structure of the T-arm is represented as an ensembleof 10 time-clusters. In all of these, U54 is completely exposed. The flexibility of the T-loop in solution in conjunction with extensive binding studies of RUMT with the TC-loop and tRNA suggest that the specificity of the RUMT/tRNA recognition is associated withtRNA tertiary structure elements. For the methylation, RUMT would simply have to breakthe tertiary interactions between the D- and T-loops, leading to a melting of the T-armstructure and making U54 available for methylation.     

Protein solution structure calculations in solution: Solvated molecular dynamics refinement of calbindin D9k

The three-dimensional solution structures of proteins determinedwith NMR-derived constraints are almost always calculated in vacuo. Thesolution structure of (Ca²⁺)_{_2}-calbindinD_9k has been redetermined by new restrained molecular dynamics(MD) calculations that include Ca²⁺ ions and explicit solventmolecules. Four parallel sets of MD refinements were run to provide accuratecomparisons of structures produced in vacuo, in vacuo withCa²⁺ ions, and with two different protocols in a solvent bathwith Ca²⁺ ions. The structural ensembles were analyzed interms of structural definition, molecular energies, packing density,solvent-accessible surface, hydrogen bonds, and the coordination of calciumions in the two binding loops. Refinement including Ca²⁺ ionsand explicit solvent results in significant improvements in the precisionand accuracy of the structure, particularly in the binding loops. Theseresults are consistent with results previously obtained in free MDsimulations of proteins in solution and show that the rMD refinedNMR-derived solution structures of proteins, especially metalloproteins, canbe significantly improved by these strategies.     

Direct NOE refinement of biomolecular structures using 2D NMR data

Summary A set of computer programs called DINOSAUR has been developed, which allows the refinement of biomolecular structures directly from 2D NOE intensities. The NOE restraining potential implemented emphasises the weak intensities corresponding to larger distances which are more likely to determine the three-dimensional structure. An approximation based on a two-spin approximation is proposed for the gradient of the NOE intensities instead of the exact solution which is extremely time-consuming. The DINOSAUR routines have been implemented in various refinement programs (Distance bound Driven Dynamics, Molecular Dynamics and Energy Minimisation) and tested on an eight-residue model peptide.