Modeling of solvent flow effects in enzyme catalysis under physiological conditions

A stochastic model for the dynamics of enzymatic catalysis in explicit, effective solvents under physiological conditions is presented. Analytically-computed first passage time densities of a diffusing particle in a spherical shell with absorbing boundaries are combined with densities obtained from explicit simulation to obtain the overall probability density for the total reaction cycle time of the enzymatic system. The method is used to investigate the catalytic transfer of a phosphoryl group in a phosphoglycerate kinase-ADP-bis phosphoglycerate system, one of the steps of glycolysis. The direct simulation of the enzyme-substrate binding and reaction is carried out using an elastic network model for the protein, and the solvent motions are described by multiparticle collision dynamics which incorporates hydrodynamic flow effects. Systems where solvent-enzyme coupling occurs through explicit intermolecular interactions, as well as systems where this coupling is taken into account by including the protein and substrate in the multiparticle collision step, are investigated and compared with simulations where hydrodynamic coupling is absent. It is demonstrated that the flow of solvent particles around the enzyme facilitates the large-scale hinge motion of the enzyme with bound substrates, and has a significant impact on the shape of the probability densities and average time scales of substrate binding for substrates near the enzyme, the closure of the enzyme after binding, and the overall time of completion of the cycle.     

Conformation and diffusion behavior of ring polymers in solution: a comparison between molecular dynamics, multiparticle collision dynamics, and lattice Boltzmann simulations

We have studied the effect of chain topology on the structural properties and diffusion of polymers in a dilute solution in a good solvent. Specifically, we have used three different simulation techniques to compare the chain size and diffusion coefficient of linear and ring polymers in solution. The polymer chain is modeled using a bead-spring representation. The solvent is modeled using three different techniques: molecular dynamics (MD) simulations with a particulate solvent in which hydrodynamic interactions are accounted through the intermolecular interactions, multiparticle collision dynamics (MPCD) with a point particle solvent which has stochastic interactions with the polymer, and the lattice Boltzmann method in which the polymer chains are coupled to the lattice fluid through friction. Our results show that the three methods give quantitatively similar results for the effect of chain topology on the conformation and diffusion behavior of the polymer chain in a good solvent. The ratio of diffusivities of ring and linear polymers is observed to be close to that predicted by perturbation calculations based on the Kirkwood hydrodynamic theory.     

Role of protein dynamics in reaction rate enhancement by enzymes

An integrated view of protein structure, dynamics, and function is emerging, where proteins are considered as dynamically active assemblies and internal motions are closely linked to function such as enzyme catalysis. Further, the motion of solvent bound to external regions of protein impacts internal motions and, therefore, protein function. Recently, we discovered a network of protein vibrations in enzyme cyclophilin A, coupled to its catalytic activity of peptidyl-prolyl cis-trans isomerization. Detailed studies suggest that this network, extending from surface regions to active site, is a conserved part of enzyme structure and has a role in promoting catalysis. In this report, theoretical investigations of concerted conformational fluctuations occurring on microsecond and longer time scales within the discovered network are presented. Using a new technique, kinetic energy was added to protein vibrational modes corresponding to conformational fluctuations in the network. The results reveal that protein dynamics promotes catalysis by altering transition state barrier crossing behavior of reaction trajectories. An increase in transmission coefficient and number of productive trajectories with increasing amounts of kinetic energy in vibrational modes is observed. Variations in active site enzyme-substrate interactions near transition state are found to be correlated with barrier recrossings. Simulations also showed that energy transferred from first solvation shell to surface residues impacts catalysis through network fluctuations. The detailed characterization of network presented here indicates that protein dynamics plays a role in rate enhancement by enzymes. Therefore, coupled networks in enzymes have wide implications in understanding allostericity and cooperative effects, as well as protein engineering and drug design.     

Millisecond dynamics in glutaredoxin during catalytic turnover is dependent on substrate binding and absent in the resting states

