Accurate Prediction of Amide Exchange in the Fast Limit Reveals Thrombin Allostery.

Amide hydrogen/deuterium exchange mass spectrometry (HDXMS) of proteins has become extremely popular for identifying ligand-binding sites, protein-protein interactions, intrinsic disorder, and allosteric changes upon protein modification. Such phenomena are revealed when amide exchange is measured in the fast limit, that is, within a few minutes of exchange in deuterated buffer. The HDXMS data have a resolution of the length of peptides and are difficult to interpret because many different phenomena lead to changes in hydrogen/deuterium exchange. We present a quantitative analysis of accelerated molecular dynamics simulations that provides impressive agreement with peptide-length HDXMS data. Comparative analysis of thrombin and a single-point mutant reveals that the simulation analysis can distinguish the subtle differences in exchange due to mutation. In addition, the results provide a deeper understanding of the underlying changes in dynamics revealed by the HDXMS that extend far from the site of mutation.

Dynamic Consequences of Mutation of Tryptophan 215 in Thrombin.

Thrombin normally cleaves fibrinogen to promote coagulation; however, binding of thrombomodulin to thrombin switches the specificity of thrombin toward protein C, triggering the anticoagulation pathway. The W215A thrombin mutant was reported to have decreased activity toward fibrinogen without significant loss of activity toward protein C. To understand how mutation of Trp215 may alter thrombin specificity, hydrogen-deuterium exchange experiments (HDXMS), accelerated molecular dynamics (AMD) simulations, and activity assays were carried out to compare the dynamics of Trp215 mutants with those of wild type (WT) thrombin. Variation in NaCl concentration had no detectable effect on the sodium-binding (220sCT) loop, but appeared to affect other surface loops. Trp215 mutants showed significant increases in amide exchange in the 170sCT loop consistent with a loss of H-bonding in this loop identified by the AMD simulations. The W215A thrombin showed increased amide exchange in the 220sCT loop and in the N-terminus of the heavy chain. The AMD simulations showed that a transient conformation of the W215A thrombin has a distorted catalytic triad. HDXMS experiments revealed that mutation of Phe227, which engages in a π-stacking interaction with Trp215, also caused significantly increased amide exchange in the 170sCT loop. Activity assays showed that only the F227V mutant had wild type catalytic activity, whereas all other mutants showed markedly lower activity. Taken together, the results explain the reduced pro-coagulant activity of the W215A mutant and demonstrate the allosteric connection between Trp215, the sodium-binding loop, and the active site.

NMR reveals a dynamic allosteric pathway in thrombin.

Although serine proteases are found ubiquitously in both eukaryotes and prokaryotes, and they comprise the largest of all of the peptidase families, their dynamic motions remain obscure. The backbone dynamics of the coagulation serine protease, apo-thrombin (S195M-thrombin), were compared to the substrate-bound form (PPACK-thrombin). R1, R2, 15N-{1H}NOEs, and relaxation dispersion NMR experiments were measured to capture motions across the ps to ms timescale. The ps-ns motions were not significantly altered upon substrate binding. The relaxation dispersion data revealed that apo-thrombin is highly dynamic, with µs-ms motions throughout the molecule. The region around the N-terminus of the heavy chain, the Na+-binding loop, and the 170 s loop, all of which are implicated in allosteric coupling between effector binding sites and the active site, were dynamic primarily in the apo-form. Most of the loops surrounding the active site become more ordered upon PPACK-binding, but residues in the N-terminal part of the heavy chain, the γ-loop, and anion-binding exosite 1, the main allosteric binding site, retain µs-ms motions. These residues form a dynamic allosteric pathway connecting the active site to the main allosteric site that remains in the substrate-bound form.

Trapping the dynamic acyl carrier protein in fatty acid biosynthesis.

