Interdomain dynamics and ligand binding: molecular dynamics simulations of glutamine binding protein

Periplasmic binding proteins from Gram-negative bacteria possess a common architecture, comprised of two domains linked by a hinge region, a fold which they share with the neurotransmitter-binding domains of ionotropic glutamate receptors (GluRs). Glutamine-binding protein (GlnBP) is one such protein, whose crystal structure has been solved in both open and closed forms. Multi-nanosecond molecular dynamics simulations have been used to explore motions about the hinge region and how they are altered by ligand binding. Glutamine binding is seen to significantly reduce inter-domain motions about the hinge region. Essential dynamics analysis of inter-domain motion revealed the presence of both hinge-bending and twisting motions, as has been reported for a related sugar-binding protein. Significantly, the influence of the ligand on GlnBP dynamics is similar to that previously observed in simulations of rat glutamate receptor (GluR2) ligand-binding domain. The essential dynamics analysis of GlnBP also revealed a third class of motion which suggests a mechanism for signal transmission in GluRs.     

A comparison of structural and dynamic properties of different simulation methods applied to SH3.

The dynamic and static properties of molecular dynamics simulations using various methods for treating solvent were compared. The SH3 protein domain was chosen as a test case because of its small size and high surface-to-volume ratio. The simulations were analyzed in structural terms by examining crystal packing, distribution of polar residues, and conservation of secondary structure. In addition, the "essential dynamics" method was applied to compare each of the molecular dynamics trajectories with a full solvent simulation. This method proved to be a powerful tool for the comparison of large concerted atomic motions in SH3. It identified methods of simulation that yielded significantly different dynamic properties compared to the full solvent simulation. Simulating SH3 using the stochastic dynamics algorithm with a vacuum (reduced charge) force field produced properties close to those of the full solvent simulation. The application of a recently described solvation term did not improve the dynamic properties. The large concerted atomic motions in the full solvent simulation as revealed by the essential dynamics method were analyzed for possible biological implications. Two loops, which have been shown to be involved in ligand binding, were seen to move in concert to open and close the ligand-binding site.     

Structure and internal mobility of proteins: a molecular dynamics study of hen egg white lysozyme

The structure and internal motions of the protein hen egg white lysozyme are studied by analysis of simulation and experimental data. A molecular dynamics simulation and an energy minimization of the protein in vacuum have been made and the results compared with high-resolution structures and temperature factors of hen egg white lysozyme in two different crystal forms and of the homologous protein human lysozyme. The structures obtained from molecular dynamics and energy minimization have root-mean-square deviations for backbone atoms of 2.3 Å and 1.1–1.3 Å, respectively, relative to the crystal structures; the different crystal structures have root-mean-square deviations of 0.73–0.81 Å for the backbone atoms. In comparing the backbone dihedral angles, the difference between the dynamics and the crystal structure on which it is based is the same as that between any two crystal structures. The internal fluctuations of atomic positions calculated from the molecular dynamics trajectory agree well with the temperature factors from the three structures. Simulation and crystal results both show that there are large motions for residues involved in exposed turns of the backbone chain, relatively smaller motions for residues involved in the middle of helices or β-sheet structures, and relatively small motions of residues near disulfide bridges. Also, both the simulation and crystal data show that side-chain atoms have larger fluctuations than main-chain atoms. Moreover, the regions that have large deviations among the x-ray crystal structures, which indicates flexibility, are found to have large fluctuations in the simulation.     

Insight derived from molecular dynamics simulation into dynamics and molecular motions of cuticle-degrading serine protease Ver112

