Solvent fluctuations in hydrophobic cavity-ligand binding kinetics.

Water plays a crucial part in virtually all protein-ligand binding processes in and out of equilibrium. Here, we investigate the role of water in the binding kinetics of a ligand to a prototypical hydrophobic pocket by explicit-water molecular dynamics (MD) simulations and implicit diffusional approaches. The concave pocket in the unbound state exhibits wet/dry hydration oscillations whose magnitude and time scale are significantly amplified by the approaching ligand. In turn, the ligand's stochastic motion intimately couples to the slow hydration fluctuations, leading to a sixfold-enhanced friction in the vicinity of the pocket entrance. The increased friction considerably decelerates association in the otherwise barrierless system, indicating the importance of molecular-scale hydrodynamic effects in cavity-ligand binding arising due to capillary fluctuations. We derive and analyze the diffusivity profile and show that the mean first passage time distribution from the MD simulation can be accurately reproduced by a standard Brownian dynamics simulation if the appropriate position-dependent friction profile is included. However, long-time decays in the water-ligand (random) force autocorrelation demonstrate violation of the Markovian assumption, challenging standard diffusive approaches for rate prediction. Remarkably, the static friction profile derived from the force correlations strongly resembles the profile derived on the Markovian assumption apart from a simple shift in space, which can be rationalized by a time-space retardation in the ligand's downhill dynamics toward the pocket. The observed spatiotemporal hydrodynamic coupling may be of biological importance providing the time needed for conformational receptor-ligand adjustments, typical of the induced-fit paradigm.

LSD1/CoREST is an allosteric nanoscale clamp regulated by H3-histone-tail molecular recognition.

The complex of lysine-specific demethylase-1 (LSD1/KDM1A) with its corepressor protein CoREST is an exceptionally relevant target for epigenetic drugs. Here, we provide insight into the local and global changes of LSD1/CoREST conformational dynamics that occur upon H3 binding on the basis of a total cumulative time of one microsecond molecular dynamics simulation. The LSD1/CoREST complex functions as an allosteric nanoscale-binding clamp, which is regulated by substrate binding. In the unbound state, LSD1/CoREST reversibly visits clamp states that are more open or significantly more closed compared with the available X-ray crystal structures. The Lys triad of residues Lys355, Lys357, and Lys359 gates the entrance of the H3 pocket. H3 binding shifts the pocket breathing dynamics toward open, higher-volume states while reducing the overall flexibility of the LSD1/CoREST nanoscale clamp. We show that the H3 pocket is an allosteric site for the regulation of the rotation of the amino oxidase domain with respect to the Tower domain. The allosteric mechanism relies on the specific reduction of nanoscale domain rotation upon local H3-tail binding. Instead, clamp opening/closing motions that do not involve domain rotation only reduce in amplitude yet are dominant in the bound state. Overall, our data suggest that the H3 binding pocket is a central target site to (i) switch off LSD1 amino oxidase activity, thus H3-tail demethylation; (ii) block the competitive binding of transcription factors; and (iii) prevent chromatin anchoring to LSD1/CoREST. This study underscores the importance of receptor flexibility for future epigenetic drug discovery.

Molecular recognition and ligand association.

We review recent developments in our understanding of molecular recognition and ligand association, focusing on two major viewpoints: (a) studies that highlight new physical insight into the molecular recognition process and the driving forces determining thermodynamic signatures of binding and (b) recent methodological advances in applications to protein-ligand binding. In particular, we highlight the challenges posed by compensating enthalpic and entropic terms, competing solute and solvent contributions, and the relevance of complex configurational ensembles comprising multiple protein, ligand, and solvent intermediate states. As more complete physics is taken into account, computational approaches increase their ability to complement experimental measurements, by providing a microscopic, dynamic view of ensemble-averaged experimental observables. Physics-based approaches are increasingly expanding their power in pharmacology applications.

Expanding the druggable space of the LSD1/CoREST epigenetic target: new potential binding regions for drug-like molecules, peptides, protein partners, and chromatin.

