Molecular simulation uncovers the conformational space of the λ Cro dimer in solution

The significant variation among solved structures of the λ Cro dimer suggests its flexibility. However, contacts in the crystal lattice could have stabilized a conformation which is unrepresentative of its dominant solution form. Here we report on the conformational space of the Cro dimer in solution using replica exchange molecular dynamics in explicit solvent. The simulated ensemble shows remarkable correlation with available x-ray structures. Network analysis and a free energy surface reveal the predominance of closed and semi-open dimers, with a modest barrier separating these two states. The fully open conformation lies higher in free energy, indicating that it requires stabilization by DNA or crystal contacts. Most NMR models are found to be unstable conformations in solution. Intersubunit salt bridging between Arg⁴ and Glu⁵³ during simulation stabilizes closed conformations. Because a semi-open state is among the low-energy conformations sampled in simulation, we propose that Cro-DNA binding may not entail a large conformational change relative to the dominant dimer forms in solution.     

Effects of PRE and POST therapy drug-pressure selected mutations on HIV-1 protease conformational sampling

Conformational sampling of pre- and post-therapy subtype B HIV-1 protease sequences derived from a pediatric subject infected via maternal transmission with HIV-1 were characterized by double electron–electron resonance spectroscopy. The conformational ensemble of the PRE construct resembles native-like inhibitor bound states. In contrast, the POST construct, which contains accumulated drug-pressure selected mutations, has a predominantly semi-open conformational ensemble, with increased populations of open-like states. The single point mutant L63P, which is contained in PRE and POST, has decreased dynamics, particularly in the flap region, and also displays a closed-like conformation of inhibitor-bound states. These findings support our hypothesis that secondary mutations accumulate in HIV-1 protease to shift conformational sampling to stabilize open-like conformations, while maintaining the predominant semi-open conformation for activity.     

Evidence of Functional Protein Dynamics from X-Ray Crystallographic Ensembles

It is widely recognized that representing a protein as a single static conformation is inadequate to describe the dynamics essential to the performance of its biological function. We contrast the amino acid displacements below and above the protein dynamical transition temperature, T_D∼215K, of hen egg white lysozyme using X-ray crystallography ensembles that are analyzed by molecular dynamics simulations as a function of temperature. We show that measuring structural variations across an ensemble of X-ray derived models captures the activation of conformational states that are of functional importance just above T_D, and they remain virtually identical to structural motions measured at 300K. Our results highlight the ability to observe functional structural variations across an ensemble of X-ray crystallographic data, and that residue fluctuations measured in MD simulations at room temperature are in quantitative agreement with the experimental observable.     

Highly Disordered Amyloid-β Monomer Probed by Single-Molecule FRET and MD Simulation

Monomers of amyloid-β (Aβ) protein are known to be disordered, but there is considerable controversy over the existence of residual or transient conformations that can potentially promote oligomerization and fibril formation. We employed single-molecule Förster resonance energy transfer (FRET) spectroscopy with site-specific dye labeling using an unnatural amino acid and molecular dynamics simulations to investigate conformations and dynamics of Aβ isoforms with 40 (Aβ40) and 42 residues (Aβ42). The FRET efficiency distributions of both proteins measured in phosphate-buffered saline at room temperature show a single peak with very similar FRET efficiencies, indicating there is apparently only one state. 2D FRET efficiency-donor lifetime analysis reveals, however, that there is a broad distribution of rapidly interconverting conformations. Using nanosecond fluorescence correlation spectroscopy, we measured the timescale of the fluctuations between these conformations to be ～35 ns, similar to that of disordered proteins. These results suggest that both Aβ40 and Aβ42 populate an ensemble of rapidly reconfiguring unfolded states, with no long-lived conformational state distinguishable from that of the disordered ensemble. To gain molecular-level insights into these observations, we performed molecular dynamics simulations with a force field optimized to describe disordered proteins. We find, as in experiments, that both peptides populate configurations consistent with random polymer chains, with the vast majority of conformations lacking significant secondary structure, giving rise to very similar ensemble-averaged FRET efficiencies.     

Phosphorylated Phospholamban Stabilizes a Compact Conformation of?the Cardiac Calcium-ATPase

