Brownian and essential dynamics studies of the HIV-1 integrase catalytic domain

The three-dimensional structure of the active site region of the enzyme HIV-1 integrase is not unambiguously known. This region includes a flexible peptide loop that cannot be well resolved in crystallographic determinations. Here we present two different computational approaches with different levels of resolution and on different time-scales to understand this flexibility and to analyze the dynamics of this part of the protein. We have used molecular dynamics simulations with an atomic model to simulate the region in a realistic and reliable way for 1 ns. It is found that parts of the loop wind up after 300 ps to extend an existing helix. This indicates that the helix is longer than in the earlier crystal structures that were used as basis for this study. Very recent crystal data confirms this finding, underlining the predictive value of accurate MD simulations. Essential dynamics analysis of the MD trajectory yields an anharmonic motion of this loop. We have supplemented the MD data with a much lower resolution Brownian dynamics simulation of 600 ns length. It provides ideas about the slow-motion dynamics of the loop. It is found that the loop explores a conformational space much larger than in the MD trajectory, leading to a "gating"-like motion with respect to the active site.     

Conformational ensemble of an intrinsically flexible loop in mitochondrial import protein Tim21 studied by modeling and molecular dynamics simulations

BackgroundTim21, a subunit of a highly dynamic translocase of the inner mitochondrial membrane (TIM23) complex, translocates proteins by interacting with subunits in the translocase of the outer membrane (TOM) complex and Tim23 channel in the TIM23 complex. A loop segment in Tim21, which is in close proximity of the binding site of Tim23, has different conformations in X-ray, NMR and new crystal contact-free space (CCFS) structures. MD simulations can provide information on the structure and dynamics of the loop in solution.MethodsThe conformational ensemble of the loop was characterized using loop modeling and molecular dynamics (MD) simulations.ResultsMD simulations confirmed mobility of the loop. Multidimensional scaling and clustering were used to characterize the dynamic conformational ensemble of the loop. Free energy landscape showed that the CCFS crystal structure occupied a low energy region as compared to the conventional X-ray crystal structure. Analysis of crystal packing indicates that the CCFS provides larger conformational space for the motions of the loop.ConclusionsOur work reported the conformational ensemble of the loop in solution, which is in agreement with the structure obtained from CCFS approach. The combination of the experimental techniques and computational methods is beneficial for studying highly flexible regions of proteins.General significanceComputational methods, such as loop modeling and MD simulations, have proved to be useful for studying conformational flexibility of proteins. These methods in integration with experimental techniques such as CCFS has the potential to transform the studies on flexible regions of proteins.     

Alpha7 Helix Plays an Important Role in the Conformational Stability of PTP1B

Abstract

The C-terminus of Protein Tyrosine Phosphatase 1B (PTP1B) includes an α-helix (α7), which forms an allosteric binding site 20 Å away from the active site. This helix is specific to PTP1B and its truncation decreases the catalytic activity significantly. Here, molecular dynamics (MD) simulations in the presence and absence of α7 were performed to investigate the role played by α7. The highly mobile α7 was found to maintain its contacts with loop 11 (L11)- α3 helix throughout the simulations. The interactions of Tyr152 on L11, Tyr176, Thr177 on the catalytically important WPD loop and Ser190 on α3 are important for the conformational stability and the concerted motions of the regions surrounding the WPD loop. In the absence of α7, L11 and WPD loop move away from their crystal structure conformations, resulting in the loss of the interactions in this region, and a decrease in the residue displacement correlations in the vicinity of WPD loop. Therefore, we suggest that one of the functionally important roles of α7 may be to limit the L11 and α3 motions, and, facilitate the WPD loop motions. Truncation of α7 in PTP1B is found to affect distant regions as well, such as the substrate recognition site and the phosphate binding-loop (P-loop), changing the conformations of these regions significantly. Our results show that the PTP1B specific α7 is important for the conformation and dynamics of the WPD loop, and also may play a role in ligand binding.     

