Investigation of Structural Dynamics of Enzymes and Protonation States of Substrates Using Computational Tools.

This review discusses the use of molecular modeling tools, together with existing experimental findings, to provide a complete atomic-level description of enzyme dynamics and function. We focus on functionally relevant conformational dynamics of enzymes and the protonation states of substrates. The conformational fluctuations of enzymes usually play a crucial role in substrate recognition and catalysis. Protein dynamics can be altered by a tiny change in a molecular system such as different protonation states of various intermediates or by a significant perturbation such as a ligand association. Here we review recent advances in applying atomistic molecular dynamics (MD) simulations to investigate allosteric and network regulation of tryptophan synthase (TRPS) and protonation states of its intermediates and catalysis. In addition, we review studies using quantum mechanics/molecular mechanics (QM/MM) methods to investigate the protonation states of catalytic residues of β-Ketoacyl ACP synthase I (KasA). We also discuss modeling of large-scale protein motions for HIV-1 protease with coarse-grained Brownian dynamics (BD) simulations

Mechanism of PhosphoThreonine/Serine Recognition and Specificity for Modular Domains from All-atom Molecular Dynamics

<p>Abstract</p> <p>Background</p> <p>Phosphopeptide-binding domains mediate many vital cellular processes such as signal transduction and protein recognition. We studied three well-known domains important for signal transduction: BRCT repeats, WW domain and forkhead-associated (FHA) domain. The first two recognize both phosphothreonine (pThr) and phosphoserine (pSer) residues, but FHA has high specificity for pThr residues. Here we used molecular dynamics (MD) simulations to reveal how FHA exclusively chooses pThr and how BRCT and WW recognize both pThr/pSer. The work also investigated the energies and thermodynamic information of intermolecular interactions.</p> <p>Results</p> <p>Simulations carried out included wide-type and mutated systems. Through analysis of MD simulations, we found that the conserved His residue defines dual loops feature of the FHA domain, which creates a small cavity reserved for only the methyl group of pThr. These well-organized loop interactions directly response to the pThr binding selectivity, while single loop (the 2nd phosphobinding site of FHA) or in combination with α-helix (BRCT repeats) or β-sheet (WW domain) fail to differentiate pThr/pSer.</p> <p>Conclusions</p> <p>Understanding the domain pre-organizations constructed by conserved residues and the driving force of domain-phosphopeptide recognition provides structural insight into pThr specific binding, which also helps in engineering proteins and designing peptide inhibitors.</p

Insights from Free-Energy Calculations: Protein Conformational Equilibrium, Driving Forces, and Ligand-Binding Modes

AbstractAccurate free-energy calculations provide mechanistic insights into molecular recognition and conformational equilibrium. In this work, we performed free-energy calculations to study the thermodynamic properties of different states of molecular systems in their equilibrium basin, and obtained accurate absolute binding free-energy calculations for protein-ligand binding using a newly developed M2 algorithm. We used a range of Asp-Phe-Gly (DFG)-in/out p38α mitogen-activated protein kinase inhibitors as our test cases. We also focused on the flexible DFG motif, which is closely connected to kinase activation and inhibitor binding. Our calculations explain the coexistence of DFG-in and DFG-out states of the loop and reveal different components (e.g., configurational entropy and enthalpy) that stabilize the apo p38α conformations. To study novel ligand-binding modes and the key driving forces behind them, we computed the absolute binding free energies of 30 p38α inhibitors, including analogs with unavailable experimental structures. The calculations revealed multiple stable, complex conformations and changes in p38α and inhibitor conformations, as well as balance in several energetic terms and configurational entropy loss. The results provide relevant physics that can aid in designing inhibitors and understanding protein conformational equilibrium. Our approach is fast for use with proteins that contain flexible regions for structure-based drug design

ARGONAUTE PIWI domain and microRNA duplex structure regulate small RNA sorting in Arabidopsis.

Small RNAs (sRNAs) are loaded into ARGONAUTE (AGO) proteins to induce gene silencing. In plants, the 5'-terminal nucleotide is important for sRNA sorting into different AGOs. Here we show that microRNA (miRNA) duplex structure also contributes to miRNA sorting. Base pairing at the 15th nucleotide of a miRNA duplex is important for miRNA sorting in both Arabidopsis AGO1 and AGO2. AGO2 favours miRNA duplexes with no middle mismatches, whereas AGO1 tolerates, or prefers, duplexes with central mismatches. AGO structure modelling and mutational analyses reveal that the QF-V motif within the conserved PIWI domain contributes to recognition of base pairing at the 15th nucleotide of a duplex, while the DDDE catalytic core of AtAGO2 is important for recognition of the central nucleotides. Finally, we rescued the adaxialized phenotype of ago1-12, which is largely due to miR165 loss-of-function, by changing miR165 duplex structure which we predict redirects it to AGO2

Ligand Binding Pathways and Conformational Transitions of the HIV Protease

