Chia-en A. Chang

Evaluating the Accuracy of the Quasiharmonic Approximation

The quasiharmonic approximation (QH) allows the configurational entropy of a molecule to be estimated on the basis of a molecular dynamics simulation, through construction of a Gaussian probability distribution of conformations with variances equal to those provided by the simulation. At its introduction in 1981, the QH method was successfully applied to simple molecular systems with only one highly occupied energy well, and fluctuations were analyzed in a system of internal bond-angle-torsion coordinates. However, more recent studies have applied the QH method to complex biomolecular systems and have relied upon Cartesian coordinates. The present study evaluates the accuracy of the QH method through comparisons with more detailed methods. The chief findings are that the QH method can markedly overestimate the configurational entropy for systems with multiple occupied energy wells and that such errors tend to be magnified by the use of Cartesian coordinates instead of bond-angle-torsion coordinates.

Structural basis for selective activation of ABA receptors.

Changing environmental conditions and lessening fresh water supplies have sparked intense interest in understanding and manipulating abscisic acid (ABA) signaling, which controls adaptive responses to drought and other abiotic stressors. We recently discovered a selective ABA agonist, pyrabactin, and used it to discover its primary target PYR1, the founding member of the PYR/PYL family of soluble ABA receptors. To understand pyrabactins selectivity, we have taken a combined structural, chemical and genetic approach. We show that subtle differences between receptor binding pockets control ligand orientation between productive and nonproductive modes. Nonproductive binding occurs without gate closure and prevents receptor activation. Observations in solution show that these orientations are in rapid equilibrium that can be shifted by mutations to control maximal agonist activity. Our results provide a robust framework for the design of new agonists and reveal a new mechanism for agonist selectivity.

Entropic contributions and the influence of the hydrophobic environment in promiscuous protein–protein association

The mechanisms by which a promiscuous protein can strongly interact with several different proteins using the same binding interface are not completely understood. An example is protein kinase A (PKA), which uses a single face on its docking/dimerization domain to interact with multiple A-kinase anchoring proteins (AKAP) that localize it to different parts of the cell. In the current study, the configurational entropy contributions to the binding between the AKAP protein HT31 with the D/D domain of RII α-regulatory subunit of PKA were examined. The results show that the majority of configurational entropy loss for the interaction was due to decreased fluctuations within rotamer states of the side chains. The result is in contrast to the widely held approximation that the decrease in the number of rotamer states available to the side chains forms the major component. Further analysis showed that there was a direct linear relationship between total configurational entropy and the number of favorable, alternative contacts available within hydrophobic environments. The hydrophobic binding pocket of the D/D domain provides alternative contact points for the side chains of AKAP peptides that allow them to adopt different binding conformations. The increase in binding conformations provides an increase in binding entropy and hence binding affinity. We infer that a general strategy for a promiscuous protein is to provide alternative contact points at its interface to increase binding affinity while the plasticity required for binding to multiple partners is retained. Implications are discussed for understanding and treating diseases in which promiscuous protein interactions are used.

Binding Pathways of Ligands to HIV-1 Protease: Coarse-grained and Atomistic Simulations

Multiscale simulations (coarse‐grained Brownian dynamics simulations and all‐atom molecular dynamics simulations in implicit solvent) were applied to reveal the binding processes of ligands as they enter the binding site of the HIV‐1 protease. The initial structures used for the molecular dynamics simulations were generated based on the Brownian dynamics trajectories, and this is the first molecular dynamics simulation of modeling the association of a ligand with the protease. We found that a protease substrate successfully binds to the protein when the flaps are fully open. Surprisingly, a smaller cyclic urea inhibitor (XK263) can reach the binding site when the flaps are not fully open. However, if the flaps are nearly closed, the inhibitor must rearrange or binding can fail because the inhibitor cannot attain proper conformations to enter the binding site. Both the peptide substrate and XK263 can also affect the proteins internal motion, which may help the flaps to open. Simulations allow us to efficiently study the ligand binding processes and may help those who study drug discovery to find optimal association pathways and to design those ligands with the best binding kinetics.

