Comparative Molecular Dynamics Analysis of RNase-S Complex Formation.

Limited proteolysis of RNase-A yields a short N-terminal S-peptide segment and the larger S-protein. Binding of S-peptide to S-protein results in the formation of an enzymatically active RNase-S protein. S-peptide undergoes a transition from intrinsic disorder to an ordered helical state upon association with S-protein to form RNase-S and is an excellent model system to study coupled folding and binding. To better understand the dynamics of the RNases-S complex and its isolated partners, comparative molecular dynamics simulations have been performed. In agreement with experiment, we find significant conformational fluctuations of the isolated S-peptide compatible with a disordered regime and only little residual helical structure. In the RNase-S complex, the N-terminal helix of S-peptide unfolds and refolds repeatedly on the microsecond timescale, indicating that the α-helical structure is only part of the equilibrium regime for these residues whereas the C-terminal residues are confined to the helical conformation that is found in the x-ray structure. This is also in line with systematic, in silico Alanine scanning free-energy simulations, which indicate that the major contribution to complex stability emerges from the C-terminal helical turn, consisting of residues 8-13 in S-peptide whereas the N-terminal S-peptide residues 1-7 make only minor contributions. Comparative simulations of S-protein in the presence and absence of S-peptide reveal that the isolated S-protein is significantly more flexible than in the complex, and undergoes a global pincerlike conformational change that narrows the S-peptide binding cleft. The narrowed binding cleft adds a barrier for complex formation likely influencing the binding kinetics. This conformational change is reversed by S-peptide association, which also stabilizes conformational fluctuations in S-protein. Such global motions associated with binding are also likely to play a role for other coupled peptide folding and binding processes at peptide binding regions on protein surfaces.

Multiscale Simulation of Receptor-Drug Association Kinetics: Application to Neuraminidase Inhibitors.

A detailed understanding of the drug-receptor association process is of fundamental importance for drug design. Due to the long time scales of typical binding kinetics, the atomistic simulation of the ligand traveling from bulk solution into the binding site is still computationally challenging. In this work, we apply a multiscale approach of combined Molecular Dynamics (MD) and Brownian Dynamics (BD) simulations to investigate association pathway ensembles for the two prominent H1N1 neuraminidase inhibitors oseltamivir and zanamivir. Including knowledge of the approximate binding site location allows for the selective confinement of detailed but expensive MD simulations and application of less demanding BD simulations for the diffusion controlled part of the association pathway. We evaluate a binding criterion based on the residence time of the inhibitor in the binding pocket and compare it to geometric criteria that require prior knowledge about the binding mechanism. The method ranks the association rates of both inhibitors in qualitative agreement with experiment and yields reasonable absolute values depending, however, on the reaction criteria. The simulated association pathway ensembles reveal that, first, ligands are oriented in the electrostatic field of the receptor. Subsequently, a salt bridge is formed between the inhibitor's carboxyl group and neuraminidase residue Arg368, followed by adopting the native binding mode. Unexpectedly, despite oseltamivir's higher overall association rate, the rate into the intermediate salt-bridge state was found to be higher for zanamivir. The present methodology is intrinsically parallelizable and, although computationally demanding, allows systematic binding rate calculation on selected sets of potential drug molecules.

Covalent dye attachment influences the dynamics and conformational properties of flexible peptides.

Fluorescence spectroscopy techniques like Förster resonance energy transfer (FRET) and fluorescence correlation spectroscopy (FCS) have become important tools for the in vitro and in vivo investigation of conformational dynamics in biomolecules. These methods rely on the distance-dependent quenching of the fluorescence signal of a donor fluorophore either by a fluorescent acceptor fluorophore (FRET) or a non-fluorescent quencher, as used in FCS with photoinduced electron transfer (PET). The attachment of fluorophores to the molecule of interest can potentially alter the molecular properties and may affect the relevant conformational states and dynamics especially of flexible biomolecules like intrinsically disordered proteins (IDP). Using the intrinsically disordered S-peptide as a model system, we investigate the impact of terminal fluorescence labeling on the molecular properties. We perform extensive molecular dynamics simulations on the labeled and unlabeled peptide and compare the results with in vitro PET-FCS measurements. Experimental and simulated timescales of end-to-end fluctuations were found in excellent agreement. Comparison between simulations with and without labels reveal that the π-stacking interaction between the fluorophore labels traps the conformation of S-peptide in a single dominant state, while the unlabeled peptide undergoes continuous conformational rearrangements. Furthermore, we find that the open to closed transition rate of S-peptide is decreased by at least one order of magnitude by the fluorophore attachment. Our approach combining experimental and in silico methods provides a benchmark for the simulations and reveals the significant effect that fluorescence labeling can have on the conformational dynamics of small biomolecules, at least for inherently flexible short peptides. The presented protocol is not only useful for comparing PET-FCS experiments with simulation results but provides a strategy to minimize the influence on molecular properties when chosing labeling positions for fluorescence experiments.

Molecular Dynamics Analysis of 4E-BP2 Protein Fold Stabilization Induced by Phosphorylation.