Cuticle-degrading serine protease Ver112, which derived from a nematophagous fungus Lecanicillium psalliotae, has been exhibited to have high cuticle-degrading and nematicidal activities. We have performed molecular dynamics (MD) simulation based on the crystal structure of Ver112 to investigate its dynamic properties and large-scale concerted motions. The results indicate that the structural core of Ver112 shows a small fluctuation amplitude, whereas the substrate binding sites, and the regions close to and opposite the substrate binding sites experience significant conformational fluctuations. The large concerted motions obtained from essential dynamics (ED) analysis of MD trajectory can lead to open or close of the substrate binding sites, which are proposed to be linked to the functional properties of Ver112, such as substrate binding, orientation, catalytic, and release. The significant motion in the loop regions that is located opposite the binding sites are considered to play an important role in modulating the dynamics of the substrate binding sites. Furthermore, the bottom of free energy landscape (FEL) of Ver112 are rugged, which is mainly caused by the fluctuations of substrate binding regions and loops located opposite the binding site. In addition, the mechanism underlying the high flexibility and catalytic activity of Ver112 was also discussed. Our simulation study complements the biochemical and structural studies, and provides insight into the dynamics-function relationship of cuticle-degrading serine protease Ver112.     

Investigation of the mechanism of domain closure in citrate synthase by molecular dynamics simulation

Six, 2 ns molecular dynamics simulations have been performed on the homodimeric enzyme citrate synthase. In three, both monomers were started from the open, unliganded X-ray conformation. In the remaining three, both monomers started from a closed, liganded X-ray conformation, with the ligands removed. Projecting the motion from the simulations onto the experimental domain motion revealed that the free-energy profile is rather flat around the open conformation, with steep sides. The most closed conformations correspond to hinge-bending angles of 12-14 compared to the 20 degrees that occurs upon the binding of oxaloacetate. It is also found that the open, unliganded X-ray conformation is situated at the edge of the steep rise in free energy, although conformations that are about 5 degrees more open were sampled. A rigid-body essential dynamics analysis of the combined open trajectories has shown that domain motions in the direction of the closed X-ray conformation are compatible with the natural domain motion of the unliganded protein, which has just two main degrees of freedom. The simulations starting from the closed conformation suggest a free-energy profile with a small barrier in going from the closed to open conformation. A combined essential dynamics and hinge-bending analysis of a trajectory that spontaneously converts from the closed to open state shows an almost exact correspondence to the experimental transition that occurs upon ligand binding. The simulations support the conclusion from an earlier analysis of the experimental transition that the beta-hairpin acts as a mechanical hinge by attaching the small domain to the large domain through a conserved main-chain hydrogen bond and salt-bridges, and allowing rotation to occur via its two flexible termini. The results point to a mechanism of domain closure in citrate synthase that has analogy to the process of closing a door.     

Functional conformational changes of endo-1,4-xylanase II from Trichoderma reesei: A molecular dynamics study

Recent crystallographic studies have revealed a range of structural changes in the three-dimensional structure of endo-1,4-xylanase (XYNII) from Trichoderma reesei. The observed conformational changes can be described as snapshots of an open-close movement of the active site of XYNII. These structures were further analyzed in this study. In addition, a total of four 1 ns molecular dynamics (MD) simulations were performed representing different states of the enzyme. A comparison of the global and local changes found in the X-ray structures and the MD runs suggested that the simulations reproduced a similar kind of active site opening and closing as predicted by the crystal structures. The open-close movement was characterized by the use of distance difference matrixes and the Hingefind program (Wriggers and Schulten, Proteins 29:1–14, 1997) to be a ‘hinge-bending’ motion involving two large rigidly-moving regions and an extended hinge. This conformational feature is probably inherent to this molecular architecture and probably plays a role in the function of XYNII. Proteins 31:434–444, 1998. © 1998 Wiley-Liss, Inc.     

Direct Observation of T4 Lysozyme Hinge-Bending Motion by Fluorescence Correlation Spectroscopy

