A fundamental problem in computational biophysics is to deduce the function of a protein from the structure. Many biological macromolecules such as enzymes, molecular motors or membrane transport proteins perform their function by cycling between multiple conformational states. Understanding such…
A fundamental problem in computational biophysics is to deduce the function of a protein from the structure. Many biological macromolecules such as enzymes, molecular motors or membrane transport proteins perform their function by cycling between multiple conformational states. Understanding such conformational transitions, which typically occur on the millisecond to second time scale, is central to understanding protein function. Molecular dynamics (MD) computer simulations have become an important tool to connect molecular structure to function, but equilibrium MD simulations are rarely able to sample on time scales longer than a few microseconds – orders of magnitudes shorter than the time scales of interest. A range of different simulation methods have been proposed to overcome this time-scale limitation. These include calculations of the free energy landscape and path sampling methods to directly sample transitions between known conformations. All these methods solve the problem to sample infrequently occupied but important regions of configuration space. Many path-sampling algorithms have been applied to the closed – open transition of the enzyme adenylate kinase (AdK), which undergoes a large, clamshell-like conformational transition between an open and a closed state. Here we review approaches to sample macromolecular transitions through the lens of AdK. We focus our main discussion on the current state of knowledge – both from simulations and experiments – about the transition pathways of ligand-free AdK, its energy landscape, transition rates and interactions with substrates. We conclude with a comparison of the discussed approaches with a view towards quantitative evaluation of path-sampling methods.
Date Created
The date the item was original created (prior to any relationship with the ASU Digital Repositories.)
The continuous random network (CRN) model of network glasses is widely accepted as a model for materials such as vitreous silica and amorphous silicon. Although it
has been more than eighty years since the proposal of the CRN, there has not…
The continuous random network (CRN) model of network glasses is widely accepted as a model for materials such as vitreous silica and amorphous silicon. Although it
has been more than eighty years since the proposal of the CRN, there has not been conclusive experimental evidence of the structure of glasses and amorphous
materials. This has now changed with the advent of two-dimensional amorphous materials. Now, not only the distribution of rings but the actual atomic ring
structure can be imaged in real space, allowing for greater charicterization of these types of networks. This dissertation reports the first work done
on the modelling of amorphous graphene and vitreous silica bilayers. Models of amorphous graphene have been created using a Monte Carlo bond-switching method
and MD method. Vitreous silica bilayers have been constructed using models of amorphous graphene and the ring statistics of silica bilayers has been studied.
Date Created
The date the item was original created (prior to any relationship with the ASU Digital Repositories.)
Transmembrane proteins are responsible for transporting ions and small molecules across the hydrophobic region of the cell membrane. We are reviewing the evidence for regulation of these transport processes by interactions with the lipids of the membrane. We focus on…
Transmembrane proteins are responsible for transporting ions and small molecules across the hydrophobic region of the cell membrane. We are reviewing the evidence for regulation of these transport processes by interactions with the lipids of the membrane. We focus on ion channels, including potassium channels, mechanosensitive and pentameric ligand gated ion channels, and active transporters, including pumps, sodium or proton driven secondary transporters and ABC transporters. For ion channels it has been convincingly shown that specific lipid-protein interactions can directly affect their function. In some cases, a combined approach of molecular and structural biology together with computer simulations has revealed the molecular mechanisms. There are also many transporters whose activity depends on lipids but understanding of the molecular mechanisms is only beginning.
Date Created
The date the item was original created (prior to any relationship with the ASU Digital Repositories.)
Molecular dynamics simulations were used to study properties of water at the interface with nanometer-size solutes. We simulated nonpolar attractive Kihara cavities given by a Lennard-Jones potential shifted by a core radius. The dipolar response of the hydration layer to…
Molecular dynamics simulations were used to study properties of water at the interface with nanometer-size solutes. We simulated nonpolar attractive Kihara cavities given by a Lennard-Jones potential shifted by a core radius. The dipolar response of the hydration layer to a uniform electric field substantially exceeds that of the bulk. For strongly attractive solutes, the collective dynamics of the hydration layer become slow compared to bulk water, as the solute size is increased. The statistics of electric field fluctuations at the solute center are Gaussian and tend toward the dielectric continuum limit with increasing solute size. A dipolar probe placed at the center of the solute is sensitive neither to the polarity excess nor to the slowed dynamics of the hydration layer. A point dipole was introduced close to the solute-water interface to further study the statistics of electric field fluctuations generated by the water. For small dipole magnitudes, the free energy surface is single-welled, with approximately Gaussian statistics. When the dipole is increased, the free energy surface becomes double-welled, before landing in an excited state, characterized again by a single-welled surface. The intermediate region is fairly broad and is characterized by electrostatic fluctuations significantly in excess of the prediction of linear response. We simulated a solute having the geometry of C180 fullerene, with dipoles introduced on each carbon. For small dipole moments, the solvent response follows the results seen for a single dipole; but for larger dipole magnitudes, the fluctuations of the solute-solvent energy pass through a second maximum. The juxtaposition of the two transitions leads to an approximately cubic scaling of the chemical potential with the dipole strengh. Umbrella sampling techniques were used to generate free energy surfaces of the electric potential fluctuations at the heme iron in Cytochrome B562. The results were unfortunately inconclusive, as the ionic background was not effectively represented in the finite-size system.
Date Created
The date the item was original created (prior to any relationship with the ASU Digital Repositories.)