Exploration of disorder in protein structures by X-ray restrained molecular dynamics

Conformational disorder in crystal structures of ribonuclease-A and crambin is studied by including two independent structures in least-squares optimizations against X-ray data. The optimizations are carried out by X-ray restrained molecular dynamics (simulated annealing refinement) and by conventional least-squares optimization. Starting from two identical structures, the optimizations against X-ray data lead to significant deviations between the two, with rms backbone displacements of 0.45 A for refinement of ribonuclease at 1.53 A resolution, and 0.31 A for crambin at 0.945 A. More than 15 independent X-ray restrained molecular dynamics runs have been carried out for ribonuclease, and the displacements between the resulting structures are highly reproducible for most atoms. These include residues with two or more conformations with significant dihedral angle differences and alternative hydrogen bonding, as well as groups of residues that undergo displacements that are suggestive of rigid-body librations. The crystallographic R-values obtained are approximately 13%, as compared to 15.3% for a comparable refinement with a single structure. Least-squares optimization without an intervening restrained molecular dynamics stage is sufficient to reproduce most of the observed displacements. Similar results are obtained for crambin, where the higher resolution of the X-ray data allows for refinement of unconstrained individual anisotropic temperature factors. These are shown to be correlated with the displacements in the two-structure refinements.     

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.     

Automatic identification of discrete substates in proteins: Singular value decomposition analysis of time-averaged crystallographic refinements

The singular value decomposition (SVD) provides a method for decomposing a molecular dynamics trajectory into fundamental modes of atomic motion. The right singular vectors are projections of the protein conformations onto these modes showing the protein motion in a generalized low-dimensional basis. Statistical analysis of the right singular vectors can be used to classify discrete configurational substates in the protein. The configuration space portraits formed from the right singular vectors can also be used to visualize complex high-dimensional motion and to examine the extent of configuration space sampling by the simulation. © 1995 Wiley-Liss, Inc.     

The crystal structure of free human hypoxanthine-guanine phosphoribosyltransferase reveals extensive conformational plasticity throughout the catalytic cycle

Human hypoxanthine-guanine phosphoribosyltransferase (HGPRT) catalyses the synthesis of the purine nucleoside monophosphates, IMP and GMP, by the addition of a 6-oxopurine base, either hypoxanthine or guanine, to the 1-beta-position of 5-phospho-alpha-d-ribosyl-1-pyrophosphate (PRib-PP). The mechanism is sequential, with PRib-PP binding to the free enzyme prior to the base. After the covalent reaction, pyrophosphate is released followed by the nucleoside monophosphate. A number of snapshots of the structure of this enzyme along the reaction pathway have been captured. These include the structure in the presence of the inactive purine base analogue, 7-hydroxy [4,3-d] pyrazolo pyrimidine (HPP) and PRib-PP.Mg2+, and in complex with IMP or GMP. The third structure is that of the immucillinHP.Mg(2+).PP(i) complex, a transition-state analogue. Here, the first crystal structure of free human HGPRT is reported to 1.9A resolution, showing that significant conformational changes have to occur for the substrate(s) to bind and for catalysis to proceed. Included in these changes are relative movement of subunits within the tetramer, rotation and extension of an active-site alpha-helix (D137-D153), reorientation of key active-site residues K68, D137 and K165, and the rearrangement of three active-site loops (100-128, 165-173 and 186-196). Toxoplasma gondii HGXPRT is the only other 6-oxopurine phosphoribosyltransferase structure solved in the absence of ligands. Comparison of this structure with human HGPRT reveals significant differences in the two active sites, including the structure of the flexible loop containing K68 (human) or K79 (T.gondii).     

Direct structure refinement of high molecular weight proteins against residual dipolar couplings and carbonyl chemical shift changes upon alignment: an application to maltose binding protein

The global fold of maltose binding protein in complex with -cyclodextrin has been determined using a CNS-based torsion angle molecular dynamics protocol involving direct refinement against dipolar couplings and carbonyl chemical shift changes that occur upon alignment. The shift changes have been included as structural restraints using a new module, CANI, that has been incorporated into CNS. Force constants and timesteps have been determined that are particularly effective in structure refinement applications involving high molecular weight proteins with small to moderate numbers of NOE restraints. Solution structures of the N- and C-domains of MBP calculated with this new protocol are within 2 Å of the X-ray conformation.     

Molecular dynamics simulation of metal free structure of Lmb,a laminin-binding adhesin of Streptococcus agalactiae: metal removal and its structural implications

Metal-binding receptors are one of the extracellular components of ATP-binding cassette transporters that are essential for regulation of metal homeostasis in bacteria. Laminin-binding adhesin (Lmb) of Streptococcus agalactiae falls under this class of solute binding proteins. It binds to zinc with a high affinity. Crystal structure of Lmb solved previously by our group reveals that the zinc is tetrahedrally coordinated by three histidines and a glutamate at the interdomain cleft. Lmb contains a long disordered loop close to the metal-binding site whose precise function is unknown. Several experimental attempts to produce apo-Lmb failed and this prompted us to carry out in silico studies to analyse the structural importance of the metal in Lmb. Here, we present the results of the molecular dynamics (MD) simulation studies of native, apo-(metal removed) and the long loop truncated Lmb models along with a homologous protein, TroA from Treponema pallidum that was taken up for validating the MD results of Lmb. Absence of a metal results in significant structural changes in Lmb, particularly at the metal-binding pocket and with the long loop, although the overall fold is retained. This study thus revealed that the Lmb can exist in different conformational states with subtle differences in the overall fold based on the presence or absence of the metal. This could be functionally important for a putative metal uptake and release and also for the adhesive function of Lmb in recognizing laminin, which contains a high number of zinc finger motifs.     

Backbone conformational flexibility of the lipid modified membrane anchor of the human N-Ras protein investigated by solid-state NMR and molecular dynamics simulation

The lipid modified human N-Ras protein, implicated in human cancer development, is of particular interest due to its membrane anchor that determines the activity and subcellular location of the protein. Previous solid-state NMR investigations indicated that this membrane anchor is highly dynamic, which may be indicative of backbone conformational flexibility. This article aims to address if a dynamic exchange between three structural models exist that had been determined previously. We applied a combination of solid-state nuclear magnetic resonance (NMR) methods and replica exchange molecular dynamics (MD) simulations using a Ras peptide that represents the terminal seven amino acids of the human N-Ras protein. Analysis of correlations between the conformations of individual amino acids revealed that Cys 181 and Met 182 undergo collective conformational exchange. Two major structures constituting about 60% of all conformations could be identified. The two conformations found in the simulation are in rapid exchange, which gives rise to low backbone order parameters and nuclear spin relaxation as measured by experimental NMR methods. These parameters were also determined from two 300 ns conventional MD simulations, providing very good agreement with the experimental data.