Conformational dynamics is important for enzyme function. Which motions of enzymes determine catalytic efficiency and whether the same motions are important for all enzymes, however, are not well understood. Here we address conformational dynamics in glutaredoxin during catalytic turnover with a combination of NMR magnetization transfer, R(2) relaxation dispersion, and ligand titration experiments. Glutaredoxins catalyze a glutathione exchange reaction, forming a stable glutathinoylated enzyme intermediate. The equilibrium between the reduced state and the glutathionylated state was biochemically tuned to exchange on the millisecond time scale. The conformational changes of the protein backbone during catalysis were followed by (15)N nuclear spin relaxation dispersion experiments. A conformational transition that is well described by a two-state process with an exchange rate corresponding to the glutathione exchange rate was observed for 23 residues. Binding of reduced glutathione resulted in competitive inhibition of the reduced enzyme having kinetics similar to that of the reaction. This observation couples the motions observed during catalysis directly to substrate binding. Backbone motions on the time scale of catalytic turnover were not observed for the enzyme in the resting states, implying that alternative conformers do not accumulate to significant concentrations. These results infer that the turnover rate in glutaredoxin is governed by formation of a productive enzyme-substrate encounter complex, and that catalysis proceeds by an induced fit mechanism rather than by conformer selection driven by intrinsic conformational dynamics.     

Mesoscopic description of solvent effects on polymer dynamics

Solvent effects on polymer dynamics and structure are investigated using a mesoscopic solvent model that accounts for hydrodynamic interactions among the polymer beads. The simulation method combines molecular dynamics of the polymer chain, interacting with the solvent molecules through intermolecular forces, with mesoscopic multiparticle collision dynamics for the solvent molecules. Changes in the intermolecular forces between the polymer beads and mesoscopic solvent molecules are used to vary the solvent conditions from those for good to poor solvents. Polymer collapse and expansion dynamics following changes in solvent conditions are studied for homopolymer and block copolymer solutions. The frictional properties of polymers are also investigated.     

Two-particle friction in a mesoscopic solvent

The effects of hydrodynamic interactions on the friction tensors for two particles in solution are studied. The particles have linear dimensions on nanometer scales and are either simple spherical particles interacting with the solvent through repulsive Lennard-Jones forces or are composite cluster particles whose atomic components interact with the solvent through repulsive Lennard-Jones forces. The solvent dynamics is modeled at a mesoscopic level through multiparticle collisions that conserve mass, momentum, and energy. The dependence of the two-particle relative friction tensors on the interparticle separation indicates the importance of hydrodynamic interactions for these nanoparticles.     

Diffusion behavior of polyelectrolytes in dilute solution: coupling effects of hydrodynamic and Coulomb interactions

The diffusion behavior of polyelectrolytes in dilute salt-free solution is studied through a hybrid mesoscale simulation technique that combines the molecular dynamics method and the multiparticle collision dynamics approach. To elucidate the effects of hydrodynamic interactions (HI), we compare results for hydrodynamic and random solvents. When HI are taken into account, we find that the chain diffusivity decreases initially and then increases gradually with the increasing strength of the Coulomb interaction. By contrast, when HI are switched off, the electrostatic-dependent diffusivity shows three distinct regions, and a plateau of approximately constant diffusivity manifests between two decreasing regions. The findings reveal that the dynamics of polyelectrolytes in dilute solution depend on the coupling effects of hydrodynamic and Coulomb interactions, and that these dynamics can be understood by considering the conformational changes of chains, the counterion condensation, and the dynamics of counterions.     

Hydration structure of human lysozyme investigated by molecular dynamics simulation and cryogenic X-ray crystal structure analyses: on the correlation between crystal water sites,solvent density,and solvent dipole

The hydration structure of human lysozyme was studied with cryogenic X-ray diffraction experiment and molecular dynamics simulations. The crystal structure analysis at a resolution of 1.4 A provided 405 crystal water molecules around the enzyme. In the simulations at 300 K, the crystal structure was immersed in explicit water molecules. We examined correlations between crystal water sites and two physical quantities calculated from the 1-ns simulation trajectories: the solvent density reflecting the time-averaged distribution of water molecules, and the solvent dipole measuring the orientational ordering of water molecules around the enzyme. The local high solvent density sites were consistent with the crystal water sites, and better correlation was observed around surface residues with smaller conformational fluctuations during the simulations. Solvent dipoles around those sites exhibited coherent and persistent ordering, indicating that the hydration water molecules at the crystal water sites were highly oriented through the interactions with hydrophilic residues. Those water molecules restrained the orientational motions of adjoining water molecules and induced a solvent dipole field, which was persistent during the simulations around the enzyme. The coherent ordering was particularly prominent in and around the active site cleft of the enzyme. Because the ordering was significant up to the third to fourth solvent layer region from the enzyme surface, the coherently ordered solvent dipoles likely contributed to the molecular recognition of the enzyme in a long-distance range. The present work may provide a new approach combining computational and the experimental studies to understand protein hydration.     