Acyl carrier protein (ACP) transports the growing fatty acid chain between enzymatic domains of fatty acid synthase (FAS) during biosynthesis. Because FAS enzymes operate on ACP-bound acyl groups, ACP must stabilize and transport the growing lipid chain. ACPs have a central role in transporting starting materials and intermediates throughout the fatty acid biosynthetic pathway. The transient nature of ACP-enzyme interactions impose major obstacles to obtaining high-resolution structural information about fatty acid biosynthesis, and a new strategy is required to study protein-protein interactions effectively. Here we describe the application of a mechanism-based probe that allows active site-selective covalent crosslinking of AcpP to FabA, the Escherichia coli ACP and fatty acid 3-hydroxyacyl-ACP dehydratase, respectively. We report the 1.9 Å crystal structure of the crosslinked AcpP-FabA complex as a homodimer in which AcpP exhibits two different conformations, representing probable snapshots of ACP in action: the 4'-phosphopantetheine group of AcpP first binds an arginine-rich groove of FabA, then an AcpP helical conformational change locks AcpP and FabA in place. Residues at the interface of AcpP and FabA are identified and validated by solution nuclear magnetic resonance techniques, including chemical shift perturbations and residual dipolar coupling measurements. These not only support our interpretation of the crystal structures but also provide an animated view of ACP in action during fatty acid dehydration. These techniques, in combination with molecular dynamics simulations, show for the first time that FabA extrudes the sequestered acyl chain from the ACP binding pocket before dehydration by repositioning helix III. Extensive sequence conservation among carrier proteins suggests that the mechanistic insights gleaned from our studies may be broadly applicable to fatty acid, polyketide and non-ribosomal biosynthesis. Here the foundation is laid for defining the dynamic action of carrier-protein activity in primary and secondary metabolism, providing insight into pathways that can have major roles in the treatment of cancer, obesity and infectious disease.

Correlated motions and residual frustration in thrombin.

Thrombin is the central protease in the cascade of blood coagulation proteases. The structure of thrombin consists of a double ß-barrel core surrounded by connecting loops and helices. Compared to chymotrypsin, thrombin has more extended loops that are thought to have arisen from insertions in the serine protease that evolved to impart greater specificity. Previous experiments showed thermodynamic coupling between ligand binding at the active site and distal exosites. We present a combined approach of molecular dynamics (MD), accelerated molecular dynamics (AMD), and analysis of the residual local frustration of apo-thrombin and active-site-bound (PPACK-thrombin). Community analysis of the MD ensembles identified changes upon active site occupation in groups of residues linked through correlated motions and physical contacts. AMD simulations, calibrated on measured residual dipolar couplings, reveal that upon active site ligation, correlated loop motions are quenched, but new ones connecting the active site with distal sites where allosteric regulators bind emerge. Residual local frustration analysis reveals a striking correlation between frustrated contacts and regions undergoing slow time scale dynamics. The results elucidate a motional network that probably evolved through retention of frustrated contacts to provide facile conversion between ensembles of states.

Allosteric networks in thrombin distinguish procoagulant vs. anticoagulant activities.

The serine protease α-thrombin is a dual-action protein that mediates the blood-clotting cascade. Thrombin alone is a procoagulant, cleaving fibrinogen to make the fibrin clot, but the thrombin-thrombomodulin (TM) complex initiates the anticoagulant pathway by cleaving protein C. A TM fragment consisting of only the fifth and sixth EGF-like domains (TM56) is sufficient to bind thrombin, but the presence of the fourth EGF-like domain (TM456) is critical to induce the anticoagulant activity of thrombin. Crystallography of the thrombin-TM456 complex revealed no significant structural changes in thrombin, suggesting that TM4 may only provide a scaffold for optimal alignment of protein C for its cleavage by thrombin. However, a variety of experimental data have suggested that the presence of TM4 may affect the dynamic properties of the active site loops. In the present work, we have used both conventional and accelerated molecular dynamics simulation to study the structural dynamic properties of thrombin, thrombin:TM56, and thrombin:TM456 across a broad range of time scales. Two distinct yet interrelated allosteric pathways are identified that mediate both the pro- and anticoagulant activities of thrombin. One allosteric pathway, which is present in both thrombin:TM56 and thrombin:TM456, directly links the TM5 domain to the thrombin active site. The other allosteric pathway, which is only present on slow time scales in the presence of the TM4 domain, involves an extended network of correlated motions linking the TM4 and TM5 domains and the active site loops of thrombin.

