All-atom simulations of structures and energetics of c-di-GMP-bound and free PleD

Frequency response of a protein to local conformational perturbations

Signals created by local perturbations are known to propagate long distances through proteins via backbone connectivity and nonbonded interactions. In the current study, signal propagation from the flexible ligand binding loop to the rest of Protein Tyrosine Phosphatase 1B (PTP1B) was investigated using frequency response techniques. Using restrained Targeted Molecular Dynamics (TMD) potential on WPD and R loops, PTP1B was driven between its crystal structure conformations at different frequencies. Propagation of the local perturbation signal was manifested via peaks at the fundamental frequency and upper harmonics of 1/f distributed spectral density of atomic variables, such as Cα atoms, dihedral angles, or polar interaction distances. Frequency of perturbation was adjusted high enough (simulation length >∼10×period of a perturbation cycle) not to be clouded by random diffusional fluctuations, and low enough (<∼0.8 ns(-1)) not to attenuate the propagating signal and enhance the contribution of the side-chains to the dissipation of the signals. Employing Discrete Fourier Transform (DFT) to TMD simulation trajectories of 16 cycles of conformational transitions at periods of 1.2 to 5 ns yielded Cα displacements consistent with those obtained from crystal structures. Identification of the perturbed atomic variables by statistical t-tests on log-log scale spectral densities revealed the extent of signal propagation in PTP1B, while phase angles of the filtered trajectories at the fundamental frequency were used to cluster collectively fluctuating elements. Hydrophobic interactions were found to have a higher contribution to signal transduction between side-chains compared to the role of polar interactions. Most of in-phase fluctuating residues on the signaling pathway were found to have high identity among PTP domains, and located over a wide region of PTP1B including the allosteric site. Due to its simplicity and efficiency, the suggested technique may find wide applications in identification of signaling pathways of different proteins.

Coupled motions during dynamics reveal a tunnel toward the active site regulated by the N-terminal α-helix in an acylaminoacyl peptidase

Acylaminoacyl peptidase (AAP) subfamily belongs to the prolyl oligopeptidase (POP) family of serine-proteases. There is a great interest in the definition of molecular mechanisms related to the activity and substrate recognition of these complex multi-domain enzymes. The active site relies at the interface between the C-terminal catalytic domain and the β-propeller domain, whose N-terminal region acts as a bridge to the hydrolase domain. In AAP, the N-terminal extension is characterized by a structurally conserved α1-helix, which is known to affect thermal stability and thermal dependence of the catalytic activity. In the present contribution, results from hundreds nanosecond all-atom molecular dynamics simulations, along with analyses of the networks of cross-correlated motions of a member of the AAP subfamily are discussed. The MD investigation identifies a tunnel that from the surrounding of the N-terminal α1-helix bring to the catalytic site. This cavity seems to be regulated by conformational changes of the α1-helix itself during the dynamics. The evidence here provided can be a useful guide for a better understanding of the mechanistic aspects related to AAP activity, but also for drug design purposes.

Stability and dynamics of cyclic diguanylic acid in solution

Cyclic diguanylic acid (CDG) is a ubiquitous messenger involved in bacterial signaling networks. Despite its central role in motility, biofilm formation, virulence, and flagellum development, fundamental properties such as its aggregation state are still poorly understood. Here the dynamics and stability of metal-free and Mg(2+)-bound CDG are characterized. Atomistic simulations establish that the CDG dimer is slightly favored (by -5 kcal mol(-1)) over its dissociated form (2 CDG), while the Mg(2+) ion coordinated in the X-ray structure readily dissociates from (CDG)(2) in solution and prefers water coordination. As a ligand in a protein, CDG binds both as a U-shaped and a quasilinear monomer. The current results indicate that the energy difference between these two conformations is only a few kilocalories per mole, which explains the facile adaptation to different protein environments. This, together with the slight preference of (CDG)(2) over 2 CDG suggests that (CDG)(2) binding to a protein does probably not occur via sequential binding of two individual monomers.

Effect of ligand binding on the intraminimum dynamics of proteins

Effects of ligand binding on protein dynamics are studied via molecular dynamics (MD) simulations on two different enzymes, dihydrofolate reductase (DHFR) and triosephosphate isomerase (TIM), in their unliganded (free) and liganded states. Domain motions in MD trajectories are analyzed by collectivities and rotation angles along the principal components (PCs). DHFR in the free state has well-defined domain rotations, whereas rotations are slightly damped in the binary complex with nicotinamide adenine dinucleotide phosphate (NADPH), and remarkably distorted in the presence of NADP(+) , showing that NADP(+) is solely responsible for the loss of correlation of the domains in DHFR. Although mean square fluctuations of MD simulations in the same PC subspaces are similar for different ligation states, linear stochastic time series models show that backbone flexibility along the first five PCs is decreased upon NADPH and NADP(+) binding in subpicosecond scale. This shows that mobility of the protein along the PCs is closely related with intraminimum dynamics, and alterations in ligation states may change the intraminimum dynamics significantly. Low vibrational frequencies of the alpha-carbon atoms of DHFR are determined from the time series models of a larger number of low indexed PCs, and it is found that number of modes in the lowest frequencies is reduced upon ligand binding. A similar result is obtained for TIM in the unliganded and dihydroxyacetone phosphate bound states. We suggest that stochastic time series modeling is a promising method to be used in determining subtle perturbations in protein dynamics.

Structural and mechanistic determinants of c-di-GMP signalling

Bis-(3'-5')-cyclic dimeric GMP (c-di-GMP) is a ubiquitous second messenger that regulates cell surface-associated traits in bacteria. Components of this regulatory network include GGDEF and EAL domain-containing proteins that determine the cellular concentrations of c-di-GMP by mediating its synthesis and degradation, respectively. Crystal structure analyses in combination with functional studies have revealed the catalytic mechanisms and regulatory principles involved. Downstream, c-di-GMP is recognized by PilZ domain-containing receptors that can undergo large-scale domain rearrangements on ligand binding. Here, we review recent data on the structure and functional properties of the protein families that are involved in c-di-GMP signalling and discuss the mechanistic implications.