Search

Principal component analysis (PCA) represents a standard approach to identify collective variables {xi} = x, which can be used to construct the free energy landscape ΔG(x) of a molecular system. While PCA is routinely applied to equilibrium molecular dynamics (MD) simulations, it is less obvious as to how to extend the approach to nonequilibrium simulation techniques. This includes, e.g., the definition of the statistical averages employed in PCA as well as the relation between the equilibrium free energy landscape ΔG(x) and the energy landscapes Δ G (x) obtained from nonequilibrium MD. As an example for a nonequilibrium method, "targeted MD" is considered which employs a moving distance constraint to enforce rare transitions along some biasing coordinate s. The introduced bias can be described by a weighting function P(s), which provides a direct relation between equilibrium and nonequilibrium data, and thus establishes a well-defined way to perform PCA on nonequilibrium data. While the resulting distribution P (x) and energy Δ G ∝ ln P will not reflect the equilibrium state of the system, the nonequilibrium energy landscape Δ G (x) may directly reveal the molecular reaction mechanism. Applied to targeted MD simulations of the unfolding of decaalanine, for example, a PCA performed on backbone dihedral angles is shown to discriminate several unfolding pathways. Although the formulation is in principle exact, its practical use depends critically on the choice of the biasing coordinate s, which should account for a naturally occurring motion between two well-defined end-states of the system. [ABSTRACT FROM AUTHOR]

Vibrational energy transfer (VET) is essential for protein function. It is responsible for efficient energy dissipation in reaction sites, and has been linked to pathways of allosteric communication. While it is understood that VET occurs via backbone as well as via non-covalent contacts, little is known about the competition of these two transport channels, which determines the VET pathways. To tackle this problem, we equipped the β-hairpin fold of a tryptophan zipper with pairs of non-canonical amino acids, one serving as a VET injector and one as a VET sensor in a femtosecond pump probe experiment. Accompanying extensive non-equilibrium molecular dynamics simulations combined with a master equation analysis unravel the VET pathways. Our joint experimental/computational endeavor reveals the efficiency of backbone vs. contact transport, showing that even if cutting short backbone stretches of only 3 to 4 amino acids in a protein, hydrogen bonds are the dominant VET pathway. Vibrational energy transfer (VET) is essential for protein function as it is responsible for efficient energy dissipation in reaction sites and is linked to pathways of allosteric communication. Here authors equipped a tryptophan zipper with a VET injector and a VET sensor for femtosecond pump probe experiments to map the VET. [ABSTRACT FROM AUTHOR]

While allostery is of paramount importance for protein regulation, the underlying dynamical process of ligand (un)binding at one site, resulting time evolution of the protein structure, and change of the binding affinity at a remote site are not well understood. Here the ligand-induced conformational transition in a widely studied model system of allostery, the PDZ2 domain, is investigated by transient infrared spectroscopy accompanied by molecular dynamics simulations. To this end, an azobenzene derived photoswitch is linked to a peptide ligand in a way that its binding affinity to the PDZ2 domain changes upon switching, thus initiating an allosteric transition in the PDZ2 domain protein. The subsequent response of the protein, covering four decades of time, ranging from ∼1 ns to ∼10 µs, can be rationalized by a remodeling of its rugged free-energy landscape, with very subtle shifts in the populations of a small number of structurally well-defined states. It is proposed that structurally and dynamically driven allostery, often discussed as limiting scenarios of allosteric communication, actually go hand-in-hand, allowing the protein to adapt its free-energy landscape to incoming signals. [ABSTRACT FROM AUTHOR]

Coarse-graining of fully atomistic molecular dynamics simulations is a long-standing goal in order to allow the description of processes occurring on biologically relevant timescales. For example, the prediction of pathways, rates and rate-limiting steps in protein-ligand unbinding is crucial for modern drug discovery. To achieve the enhanced sampling, we perform dissipation-corrected targeted molecular dynamics simulations, which yield free energy and friction profiles of molecular processes under consideration. Subsequently, we use these fields to perform temperature-boosted Langevin simulations which account for the desired kinetics occurring on multisecond timescales and beyond. Adopting the dissociation of solvated sodium chloride, trypsin-benzamidine and Hsp90-inhibitor protein-ligand complexes as test problems, we reproduce rates from molecular dynamics simulation and experiments within a factor of 2–20, and dissociation constants within a factor of 1–4. Analysis of friction profiles reveals that binding and unbinding dynamics are mediated by changes of the surrounding hydration shells in all investigated systems. Protein-ligand unbinding processes are out of reach for atomistic simulations due to time-scale involved. Here the authors demonstrate an approach relying on dissipation-corrected targeted molecular dynamics that enables to provide binding and unbinding rates with a speed-up of several orders of magnitude. [ABSTRACT FROM AUTHOR]

Accessory β-subunits of the KCNE gene family modulate the function of various cation channel α-subunits by the formation of heteromultimers. Among the most dramatic changes of biophysical properties of a voltage-gated channel by KCNEs are the effects of KCNE1 on KCNQ1 channels. KCNQ1 and KCNE1 are believed to form nativeIKs channels. Here, we characterize molecular determinants of KCNE1 interaction with KCNQ1 channels by scanning mutagenesis, double mutant cycle analysis, and molecular dynamics simulations. Our findings suggest that KCNE1 binds to the outer face of the KCNQ1 channel pore domain, modifies interactions between voltage sensor, S4-S5 linker and the pore domain, leading to structural modifications of the selectivity filter and voltage sensor domain. Molecular dynamics simulations suggest a stable interaction of the KCNE1 transmembrane α-helix with the pore domain S5/S6 and part of the voltage sensor domain S4 of KCNQ1 in a putative pre-open channel state. Formation of this state may induce slow activation gating, the pivotal characteristic of native cardiac IKs channels. This new KCNQ1-KCNE1 model may become useful for dynamic modeling of disease-associated mutant IKs channels. Copyright © 2011 S. Karger AG, Basel [ABSTRACT FROM AUTHOR]