Disease-associated mutations disrupt functionally important regions of intrinsic protein disorder.

The effects of disease mutations on protein structure and function have been extensively investigated, and many predictors of the functional impact of single amino acid substitutions are publicly available. The majority of these predictors are based on protein structure and evolutionary conservation, following the assumption that disease mutations predominantly affect folded and conserved protein regions. However, the prevalence of the intrinsically disordered proteins (IDPs) and regions (IDRs) in the human proteome together with their lack of fixed structure and low sequence conservation raise a question about the impact of disease mutations in IDRs. Here, we investigate annotated missense disease mutations and show that 21.7% of them are located within such intrinsically disordered regions. We further demonstrate that 20% of disease mutations in IDRs cause local disorder-to-order transitions, which represents a 1.7-2.7 fold increase compared to annotated polymorphisms and neutral evolutionary substitutions, respectively. Secondary structure predictions show elevated rates of transition from helices and strands into loops and vice versa in the disease mutations dataset. Disease disorder-to-order mutations also influence predicted molecular recognition features (MoRFs) more often than the control mutations. The repertoire of disorder-to-order transition mutations is limited, with five most frequent mutations (R→W, R→C, E→K, R→H, R→Q) collectively accounting for 44% of all deleterious disorder-to-order transitions. As a proof of concept, we performed accelerated molecular dynamics simulations on a deleterious disorder-to-order transition mutation of tumor protein p63 and, in agreement with our predictions, observed an increased α-helical propensity of the region harboring the mutation. Our findings highlight the importance of mutations in IDRs and refine the traditional structure-centric view of disease mutations. The results of this study offer a new perspective on the role of mutations in disease, with implications for improving predictors of the functional impact of missense mutations.

Exploring the Photophysical Properties of Molecular Systems Using Excited State Accelerated ab Initio Molecular Dynamics.

In the present work, we employ excited state accelerated ab initio molecular dynamics (A-AIMD) to efficiently study the excited state energy landscape and photophysical topology of a variety of molecular systems. In particular, we focus on two important challenges for the modeling of excited electronic states: (i) the identification and characterization of conical intersections and crossing seams, in order to predict different and often competing radiationless decay mechanisms, and (ii) the description of the solvent effect on the absorption and emission spectra of chemical species in solution. In particular, using as examples the Schiff bases formaldimine and salicylidenaniline, we show that A-AIMD can be readily employed to explore the conformational space around crossing seams in molecular systems with very different photochemistry. Using acetone in water as an example, we demonstrate that the enhanced configurational space sampling may be used to accurately and efficiently describe both the prominent features and line-shapes of absorption and emission spectra.

The dynamic structure of thrombin in solution.

The backbone dynamics of human α-thrombin inhibited at the active site serine were analyzed using R(1), R(2), and heteronuclear NOE experiments, variable temperature TROSY 2D [(1)H-(15)N] correlation spectra, and R(ex) measurements. The N-terminus of the heavy chain, which is formed upon zymogen activation and inserts into the protein core, is highly ordered, as is much of the double beta-barrel core. Some of the surface loops, by contrast, remain very dynamic with order parameters as low as 0.5 indicating significant motions on the ps-ns timescale. Regions of the protein that were thought to be dynamic in the zymogen and to become rigid upon activation, in particular the γ-loop, the 180s loop, and the Na(+) binding site have order parameters below 0.8. Significant R(ex) was observed in most of the γ-loop, in regions proximal to the light chain, and in the ß-sheet core. Accelerated molecular dynamics simulations yielded a molecular ensemble consistent with measured residual dipolar couplings that revealed dynamic motions up to milliseconds. Several regions, including the light chain and two proximal loops, did not appear highly dynamic on the ps-ns timescale, but had significant motions on slower timescales.