The sarcoendoplasmic reticulum calcium ATPase (SERCA) plays a key role in cardiac calcium handling and is considered a high-value target for the treatment of heart failure. SERCA undergoes conformational changes as it harnesses the chemical energy of ATP for active transport. X-ray crystallography has provided insight into SERCA structural substates, but it is not known how well these static snapshots describe in vivo conformational dynamics. The goals of this work were to quantify the direction and magnitude of SERCA motions as the pump performs work in live cardiac myocytes, and to identify structural determinants of SERCA regulation by phospholamban. We measured intramolecular fluorescence resonance energy transfer (FRET) between fluorescent proteins fused to SERCA cytoplasmic domains. We detected four discrete structural substates for SERCA expressed in cardiac muscle cells. The relative populations of these discrete states oscillated with electrical pacing. Low FRET states were most populated in low Ca (diastole), and were indicative of an open, disordered structure for SERCA in the E2 (Ca-free) enzymatic substate. High FRET states increased with Ca (systole), suggesting rigidly closed conformations for the E1 (Ca-bound) enzymatic substates. Notably, a special compact E1 state was observed after treatment with β-adrenergic agonist or with coexpression of phosphomimetic mutants of phospholamban. The data suggest that SERCA calcium binding induces the pump to undergo a transition from an open, dynamic conformation to a closed, ordered structure. Phosphorylated phospholamban stabilizes a unique conformation of SERCA that is characterized by a compact architecture.     

Phosphorylated Phospholamban Stabilizes a Compact Conformation of the Cardiac Calcium-ATPase

The sarcoendoplasmic reticulum calcium ATPase (SERCA) plays a key role in cardiac calcium handling and is considered a high-value target for the treatment of heart failure. SERCA undergoes conformational changes as it harnesses the chemical energy of ATP for active transport. X-ray crystallography has provided insight into SERCA structural substates, but it is not known how well these static snapshots describe in vivo conformational dynamics. The goals of this work were to quantify the direction and magnitude of SERCA motions as the pump performs work in live cardiac myocytes, and to identify structural determinants of SERCA regulation by phospholamban. We measured intramolecular fluorescence resonance energy transfer (FRET) between fluorescent proteins fused to SERCA cytoplasmic domains. We detected four discrete structural substates for SERCA expressed in cardiac muscle cells. The relative populations of these discrete states oscillated with electrical pacing. Low FRET states were most populated in low Ca (diastole), and were indicative of an open, disordered structure for SERCA in the E2 (Ca-free) enzymatic substate. High FRET states increased with Ca (systole), suggesting rigidly closed conformations for the E1 (Ca-bound) enzymatic substates. Notably, a special compact E1 state was observed after treatment with β-adrenergic agonist or with coexpression of phosphomimetic mutants of phospholamban. The data suggest that SERCA calcium binding induces the pump to undergo a transition from an open, dynamic conformation to a closed, ordered structure. Phosphorylated phospholamban stabilizes a unique conformation of SERCA that is characterized by a compact architecture.     

DNP-Enhanced MAS NMR: A Tool to Snapshot Conformational Ensembles of α-Synuclein in Different States

Intrinsically disordered proteins dynamically sample a wide conformational space and therefore do not adopt a stable and defined three-dimensional conformation. The structural heterogeneity is related to their proper functioning in physiological processes. Knowledge of the conformational ensemble is crucial for a complete comprehension of this kind of proteins. We here present an approach that utilizes dynamic nuclear polarization-enhanced solid-state NMR spectroscopy of sparsely isotope-labeled proteins in frozen solution to take snapshots of the complete structural ensembles by exploiting the inhomogeneously broadened line-shapes. We investigated the intrinsically disordered protein α-synuclein (α-syn), which plays a key role in the etiology of Parkinson’s disease, in three different physiologically relevant states. For the free monomer in frozen solution we could see that the so-called “random coil conformation” consists of α-helical and β-sheet-like conformations, and that secondary chemical shifts of neighboring amino acids tend to be correlated, indicative of frequent formation of secondary structure elements. Based on these results, we could estimate the number of disordered regions in fibrillar α-syn as well as in α-syn bound to membranes in different protein-to-lipid ratios. Our approach thus provides quantitative information on the propensity to sample transient secondary structures in different functional states. Molecular dynamics simulations rationalize the results.     