Endosidin2 targets conserved exocyst complex subunit EXO70 to inhibit exocytosis

Significance The exocyst complex is a conserved protein complex that tethers the secretory vesicles to the site of membrane fusion during exocytosis, an essential cellular process that transports molecules, such as protein, to the cell surface or extracellular space. We identified a small molecule that targets the EXO70 (exocyst component of 70 kDa) subunit of the exocyst complex to inhibit exocytosis. This compound made it possible to control the dynamics of the exocytosis process in a dosage-dependent manner in different organisms and overcame the mutant lethality and genetic redundancy issues in studying mechanisms of exocyst complex regulation. Further design of molecules with higher affinity and more potent activity may make it possible to use drugs to control human diseases related to exocytosis, such as cancer and diabetes. The exocyst complex regulates the last steps of exocytosis, which is essential to organisms across kingdoms. In humans, its dysfunction is correlated with several significant diseases, such as diabetes and cancer progression. Investigation of the dynamic regulation of the evolutionarily conserved exocyst-related processes using mutants in genetically tractable organisms such as Arabidopsis thaliana is limited by the lethality or the severity of phenotypes. We discovered that the small molecule Endosidin2 (ES2) binds to the EXO70 (exocyst component of 70 kDa) subunit of the exocyst complex, resulting in inhibition of exocytosis and endosomal recycling in both plant and human cells and enhancement of plant vacuolar trafficking. An EXO70 protein with a C-terminal truncation results in dominant ES2 resistance, uncovering possible distinct regulatory roles for the N terminus of the protein. This study not only provides a valuable tool in studying exocytosis regulation but also offers a potentially new target for drugs aimed at addressing human disease.

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 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 favors miRNA duplexes with no middle mismatches, whereas AGO1 tolerates, or prefers, duplexes with central mismatches. AGO structure modeling 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.

The role of oligomerization and cooperative regulation in protein function: the case of tryptophan synthase.

The oligomerization/co-localization of protein complexes and their cooperative regulation in protein function is a key feature in many biological systems. The synergistic regulation in different subunits often enhances the functional properties of the multi-enzyme complex. The present study used molecular dynamics and Brownian dynamics simulations to study the effects of allostery, oligomerization and intermediate channeling on enhancing the protein function of tryptophan synthase (TRPS). TRPS uses a set of α/β–dimeric units to catalyze the last two steps of L-tryptophan biosynthesis, and the rate is remarkably slower in the isolated monomers. Our work shows that without their binding partner, the isolated monomers are stable and more rigid. The substrates can form fairly stable interactions with the protein in both forms when the protein reaches the final ligand–bound conformations. Our simulations also revealed that the α/β–dimeric unit stabilizes the substrate–protein conformation in the ligand binding process, which lowers the conformation transition barrier and helps the protein conformations shift from an open/inactive form to a closed/active form. Brownian dynamics simulations with a coarse-grained model illustrate how protein conformations affect substrate channeling. The results highlight the complex roles of protein oligomerization and the fine balance between rigidity and dynamics in protein function.

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.

Protonation States of the Tryptophan Synthase Internal Aldimine Active Site from Solid-State NMR Spectroscopy: Direct Observation of the Protonated Schiff Base Linkage to Pyridoxal-5′-Phosphate

The acid–base chemistry that drives catalysis in pyridoxal-5′-phosphate (PLP)-dependent enzymes has been the subject of intense interest and investigation since the initial identification of PLP’s role as a coenzyme in this extensive class of enzymes. It was first proposed over 50 years ago that the initial step in the catalytic cycle is facilitated by a protonated Schiff base form of the holoenzyme in which the linking lysine ε-imine nitrogen, which covalently binds the coenzyme, is protonated. Here we provide the first 15N NMR chemical shift measurements of such a Schiff base linkage in the resting holoenzyme form, the internal aldimine state of tryptophan synthase. Double-resonance experiments confirm the assignment of the Schiff base nitrogen, and additional 13C, 15N, and 31P chemical shift measurements of sites on the PLP coenzyme allow a detailed model of coenzyme protonation states to be established.

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.