It is important to determine the binding pathways and mechanisms of ligand molecules to target proteins to effectively design therapeutic drugs. Molecular dynamics (MD) is a promising computational tool that allows us to simulate protein–drug binding at an atomistic level. However, the gap between the time scales of current simulations and those of many drug binding processes has limited the usage of conventional MD, which has been reflected in studies of the HIV protease. Here, we have applied a robust enhanced simulation method, Gaussian accelerated molecular dynamics (GaMD), to sample binding pathways of the XK263 ligand and associated protein conformational changes in the HIV protease. During two of 10 independent GaMD simulations performed over 500–2500 ns, the ligand was observed to successfully bind to the protein active site. Although GaMD-derived free energy profiles were not fully converged because of insufficient sampling of the complex system, the simulations still allowed us to identify relatively low-energy intermediate conformational states during binding of the ligand to the HIV protease. Relative to the X-ray crystal structure, the XK263 ligand reached a minimum root-mean-square deviation (RMSD) of 2.26 Å during 2.5 μs of GaMD simulation. In comparison, the ligand RMSD reached a minimum of only ~5.73 Å during an earlier 14 μs conventional MD simulation. This work highlights the enhanced sampling power of the GaMD approach and demonstrates its wide applicability to studies of drug–receptor interactions for the HIV protease and by extension many other target proteins

Detection of concealed and buried chemicals by using multifrequency excitations

This is the published version. Copyright © 2010 American Institute of PhysicsIn this paper, we present a new type of concealed and buried chemical detection system by stimulating and enhancing spectroscopic signatures with multifrequency excitations, which includes a low frequency gradient dcelectric field, a high frequency microwave field, and higher frequency infrared (IR) radiations. Each excitation frequency plays a unique role. The microwave, which can penetrate into the underground and/or pass through the dielectric covers with low attenuation, could effectively transform its energy into the concealed and buried chemicals and increases its evaporation rate from the sample source. Subsequently, a gradient dcelectric field, generated by a Van De Graaff generator, not only serves as a vapor accelerator for efficiently expediting the transportation process of the vapor release from the concealed and buried chemicals but also acts as a vapor concentrator for increasing the chemical concentrations in the detection area, which enables the trace level chemical detection. Finally, the stimulated and enhanced vapors on the surface are detected by the IR spectroscopic fingerprints. Our theoretical and experimental results demonstrate that more than sixfold increase in detection signal can be achieved by using this proposed technology. The proposed technology can also be used for standoff detection of concealed and buried chemicals by adding the remote IR and/or thermal spectroscopic and imaging detection systems

Dynamics of the Acetylcholinesterase Tetramer

Acetylcholinesterase rapidly hydrolyzes the neurotransmitter acetylcholine in cholinergic synapses, including the neuromuscular junction. The tetramer is the most important functional form of the enzyme. Two low-resolution crystal structures have been solved. One is compact with two of its four peripheral anionic sites (PAS) sterically blocked by complementary subunits. The other is a loose tetramer with all four subunits accessible to solvent. These structures lacked the C-terminal amphipathic t-peptide (WAT domain) that interacts with the proline-rich attachment domain (PRAD). A complete tetramer model (AChEt) was built based on the structure of the PRAD/WAT complex and the compact tetramer. Normal mode analysis suggested that AChEt could exist in several conformations with subunits fluctuating relative to one another. Here, a multiscale simulation involving all-atom molecular dynamics and Cα-based coarse-grained Brownian dynamics simulations was carried out to investigate the large-scale intersubunit dynamics in AChEt. We sampled the ns-μs timescale motions and found that the tetramer indeed constitutes a dynamic assembly of monomers. The intersubunit fluctuation is correlated with the occlusion of the PAS. Such motions of the subunits “gate” ligand-protein association. The gates are open more than 80% of the time on average, which suggests a small reduction in ligand-protein binding. Despite the limitations in the starting model and approximations inherent in coarse graining, these results are consistent with experiments which suggest that binding of a substrate to the PAS is only somewhat hindered by the association of the subunits

T-Analyst: a program for efficient analysis of protein conformational changes by torsion angles

T-Analyst is a user-friendly computer program for analyzing trajectories from molecular modeling. Instead of using Cartesian coordinates for protein conformational analysis, T-Analyst is based on internal bond-angle-torsion coordinates in which internal torsion angle movements, such as side-chain rotations, can be easily detected. The program computes entropy and automatically detects and corrects angle periodicity to produce accurate rotameric states of dihedrals. It also clusters multiple conformations and detects dihedral rotations that contribute hinge-like motions. Correlated motions between selected dihedrals can also be observed from the correlation map. T-Analyst focuses on showing changes in protein flexibility between different states and selecting representative protein conformations for molecular docking studies. The program is provided with instructions and full source code in Perl

Diffusional Channeling in the Sulfate-Activating Complex: Combined Continuum Modeling and Coarse-Grained Brownian Dynamics Studies

Enzymes required for sulfur metabolism have been suggested to gain efficiency by restricted diffusion (i.e., channeling) of an intermediate APS2– between active sites. This article describes modeling of the whole channeling process by numerical solution of the Smoluchowski diffusion equation, as well as by coarse-grained Brownian dynamics. The results suggest that electrostatics plays an essential role in the APS2– channeling. Furthermore, with coarse-grained Brownian dynamics, the substrate channeling process has been studied with reactions in multiple active sites. Our simulations provide a bridge for numerical modeling with Brownian dynamics to simulate the complicated reaction and diffusion and raise important questions relating to the electrostatically mediated substrate channeling in vitro, in situ, and in vivo