Binding induced conformational changes of proteins correlate with their intrinsic fluctuations: a case study of antibodies

Background  

How antibodies recognize and bind to antigens can not be totally explained by rigid shape and electrostatic complimentarity models. Alternatively, pre-existing equilibrium hypothesis states that the native state of an antibody is not defined by a single rigid conformation but instead with an ensemble of similar conformations that co-exist at equilibrium. Antigens bind to one of the preferred conformations making this conformation more abundant shifting the equilibrium.     

On the problem of comparing protein structures: Development and applications of a new method for the assessment of structural similarities of polypeptide conformations

The development of prediction schemes and the search for evolutionary relationships amongst proteins require reliable methods for the comparison of polypeptide structures. It is shown that methods which attempt to describe structural similarities by a single value generally do not yield reliable estimates of the relatedness of two conformations. A new method is reported, called the D_k procedure, which yields a spectrum of deviations between two structures. Each particular D_k value is a measure of the similarity of the diagonals of the distance matrices of the compared conformations, k being the distance of the diagonals relative to the main diagonal.The method has the following features. (1) D_k is independent of chain length; (2) the method yields the relatedness of two conformations in terms of different structural levels; (3) D_k is a high-speed algorithm; (4) the D_k deviations of random structures from any particular conformation are predictable.The following applications are reported. (1) The success of both secondary and tertiary structure predictions are measured in terms of a reliable quality index. (2) The route of a conformation during simulation studies is followed on different structural levels, which exhibit the characteristics of the simulation. (3) The significance of hypotheses on protein folding subject to prediction schemes can be established. (4) A priori information (fixing pieces of secondary structure derived from X-ray investigations during prediction) is extractable from the predicted structures. (5) The evolutionary relatedness of two nucleotide binding proteins is established.The simplicity and speed of the D_k procedure allow its implementation even on minicomputers.     

DFTMD studies of β-cellobiose: conformational preference using implicit solvent

Previous DFT in vacuo studies on the conformational preferences for cellobiose showed that upon optimization the φ_H-anti conformations were of lower energy than the syn forms. Upon optimization using an implicit solvation method, COSMO, the syn or observed form was still not predicted to be of lower energy than the φ_H-anti form, even though optimization after addition of several explicit water molecules did show a relative energy difference favoring the syn form. In order to examine the predictive ability of COSMO on this carbohydrate, constant energy dynamics, DFTMD, simulations were carried out on low energy syn and φ_H-anti conformations with and without COSMO included during the dynamics. The resulting analysis confirmed that when COSMO is included in the dynamics, the syn conformations become energetically favored over the φ_H-anti forms suggesting that both solvent and entropy play roles in dictating the solution conformation of cellobiose. Analysis of the dynamic runs includes distributions of selected dihedral angles versus time, conformational transitions, and populations of some quasi-planar, boat, skew forms during the simulations.     

A conformation ensemble approach to protein residue-residue contact

Background

Protein residue-residue contact prediction is important for protein model generation and model evaluation. Here we develop a conformation ensemble approach to improve residue-residue contact prediction. We collect a number of structural models stemming from a variety of methods and implementations. The various models capture slightly different conformations and contain complementary information which can be pooled together to capture recurrent, and therefore more likely, residue-residue contacts.

Results

We applied our conformation ensemble approach to free modeling targets from both CASP8 and CASP9. Given a diverse ensemble of models, the method is able to achieve accuracies of. 48 for the top L/5 medium range contacts and. 36 for the top L/5 long range contacts for CASP8 targets (L being the target domain length). When applied to targets from CASP9, the accuracies of the top L/5 medium and long range contact predictions were. 34 and. 30 respectively.

Conclusions

When operating on a moderately diverse ensemble of models, the conformation ensemble approach is an effective means to identify medium and long range residue-residue contacts. An immediate benefit of the method is that when tied with a scoring scheme, it can be used to successfully rank models.     

Conformational equilibria and free energy profiles for the allosteric transition of the ribose-binding protein

