Ligand unbinding mechanisms and kinetics for T4 lysozyme mutants from τRAMD simulations.
Ontology highlight
ABSTRACT:
SUBMITTER: Nunes-Alves A
PROVIDER: S-EPMC8244441 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Drug discovery is expensive and high-risk. Its main reasons of failure are lack of efficacy and toxicity of a drug candidate. Binding affinity for the biological target has been usually considered one of the most relevant figures of merit to judge a drug candidate along with bioavailability, selectivity and metabolic properties, which could depend on off-target interactions. Nevertheless, affinity does not always satisfactorily correlate with in vivo drug efficacy. It is indeed becoming increasingly evident that the time a drug spends in contact with its target (aka residence time) can be a more reliable figure of merit. Experimental kinetic measurements are operatively limited by the cost and the time needed to synthesize compounds to be tested, to express and purify the target, and to setup the assays. We present here a simple and efficient molecular-dynamics-based computational approach to prioritize compounds according to their residence time. We devised a multiple-replica scaled molecular dynamics protocol with suitably defined harmonic restraints to accelerate the unbinding events while preserving the native fold. Ligands are ranked according to the mean observed scaled unbinding time. The approach, trivially parallel and easily implementable, was validated against experimental information available on biological systems of pharmacological relevance.
| S-EPMC4477625 | biostudies-literature
Project description:The studies presented here explore the relationship between protein packing and molecular flexibility using ligand-binding cavity mutants of T4 lysozyme. Although previously reported crystal structures of the mutants investigated show single conformations that are similar to the WT protein, site-directed spin labeling in solution reveals additional conformational substates in equilibrium exchange with a WT-like population. Remarkably, binding of ligands, including the general anesthetic halothane shifts the population to the WT-like state, consistent with a conformational selection model of ligand binding, but structural adaptation to the ligand is also apparent in one mutant. Distance mapping with double electron-electron resonance spectroscopy and the absence of ligand binding suggest that the new substates induced by the cavity-creating mutations represent alternate packing modes in which the protein fills or partially fills the cavity with side chains, including the spin label in one case; external ligands compete with the side chains for the cavity space, stabilizing the WT conformation. The results have implications for mechanisms of anesthesia, the response of proteins to hydrostatic pressure, and protein engineering.
| S-EPMC3832024 | biostudies-literature
Project description:The residence time (RT), the time for which a drug remains bound to its biological target, is a critical parameter for drug design. The prediction of this key kinetic property has been proven to be challenging and computationally demanding in the framework of atomistic simulations. In the present work, we setup and applied two distinct metadynamics protocols to estimate the RTs of muscarinic M3 receptor antagonists. In the first method, derived from the conformational flooding approach, the kinetics of unbinding is retrieved from a physics-based parameter known as the acceleration factor α (i.e., the running average over time of the potential deposited in the bound state). Such an approach is expected to recover the absolute RT value for a compound of interest. In the second method, known as the tMETA-D approach, a qualitative estimation of the RT is given by the time of simulation required to drive the ligand from the binding site to the solvent bulk. This approach has been developed to reproduce the change of experimental RTs for compounds targeting the same target. Our analysis shows that both computational protocols are able to rank compounds in agreement with their experimental RTs. Quantitative structure-kinetics relationship (SKR) models can be identified and employed to predict the impact of a chemical modification on the experimental RT once a calibration study has been performed.
| S-EPMC10170513 | biostudies-literature
Project description:Application of hydrostatic pressure shifts protein conformational equilibria in a direction to reduce the volume of the system. A current view is that the volume reduction is dominated by elimination of voids or cavities in the protein interior via cavity hydration, although an alternative mechanism wherein cavities are filled with protein side chains resulting from a structure relaxation has been suggested [López CJ, Yang Z, Altenbach C, Hubbell WL (2013) Proc Natl Acad Sci USA 110(46):E4306-E4315]. In the present study, mechanisms for elimination of cavities under high pressure are investigated in the L99A cavity mutant of T4 lysozyme and derivatives thereof using site-directed spin labeling, pressure-resolved double electron-electron resonance, and high-pressure circular dichroism spectroscopy. In the L99A mutant, the ground state is in equilibrium with an excited state of only ∼ 3% of the population in which the cavity is filled by a protein side chain [Bouvignies et al. (2011) Nature 477(7362):111-114]. The results of the present study show that in L99A the native ground state is the dominant conformation to pressures of 3 kbar, with cavity hydration apparently taking place in the range of 2-3 kbar. However, in the presence of additional mutations that lower the free energy of the excited state, pressure strongly populates the excited state, thereby eliminating the cavity with a native side chain rather than solvent. Thus, both cavity hydration and structure relaxation are mechanisms for cavity elimination under pressure, and which is dominant is determined by details of the energy landscape.
| S-EPMC4434698 | biostudies-literature
Project description:The ability to predict the mechanisms and the associated rate constants of protein-ligand unbinding is of great practical importance in drug design. In this work we demonstrate how a recently introduced metadynamics-based approach allows exploration of the unbinding pathways, estimation of the rates, and determination of the rate-limiting steps in the paradigmatic case of the trypsin-benzamidine system. Protein, ligand, and solvent are described with full atomic resolution. Using metadynamics, multiple unbinding trajectories that start with the ligand in the crystallographic binding pose and end with the ligand in the fully solvated state are generated. The unbinding rate k off is computed from the mean residence time of the ligand. Using our previously computed binding affinity we also obtain the binding rate k on. Both rates are in agreement with reported experimental values. We uncover the complex pathways of unbinding trajectories and describe the critical rate-limiting steps with unprecedented detail. Our findings illuminate the role played by the coupling between subtle protein backbone fluctuations and the solvation by water molecules that enter the binding pocket and assist in the breaking of the shielded hydrogen bonds. We expect our approach to be useful in calculating rates for general protein-ligand systems and a valid support for drug design.
| S-EPMC4321287 | biostudies-literature
Project description:A key factor influencing a drug's efficacy is its residence time in the binding pocket of the host protein. Using atomistic computer simulation to predict this residence time and the associated dissociation process is a desirable but extremely difficult task due to the long timescales involved. This gets further complicated by the presence of biophysical factors such as steric and solvation effects. In this work, we perform molecular dynamics (MD) simulations of the unbinding of a popular prototypical hydrophobic cavity-ligand system using a metadynamics-based approach that allows direct assessment of kinetic pathways and parameters. When constrained to move in an axial manner, the unbinding time is found to be on the order of 4,000 s. In accordance with previous studies, we find that the cavity must pass through a region of sharp wetting transition manifested by sudden and high fluctuations in solvent density. When we remove the steric constraints on ligand, the unbinding happens predominantly by an alternate pathway, where the unbinding becomes 20 times faster, and the sharp wetting transition instead becomes continuous. We validate the unbinding timescales from metadynamics through a Poisson analysis, and by comparison through detailed balance to binding timescale estimates from unbiased MD. This work demonstrates that enhanced sampling can be used to perform explicit solvent MD studies at timescales previously unattainable, to our knowledge, obtaining direct and reliable pictures of the underlying physiochemical factors including free energies and rate constants.
| S-EPMC4593114 | biostudies-literature
Project description:DEER (double electron-electron resonance) is a powerful pulsed ESR (electron spin resonance) technique allowing the determination of distance histograms between pairs of nitroxide spin-labels linked to a protein in a native-like solution environment. However, exploiting the huge amount of information provided by ESR/DEER histograms to refine structural models is extremely challenging. In this study, a restrained ensemble (RE) molecular dynamics (MD) simulation methodology is developed to address this issue. In RE simulation, the spin-spin distance distribution histograms calculated from a multiple-copy MD simulation are enforced, via a global ensemble-based energy restraint, to match those obtained from ESR/DEER experiments. The RE simulation is applied to 51 ESR/DEER distance histogram data from spin-labels inserted at 37 different positions in T4 lysozyme (T4L). The rotamer population distribution along the five dihedral angles connecting the nitroxide ring to the protein backbone is determined and shown to be consistent with available information from X-ray crystallography. For the purpose of structural refinement, the concept of a simplified nitroxide dummy spin-label is designed and parametrized on the basis of these all-atom RE simulations with explicit solvent. It is demonstrated that RE simulations with the dummy nitroxide spin-labels imposing the ESR/DEER experimental distance distribution data are able to systematically correct and refine a series of distorted T4L structures, while simple harmonic distance restraints are unsuccessful. This computationally efficient approach allows experimental restraints from DEER experiments to be incorporated into RE simulations for efficient structural refinement.
| S-EPMC3684008 | biostudies-literature
Project description:Ligand binding sites in proteins are often localized to deeply buried cavities, inaccessible to bulk solvent. Yet, in many cases binding of cognate ligands occurs rapidly. An intriguing system is presented by the L99A cavity mutant of T4 Lysozyme (T4L L99A) that rapidly binds benzene (~106 M-1s-1). Although the protein has long served as a model system for protein thermodynamics and crystal structures of both free and benzene-bound T4L L99A are available, the kinetic pathways by which benzene reaches its solvent-inaccessible binding cavity remain elusive. The current work, using extensive molecular dynamics simulation, achieves this by capturing the complete process of spontaneous recognition of benzene by T4L L99A at atomistic resolution. A series of multi-microsecond unbiased molecular dynamics simulation trajectories unequivocally reveal how benzene, starting in bulk solvent, diffuses to the protein and spontaneously reaches the solvent inaccessible cavity of T4L L99A. The simulated and high-resolution X-ray derived bound structures are in excellent agreement. A robust four-state Markov model, developed using cumulative 60 μs trajectories, identifies and quantifies multiple ligand binding pathways with low activation barriers. Interestingly, none of these identified binding pathways required large conformational changes for ligand access to the buried cavity. Rather, these involve transient but crucial opening of a channel to the cavity via subtle displacements in the positions of key helices (helix4/helix6, helix7/helix9) leading to rapid binding. Free energy simulations further elucidate that these channel-opening events would have been unfavorable in wild type T4L. Taken together and via integrating with results from experiments, these simulations provide unprecedented mechanistic insights into the complete ligand recognition process in a buried cavity. By illustrating the power of subtle helix movements in opening up multiple pathways for ligand access, this work offers an alternate view of ligand recognition in a solvent-inaccessible cavity, contrary to the common perception of a single dominant pathway for ligand binding.
| S-EPMC5979041 | biostudies-literature
Project description:The atomic-level mechanisms that coordinate ligand release from protein pockets are only known for a handful of proteins. Here, we report results from accelerated molecular dynamics simulations for benzene dissociation from the buried cavity of the T4 lysozyme Leu99Ala mutant (L99A). In these simulations, benzene is released through a previously characterized, sparsely populated room-temperature excited state of the mutant, explaining the coincidence for experimentally measured benzene off rate and apo protein slow-timescale NMR relaxation rates between ground and excited states. The path observed for benzene egress is a multistep ligand migration from the buried cavity to ultimate release through an opening between the F/G-, H-, and I-helices and requires a number of cooperative multiresidue and secondary-structure rearrangements within the C-terminal domain of L99A. These rearrangements are identical to those observed along the ground state to excited state transitions characterized by molecular dynamic simulations run on the Anton supercomputer. Analyses of the molecular properties of the residues lining the egress path suggest that protein surface electrostatic potential may play a role in the release mechanism. Simulations of wild-type T4 lysozyme also reveal that benzene-egress-associated dynamics in the L99A mutant are potentially exaggerations of the substrate-processivity-related dynamics of the wild type.
| S-EPMC6349996 | biostudies-literature
Project description:Protein crystallography and calorimetry were used to characterize the binding of 1,2-azaborines to model cavities in T4 lysozyme in direct comparison to their carbonaceous counterparts. In the apolar L99A cavity, affinity for Ab dropped only slightly versus benzene. In the cavity designed to accommodate a single hydrogen bond (L99A/M102Q), Gln102?O···H-N hydrogen bonding for Ab and BEtAb was observed in the crystallographic complexes. The strength of the hydrogen bonding was estimated as 0.94 and 0.64 kcal/mol for Ab and BEtAb, respectively. This work unambiguously demonstrates that 1,2-azaborines can be readily accommodated in classic aryl recognition pockets and establishes one of 1,2-azaborine's distinguishing features from its carbonaceous isostere benzene: its ability to serve as an NH hydrogen bond donor in a biological setting.
| S-EPMC5087281 | biostudies-literature