Kinetics of diffusion-controlled enzymatic reactions with charged substrates

The Debye-Hückel limiting law (DHL) has often been used to estimate rate constants of diffusion-controlled reactions under different ionic strengths. Two main approximations are adopted in DHL: one is that the solution of the linearized Poisson-Boltzmann equation for a spherical cavity is used to estimate the excess electrostatic free energy of a solution; the other is that details of electrostatic interactions of the solutes are neglected. This makes DHL applicable only at low ionic strengths and dilute solutions (very low substrate/solute concentrations). We show in this work that through numerical solution of the Poisson-Nernst-Planck equations, diffusion-reaction processes can be studied at a variety of conditions including realistically concentrated solutions, high ionic strength, and certainly with non-equilibrium charge distributions. Reaction rate coefficients for the acetylcholine-acetylcholinesterase system are predicted to strongly depend on both ionic strength and substrate concentration. In particular, they increase considerably with increase of substrate concentrations at a fixed ionic strength, which is open to experimental testing. This phenomenon is also verified on a simple model, and is expected to be general for electrostatically attracting enzyme-substrate systems

Molecular dynamics simulations and drug discovery

This review discusses the many roles atomistic computer simulations of macromolecular (for example, protein) receptors and their associated small-molecule ligands can play in drug discovery, including the identification of cryptic or allosteric binding sites, the enhancement of traditional virtual-screening methodologies, and the direct prediction of small-molecule binding energies. The limitations of current simulation methodologies, including the high computational costs and approximations of molecular forces required, are also discussed. With constant improvements in both computer power and algorithm design, the future of computer-aided drug design is promising; molecular dynamics simulations are likely to play an increasingly important role

Accessing a Hidden Conformation of the Maltose Binding Protein Using Accelerated Molecular Dynamics

Periplasmic binding proteins (PBPs) are a large family of molecular transporters that play a key role in nutrient uptake and chemotaxis in Gram-negative bacteria. All PBPs have characteristic two-domain architecture with a central interdomain ligand-binding cleft. Upon binding to their respective ligands, PBPs undergo a large conformational change that effectively closes the binding cleft. This conformational change is traditionally viewed as a ligand induced-fit process; however, the intrinsic dynamics of the protein may also be crucial for ligand recognition. Recent NMR paramagnetic relaxation enhancement (PRE) experiments have shown that the maltose binding protein (MBP) - a prototypical member of the PBP superfamily - exists in a rapidly exchanging (ns to µs regime) mixture comprising an open state (approx 95%), and a minor partially closed state (approx 5%). Here we describe accelerated MD simulations that provide a detailed picture of the transition between the open and partially closed states, and confirm the existence of a dynamical equilibrium between these two states in apo MBP. We find that a flexible part of the protein called the balancing interface motif (residues 175–184) is displaced during the transformation. Continuum electrostatic calculations indicate that the repacking of non-polar residues near the hinge region plays an important role in driving the conformational change. Oscillations between open and partially closed states create variations in the shape and size of the binding site. The study provides a detailed description of the conformational space available to ligand-free MBP, and has implications for understanding ligand recognition and allostery in related proteins

Climate Change and Trophic Response of the Antarctic Bottom Fauna

BACKGROUND: As Earth warms, temperate and subpolar marine species will increasingly shift their geographic ranges poleward. The endemic shelf fauna of Antarctica is especially vulnerable to climate-mediated biological invasions because cold temperatures currently exclude the durophagous (shell-breaking) predators that structure shallow-benthic communities elsewhere. METHODOLOGY/PRINCIPAL FINDINGS: We used the Eocene fossil record from Seymour Island, Antarctic Peninsula, to project specifically how global warming will reorganize the nearshore benthos of Antarctica. A long-term cooling trend, which began with a sharp temperature drop approximately 41 Ma (million years ago), eliminated durophagous predators-teleosts (modern bony fish), decapod crustaceans (crabs and lobsters) and almost all neoselachian elasmobranchs (modern sharks and rays)-from Antarctic nearshore waters after the Eocene. Even prior to those extinctions, durophagous predators became less active as coastal sea temperatures declined from 41 Ma to the end of the Eocene, approximately 33.5 Ma. In response, dense populations of suspension-feeding ophiuroids and crinoids abruptly appeared. Dense aggregations of brachiopods transcended the cooling event with no apparent change in predation pressure, nor were there changes in the frequency of shell-drilling predation on venerid bivalves. CONCLUSIONS/SIGNIFICANCE: Rapid warming in the Southern Ocean is now removing the physiological barriers to shell-breaking predators, and crabs are returning to the Antarctic Peninsula. Over the coming decades to centuries, we predict a rapid reversal of the Eocene trends. Increasing predation will reduce or eliminate extant dense populations of suspension-feeding echinoderms from nearshore habitats along the Peninsula while brachiopods will continue to form large populations, and the intensity of shell-drilling predation on infaunal bivalves will not change appreciably. In time the ecological effects of global warming could spread to other portions of the Antarctic coast. The differential responses of faunal components will reduce the endemic character of Antarctic subtidal communities, homogenizing them with nearshore communities at lower latitudes