Mesoscopic model for diffusion-influenced reaction dynamics

A hybrid mesoscopic multiparticle collision model is used to study diffusion-influenced reaction kinetics. The mesoscopic particle dynamics conserves mass, momentum, and energy so that hydrodynamic effects are fully taken into account. Reactive and nonreactive interactions with catalytic solute particles are described by full molecular dynamics. Results are presented for large-scale, three-dimensional simulations to study the influence of diffusion on the rate constants of the A + C <==> B + C reaction. In the limit of a dilute solution of catalytic C particles, the simulation results are compared with diffusion equation approaches for both the irreversible and reversible reaction cases. Simulation results for systems where the volume fraction phi of catalytic spheres is high are also presented, and collective interactions among reactions on catalytic spheres that introduce volume fraction dependence in the rate constants are studied.     

Rotational and translational relaxation times of quasi-elastic and rigid dumbbells elastically bound to polymer network junctions

The dynamics of a rigid rod located between fixed junctions of a polymer network is studied. Three approaches are used in the solution of this problem. The first is based on the viscoelastic model, where a rigid rod is simulated by an elastic dumbbell with a fixed average length; the second includes solution of equations of motion for projections of the rigid rod using the Lagrangian multipliers under the constraint condition; and the third involves solution of the diffusion equation in the presence of an elastic potential. The second and third approaches allow calculation of orientational relaxation times for rod projections under the action of a strong orienting field. The dependences of the relaxation times of orientational and translational motions of the rod projections on the coordinate axes and the orientational relaxation times of mean-square rod projections on the model parameters (the distances between fixed polymer network junctions, the length of the rigid rod, and the elastic coefficient characterizing the binding between the rod and the network) are found.     

Coarse-grain model for lipid bilayer self-assembly and dynamics: Multiparticle collision description of the solvent

A mesoscopic coarse-grain model for computationally efficient simulations of biomembranes is presented. It combines molecular dynamics simulations for the lipids, modeled as elastic chains of beads, with multiparticle collision dynamics for the solvent. Self-assembly of a membrane from a uniform mixture of lipids is observed. Simulations at different temperatures demonstrate that it reproduces the gel and liquid phases of lipid bilayers. Investigations of lipid diffusion in different phases reveals a crossover from subdiffusion to normal diffusion at long times. Macroscopic membrane properties, such as stretching and bending elastic moduli, are determined directly from the mesoscopic simulations. Velocity correlation functions for membrane flows are determined and analyzed.     

Elastic network normal modes provide a basis for protein structure refinement

It is well recognized that thermal motions of atoms in the protein native state, the fluctuations about the minimum of the global free energy, are well reproduced by the simple elastic network models (ENMs) such as the anisotropic network model (ANM). Elastic network models represent protein dynamics as vibrations of a network of nodes (usually represented by positions of the heavy atoms or by the C(α) atoms only for coarse-grained representations) in which the spatially close nodes are connected by harmonic springs. These models provide a reliable representation of the fluctuational dynamics of proteins and RNA, and explain various conformational changes in protein structures including those important for ligand binding. In the present paper, we study the problem of protein structure refinement by analyzing thermal motions of proteins in non-native states. We represent the conformational space close to the native state by a set of decoys generated by the I-TASSER protein structure prediction server utilizing template-free modeling. The protein substates are selected by hierarchical structure clustering. The main finding is that thermal motions for some substates, overlap significantly with the deformations necessary to reach the native state. Additionally, more mobile residues yield higher overlaps with the required deformations than do the less mobile ones. These findings suggest that structural refinement of poorly resolved protein models can be significantly enhanced by reduction of the conformational space to the motions imposed by the dominant normal modes.     