Solution NMR and computer simulation studies of active site loop motion in triosephosphate isomerase

Solution NMR spin relaxation experiments and classical MD simulations are used to study the dynamics of triosephosphate isomerase (TIM) in complex with glycerol 3-phosphate (G3P). Three regions in TIM exhibit conformational transitions on the micros-ms time scale as detected by chemical exchange broadening effects in NMR spectroscopy: residue Lys 84 on helix C, located at the dimeric interface; active site loop 6; and helix G. The results indicate that the conformational exchange process affecting the residues of loop 6 is the correlated opening and closing of the loop. Distinct processes are responsible for the chemical exchange linebroadening observed in the other regions of TIM. MD simulations confirm that motions of individual residues within the active site loop are correlated and suggest that the chemical exchange processes observed for residues in helix G arise from transitions between 3(10)- and alpha-helical structures. The results of the joint NMR and MD study provide global insight into the role of conformational dynamic processes in the function of TIM.     

Influence of the acetylcholinesterase active site protonation on omega loop and active site dynamics

Existence of alternative entrances in acetylcholinesterase (AChE) could explain the contrast between the very high AChE catalytic efficiency and the narrow and long access path to the active site revealed by X-ray crystallography. Alternative entrances could facilitate diffusion of the reaction products or at least water and ions from the active site. Previous molecular dynamics simulations identified side door and back door as the most probable alternative entrances. The simulations of non-inhibited AChE suggested that the back door opening events occur only rarely (0.8% of the time in the 10ns trajectory). Here we present a molecular dynamics simulation of non-inhibited AChE, where the back door opening appears much more often (14% of the time in the 12ns trajectory) and where the side door opening was observed quite frequently (78% of trajectory time). We also present molecular dynamics, where the back door does not open at all, or where large conformational changes of the AChE omega loop occur together with alternative passage opening events. All these differences in AChE dynamical behavior are caused by different protonation states of two glutamate residues located on bottom of the active site gorge (Glu202 and G450 in Mus musculus AChE). Our results confirm the results of previous molecular dynamics simulations, expand the view and suggest the probable reasons for the overall conformational behavior of AChE omega loop.     

Alpha7 helix plays an important role in the conformational stability of PTP1B

The C-terminus of Protein Tyrosine Phosphatase 1B (PTP1B) includes an α-helix α7), which forms an allosteric binding site 20 ? away from the active site. This helix is specific to PTP1B and its truncation decreases the catalytic activity significantly. Here, molecular dynamics (MD) simulations in the presence and absence of α7 were performed to investigate the role played by α7. The highly mobile α7 was found to maintain its contacts with loop 11 (L11)α3 helix throughout the simulations. The interactions of Tyr152 on L11, Tyr176, Thr177 on the catalytically important WPD loop and Ser190 on α3 are important for the conformational stability and the concerted motions of the regions surrounding the WPD loop. In the absence of α7, L11 and WPD loop move away from their crystal structure conformations, resulting in the loss of the interactions in this region, and a decrease in the residue displacement correlations in the vicinity of WPD loop. Therefore, we suggest that one of the functionally important roles of α7 may be to limit the L11 and α3 motions, and, facilitate the WPD loop motions. Truncation of α7 in PTP1B is found to affect distant regions as well, such as the substrate recognition site and the phosphate binding-loop (P-loop), changing the conformations of these regions significantly. Our results show that the PTP1B specific α7 is important for the conformation and dynamics of the WPD loop, and also may play a role in ligand binding.     

Functional and structural insights revealed by molecular dynamics simulations of an essential RNA editing ligase in Trypanosoma brucei