Lysine specific demethylase-1 (LSD1/KDM1A) in complex with its corepressor protein CoREST is a promising target for epigenetic drugs. No therapeutic that targets LSD1/CoREST, however, has been reported to date. Recently, extended molecular dynamics (MD) simulations indicated that LSD1/CoREST nanoscale clamp dynamics is regulated by substrate binding and highlighted key hinge points of this large-scale motion as well as the relevance of local residue dynamics. Prompted by the urgent need for new molecular probes and inhibitors to understand LSD1/CoREST interactions with small-molecules, peptides, protein partners, and chromatin, we undertake here a configurational ensemble approach to expand LSD1/CoREST druggability. The independent algorithms FTMap and SiteMap and our newly developed Druggable Site Visualizer (DSV) software tool were used to predict and inspect favorable binding sites. We find that the hinge points revealed by MD simulations at the SANT2/Tower interface, at the SWIRM/AOD interface, and at the AOD/Tower interface are new targets for the discovery of molecular probes to block association of LSD1/CoREST with chromatin or protein partners. A fourth region was also predicted from simulated configurational ensembles and was experimentally validated to have strong binding propensity. The observation that this prediction would be prevented when using only the X-ray structures available (including the X-ray structure bound to the same peptide) underscores the relevance of protein dynamics in protein interactions. A fifth region was highlighted corresponding to a small pocket on the AOD domain. This study sets the basis for future virtual screening campaigns targeting the five novel regions reported herein and for the design of LSD1/CoREST mutants to probe LSD1/CoREST binding with chromatin and various protein partners.

Substrate selection influences molecular recognition in a screen for lymphoid tyrosine phosphatase inhibitors.

Assay design is an important variable that influences the outcome of an inhibitor screen. Here, we have investigated the hypothesis that protein tyrosine phosphatase inhibitors with improved biological activity could be identified from a screen by using a biologically relevant peptide substrate, rather than traditional phosphotyrosine mimetic substrates. A 2000-member library of drugs and drug-like compounds was screened for inhibitors of lymphoid tyrosine phosphatase (LYP) by using both a peptide substrate (Ac-ARLIEDNE-pCAP-TAREG-NH2, peptide 1) and a small-molecule phosphotyrosine mimetic substrate (difluoromethyl umbelliferyl phosphate, DiFMUP). The results demonstrate that compounds that inhibited enzyme activity on the peptide substrate had greater biological activity than compounds that only inhibited enzyme activity on DiFMUP. Finally, epigallocatechin-3,5-digallate was identified as the most potent inhibitor of lymphoid tyrosine phosphatase activity to date, with an IC50 of 50 nM and significant activity in T-cells. Molecular docking simulations provided a first model for binding of this potent inhibitor to LYP; this will constitute the platform for ongoing lead optimization efforts.

LSD1/CoREST reversible opening-closing dynamics: discovery of a nanoscale clamp for chromatin and protein binding.

LSD1 associated with its corepressor protein CoREST is an exceptionally relevant target for epigenetic drugs. Hypotheses for the role of LSD1/CoREST as a multidocking site for chromatin and protein binding would require significant molecular flexibility, and LSD1/CoREST large-amplitude conformational dynamics is currently unknown. Here, molecular dynamics simulation reveals that the LSD1/CoREST complex in solution functions as a reversible nanoscale binding clamp. We show that the H3 histone tail binding pocket is a potential allosteric site for regulation of the rotation of SWIRM/SANT2 domains around the Tower domain. Thus, targeting this site and including receptor flexibility are crucial strategies for future drug discovery.

Multiple pathways guide oxygen diffusion into flavoenzyme active sites.

Dioxygen (O(2)) and other gas molecules have a fundamental role in a variety of enzymatic reactions. However, it is only poorly understood which O(2) uptake mechanism enzymes employ to promote efficient catalysis and how general this is. We investigated O(2) diffusion pathways into monooxygenase and oxidase flavoenzymes, using an integrated computational and experimental approach. Enhanced-statistics molecular dynamics simulations reveal spontaneous protein-guided O(2) diffusion from the bulk solvent to preorganized protein cavities. The predicted protein-guided diffusion paths and the importance of key cavity residues for oxygen diffusion were verified by combining site-directed mutagenesis, rapid kinetics experiments, and high-resolution X-ray structures. This study indicates that monooxygenase and oxidase flavoenzymes employ multiple funnel-shaped diffusion pathways to absorb O(2) from the solvent and direct it to the reacting C4a atom of the flavin cofactor. The difference in O(2) reactivity among dehydrogenases, monooxygenases, and oxidases ultimately resides in the fine modulation of the local environment embedding the reactive locus of the flavin.