Modulation of Kinetic Pathways of Enzyme–Substrate Interaction in a Microfluidic Channel: Nanoscopic Water Dynamics as a Switch

Enzyme-mediated catalysis is attributed to enzyme–substrate interactions, with models such as “induced fit” and “conformational selection” emphasizing the role of protein conformational transitions. The dynamic nature of the protein structure, thus, plays a crucial role in molecular recognition and substrate binding. As large-scale protein motions are coupled to water motions, hydration dynamics play a key role in protein dynamics, and hence, in enzyme catalysis. Here, microfluidic techniques and time-dependent fluorescence Stokes shift (TDFSS) measurements are employed to elucidate the role of nanoscopic water dynamics in the interaction of an enzyme, α-Chymotrypsin (CHT), with a substrate, Ala-Ala-Phe-7-amido-4-methylcoumarin (AMC) in the cationic reverse micelles of benzylhexadecyldimethylammonium chloride (BHDC/benzene) and anionic reverse micelles of sodium bis(2-ethylhexyl)sulfosuccinate (AOT/benzene). The kinetic pathways unraveled from the microfluidic setup are consistent with the “conformational selection” fit for the interaction of CHT with AMC in the cationic reverse micelles, whereas an “induced fit” mechanism is indicated for the anionic reverse micelles. In the cationic reverse micelles of BHDC, faster hydration dynamics (≈550 ps) aid the pathway of “conformational selection”, whereas in the anionic reverse micelles of AOT, the significantly slower dynamics of hydration (≈1600 ps) facilitate an “induced fit” mechanism for the formation of the final enzyme–substrate complex. The role of water dynamics in dictating the mechanism of enzyme–substrate interaction becomes further manifest in the neutral reverse micelles of Brij-30 and Triton X-100. In the former, the faster water dynamics aid the “conformational selection” pathway, whereas the significantly slower dynamics of water molecules in the latter are conducive to the “induced fit” mechanism in the enzyme–substrate interaction. Thus, nanoscopic water dynamics act as a switch in modulating the pathway of recognition of an enzyme (CHT) by the substrate (AMC) in reverse micelles.     

Characterization of the transition state of functional enzyme dynamics

Through characterization of the solvent isotope effect on protein dynamics, we have examined determinants of the rate limitation to enzyme catalysis. A global conformational change in Ribonuclease A limits the overall rate of catalytic turnover. Here we show that this motion is sensitive to solvent deuterium content; the isotope effect is 2.2, a value equivalent to the isotope effect on the catalytic rate constant. We further demonstrate that the protein motion possesses a linear proton inventory plot, indicating that a single proton is transferred in the transition state. These results provide compelling evidence for close coupling between enzyme dynamics and function and demonstrate that characterization of the transition state for protein motion in atomic detail is experimentally accessible.     

Design of chemically propelled nanodimer motors

The self-propelled motion of nanodimers fueled by a chemical reaction taking place under nonequilibrium steady state conditions is investigated. The nanodimer consists of a pair of catalytic and chemically inactive spheres, in general with different sizes, with a fixed internuclear separation. The solvent in which the dimer moves is treated at a particle-based mesoscopic level using multiparticle collision dynamics. The directed motion of the dimer can be controlled by adjusting the interaction potentials between the solvent molecules and the dimer spheres, the internuclear separation, and sphere sizes. Dimers can be designed so that the directed motion along the internuclear axis occurs in either direction and is much larger than the thermal velocity fluctuations, a condition needed for such nanodimers to perform tasks involving targeted dynamics.     

Functional concerted motions in the bovine serum retinol-binding protein

The large concerted motions in the apo/holo bovine serum retinol- binding protein were studied using molecular dynamics simulation and 'essential dynamics' analysis. Initially, concerted motions were calculated from conformational differences between various crystal structures. The dynamic behaviour of the protein in the configurational space directions, described by these concerted motions, is analysed. This reveals that the large backbone dynamics of the protein is not influenced by the presence of retinol. Study of free retinol dynamics and retinol in the retinol binding site reveals that the protein binds retinol in a favourable conformation, as opposed to what has been previously described for the bovine cellular retinol-binding protein.