Enhanced Lipid Diffusion and Mixing in Accelerated Molecular Dynamics.

Accelerated molecular dynamics (aMD) is an enhanced sampling technique that expedites conformational space sampling by reducing the barriers separating various low-energy states of a system. Here, we present the first application of the aMD method on lipid membranes. Altogether, ∼1.5 µs simulations were performed on three systems: a pure POPC bilayer, a pure DMPC bilayer, and a mixed POPC:DMPC bilayer. Overall, the aMD simulations are found to produce significant speedup in trans-gauche isomerization and lipid lateral diffusion versus those in conventional MD (cMD) simulations. Further comparison of a 70-ns aMD run and a 300-ns cMD run of the mixed POPC:DMPC bilayer shows that the two simulations yield similar lipid mixing behaviors, with aMD generating a 2-3-fold speedup compared to cMD. Our results demonstrate that the aMD method is an efficient approach for the study of bilayer structural and dynamic properties. On the basis of simulations of the three bilayer systems, we also discuss the impact of aMD parameters on various lipid properties, which can be used as a guideline for future aMD simulations of membrane systems.

Studying functional dynamics in bio-molecules using accelerated molecular dynamics.

Many biologically important processes such as enzyme catalysis, signal transduction, ligand binding and allosteric regulation occur on the micro- to millisecond time-scale. Despite the sustained and rapid increase in available computational power and the development of efficient simulation algorithms, molecular dynamics (MD) simulations of proteins and bio-machines are generally limited to time-scales of tens to hundreds of nano-seconds. In this perspective article we present a comprehensive review of Accelerated Molecular Dynamics (AMD), an extended biased potential molecular dynamics approach that allows for the efficient study of bio-molecular systems up to time-scales several orders of magnitude greater than those accessible using standard classical MD methods, whilst still maintaining a fully atomistic representation of the system. Compared to many other approaches, AMD affords efficient enhanced conformational space sampling without any a priori understanding of the underlying free energy surface, nor does it require the specific prior definition of a reaction coordinate or a set of collective variables. Successful applications of the AMD method, including the study of slow time-scale functional dynamics in folded proteins and the conformational behavior of natively unstructured proteins are discussed and an outline of the different variants and extensions to the standard AMD approach is presented.

Accelerating chemical reactions: exploring reactive free-energy surfaces using accelerated ab initio molecular dynamics.

A biased potential molecular dynamics simulation approach, accelerated molecular dynamics (AMD), has been implemented in the framework of ab initio molecular dynamics for the study of chemical reactions. Using two examples, the double proton transfer reaction in formic acid dimer and the hypothetical adiabatic ring opening and subsequent rearrangement reactions in methylenecyclopropane, it is demonstrated that ab initio AMD can be readily employed to efficiently explore the reactive potential energy surface, allowing the prediction of chemical reactions and the identification of metastable states. An adaptive variant of the AMD method is developed, which additionally affords an accurate representation of both the free-energy surface and the mechanism associated with the chemical reaction of interest and can also provide an estimate of the reaction rate.

Accessing a hidden conformation of the maltose binding protein using accelerated molecular dynamics.

Periplasmic binding proteins (PBPs) are a large family of molecular transporters that play a key role in nutrient uptake and chemotaxis in Gram-negative bacteria. All PBPs have characteristic two-domain architecture with a central interdomain ligand-binding cleft. Upon binding to their respective ligands, PBPs undergo a large conformational change that effectively closes the binding cleft. This conformational change is traditionally viewed as a ligand induced-fit process; however, the intrinsic dynamics of the protein may also be crucial for ligand recognition. Recent NMR paramagnetic relaxation enhancement (PRE) experiments have shown that the maltose binding protein (MBP) - a prototypical member of the PBP superfamily - exists in a rapidly exchanging (ns to µs regime) mixture comprising an open state (approx 95%), and a minor partially closed state (approx 5%). Here we describe accelerated MD simulations that provide a detailed picture of the transition between the open and partially closed states, and confirm the existence of a dynamical equilibrium between these two states in apo MBP. We find that a flexible part of the protein called the balancing interface motif (residues 175-184) is displaced during the transformation. Continuum electrostatic calculations indicate that the repacking of non-polar residues near the hinge region plays an important role in driving the conformational change. Oscillations between open and partially closed states create variations in the shape and size of the binding site. The study provides a detailed description of the conformational space available to ligand-free MBP, and has implications for understanding ligand recognition and allostery in related proteins.