Predictive power of molecular dynamics receptor structures in virtual screening.

Molecular dynamics (MD) simulation is a well-established method for understanding protein dynamics. Conformations from unrestrained MD simulations have yet to be assessed for blind virtual screening (VS) by docking. This study presents a critical analysis of the predictive power of MD snapshots to this regard, evaluating two well-characterized systems of varying flexibility in ligand-bound and unbound configurations. Results from such VS predictions are discussed with respect to experimentally determined structures. In all cases, MD simulations provide snapshots that improve VS predictive power over known crystal structures, possibly due to sampling more relevant receptor conformations. Additionally, MD can move conformations previously not amenable to docking into the predictive range.

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

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

Water in cavity-ligand recognition.

We use explicit solvent molecular dynamics simulations to estimate free energy, enthalpy, and entropy changes along the cavity-ligand association coordinate for a set of seven model systems with varying physicochemical properties. Owing to the simplicity of the considered systems we can directly investigate the role of water thermodynamics in molecular recognition. A broad range of thermodynamic signatures is found in which water (rather than cavity or ligand) enthalpic or entropic contributions appear to drive cavity-ligand binding or rejection. The unprecedented, nanoscale picture of hydration thermodynamics can help the interpretation and design of protein-ligand binding experiments. Our study opens appealing perspectives to tackle the challenge of solvent entropy estimation in complex systems and for improving molecular simulation models.

Water-membrane partition thermodynamics of an amphiphilic lipopeptide: an enthalpy-driven hydrophobic effect.

To shed light on the driving force for the hydrophobic effect that partitions amphiphilic lipoproteins between water and membrane, we carried out an atomically detailed thermodynamic analysis of a triply lipid modified H-ras heptapeptide anchor (ANCH) in water and in a DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine) bilayer. Combining molecular mechanical and continuum solvent approaches with an improved technique for solute entropy calculation, we obtained an overall transfer free energy of approximately -13 kcal mol(-1). This value is in qualitative agreement with free energy changes derived from a potential of mean force calculation and indirect experimental observations. Changes in free energies of solvation and ANCH conformational reorganization are unfavorable, whereas ANCH-DMPC interactions-especially van der Waals-favor insertion. These results are consistent with an enthalpy-driven hydrophobic effect, in accord with earlier calorimetric data on the membrane partition of other amphiphiles. Furthermore, structural and entropic analysis of molecular dynamics-generated ensembles suggests that conformational selection may play a hitherto unappreciated role in membrane insertion of lipid-modified peptides and proteins.

Exploring the conserved water site and hydration of a coiled-coil trimerisation motif: a MD simulation study.

The solvent structure and dynamics around ccbeta-p, a 17-residue peptide that forms a parallel three-stranded alpha-helical coiled coil in solution, was analysed through 10 ns explicit solvent molecular dynamics (MD) simulations at 278 and 330 K. Comparison with two corresponding simulations of the monomeric form of ccbeta-p was used to investigate the changes of hydration upon coiled-coil formation. Pronounced peaks in the solvent density distribution between residues Arg8 and Glu13 of neighbouring helices show the presence of water bridges between the helices of the ccbeta-p trimer; this is in agreement with the water sites observed in X-ray crystallography experiments. Interestingly, this water site is structurally conserved in many three-stranded coiled coils and, together with the Arg and Glu residues, forms part of a motif that determines three-stranded coiled-coil formation. Our findings show that little direct correlation exists between the solvent density distribution and the temporal ordering of water around the trimeric coiled coil. The MD-calculated effective residence times of up to 40 ps show rapid exchange of surface water molecules with the bulk phase, and indicate that the solvent distribution around biomolecules requires interpretation in terms of continuous density distributions rather than in terms of discrete molecules of water. Together, our study contributes to understanding the principles of three-stranded coiled-coil formation.

(Thermo)dynamic role of receptor flexibility, entropy, and motional correlation in protein-ligand binding.