Bacteriophage T4 Lysozyme (T4L) catalyzes the hydrolysis of the peptidoglycan layer of the bacterial cell wall late in the infection cycle. It has long been postulated that equilibrium dynamics enable substrate access to the active site located at the interface between the N- and C-terminal domains. Crystal structures of WT-T4L and point mutants captured a range of conformations that differ by the hinge-bending angle between the two domains. Evidence of equilibrium between open and closed conformations in solution was gleaned from distance measurements between the two domains but the nature of the equilibrium and the timescale of the underlying motion have not been investigated. Here, we used fluorescence fluctuation spectroscopy to directly detect T4L equilibrium conformational fluctuations in solution. For this purpose, Tetramethylrhodamine probes were introduced at pairs of cysteines in regions of the molecule that undergo relative displacement upon transition from open to closed conformations. Correlation analysis of Tetramethylrhodamine intensity fluctuations reveals hinge-bending motion that changes the relative distance and orientation of the N- and C-terminal domains with ≅15 μs relaxation time. That this motion involves interconversion between open and closed conformations was further confirmed by the dampening of its amplitude upon covalent substrate trapping. In contrast to the prevalent two-state model of T4L equilibrium, molecular brightness and number of particles obtained from cumulant analysis suggest that T4L populates multiple intermediate states, consistent with the wide range of hinge-bending angles trapped in the crystal structure of T4L mutants.     

E. coli 5'-nucleotidase undergoes a hinge-bending domain rotation resembling a ball-and-socket motion

Structures of nine independent conformers of E. coli 5'-nucleotidase (5'-NT) have been analyzed using four different crystal forms. These data show that the two-domain protein undergoes an unusual 96 degrees hinge-bending domain rotation. Structures of the open and closed forms with substrates and inhibitors reveal that the substrate moves by approximately 25 A with the large domain rotation into the catalytic site. The domain motions derived from a comparison of the nine conformations agree well with motions obtained from a normal mode analysis in that all independent domain rotations are around axes that are roughly located in the plane which includes the domain centers and the hinge. Two residues, Lys355 and Gly356, form the core of the hinge region and undergo a large change of the main-chain torsion angles. The hinge-bending movement observed for 5'-nucleotidase differs markedly from a classical hinge-bending closure motion which involves an opening of the substrate or ligand-binding cleft between two domains. In contrast, the movement observed in 5'-nucleotidase resembles that of a ball-and-socket joint. The smaller C-terminal domain rotates approximately around its center such that the residues at the domain interface move in a sliding motion along the interface. Few direct interdomain contacts and a layer of water molecules between the two domains facilitate the sliding motion.     

The consistency of large concerted motions in proteins in molecular dynamics simulations.

A detailed investigation is presented into the effect of limited sampling time and small changes in the force field on molecular dynamics simulations of a protein. Thirteen independent simulations of the B1 IgG-binding domain of streptococcal protein G were performed, with small changes in the simulation parameters in each simulation. Parameters studied included temperature, bond constraints, cut-off radius for electrostatic interactions, and initial placement of hydrogen atoms. The essential dynamics technique was used to reveal dynamic differences between the simulations. Similar essential dynamics properties were found for all simulations, indicating that the large concerted motions found in the simulations are not particularly sensitive to small changes in the force field. A thorough investigation into the stability of the essential dynamics properties as derived from a molecular dynamics simulation of a few hundred picoseconds is provided. Although the definition of the essential modes of motion has not fully converged in these short simulations, the subspace in which these modes are confined is found to be reproducible.     

Molecular dynamics of hinge-bending motion of IgG vanishing with hydrolysis by papain

We have performed dielectric relaxation measurements via a time domain reflectometry (TDR) method to study dynamic behaviors of the segmental flexibility of immunoglobulin G (IgG) in aqueous solution without antigen binding. In general, an intermediate relaxation process due to bound water is observed around 100 MHz at 25 degrees C for common proteins between two relaxation processes due to overall rotation and reorientation of free water. However, the intermediate process observed around 6 MHz for IgG was due to both bound water and hinge-bending motion. The apparent activation energy of 33 kJ/mol was larger than 27 kJ/mol for only bound water, and the relaxation strength was about five times as large as expected for bound water. The shape of the relaxation curve was very broad and asymmetric. These characteristic differences arising from the hinge-bending motion of IgG disappeared for fragments decomposed from IgG hydrolyzed by papain, since the hinge-bending motion did not exist in this case. We have separated the relaxation processes due to hinge-bending motion and bound water for IgG and obtained the Fab-Fab angle of IgG as about 130 degrees by Kirkwood's correlation parameter and the activation energy of 34 kJ/mol for hinge-bending motion.     

