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.

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.

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.

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.

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.

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.

Functional dynamics of the folded ankyrin repeats of I kappa B alpha revealed by nuclear magnetic resonance.

Inhibition of nuclear factor kappaB (NF-kappaB) is mainly accomplished by IkappaB alpha, which consists of a signal response sequence at the N-terminus, a six-ankyrin repeat domain (ARD) that binds NF-kappaB, and a C-terminal PEST sequence. Previous studies with the ARD revealed that the fifth and sixth repeats are only partially folded in the absence of NF-kappaB. Here we report NMR studies of a truncated version of IkappaB alpha, containing only the first four ankyrin repeats, IkappaB alpha(67-206). This four-repeat segment is well-structured in the free state, enabling full resonance assignments to be made. H-D exchange, backbone dynamics, and residual dipolar coupling (RDC) experiments reveal regions of flexibility. In addition, regions consistent with the presence of micro- to millisecond motions occur periodically throughout the repeat structure. Comparison of the RDCs with the crystal structure gave only moderate agreement, but an ensemble of structures generated by accelerated molecular dynamics gave much better agreement with the measured RDCs. The regions showing flexibility correspond to those implicated in entropic compensation for the loss of flexibility in ankyrin repeats 5 and 6 upon binding to NF-kappaB. The regions showing micro- to millisecond motions in the free protein are the ends of the beta-hairpins that directly interact with NF-kappaB in the complex.

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.

Structural dynamics of protein backbone phi angles: extended molecular dynamics simulations versus experimental (3) J scalar couplings.

(3)J scalar couplings report on the conformational averaging of backbone phi angles in peptides and proteins, and therefore represent a potentially powerful tool for studying the details of both structure and dynamics in solution. We have compared an extensive experimental dataset with J-couplings predicted from unrestrained molecular dynamics simulation using enhanced sampling available from accelerated molecular dynamics or using long timescale trajectories (200 ns). The dynamic fluctuations predicted to be present along the backbone, in agreement with residual dipolar coupling analysis, are compatible with the experimental (3)J scalar couplings providing a slightly better reproduction of these experimental parameters than a high-resolution static structure.

Structural biology by NMR: structure, dynamics, and interactions.

The function of bio-macromolecules is determined by both their 3D structure and conformational dynamics. These molecules are inherently flexible systems displaying a broad range of dynamics on time-scales from picoseconds to seconds. Nuclear Magnetic Resonance (NMR) spectroscopy has emerged as the method of choice for studying both protein structure and dynamics in solution. Typically, NMR experiments are sensitive both to structural features and to dynamics, and hence the measured data contain information on both. Despite major progress in both experimental approaches and computational methods, obtaining a consistent view of structure and dynamics from experimental NMR data remains a challenge. Molecular dynamics simulations have emerged as an indispensable tool in the analysis of NMR data.

Characterization of protein dynamics from residual dipolar couplings using the three dimensional Gaussian axial fluctuation model.

Residual dipolar couplings are potentially very powerful probes of slower protein motions, providing access to dynamic events occurring on functionally important timescales up to the millisecond. One recent approach uses the three dimensional Gaussian Axial Fluctuation model (3D GAF) to determine the major directional modes and associated amplitudes of motions along the peptide chain. In this study we have used standard and accelerated molecular dynamics simulations to determine the accuracy of 3D GAF-based approaches in characterizing the nature and extent of local molecular motions. We compare modes determined directly from the trajectories with motional parameterization derived from RDCs simulated from the same trajectories. Three approaches are tested, that either suppose a known three-dimensional structure, simultaneously determine backbone structure and dynamics, or determine dynamic modes in the absence of a structural model. The results demonstrate the robustness of the 3D GAF analysis even in the presence of large-scale motions, and illustrate the remarkably quantitative nature of the extracted amplitudes. These observations suggest that the approach can be generally used for the study of functionally interesting biomolecular motions.

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.

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.

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.

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.