The binding of 2-amino-5-methylthiazole to the W191G cavity mutant of cytochrome c peroxidase is an ideal test case to investigate the entropic contribution to the binding free energy due to changes in receptor flexibility. The dynamic and thermodynamic role of receptor flexibility are studied by 50 ns-long explicit-solvent molecular dynamics simulations of three separate receptor ensembles: W191G binding a K(+) ion, W191G-2a5mt complex with a closed 190-195 gating loop, and apo with an open loop. We employ a method recently proposed to estimate accurate absolute single-molecule configurational entropies and their differences for systems undergoing conformational transitions. We find that receptor flexibility plays a generally underestimated role in protein-ligand binding (thermo)dynamics and that changes of receptor motional correlation determine such large entropy contributions.

E9-Im9 colicin DNase-immunity protein biomolecular association in water: a multiple-copy and accelerated molecular dynamics simulation study.

Protein-protein transient and dynamic interactions underlie all biological processes. The molecular dynamics (MD) of the E9 colicin DNase protein, its Im9 inhibitor protein, and their E9-Im9 recognition complex are investigated by combining multiple-copy (MC) MD and accelerated MD (aMD) explicit-solvent simulation approaches, after validation with crystalline-phase and solution experiments. Im9 shows higher flexibility than its E9 counterpart. Im9 displays a significant reduction of backbone flexibility and a remarkable increase in motional correlation upon E9 association. Im9 loops 23-31 and 54-64 open with respect to the E9-Im9 X-ray structure and show high conformational diversity. Upon association a large fraction (approximately 20 nm2) of E9 and Im9 protein surfaces become inaccessible to water. Numerous salt bridges transiently occurring throughout our six 50 ns long MC-MD simulations are not present in the X-ray model. Among these Im9 Glu31-E9 Arg96 and Im9 Glu41-Lys89 involve interface interactions. Through the use of 10 ns of Im9 aMD simulation, we reconcile the largest thermodynamic impact measured for Asp51Ala mutation with Im9 structure and dynamics. Lys57 acts as an essential molecular switch to shift Im9 surface loop towards an ideal configuration for E9 inhibition. This is achieved by switching Asp60-Lys57 and Asp62-Lys57 hydrogen bonds to Asp51-Lys57 salt bridge. E9-Im9 recognition involves shifts of conformational distributions, reorganization of intramolecular hydrogen bond patterns, and formation of new inter- and intramolecular interactions. The description of key transient biological interactions can be significantly enriched by the dynamic and atomic-level information provided by computer simulations.

Configurational entropy elucidates the role of salt-bridge networks in protein thermostability.

Detailed knowledge of how networks of surface salt bridges contribute to protein thermal stability is essential not only to understand protein structure and function but also to design thermostable proteins for industrial applications. Experimental studies investigating thermodynamic stability through measurements of free energy associated with mutational alterations in proteins provide only macroscopic evidence regarding the structure of salt-bridge networks and assessment of their contribution to protein stability. Using explicit-solvent molecular dynamics simulations to provide insight on the atomic scale, we investigate here the structural stability, defined in terms of root-mean-square fluctuations, of a short polypeptide designed to fold into a stable trimeric coiled coil with a well-packed hydrophobic core and an optimal number of intra- and interhelical surface salt bridges. We find that the increase of configurational entropy of the backbone and side-chain atoms and decreased pair correlations of these with increased temperature are consistent with nearly constant atom-positional root-mean-square fluctuations, increased salt-bridge occupancies, and stronger electrostatic interactions in the coiled coil. Thus, our study of the coiled coil suggests a mechanism in which well-designed salt-bridge networks could accommodate stochastically the disorder of increased thermal motion to produce thermostability.

Biomolecular modeling: Goals, problems, perspectives.

Computation based on molecular models is playing an increasingly important role in biology, biological chemistry, and biophysics. Since only a very limited number of properties of biomolecular systems is actually accessible to measurement by experimental means, computer simulation can complement experiment by providing not only averages, but also distributions and time series of any definable quantity, for example, conformational distributions or interactions between parts of systems. Present day biomolecular modeling is limited in its application by four main problems: 1) the force-field problem, 2) the search (sampling) problem, 3) the ensemble (sampling) problem, and 4) the experimental problem. These four problems are discussed and illustrated by practical examples. Perspectives are also outlined for pushing forward the limitations of biomolecular modeling.

Comparison of atomic-level and coarse-grained models for liquid hydrocarbons from molecular dynamics configurational entropy estimates.