On the analysis of sedimentation velocity in the study of protein complexes

Sedimentation velocity analytical ultracentrifugation has experienced a significant transformation, precipitated by the possibility of efficiently fitting Lamm equation solutions to the experimental data. The precision of this approach depends on the ability to account for the imperfections of the experiment, both regarding the sample and the instrument. In the present work, we explore in more detail the relationship between the sedimentation process, its detection, and the model used in the mathematical data analysis. We focus on configurations that produce steep and fast-moving sedimentation boundaries, such as frequently encountered when studying large multi-protein complexes. First, as a computational tool facilitating the analysis of heterogeneous samples, we introduce the strategy of partial boundary modeling. It can simplify the modeling by restricting the direct boundary analysis to species with sedimentation coefficients in a predefined range. Next, we examine factors related to the experimental detection, including the magnitude of optical aberrations generated by out-of-focus solution columns at high protein concentrations, the relationship between the experimentally recorded signature of the meniscus and the meniscus parameter in the data analysis, and the consequences of the limited radial and temporal resolution of the absorbance optical scanning system. Surprisingly, we find that large errors can be caused by the finite scanning speed of the commercial absorbance optics, exceeding the statistical errors in the measured sedimentation coefficients by more than an order of magnitude. We describe how these effects can be computationally accounted for in SEDFIT and SEDPHAT

Computational Identification of Uncharacterized Cruzain Binding Sites

Chagas disease, caused by the unicellular parasite Trypanosoma cruzi, claims 50,000 lives annually and is the leading cause of infectious myocarditis in the world. As current antichagastic therapies like nifurtimox and benznidazole are highly toxic, ineffective at parasite eradication, and subject to increasing resistance, novel therapeutics are urgently needed. Cruzain, the major cysteine protease of Trypanosoma cruzi, is one attractive drug target. In the current work, molecular dynamics simulations and a sequence alignment of a non-redundant, unbiased set of peptidase C1 family members are used to identify uncharacterized cruzain binding sites. The two sites identified may serve as targets for future pharmacological intervention

Understanding biomolecular motion, recognition, and allostery by use of conformational ensembles

We review the role conformational ensembles can play in the analysis of biomolecular dynamics, molecular recognition, and allostery. We introduce currently available methods for generating ensembles of biomolecules and illustrate their application with relevant examples from the literature. We show how, for binding, conformational ensembles provide a way of distinguishing the competing models of induced fit and conformational selection. For allostery we review the classic models and show how conformational ensembles can play a role in unravelling the intricate pathways of communication that enable allostery to occur. Finally, we discuss the limitations of conformational ensembles and highlight some potential applications for the future

Single-Molecule Force Spectroscopy: Experiments, Analysis, and Simulations

International audienceThe mechanical properties of cells and of subcellular components are important to obtain a mechanistic molecular understanding of biological processes. The quantification of mechanical resistance of cells and biomolecules using biophysical methods matured thanks to the development of nanotechnologies such as optical and magnetic tweezers, the biomembrane force probe and atomic force microscopy (AFM). The quantitative nature of force spectroscopy measurements has converted AFM into a valuable tool in biophysics. Force spectroscopy allows the determination of the forces required to unfold protein domains and to disrupt individual receptor/ligand bonds. Molecular simulation as a computational microscope allows investigation of similar biological processes with an atomistic detail. In this chapter, we first provide a step-by-step protocol of force spectroscopy including sample preparation, measurement and analysis of force spectroscopy using AFM and its interpretation in terms of available theories. Next, we present the background for molecular dynamics (MD) simulations focusing on steered molecular dynamics (SMD) and the importance of bridging of computational tools with experimental technique

Molecular dynamics simulation approach to investigate dynamic behaviour of system through the application of newtonian mechanics

Molecular dynamics simulations have been successfully incorporated and evolved into a mature technique within a variety of pharmaceutical research programs to study the complex biological and chemical systems. Broadly used in modern drug design, molecular docking methods can be used effectively to understand the macromolecular structure-to-function relationships and ligand conformations adopted within the binding sites of macromolecular targets. Information gathered about the dynamic properties of ligand–receptor binding such as free energy by evaluating critical phenomena involved in the intermolecular recognition process. These results can be employed to shift the usual paradigm of structural bioinformatics from studying single structures to analyse conformational ensembles. Today, as a variety of docking algorithms are available, an understanding of advantages and limitations of each method is of fundamental importance in the development of effective strategies and the generation of relevant results. The purpose of this chapter is to examine the current molecular docking strategies used in drug discovery and medicinal chemistry, exploring the advancements in the field and role played by integration of structure-and ligand-based methods