On the Use of Accelerated Molecular Dynamics to Enhance Configurational Sampling in Ab Initio Simulations.

We have implemented the accelerated molecular dynamics approach (Hamelberg, D.; Mongan, J.; McCammon, J. A. J. Chem. Phys. 2004, 120 (24), 11919) in the framework of ab initio MD (AIMD). Using three simple examples, we demonstrate that accelerated AIMD (A-AIMD) can be used to accelerate solvent relaxation in AIMD simulations and facilitate the detection of reaction coordinates: (i) We show, for one cyclohexane molecule in the gas phase, that the method can be used to accelerate the rate of the chair-to-chair interconversion by a factor of ∼1 × 10(5), while allowing for the reconstruction of the correct canonical distribution of low-energy states; (ii) We then show, for a water box of 64 H(2)O molecules, that A-AIMD can also be used in the condensed phase to accelerate the sampling of water conformations, without affecting the structural properties of the solvent; and (iii) The method is then used to compute the potential of mean force (PMF) for the dissociation of Na-Cl in water, accelerating the convergence by a factor of ∼3-4 compared to conventional AIMD simulations.(2) These results suggest that A-AIMD is a useful addition to existing methods for enhanced conformational and phase-space sampling in solution. While the method does not make the use of collective variables superfluous, it also does not require the user to define a set of collective variables that can capture all the low-energy minima on the potential energy surface. This property may prove very useful when dealing with highly complex multidimensional systems that require a quantum mechanical treatment.

Adaptive Accelerated Molecular Dynamics (Ad-AMD) Revealing the Molecular Plasticity of P450cam.

An extended accelerated molecular dynamics (AMD) methodology called adaptive AMD is presented. Adaptive AMD (Ad-AMD) is an efficient and robust conformational space sampling algorithm that is particularly-well suited to proteins with highly structured potential energy surfaces exhibiting complex, large-scale collective conformational transitions. Ad-AMD simulations of substrate-free P450cam reveal that this system exists in equilibrium between a fully and partially open conformational state. The mechanism for substrate binding depends on the size of the ligand. Larger ligands enter the P450cam binding pocket, and the resulting substrate-bound system is trapped in an open conformation via a population shift mechanism. Small ligands, which fully enter the binding pocket, cause an induced-fit mechanism, resulting in the formation of an energetically stable closed conformational state. These results are corroborated by recent experimental studies and potentially provide detailed insight into the functional dynamics and conformational behavior of the entire cytochrome-P450 superfamily.

Enhanced conformational space sampling improves the prediction of chemical shifts in proteins.

A biased-potential molecular dynamics simulation method, accelerated molecular dynamics (AMD), was combined with the chemical shift prediction algorithm SHIFTX to calculate (1)H(N), (15)N, (13)Calpha, (13)Cbeta, and (13)C' chemical shifts of the ankyrin repeat protein IkappaBalpha (residues 67-206), the primary inhibitor of nuclear factor kappa-B (NF-kappaB). Free-energy-weighted molecular ensembles were generated over a range of acceleration levels, affording systematic enhancement of the conformational space sampling of the protein. We have found that the predicted chemical shifts, particularly for the (15)N, (13)Calpha, and (13)Cbeta nuclei, improve substantially with enhanced conformational space sampling up to an optimal acceleration level. Significant improvement in the predicted chemical shift data coincides with those regions of the protein that exhibit backbone dynamics on longer time scales. Interestingly, the optimal acceleration level for reproduction of the chemical shift data has previously been shown to best reproduce the experimental residual dipolar coupling (RDC) data for this system, as both chemical shift data and RDCs report on an ensemble and time average in the millisecond range.

