Spontaneous quaternary and tertiary T-R transitions of human hemoglobin in molecular dynamics simulation.

We present molecular dynamics simulations of unliganded human hemoglobin (Hb) A under physiological conditions, starting from the R, R2, and T state. The simulations were carried out with protonated and deprotonated HC3 histidines His(beta)146, and they sum up to a total length of 5.6 micros.

Ligand docking and binding site analysis with PyMOL and Autodock/Vina.

Docking of small molecule compounds into the binding site of a receptor and estimating the binding affinity of the complex is an important part of the structure-based drug design process. For a thorough understanding of the structural principles that determine the strength of a protein/ligand complex both, an accurate and fast docking protocol and the ability to visualize binding geometries and interactions are mandatory. Here we present an interface between the popular molecular graphics system PyMOL and the molecular docking suites Autodock and Vina and demonstrate how the combination of docking and visualization can aid structure-based drug design efforts.

The anti-protozoal drug pentamidine blocks KIR2.x-mediated inward rectifier current by entering the cytoplasmic pore region of the channel.

Background And Purpose: Pentamidine is a drug used in treatment of protozoal infections. Pentamidine treatment may cause sudden cardiac death by provoking cardiac arrhythmias associated with QTc prolongation and U-wave alterations. This proarrhythmic effect was linked to inhibition of hERG trafficking, but not to acute block of ion channels contributing to the action potential.

Toward a consensus model of the HERG potassium channel.

Malfunction of hERG potassium channels, due to inherited mutations or inhibition by drugs, can cause long QT syndrome, which can lead to life-threatening arrhythmias. A three-dimensional structure of hERG is a prerequisite to understand the molecular basis of hERG malfunction. To achieve a consensus model, we carried out an extensive analysis of hERG models based on various alignments of helix S5.

Conformational transitions upon ligand binding: holo-structure prediction from apo conformations.

Biological function of proteins is frequently associated with the formation of complexes with small-molecule ligands. Experimental structure determination of such complexes at atomic resolution, however, can be time-consuming and costly. Computational methods for structure prediction of protein/ligand complexes, particularly docking, are as yet restricted by their limited consideration of receptor flexibility, rendering them not applicable for predicting protein/ligand complexes if large conformational changes of the receptor upon ligand binding are involved.

A plate reader-based method for cell water permeability measurement.

Cell volume and water permeability measurements in cultured mammalian cells are typically conducted under a light microscope. Many of the employed approaches are time consuming and not applicable to a study of confluent epithelial cell monolayers. We present here an adaptation of a calcein-quenching-based approach for a plate reader.

Detection of functional modes in protein dynamics.

Proteins frequently accomplish their biological function by collective atomic motions. Yet the identification of collective motions related to a specific protein function from, e.g.

Secondary structure propensities in peptide folding simulations: a systematic comparison of molecular mechanics interaction schemes.

We present a systematic study directed toward the secondary structure propensity and sampling behavior in peptide folding simulations with eight different molecular dynamics force-field variants in explicit solvent. We report on the combinational result of force field, water model, and electrostatic interaction schemes and compare to available experimental characterization of five studied model peptides in terms of reproduced structure and dynamics. The total simulation time exceeded 18 mus and included simulations that started from both folded and extended conformations.

Crystal structure of a yeast aquaporin at 1.15 angstrom reveals a novel gating mechanism.

Aquaporins are transmembrane proteins that facilitate the flow of water through cellular membranes. An unusual characteristic of yeast aquaporins is that they frequently contain an extended N terminus of unknown function. Here we present the X-ray structure of the yeast aquaporin Aqy1 from Pichia pastoris at 1.

The thermodynamic influence of trapped water molecules on a protein-ligand interaction.

Water molecules doing time: Atomic-resolution crystal structures of the PPIase domain of cyclophilin G, alone and in complex with cyclosporin A, and together with MD simulations and calorimetry, reveal how trapped water molecules influence the thermodynamic profile of a protein-ligand interaction.

Determinants of water permeability through nanoscopic hydrophilic channels.

Naturally occurring pores show a variety of polarities and sizes that are presumably directly linked to their biological function. Many biological channels are selective toward permeants similar or smaller in size than water molecules, and therefore their pores operate in the regime of single-file water pores. Intrinsic factors affecting water permeability through such pores include the channel-membrane match, the structural stability of the channel, the channel geometry and channel-water affinity.

Dynamics and energetics of permeation through aquaporins. What do we learn from molecular dynamics simulations?

Aquaporins (AQPs) are a family of integral membrane proteins, which facilitate the rapid and yet highly selective flux of water and other small solutes across biological membranes. Molecular dynamics (MD) simulations contributed substantially to the understanding of the molecular mechanisms that underlie this remarkable efficiency and selectivity of aquaporin channels. This chapter reviews the current state of MD simulations of aquaporins and related aquaglyceroporins as well as the insights these simulations have provided.