Interaction of a <Emphasis Type="Italic">Salmonella enteritidis</Emphasis> O-antigen octasaccharide with the phage P22 tailspike protein by NMR spectroscopy and docking studies

The tailspike protein P22 recognizes an octasaccharide derived from the O-antigen polysaccharide of Salmonella enteritidis in a shallow groove and molecular docking successfully identifies this binding region on the protein surface. Analysis by 2D ¹H,¹H-T-ROESY and transferred NOESY NMR experiments indicate that the bound octasaccharide ligand has a conformation similar to that observed in solution. The results from a saturation transfer difference NMR experiment show that a large number of protons in the octasaccharide are in close contact with the protein as a result of binding. A comparison of the crystal structure of the complex and a molecular dynamics simulation of the octasaccharide with explicit water molecules suggest that only minor conformational changes are needed upon binding to the tailspike protein.     

Investigation of Domain Motions in Bacteriophage T4 Lysozyme

Abstract

Hinge-bending in T4 lysozyme has been inferred from single amino acid mutant crystalline allomorphs by Matthews and coworkers. This raises an important question: are the different conformers in the unit cell artifacts of crystal packing forces, or do they represent different solution state structures? The objective of this theoretical study is to determine whether domain motions and hinge-bending could be simulated in T4 lysozyme using molecular dynamics. An analysis of a 400 ps molecular dynamics simulation of the 164 amino acid enzyme T4 lysozyme is presented. Molecular dynamics calculations were computed using the Discover software package (Biosym Technologies). All hydrogen atoms were modeled explicitly with the inclusion of all 152 crystallographic waters at a temperature of 300 K. The native T4 lysozyme molecular dynamics simulation demonstrated hinge-bending in the protein. Relative domain motions between the N-terminal and C-terminal domains were evident. The enzyme hinge bending sites resulted from small changes in backbone atom conformations over several residues rather than rotation about a single bound. Two hinge loci were found in the simulation. One locus comprises residues 8–14 near the C-terminal of the A helix; the other site, residues 77–83 near the C-terminal of the C helix. Comparison of several snapshot structures from the dynamics trajectory clearly illustrates domain motions between the two lysozyme lobes. Time correlated atomic motions in the protein were analyzed using a dynamical cross-correlation map. We found a high degree of correlated atomic motions in each of the domains and, to a lesser extent, anticorrelated motions between the two domains. We also found that the hairpin loop in the N-terminal lobe (residues 19–24) acted as a mobile ‘flap’ and exhibited highly correlated dynamic motions across the cleft of the active site, especially with residue 142.     

A diverse view of protein dynamics from NMR studies of HIV-1 protease flaps

Internal motion in proteins fulfills a multitude of roles in biological processes. NMR spectroscopy has been applied to elucidate protein dynamics at the atomic level, albeit at a low resolution, and is often complemented by molecular dynamics simulation. However, it is critical to justify the consistency between simulation results and conclusions often drawn from multiple experiments in which uncertainties arising from assumed motional models may not be explicitly evaluated. To understand the role of the flaps of HIV-1 protease dimer in substrate recognition and protease function, many molecular dynamics simulations have been performed. The simulations have resulted in various proposed models of the flap dynamics, some of which are more consistent than others with our working model previously derived from experiments. However, using the working model to discriminate among the simulation results is not straightforward because the working model was derived from a combination of NMR experiments and crystal structure data. In this study, we use the NMR chemical shifts and relaxation data of the protease "monomer" rather than structural data to narrow down the possible conformations of the flaps of the "dimer". For the first time, we show that the tips of the flaps in the unliganded protease dimer interact with each other in solution. Accordingly, we discuss the consistency of the simulations with the model derived from all experimental data.     

