Constraint methods for determining pathways and free energy of activated processes

Activated processes from chemical reactions up to conformational transitions of large biomolecules are hampered by barriers which are overcome only by the input of some free energy of activation. Hence, the characteristic and rate-determining barrier regions are not sufficiently sampled by usual simulation techniques. Constraints on a reaction coordinate r have turned out to be a suitable means to explore difficult pathways without changing potential function, energy or temperature. For a dense sequence of values of r, the corresponding sequence of simulations provides a pathway for the process. As only one coordinate among thousands is fixed during each simulation, the pathway essentially reflects the system's internal dynamics. From mean forces the free energy profile can be calculated to obtain reaction rates and insight in the reaction mechanism. In the last decade, theoretical tools and computing capacity have been developed to a degree where simulations give impressive qualitative insight in the processes at quantitative agreement with experiments. Here, we give an introduction to reaction pathways and coordinates, and develop the theory of free energy as the potential of mean force. We clarify the connection between mean force and constraint force which is the central quantity evaluated, and discuss the mass metric tensor correction. Well-behaved coordinates without tensor correction are considered. We discuss the theoretical background and practical implementation on the example of the reaction coordinate of targeted molecular dynamics simulation. Finally, we compare applications of constraint methods and other techniques developed for the same purpose, and discuss the limits of the approach

Structure and mechanism of the reversible photoswitch of a fluorescent protein

Proteins that can be reversibly photoswitched between a fluorescent and a nonfluorescent state bear enormous potential in diverse fields, such as data storage, in vivo protein tracking, and subdiffraction resolution light microscopy. However, these proteins could hitherto not live up to their full potential because the molecular switching mechanism is not resolved. Here, we clarify the molecular photoswitching mechanism of asFP595, a green fluorescent protein (GFP)-like protein that can be transferred from a nonfluorescent "off" to a fluorescent "on" state and back again, by green and blue light, respectively. To this end, we establish reversible photoswitching of fluorescence in whole protein crystals and show that the switching kinetics in the crystal is identical with that in solution. Subsequent x-ray analysis demonstrated that upon the absorption of a green photon, the chromophore isomerizes from a trans (off) to a cis (on) state. Molecular dynamics calculations suggest that isomerization occurs through a bottom hula twist mechanism with concomitant rotation of both bonds of the chromophoric methine ring bridge. This insight into the switching mechanism should facilitate the targeted design of photoswitchable proteins. Reversible photoswitching of the protein chromophore system within intact crystals also constitutes a step toward the use of fluorescent proteins in three-dimensional data recording

Mesodynamics in the SARS nucleocapsid measured by NMR field cycling

Protein motions on all timescales faster than molecular tumbling are encoded in the spectral density. The dissection of complex protein dynamics is typically performed using relaxation rates determined at high and ultra-high field. Here we expand this range of the spectral density to low fields through field cycling using the nucleocapsid protein of the SARS coronavirus as a model system. The field-cycling approach enables site-specific measurements of R1 at low fields with the sensitivity and resolution of a high-field magnet. These data, together with high-field relaxation and heteronuclear NOE, provide evidence for correlated rigid-body motions of the entire β-hairpin, and corresponding motions of adjacent loops with a time constant of 0.8 ns (mesodynamics). MD simulations substantiate these findings and provide direct verification of the time scale and collective nature of these motions

Structural Basis for Polyproline-Mediated Ribosome Stalling and Rescue by the Translation Elongation Factor EF-P

Ribosomes synthesizing proteins containing consecutive proline residues become stalled and require rescue via the action of uniquely modified translation elongation factors, EF-P in bacteria, or archaeal/eukaryotic a/eIF5A. To date, no structures exist of EF-P or eIF5A in complex with translating ribosomes stalled at polyproline stretches, and thus structural insight into how EF-P/eIF5A rescue these arrested ribosomes has been lacking. Here we present cryo-EM structures of ribosomes stalled on proline stretches, without and with modified EF-P. The structures suggest that the favored conformation of the polyproline-containing nascent chain is incompatible with the peptide exit tunnel of the ribosome and leads to destabilization of the peptidyltRNA. Binding of EF-P stabilizes the P-site tRNA, particularly via interactions between its modification and the CCA end, thereby enforcing an alternative conformation of the polyproline-containing nascent chain, which allows a favorable substrate geometry for peptide bond formation

Strength of Hydrogen Bond Network Takes Crucial Roles in the Dissociation Process of Inhibitors from the HIV-1 Protease Binding Pocket

To understand the underlying mechanisms of significant differences in dissociation rate constant among different inhibitors for HIV-1 protease, we performed steered molecular dynamics (SMD) simulations to analyze the entire dissociation processes of inhibitors from the binding pocket of protease at atomistic details. We found that the strength of hydrogen bond network between inhibitor and the protease takes crucial roles in the dissociation process. We showed that the hydrogen bond network in the cyclic urea inhibitors AHA001/XK263 is less stable than that of the approved inhibitor ABT538 because of their large differences in the structures of the networks. In the cyclic urea inhibitor bound complex, the hydrogen bonds often distribute at the flap tips and the active site. In contrast, there are additional accessorial hydrogen bonds formed at the lateral sides of the flaps and the active site in the ABT538 bound complex, which take crucial roles in stabilizing the hydrogen bond network. In addition, the water molecule W301 also plays important roles in stabilizing the hydrogen bond network through its flexible movement by acting as a collision buffer and helping the rebinding of hydrogen bonds at the flap tips. Because of its high stability, the hydrogen bond network of ABT538 complex can work together with the hydrophobic clusters to resist the dissociation, resulting in much lower dissociation rate constant than those of cyclic urea inhibitor complexes. This study may provide useful guidelines for design of novel potent inhibitors with optimized interactions

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

Accurate multiple time step in biased molecular simulations

Many recently introduced enhanced sampling techniques are based on biasing coarse descriptors (collective variables) of a molecular system on the fly. Sometimes the calculation of such collective variables is expensive and becomes a bottleneck in molecular dynamics simulations. An algorithm to treat smooth biasing forces within a multiple time step framework is here discussed. The implementation is simple and allows a speed up when expensive collective variables are employed. The gain can be substantial when using massively parallel or GPU-based molecular dynamics software. Moreover, a theoretical framework to assess the sampling accuracy is introduced, which can be used to assess the choice of the integration time step in both single and multiple time step biased simulations

A Steered Molecular Dynamics Study of Binding and Translocation Processes in the GABA Transporter

The entire substrate translocation pathway in the human GABA transporter (GAT-1) was explored for the endogenous substrate GABA and the anti-convulsive drug tiagabine. Following a steered molecular dynamics (SMD) approach, in which a harmonic restraining potential is applied to the ligand, dissociation and re-association of ligands were simulated revealing events leading to substrate (GABA) translocation and inhibitor (tiagabine) mechanism of action. We succeeded in turning the transporter from the outward facing occluded to the open-to-out conformation, and also to reorient the transporter to the open-to-in conformation. The simulations are validated by literature data and provide a substrate pathway fingerprint in terms of which, how, and in which sequence specific residues are interacted with. They reveal the essential functional roles of specific residues, e.g. the role of charged residues in the extracellular vestibule including two lysines (K76 (TM1) and K448 (TM10)) and a TM6-triad (D281, E283, and D287) in attracting and relocating substrates towards the secondary/interim substrate-binding site (S2). Likewise, E101 is highlighted as essential for the relocation of the substrate from the primary substrate-binding site (S1) towards the cytoplasm