Domain motions of hyaluronan lyase underlying processive hyaluronan translocation.

Hyaluronan lyase (Hyal) is a surface enzyme occurring in many bacterial organisms including members of Streptococcus species. Streptococcal Hyal primarily degrades hyaluronan-substrate (HA) of the extracellular matrix. This degradation appears to facilitate the spread of this bacterium throughout host tissues.

tCONCOORD-GUI: visually supported conformational sampling of bioactive molecules.

Conformational flexibility of bioactive molecules poses a major challenge to computational biology. tCONCOORD generates structure ensembles based on geometrical considerations and has been successfully applied to predict protein conformational flexibility and essential degrees of freedom. We have now developed a graphical user interface (GUI) for tCONCOORD, which substantially facilitates the simulation setup and provides valuable insights into the structure analysis and constraint definition process in tCONCOORD.

The atomistic mechanism of conformational transition in adenylate kinase: a TEE-REX molecular dynamics study.

We report on an atomistic molecular dynamics simulation of the complete conformational transition of Escherichia coli adenylate kinase (ADK) using the recently developed TEE-REX algorithm. Two phases characterize the transition pathway of ADK, which folds into the domains CORE and LID and the AMP binding domain AMPbd. Starting from the closed conformation, half-opening of the AMPbd precedes a partially correlated opening of the LID and AMPbd, defining the second phase.

Recognition dynamics up to microseconds revealed from an RDC-derived ubiquitin ensemble in solution.

Protein dynamics are essential for protein function, and yet it has been challenging to access the underlying atomic motions in solution on nanosecond-to-microsecond time scales. We present a structural ensemble of ubiquitin, refined against residual dipolar couplings (RDCs), comprising solution dynamics up to microseconds. The ensemble covers the complete structural heterogeneity observed in 46 ubiquitin crystal structures, most of which are complexes with other proteins.

Self-consistent residual dipolar coupling based model-free analysis for the robust determination of nanosecond to microsecond protein dynamics.

Residual dipolar couplings (RDCs) provide information about the dynamic average orientation of inter-nuclear vectors and amplitudes of motion up to milliseconds. They complement relaxation methods, especially on a time-scale window that we have called supra-tau(c) (tau(c) < supra-tau(c) < 50 micros). Here we present a robust approach called Self-Consistent RDC-based Model-free analysis (SCRM) that delivers RDC-based order parameters-independent of the details of the structure used for alignment tensor calculation-as well as the dynamic average orientation of the inter-nuclear vectors in the protein structure in a self-consistent manner.

Not only enthalpy: large entropy contribution to ion permeation barriers in single-file channels.

The effect of channel length on the barrier for potassium ion permeation through single-file channels has been studied by means of all-atom molecular dynamics simulations. Using series of peptidic gramicidin-like and simplified ring-structured channels, both embedded in model membranes, we obtained two distinct types of behavior: saturation of the central free energy barriers for peptidic channels and a linear increase in simplified ring-structured channels with increasing channel length. The saturation of the central free energy barrier for the peptidic channels occurs at relatively short lengths, and it is correlated with the desolvation from the bulk water.

Kinetics, statistics, and energetics of lipid membrane electroporation studied by molecular dynamics simulations.

Membrane electroporation is the method to directly transfer bioactive substances such as drugs and genes into living cells, as well as preceding electrofusion. Although much information on the microscopic mechanism has been obtained both from experiment and simulation, the existence and nature of possible intermediates is still unclear. To elucidate intermediates of electropore formation by direct comparison with measured prepore formation kinetics, we have carried out 49 atomistic electroporation simulations on a palmitoyl-oleoyl-phosphatidylcholine bilayer for electric field strengths between 0.

The molecular mechanism of toxin-induced conformational changes in a potassium channel: relation to C-type inactivation.

Recently, a solid-state NMR study revealed that scorpion toxin binding leads to conformational changes in the selectivity filter of potassium channels. The exact nature of the conformational changes, however, remained elusive. We carried out all-atom molecular dynamics simulations that enabled us to cover the complete pathway of toxin approach and binding, and we validated our simulation results by using solid-state NMR data and electrophysiological measurements.

Normal modes and essential dynamics.

Normal mode analysis and essential dynamics analysis are powerful methods used for the analysis of collective motions in biomolecules. Their application has led to an appreciation of the importance of protein dynamics in function and the relationship between structure and dynamical behavior. In this chapter, the methods and their implementation are introduced and recent developments such as elastic networks and advanced sampling techniques are described.