Flexibility in a drug transport accessory protein: molecular dynamics simulations of MexA

Drug resistance in gram-negative bacteria may be conferred via efflux through a tripartite complex of an inner membrane pump, an outer membrane pore, and a periplasmic adaptor protein. These are AcrB, TolC, and AcrA, respectively, in Escherichia coli. In Pseudomonas aerugonisa, their homologs are MexB, OprM, and MexA. Defining the interdomain dynamics of the adaptor protein is essential to understanding the mechanism of complex formation. Extended (25 ns) molecular dynamics simulations of MexA have been performed to determine such interdomain dynamics. Analysis of conformational drift demonstrates substantial motions of the three domains of MexA relative to one another. Principal components analysis reveals a hinge-bending motion and rotation of the alpha-helical hairpin relative to the other domains to be the two dominant motions. These two motions provide an element of considerable flexibility which is likely to be exploited in the adaptor function of MexA.     

Folding funnels and conformational transitions via hinge-bending motions

In this article we focus on presenting a broad range of examples illustrating low-energy transitions via hinge-bending motions. The examples are divided according to the type of hinge-bending involved; namely, motions involving fragments of the protein chains, hinge-bending motions involving protein domains, and hinge-bending motions between the covalently unconnected subunits. We further make a distinction between allosterically and nonallosterically regulated proteins. These transitions are discussed within the general framework of folding and binding funnels. We propose that the conformers manifesting such swiveling motions are not the outcome of “induced fit” binding mechanism; instead, molecules exist in an ensemble of conformations that are in equilibrium in solution. These ensembles, which populate the bottoms of the funnels,a priori contain both the “open” and the “closed” conformational isomers. Furthermore, we argue that there are no fundamental differences among the physical principles behind the folding and binding funnels. Hence, there is no basic difference between funnels depicting ensembles of conformers of single molecules with fragment, or domain motions, as compared to subunits in multimeric quaternary structures, also showing such conformational transitions. The difference relates only to the size and complexity of the system. The larger the system, the more complex its corresponding fused funnel(s). In particular, funnels associated with allosterically regulated proteins are expected to be more complicated, because allostery is frequently involved with movements between subunits, and consequently is often observed in multichain and multimolecular complexes. This review centers on the critical role played by flexibility and conformational fluctuations in enzyme activity. Internal motions that extend over different time scales and with different amplitudes are known to be essential for the catalytic cycle. The conformational change observed in enzyme-substrate complexes as compared to the unbound enzyme state, and in particular the hinge-bending motions observed in enzymes with two domains, have a substantial effect on the enzymatic catalytic activity. The examples we review span the lipolytic enzymes that are particularly interesting, owing to their activation at the water-oil interface; an allosterically controlled dehydrogenase (lactate dehydrogenase); a DNA methyltransferase, with a covalently-bound intermediate; large-scale flexible loop motions in a glycolytic enzyme (TIM); domain motion in PGK, an enzyme which is essential in most cells, both for ATP generation in aerobes and for fermentation in anaerobes; adenylate kinase, showing large conformational changes, owing to their need to shield their catalytic centers from water; a calcium-binding protein (calmodulin), involved in a wide range of cellular calcium-dependent signaling; diphtheria toxin, whose large domain motion has been shown to yield “domain swapping” the hexameric glutamate dehydrogenase, which has been studied both in a thermophile and in a mesophile; an allosteric enzyme, showing subunit motion between the R and the T states (aspartate transcarbamoylase), and the historically well-studied lac represoor. Nonallosteric subunit transitions are also addressed with some examples (aspartate receptor andBamHI endonuclease). Hence, using this enzyme-catalysis-centered discussion, we address energy funnel landscapes of large-scale conformational transitions, rather than the faster, quasi-harmonic, thermal fluctuations.     