RNA editing ligase 1 (TbREL1) is required for the survival of both the insect and bloodstream forms of Trypanosoma brucei, the parasite responsible for the devastating tropical disease African sleeping sickness. The type of RNA editing that TbREL1 is involved in is unique to the trypanosomes, and no close human homolog is known to exist. In addition, the high-resolution crystal structure revealed several unique features of the active site, making this enzyme a promising target for structure-based drug design. In this work, two 20 ns atomistic molecular dynamics (MD) simulations are employed to investigate the dynamics of TbREL1, both with and without the ATP substrate present. The flexibility of the active site, dynamics of conserved residues and crystallized water molecules, and the interactions between TbREL1 and the ATP substrate are investigated and discussed in the context of TbREL1's function. Differences in local and global motion upon ATP binding suggest that two peripheral loops, unique to the trypanosomes, may be involved in interdomain signaling events. Notably, a significant structural rearrangement of the enzyme's active site occurs during the apo simulations, opening an additional cavity adjacent to the ATP binding site that could be exploited in the development of effective inhibitors directed against this protozoan parasite. Finally, ensemble averaged electrostatics calculations over the MD simulations reveal a novel putative RNA binding site, a discovery that has previously eluded scientists. Ultimately, we use the insights gained through the MD simulations to make several predictions and recommendations, which we anticipate will help direct future experimental studies and structure-based drug discovery efforts against this vital enzyme.     

On the Permeation by Dioxygen of Urate Oxidase from Aspergillus flavus in Complex with Xanthine Anion: Dioxygen Pathways and a Portrait of the Enzyme Cavities from Molecular Dynamics Simulations in Water Solution

This work describes molecular dynamics (MD) simulations in aqueous media for the complex of the homotetrameric urate oxidase (UOX) from Aspergillus flavus with xanthine anion ( 5 ) in the presence of dioxygen (O₂). After 196.6 ns of trajectory from unrestrained MD, a O₂ molecule was observed leaving the bulk solvent to penetrate the enzyme between two subunits, A/C. From here, the same O₂ molecule was observed migrating, across subunit C, to the hydrophobic cavity that shares residue V227 with the active site. The latter was finally attained, after 378.3 ns of trajectory, with O₂ at a bonding distance from 5 . The reverse same O₂ pathway, from 5 to the bulk solvent, was observed as preferred pathway under random acceleration MD (RAMD), where an external, randomly oriented force was acting on O₂. Both MD and RAMD simulations revealed several cavities populated by O₂ during its migration from the bulk solvent to the active site or backwards. Paying attention to the last hydrophobic cavity that apparently serves as O₂ reservoir for the active site, it was noticed that its volume undergoes ample fluctuations during the MD simulation, as expected from the thermal motion of a flexible protein, independently from the particular subunit and no matter whether the cavity is filled or not by O₂.     

Domain motions of the Mip protein from Legionella pneumophila

The homodimeric 45.6 kDa (total mass) Mip protein, a virulence factor from Legionella pneumophila, was investigated with solution NMR spectroscopy and molecular dynamics (MD) simulations. Two Mip monomers are dimerized via an N-terminal helix bundle that is connected via a long alpha-helix to a C-terminal FKBP domain in each subunit. More than 85% of the amino acids were identified in triple-resonance NMR spectra. (15)N relaxation analysis showed a bimodal distribution of R(1)/R(2) values, with the lower ratio in the N-terminal domain. Relaxation dispersion measurements confirmed that these reduced ratios did not originate from conformational exchange. Thus, two different correlation times (tau(c)) can be deduced, reflecting partly uncoupled motions of both domains. Relaxation data of a Mip(77)(-)(213) monomer mutant were similar to those observed in the dimer, corroborating that the FKBP domain, including part of the connecting helix, behaves as one dynamic entity. MD simulations (18 ns) of the Mip dimer also yielded two different correlation times for the two domains and thus confirm the independence of the domain motions. Principal component analysis of the dihedral space covariance matrix calculated from the MD trajectory suggests a flexible region in the long connecting helix that acts as a hinge between the two domains. Such motion provides a possible explanation of how Mip can bind to complex molecular components of the extracellular matrix and mediate alveolar damage and bacterial spread in the lung.     