Protein phosphorylation can affect the interaction with partner proteins but can also induce conformational transitions. In case of the eukaryotic translation initiation factor 4E-binding protein 2 (4E-BP2) threonine (Thr) phosphorylation at two turn motifs results in transition from a disordered to a folded structure. In order to elucidate the stabilizing mechanism we employed comparative molecular dynamics (MD) free energy simulations on the turn motifs indicating that Thr-phosphorylation favors a folded whereas dephosphorylation or substitution by Glu residues destabilizes the turn structure. In multiple unrestrained MD simulations at elevated temperature of the 4E-BP2 domain only the double phosphorylated variant remained close to the folded structure in agreement with experiment. Three surface Arg residues were identified as additional key elements for the tertiary structure stabilization of the whole phosphorylated domain. In addition to the local turn structure double phosphorylation also leads to an overall electrostatic stabilization of the folded form compared to wild type and other investigated variants of 4E-BP2. The principles of phosphorylation mediated fold stabilization identified in the present study may also be helpful for identifying other structural motifs that can be affected by phosphorylation or provide a route to design such motifs.

Mechanism of pKID/KIX Association Studied by Molecular Dynamics Free Energy Simulations.

The phosphorylated kinase-inducible domain (pKID) associates with the kinase interacting domain (KIX) via a coupled folding and binding mechanism. The pKID domain is intrinsically disordered when unbound and upon phosphorylation at Ser133 binds to the KIX domain adopting a well-defined kinked two-helix structure. In order to identify putative hot spot residues of binding that could serve as an initial stable anchor, we performed in silico alanine scanning free energy simulations. The simulations indicate that charged residues including the phosphorylated central Ser133 of pKID make significant contributions to binding. However, these are of slightly smaller magnitude compared to several hydrophobic side chains not defining a single dominant binding hot spot. Both continuous molecular dynamics (MD) simulations and free energy analysis demonstrate that phosphorylation significantly stabilizes the central kinked motif around Ser133 of pKID and shifts the conformational equilibrium toward the bound conformation already in the absence of KIX. This result supports a view that pKID/KIX association follows in part a conformational selection process. During a 1.5 µs explicit solvent MD simulation, folding of pKID on the surface of KIX was observed after an initial contact at the bound position of the phosphorylation site was enforced following a sequential process of αA helix association and a stepwise association and folding of the second αB helix compatible with available experimental results.

Adenylylation of Tyr77 stabilizes Rab1b GTPase in an active state: A molecular dynamics simulation analysis.

The pathogenic pathway of Legionella pneumophila exploits the intercellular vesicle transport system via the posttranslational attachment of adenosine monophosphate (AMP) to the Tyr77 sidechain of human Ras like GTPase Rab1b. The modification, termed adenylylation, is performed by the bacterial enzyme DrrA/SidM, however the effect on conformational properties of the molecular switch mechanism of Rab1b remained unresolved. In this study we find that the adenylylation of Tyr77 stabilizes the active Rab1b state by locking the switch in the active signaling conformation independent of bound GTP or GDP and that electrostatic interactions due to the additional negative charge in the switch region make significant contributions. The stacking interaction between adenine and Phe45 however, seems to have only minor influence on this stabilisation. The results may also have implications for the mechanistic understanding of conformational switching in other signaling proteins.

Protein-ligand docking using hamiltonian replica exchange simulations with soft core potentials.

Molecular dynamics (MD) simulations in explicit solvent allow studying receptor-ligand binding processes including full flexibility of the binding partners and an explicit inclusion of solvation effects. However, in MD simulations, the search for an optimal ligand-receptor complex geometry is frequently trapped in locally stable non-native binding geometries. A Hamiltonian replica-exchange (H-REMD)-based protocol has been designed to enhance the sampling of putative ligand-receptor complexes. It is based on softening nonbonded ligand-receptor interactions along the replicas and one reference replica under the control of the original force field. The efficiency of the method has been evaluated on two receptor-ligand systems and one protein-peptide complex. Starting from misplaced initial docking geometries, the H-REMD method reached in each case the known binding geometry significantly faster than a standard MD simulation. The approach could also be useful to identify and evaluate alternative binding geometries in a given binding region with small relative differences in binding free energy.

Role of tyrosine hot-spot residues at the interface of colicin E9 and immunity protein 9: a comparative free energy simulation study.

The endonuclease activity of the bacterial colicin 9 enzyme is controlled by the specific and high-affinity binding of immunity protein 9 (Im9). Molecular dynamics simulation studies in explicit solvent were used to investigate the free energy change associated with the mutation of two hot-spot interface residues [tyrosine (Tyr): Tyr54 and Tyr55] of Im9 to Ala. In addition, the effect of several other mutations (Leu33Ala, Leu52Ala, Val34Ala, Val37Ala, Ser48Ala, and Ile53Ala) with smaller influence on binding affinity was also studied. Good qualitative agreement of calculated free energy changes and experimental data on binding affinity of the mutations was observed. The simulation studies can help to elucidate the molecular details on how the mutations influence protein-protein binding affinity. The role of solvent and conformational flexibility of the partner proteins was studied by comparing the results in the presence or absence of solvent and with or without positional restraints. Restriction of the conformational mobility of protein partners resulted in significant changes of the calculated free energies but of similar magnitude for isolated Im9 and for the complex and therefore in only modest changes of binding free energy differences. Although the overall binding free energy change was similar for the two Tyr-Ala mutations, the physical origin appeared to be different with solvation changes contributing significantly to the Tyr55Ala mutation and to a loss of direct protein-protein interactions dominating the free energy change due to the Tyr54Ala mutation.