Dynamics of Fusarium solani cutinase investigated through structural comparison among different crystal forms of its variants

In characterizing mutants and covalently inhibited complexes of Fusarium solani cutinase, which is a 197-residue lipolytic enzyme, 34 variant structures, crystallizing in 8 different crystal forms, have been determined, mostly at high resolution. Taking advantage of this considerable body of information, a structural comparative analysis was carried out to investigate the dynamics of cutinase. Surface loops were identified as the major flexible protein regions, particularly those forming the active-site groove, whereas the elements constituting the protein scaffold were found to retain the same conformation in all the cutinase variants studied. Flexibility turned out to be correlated with thermal motion. With a given crystal packing environment, a high flexibility turned out to be correlated with a low involvement in crystal packing contacts. The high degree of crystal polymorphism, which allowed different conformations with similar energy to be detected, made it possible to identify motions which would have remained unidentified if only a single crystal form had been available. Fairly good agreement was found to exist between the data obtained from the structural comparison and those from a molecular dynamics (MD) simulation carried out on the native enzyme. The crystallographic approach used in this study turned out to be a suitable tool for investigating cutinase dynamics. Because of the availability of a set of closely related proteins in different crystal environments, the intrinsic drawback of a crystallographic approach was bypassed. By combining several static pictures, the dynamics of the protein could be monitored much more realistically than what can be achieved on the basis of static pictures alone. Proteins 26:442–458 © 1996 Wiley-Liss, Inc.     

On the similarity of properties in solution or in the crystalline state: A molecular dynamics study of hen lysozyme

As protein crystals generally possess a high water content, it is assumed that the behaviour of a protein in solution and in crystal environment is very similar. This assumption can be investigated by molecular dynamics (MD) simulation of proteins in the different environments. Two 2ns simulations of hen egg white lysozyme (HEWL) in crystal and solution environment are compared to one another and to experimental data derived from both X-ray and NMR experiments, such as crystallographic B-factors, NOE atom–atom distance bounds, ³J_{H N}-coupling constants, and ¹H-¹⁵N bond vector order parameters. Both MD simulations give very similar results. The crystal simulation reproduces X-ray and NMR data slightly better than the solution simulation.     

The Protein-protein Interactions Between SMPI and Thermolysin Studied by Molecular Dynamics and MM/PBSA Calculations

Abstract

Thermolysin is a zinc-metalloendopeptidase secreted by the gram-positive thermophilic bacterium Bacillus thermoproteolyticus. Thermolysin belongs to the gluzinicin family of enzymes, which is selectively inhibited by Steptomyces metalloproteinase inhibitor (SMPI). Very little is known about the interaction between SMPI and thermolysin. Knowledge about the protein-protein interactions is very important for designing new thermolysin inhibitors with possible industrial or pharmaceutical applications. In the present study, two binding modes between SMPI and thermolysin were studied by 2300 picoseconds (ps) of comparative molecular dynamics (MD) simulations and calculation of the free energy of binding using the molecular mechanics-Poisson-Boltmann surface area (MM/PBSA) method. One of the positions, the ‘horizontal arrow head docking’ (HAHD) was similar to the previously proposed binding mode by Tate et al. (Tate, S., Ohno, A., Seeram, S. S., Hiraga, K., Oda, K., and Kainosho, M. J. Mol. Biol. 282, 435–446 (1998)). The other position, the ‘vertical arrow head docking’ (VAHD) was obtained by a manual docking guided by the shape and charge distribution of SMPI and the binding pocket of thermolysin. The calculations showed that SMPI had stronger interactions with thermolysin in the VAHD than in the HAHD complex, and the VAHD complex was considered more realistic than the HAHD complex. SMPI interacted with thermolysin not only at the active site but had auxiliary binding sites contributing to proper interactions. The VAHD complex can be used for designing small molecule inhibitors mimicking the SMPI-thermolysin binding interfaces.