Molecular dynamics simulations of DNA and a protein-DNA complex including solvent

The results of a recent nanosecond (ns) molecular dynamics (MD) simulation of the d(CGCGAATTCGCG) double helix in water and a 100 ps MD study of the repressor-operator complex are described. The DNA simulations are analyzed in terms of the structural dynamics, fluctuations in the groove width and bending of the helical axis. The results indicate that the ns dynamical trajectory progresses through a series of three substates of B form DNA, with lifetimes of the order of hundreds of picoseconds (ps). An incipient dynamical equilibrium is evident. A comparison of the calculated axis bending with that observed in corresponding crystal structure data is presented. Simulation of the DNA in complex with the protein and that of the free DNA in solution, starting from the crystal conformation, reveal the dynamical changes that occur on complex formation.     

Effects of protonation state of Asp181 and position of active site water molecules on the conformation of PTP1B

In protein tyrosine phosphatase 1B (PTP1B), the flexible WPD loop adopts a closed conformation (WPD_closed) in the active state of PTP1B, bringing the catalytic Asp181 close to the active site pocket, while WPD loop is in an open conformation (WPD_open) in the inactive state. Previous studies showed that Asp181 may be protonated at physiological pH, and ordered water molecules exist in the active site. In the current study, molecular dynamics simulations are employed at different Asp181 protonation states and initial positions of active site water molecules, and compared with the existing crystallographic data of PTP1B. In WPD_closed conformation, the active site is found to maintain its conformation only in the protonated state of Asp181 in both free and liganded states, while Asp181 is likely to be deprotonated in WPD_open conformation. When the active site water molecule network that is a part of the free WPD_closed crystal structure is disrupted, intermediate WPD loop conformations, similar to that in the PTPRR crystal structure, are sampled in the MD simulations. In liganded PTP1B, one active site water molecule is found to be important for facilitating the orientation of Cys215 and the phosphate ion, thus may play a role in the reaction. In conclusion, conformational stability of WPD loop, and possibly catalytic activity of PTP1B, is significantly affected by the protonation state of Asp181 and position of active site water molecules, showing that these aspects should be taken into consideration both in MD simulations and inhibitor design. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Comparison of simulated and experimentally determined dynamics for a variant of the Lacl DNA-binding domain, Nlac-P.

Recent advances in the experimentally determined structures and dynamics of the domains within LacI provide a rare context for evaluating dynamics calculations. A 1500-ps trajectory was simulated for a variant of the LacI DNA-binding domain, which consists of the first three helices in LacI and the hinge helix of the homologous PurR. Order parameters derived from dynamics simulations are compared to those obtained for the LacI DNA-binding domain with 15N relaxation NMR spectroscopy (Slijper et al., 1997. Biochemistry. 36:249-254). The MD simulations suggest that the unstructured loop between helices II and III does not exist in a discrete state under the conditions of no salt and neutral pH, but occupies a continuum of states between the DNA-bound and free structures. Simulations also indicate that the unstructured region between helix III and the hinge helix is very mobile, rendering motions of the hinge helix essentially independent of the rest of the protein. Finally, the alpha-helical hydrogen bonds in the hinge helix are broken after 1250 ps, perhaps as a prelude to helix unfolding.     

Molecular insight into pseudolysin inhibition using the MM-PBSA and LIE methods