Molecular liquids can be modeled at different levels of spatial resolution. In atomic-level (AL) models, all (heavy) atoms can be explicitly simulated. In coarse-grained (CG) models, particles (beads) that represent groups of covalently bound atoms are used as elementary units. Ideally, a CG model should reproduce the thermodynamic and structural properties of the corresponding AL model after mapping to the lower-resolution scale. In the present work, two such models are investigated: (i) the classical GROMOS atomic-level model; (ii) a CG model recently proposed by Marrink et al., which maps approximately four non-hydrogen atoms to one bead [J. Phys. Chem. B 2004, 108, 750]. The study is restricted to n-alkanes whose aliphatic fragments are abundantly found in lipids of biological interest. Additionally, cis-9-octadecene is included, as a template chain of the lipid dioleoylphosphatidylcholine (DOPC). The two representations of molecules in the liquid phase are compared in terms of average molecular structures, extent of configurational space sampled, and single-molecule entropies. An approximate method is used to estimate the rotational contributions to the absolute configurational entropy. Good correspondence between the AL and CG representations is found. The loss in configurational entropy due to the reduction in degrees of freedom upon coarse-graining of the model is estimated.

Configurational entropies of lipids in pure and mixed bilayers from atomic-level and coarse-grained molecular dynamics simulations.

Single-chain and single-fragment configurational entropies of lipid tails in hydrated lipid bilayers are evaluated from molecular dynamics simulations using the quasi-harmonic approximation. The entropy distribution along individual acyl tails is obtained and compared to that of corresponding hydrocarbon chains in the liquid phase. We consider pure dipalmitoylphosphatidylcholine and mixed dioleoylphosphatidylcholine/dioleoylphosphatidylethanolamine bilayers. The systems are modeled at different levels of spatial resolution: In an atomic-level (AL) model all (heavy) atoms are explicitly simulated; in a coarse-grained (CG) model particles (beads) representing groups of covalently bound atoms are used, which map approximately four non-hydrogen atoms to one interaction site. Single-chain and single-fragment entropies and correlations between the motions of (single) acyl chains are compared. A good correspondence is found between the flexibility of the AL and CG models. The loss in configurational entropy due to the reduction in the number of degrees of freedom upon coarse-graining of the model is estimated. The CG model shows about 4 times faster convergence of the chain entropies than the more detailed AL model. Corrections to the quasi-harmonic entropy estimates were found to be small for the CG model. For the AL model, the correction due to mode anharmonicities is small, but the correction due to pairwise (supralinear) mode correlations is sizable.

A marine analgesic peptide, Contulakin-G, and neurotensin are distinct agonists for neurotensin receptors: uncovering structural determinants of desensitization properties.

Neurotensin receptors have been studied as molecular targets for the treatment of pain, schizophrenia, addiction, or cancer. Neurotensin (NT) and Contulakin-G, a glycopeptide isolated from a predatory cone snail Conus geographus, share a sequence similarity at the C-terminus, which is critical for activation of neurotensin receptors. Both peptides are potent analgesics, although affinity and agonist potency of Contulakin-G toward neurotensin receptors are significantly lower, as compared to those for NT. In this work, we show that the weaker agonist properties of Contulakin-G result in inducing significantly less desensitization of neurotensin receptors and preserving their cell-surface density. Structure-activity relationship (SAR) studies suggested that both glycosylation and charged amino acid residues in Contulakin-G or NT played important roles in desensitizing neurotensin receptors. Computational modeling studies of human neurotensin receptor NTS1 and Contulakin-G confirmed the role of glycosylation in weakening interactions with the receptors. Based on available SAR data, we designed, synthesized, and characterized an analog of Contulakin-G in which the glycosylated amino acid residue, Gal-GalNAc-Thr10, was replaced by memantine-Glu10 residue. This analog exhibited comparable agonist potency and weaker desensitization properties as compared to that of Contulakin-G, while producing analgesia in the animal model of acute pain following systemic administration. We discuss our study in the context of feasibility and safety of developing NT therapeutic agents with improved penetration across the blood-brain barrier. Our work supports engineering peptide-based agonists with diverse abilities to desensitize G-protein coupled receptors and further emphasizes opportunities for conotoxins as novel pharmacological tools and drug candidates.