Cobalt 1,3-Diisopropyl-1H-imidazol-2-ylidene Complexes: Synthesis, Solid-State Structures, and Quantum Chemistry Calculations.

The reaction of (η(5)-C(5)H(5))Co(PPh(3))(2) with 1,3-bis(isopropyl)imidazol-2-ylidene (Im(i)Pr(2)) leads to the formation of (η(5)-C(5)H(5))Co(PPh(3))(Im(i)Pr(2)). In a similar fashion (η(5)-C(5)H(5))Co(Im(i)Pr(2))(CO) is formed from (η(5)-C(5)H(5))Co(CO)(2). The barrier to rotation about the Co-C(carbene) bond in the latter complex has been determined by variable-temperature (1)H NMR spectroscopy (13.6 kcal/mol) and by computation (13.3 kcal/mol). The structural and dynamic properties of (η(5)-C(5)H(5))Co(ImMe(2))(CO) (ImMe(2) = 1,3-dimethylimidazol-2-ylidene) and (η(5)-C(5)H(5))Co(ImAr(2))(CO) (ImAr(2) = 1,3-dimesityl-2-ylidene) were examined by quantum chemistry calculations.

Toward a unified representation of protein structural dynamics in solution.

An atomic resolution description of protein flexibility is essential for understanding the role that structural dynamics play in biological processes. Despite the unique dependence of nuclear magnetic resonance (NMR) to motional averaging on different time scales, NMR-based protein structure determination often ignores the presence of dynamics, representing rapidly exchanging conformational equilibria in terms of a single static structure. In this study, we use the rich dynamic information encoded in experimental NMR parameters to develop a molecular and statistical mechanical characterization of the conformational behavior of proteins in solution. Critically, and in contrast to previously proposed techniques, we do not use empirical energy terms to restrain a conformational search, a procedure that can strongly perturb simulated dynamics in a nonpredictable way. Rather, we use accelerated molecular dynamic simulation to gradually increase the level of conformational sampling and to identify the appropriate level of sampling via direct comparison of unrestrained simulation with experimental data. This constraint-free approach thereby provides an atomic resolution free-energy weighted Boltzmann description of protein dynamics occurring on time scales over many orders of magnitude in the protein ubiquitin.

Mechanism of glycan receptor recognition and specificity switch for avian, swine, and human adapted influenza virus hemagglutinins: a molecular dynamics perspective.

Hemagglutinins (HA's) from duck, swine, and human influenza viruses have previously been shown to prefer avian and human glycan receptor analogues with distinct topological profiles, pentasaccharides LSTa (alpha-2,3 linkage) and LSTc (alpha-2,6 linkage), in comparative molecular dynamics studies. On the basis of detailed analyses of the dynamic motions of the receptor binding domains (RBDs) and interaction energy profiles with individual glycan residues, we have identified approximately 30 residue positions in the RBD that present distinct profiles with the receptor analogues. Glycan binding constrained the conformational space sampling by the HA. Electrostatic steering appeared to play a key role in glycan binding specificity. The complex dynamic behaviors of the major SSE and trimeric interfaces with or without bound glycans suggested that networks of interactions might account for species specificity in these low affinity and high avidity (multivalent) interactions between different HA and glycans. Contact frequency, energetic decomposition, and H-bond analyses revealed species-specific differences in HA-glycan interaction profiles, not readily discernible from crystal structures alone. Interaction energy profiles indicated that mutation events at the set of residues such as 145, 156, 158, and 222 would favor human or avian receptor analogues, often through interactions with distal asialo-residues. These results correlate well with existing experimental evidence, and suggest new opportunities for simulation-based vaccine and drug development.