Pseudolysin, the extracellullar elastase of Pseudomonas aeruginosa (EC: 3.4.24.26) plays an important role in the pathogenesis of P. aeruginosa infections. In the present study, molecular dynamics simulations and theoretical affinity predictions were used to gain molecular insight into pseudolysin inhibition. Four low molecular weight inhibitors were docked at their putative binding sites and molecular dynamics (MD) simulations were performed for 5.0 ns, and the free energy of binding was calculated by the linear interaction energy method. The number and the contact surface area of stabilizing hydrophobic, aromatic, and hydrogen bonding interactions appears to reflect the affinity differences between the inhibitors. The proteinaceous inhibitor, Streptomyces metalloproteinase inhibitor (SMPI) was docked in three different binding positions and MD simulations were performed for 3.0 ns. The MD trajectories were used for molecular mechanics-Poisson-Boltzmann surface area analysis of the three binding positions. Computational alanine scanning of the average pseudolysin-SMPI complexes after MD revealed residues at the pseudolysin-SMPI interface giving the main contribution to the free energy of binding. The calculations indicated that SMPI interacts with pseudolysin via the rigid active site loop, but that also contact sites outside this loop contribute significantly to the free energy of association.     

Molecular Dynamic Simulations of the N-Terminal Receiver Domain of NtrC Reveal Intrinsic Conformational Flexibility in the Inactive State

Abstract

The N-terminal receiver domain of NtrC is the molecular switch in the two-component signal transduction. It is the first protein where structures of both the active (phosphyroylated) and inactive (unphosphyroylated) states are determined experimentally. Phosphorylation of the NtrC at the active site induces large structural change. NMR experiments suggested that the wild type unphosphorylated NtrC adopts both the active and the inactive conformations and the phosphorylation stabilizes the active conformations. We applied free (unconstrained) molecular dynamic (MD) simulation to examine the intrinsic flexibilities and stabilities of the NtrC receiver domain in both the active and inactive conformations. Molecular dynamic simulations showed that the inactive state of NtrC receiver domain is more flexible than the active state. There were large movements in helix 4 and loop β3-α3 which coincide with major structural differences between the inactive and active states. We observed large root-mean-square deviations from the initial starting structure and the large root-mean-square fluctuations during MD simulation for the inactive state. We then investigated the activation pathway with Targeted MD simulation. We show that the intrinsic flexibility in the loop β3-α3 plays an important role in triggering the conformational change. Phosphorylation at the active site may serve to stabilize the conformational change. These results together suggest that the unphosphorylated NtrC receiver domain could be involved in a conformational equilibrium between two different states.     

Molecular dynamic simulations of the N-terminal receiver domain of NtrC reveal intrinsic conformational flexibility in the inactive state

The N-terminal receiver domain of NtrC is the molecular switch in the two-component signal transduction. It is the first protein where structures of both the active (phosphyroylated) and inactive (unphosphyroylated) states are determined experimentally. Phosphorylation of the NtrC at the active site induces large structural change. NMR experiments suggested that the wild type unphosphorylated NtrC adopts both the active and the inactive conformations and the phosphorylation stabilizes the active conformations. We applied free (unconstrained) molecular dynamic (MD) simulation to examine the intrinsic flexibilities and stabilities of the NtrC receiver domain in both the active and inactive conformations. Molecular dynamic simulations showed that the inactive state of NtrC receiver domain is more flexible than the active state. There were large movements in helix 4 and loop beta3-alpha3 which coincide with major structural differences between the inactive and active states. We observed large root-mean-square deviations from the initial starting structure and the large root-mean-square fluctuations during MD simulation for the inactive state. We then investigated the activation pathway with Targeted MD simulation. We show that the intrinsic flexibility in the loop beta3-alpha3 plays an important role in triggering the conformational change. Phosphorylation at the active site may serve to stabilize the conformational change. These results together suggest that the unphosphorylated NtrC receiver domain could be involved in a conformational equilibrium between two different states.     

Dynamics of the WPD loop of the Yersinia protein tyrosine phosphatase