The ribose-binding protein (RBP) is a sugar-binding bacterial periplasmic protein whose function is associated with a large allosteric conformational change from an open to a closed conformation upon binding to ribose. The crystal structures of RBP in open and closed conformations have been solved. It has been hypothesized that the open and closed conformations exist in a dynamic equilibrium in solution, and that sugar binding shifts the population from open conformations to closed conformations. Here, we study by computer simulations the thermodynamic changes that accompany this conformational change, and model the structural changes that accompany the allosteric transition, using umbrella sampling molecular dynamics and the weighted histogram analysis method. The open state is comprised of a diverse ensemble of conformations; the open ribose-free X-ray crystal conformations being representative of this ensemble. The unligated open form of RBP is stabilized by conformational entropy. The simulations predict detectable populations of closed ribose-free conformations in solution. Additional interdomain hydrogen bonds stabilize this state. The predicted shift in equilibrium from the open to the closed state on binding to ribose is in agreement with experiments. This is driven by the energetic stabilization of the closed conformation due to ribose-protein interactions. We also observe a significant population of a hitherto unobserved ribose-bound partially open state. We believe that this state is the one that has been suggested to play a role in the transfer of ribose to the membrane-bound permease complex.     

An integrative approach combining ion mobility mass spectrometry,X-ray crystallography,and nuclear magnetic resonance spectroscopy to study the conformational dynamics of α1-antitrypsin upon ligand binding

Native mass spectrometry (MS) methods permit the study of multiple protein species within solution equilibria, whereas ion mobility (IM)-MS can report on conformational behavior of specific states. We used IM-MS to study a conformationally labile protein (α₁-antitrypsin) that undergoes pathological polymerization in the context of point mutations. The folded, native state of the Z-variant remains highly polymerogenic in physiological conditions despite only minor thermodynamic destabilization relative to the wild-type variant. Various data implicate kinetic instability (conformational lability within a native state ensemble) as the basis of Z α₁-antitrypsin polymerogenicity. We show the ability of IM-MS to track such disease-relevant conformational behavior in detail by studying the effects of peptide binding on α₁-antitrypsin conformation and dynamics. IM-MS is, therefore, an ideal platform for the screening of compounds that result in therapeutically beneficial kinetic stabilization of native α₁-antitrypsin. Our findings are confirmed with high-resolution X-ray crystallographic and nuclear magnetic resonance spectroscopic studies of the same event, which together dissect structural changes from dynamic effects caused by peptide binding at a residue-specific level. IM-MS methods, therefore, have great potential for further study of biologically relevant thermodynamic and kinetic instability of proteins and provide rapid and multidimensional characterization of ligand interactions of therapeutic interest.PDB Code(s): 4PYW     

Analysis of the bacterial luciferase mobile loop by replica-exchange molecular dynamics

Bacterial luciferase contains an extended 29-residue mobile loop. Movements of this loop are governed by binding of either flavin mononucleotide (FMNH₂) or polyvalent anions. To understand this process, loop dynamics were investigated using replica-exchange molecular dynamics that yielded conformational ensembles in either the presence or absence of FMNH₂. The resulting data were analyzed using clustering and network analysis. We observed the closed conformations that are visited only in the simulations with the ligand. Yet the mobile loop is intrinsically flexible, and FMNH₂ binding modifies the relative populations of conformations. This model provides unique information regarding the function of a crystallographically disordered segment of the loop near the binding site. Structures at or near the fringe of this network were compatible with flavin binding or release. Finally, we demonstrate that the crystallographically observed conformation of the mobile loop bound to oxidized flavin was influenced by crystal packing. Thus, our study has revealed what we believe are novel conformations of the mobile loop and additional context for experimentally determined structures.     

Preferred conformations of a linear RGD tripeptide.

Conformational study of RGD tripeptides in the nonhydrated and hydrated states was carried out using an empirical potential function ECEPP/3 and the hydration shell model in order to investigate preferred conformations and factors responsible for their stability. RGD tripeptides in the nonhydrated and hydrated states can be interpreted as existing as an ensemble of feasible conformations rather than as a single dominant conformation from the analysis of distributions of backbone conformations, hydrogen bonds and beta-turns. The different distributions of conformations for the neutral and zwitterionic RGD tripeptides in both states may indicate that the conformation of the RGD tripeptide is liable to depend on solvent polarity and pH values. beta-Turn populations for the neutral tripeptide in both states are reasonably consistent with NMR measurements on linear RGD-containing peptides. The degradation of RGD tripeptide seems to be attributed mainly to the hydrogen bonds between the Asp side-chain and the backbone of Asp residue or C-terminal NHMe group, rather than to the flexible backbones of Gly and Asp residues.     