The bacterial protein tyrosine phosphatase YopH is an essential virulence determinant in Yersinia spp., causing gastrointestinal diseases and the plague. Like eukaryotic PTPases, YopH catalyzes the hydrolysis of the phosphate moiety of phosphotyrosine within a highly conserved binding pocket, which is also characterized by the closure of the so-called "WPD loop" upon ligand binding. In this study, we investigate the conformational changes and dynamics of the WPD loop by molecular dynamics simulations. Consistent with experimental observations, our simulations show that the WPD loop of YopH is intrinsically flexible and fluctuates between the open and closed conformation with a frequency of approximately 4 ns for the apo, native protein. The region of helix alpha4 spanning loop 384-392, which has been revealed experimentally as a second substrate-binding site in YopH, is found to be highly associated with the WPD loop, stabilizing it in the closed, active conformation, and providing a structural basis for the cooperation of the second-substrate binding site in substrate recognition. Loop L4 (residues 323-327) is shown to be involved in a parallel, correlated motion mode with the WPD loop that contributes the stabilization of a more extended open conformation. In addition, we have simulated the loop reopening in the ligand-bound protein complex by applying the locally enhanced sampling method. Finally, the dynamic behavior of the WPD loop for the C403S mutant differs from the wild-type YopH remarkably. These results shed light on the role of the WPD loop in PTPase-mediated catalysis, and are useful in structure-based design for novel, selective YopH inhibitors as antibacterial drugs.     

A molecular dynamics study of the ligand release path in yeast cytosine deaminase

Yeast cytosine deaminase, a zinc metalloenzyme, catalyzes the deamination of cytosine to uracil. Experimental and computational evidence indicates that the rate-limiting step is product release, instead of the chemical reaction step. In this work, we use molecular dynamics to suggest ligand exit paths. Simulation at 300 K shows that the active site is well protected by the C-terminal helix (residues 150-158) and F-114 loop (residues 111-117) and that on the molecular dynamics timescale water does not flow in or out of the active site. In contrast, simulation at 320 K shows a significant increase in flexibility of the C-terminal helix and F-114 loop. The motions of these two regions at 320 K open the active site and permit water molecules to diffuse into and out of the active site through two paths with one much more favored than the other. Cytosine is pushed out of the active site by a restraint method in two directions specified by these two paths. In path 1 the required motion of the protein is local-involving only the C-terminal helix and F-114 loop-and two residues, F-114 and I-156, are identified that have to be moved away to let cytosine out; whereas in path 2, the protein has to rearrange itself much more extensively, and the changes are also much larger compared to the path 1 simulation.     

Ribosomal RNA Kink-turn Motif—A Flexible Molecular Hinge

Abstract

Ribosomal RNA K-turn motifs are asymmetric internal loops characterized by a sharp bend in the phosphodiester backbone resulting in “V” shaped structures, recurrently observed in ribosomes and showing a high degree of sequence conservation. We have carried out extended explicit solvent molecular dynamics simulations of selected K-turns, in order to investigate their intrinsic structural and dynamical properties. The simulations reveal an unprecedented dynamical flexibility of the K-turns around their X-ray geometries. The K-turns sample, on the nanosecond timescale, different conformational substates. The overall behavior of the simulations suggests that the sampled geometries are essentially isoenergetic and separated by minimal energy barriers. The nanosecond dynamics of isolated K-turns can be qualitatively considered as motion of two rigid helix stems controlled by a very flexible internal loop which then leads to substantial hinge-like motions between the two stems. This internal dynamics of K-turns is strikingly different for example from the bacterial 5S rRNA Loop E motif or BWYV frameshifting pseudoknot which appear to be rigid in the same type of simulations. Bistability and flexibility of K-turns was also suggested by several recent biochemical studies. Although the results of MD simulations should be considered as a qualitative picture of the K-turn dynamics due to force field and sampling limitations, the main advantage of the MD technique is its ability to investigate the region close to K-turn riboso- mal-like geometries. This part of the conformational space is not well characterized by the solution experiments due to large-scale conformational changes seen in the experiments. We suggest that K-turns are well suited to act as flexible structural elements of ribosomal RNA. They can for example be involved in mediation of large-scale motions or they can allow a smooth assembling of the other parts of the ribosome.