A flexible-protein molecular docking study of the binding of ruthenium complex compounds to PIM1, GSK-3β, and CDK2/Cyclin A protein kinases

We employ ensemble docking simulations to characterize the interactions of two enantiomeric forms of a Ru-complex compound (1-R and 1-S) with three protein kinases, namely PIM1, GSK-3β, and CDK2/cyclin A. We show that our ensemble docking computational protocol adequately models the structural features of these interactions and discriminates between competing conformational clusters of ligand-bound protein structures. Using the determined X-ray crystal structure of PIM1 complexed to the compound 1-R as a control, we discuss the importance of including the protein flexibility inherent in the ensemble docking protocol, for the accuracy of the structure prediction of the bound state. A comparison of our ensemble docking results suggests that PIM1 and GSK-3β bind the two enantiomers in similar fashion, through two primary binding modes: conformation I, which is very similar to the conformation presented in the existing PIM1/compound 1-R crystal structure; conformation II, which represents a 180° flip about an axis through the NH group of the pyridocarbazole moiety, relative to conformation I. In contrast, the binding of the enantiomers to CDK2 is found to have a different structural profile including a suggested bound conformation, which lacks the conserved hydrogen bond between the kinase and the ligand (i.e., ATP, staurosporine, Ru-complex compound). The top scoring conformation of the inhibitor bound to CDK2 is not present among the top-scoring conformations of the inhibitor bound to either PIM1 or GSK-3β and vice-versa. Collectively, our results help provide atomic-level insights into inhibitor selectivity among the three kinases.

Molecular Dynamics Study of the Opening Mechanism for DNA Polymerase I

During DNA replication, DNA polymerases follow an induced fit mechanism in order to rapidly distinguish between correct and incorrect dNTP substrates. The dynamics of this process are crucial to the overall effectiveness of catalysis. Although X-ray crystal structures of DNA polymerase I with substrate dNTPs have revealed key structural states along the catalytic pathway, solution fluorescence studies indicate that those key states are populated in the absence of substrate. Herein, we report the first atomistic simulations showing the conformational changes between the closed, open, and ajar conformations of DNA polymerase I in the binary (enzyme∶DNA) state to better understand its dynamics. We have applied long time-scale, unbiased molecular dynamics to investigate the opening process of the fingers domain in the absence of substrate for B. stearothermophilis DNA polymerase in silico. These simulations are biologically and/or physiologically relevant as they shed light on the transitions between states in this important enzyme. All closed and ajar simulations successfully transitioned into the fully open conformation, which is known to be the dominant binary enzyme-DNA conformation from solution and crystallographic studies. Furthermore, we have detailed the key stages in the opening process starting from the open and ajar crystal structures, including the observation of a previously unknown key intermediate structure. Four backbone dihedrals were identified as important during the opening process, and their movements provide insight into the recognition of dNTP substrate molecules by the polymerase binary state. In addition to revealing the opening mechanism, this study also demonstrates our ability to study biological events of DNA polymerase using current computational methods without biasing the dynamics.     

Assessing the effect of conformational averaging on the measured values of observables

Experiment and computer simulation are two complementary tools to understand the dynamics and behavior of biopolymers in solution. One particular area of interest is the ensemble of conformations populated by a particular molecule in solution. For example, what fraction of a protein sample exists in its folded conformation? How often does a particular peptide form an alpha helix versus a beta hairpin? To address these questions, it is important to determine the sensitivity of a particular experiment to changes in the distribution of molecular conformations. Consequently, a general analytic formalism is proposed to determine the sensitivity of a spectroscopic observable to the underlying distribution of conformations. A particular strength of the approach is that it provides an expression for a weighted average across conformational substates that is independent of the averaging function used. The formalism is described and applied to experimental and simulated nuclear Overhauser enhancement (NOE) and ³ J-coupling data on peptides in solution.