Rotamer Dynamics: Analysis of Rotamers in Molecular Dynamics Simulations of Proteins

Given by χ torsional angles, rotamers describe the side-chain conformations of amino acid residues in a protein based on the rotational isomers (hence the word rotamer). Constructed rotamer libraries, based on either protein crystal structures or dynamics studies, are the tools for classifying rotamers (torsional angles) in a way that reflect their frequency in nature. Rotamer libraries are routinely used in structure modeling and evaluation. In this perspective article, we would like to encourage researchers to apply rotamer analyses beyond their traditional use. Molecular dynamics (MD) of proteins highlight the in silico behavior of molecules in solution and thus can identify favorable side-chain conformations. In this article, we used simple computational tools to study rotamer dynamics (RD) in MD simulations. First, we isolated each frame in the MD trajectories in separate Protein Data Bank files via the cpptraj module in AMBER. Then, we extracted torsional angles via the Bio3D module in R language. The classification of torsional angles was also done in R according to the penultimate rotamer library. RD analysis is useful for various applications such as protein folding, study of rotamer-rotamer relationship in protein-protein interaction, real-time correlation between secondary structures and rotamers, study of flexibility of side chains in binding site for molecular docking preparations, use of RD as guide in functional analysis and study of structural changes caused by mutations, providing parameters for improving coarse-grained MD accuracy and speed, and many others. Major challenges facing RD to emerge as a new scientific field involve the validation of results via easy, inexpensive wet-lab methods. This realm is yet to be explored.     

Lipid Configurations from Molecular Dynamics Simulations

The extent to which current force fields faithfully reproduce conformational properties of lipids in bilayer membranes, and whether these reflect the structural principles established for phospholipids in bilayer crystals, are central to biomembrane simulations. We determine the distribution of dihedral angles in palmitoyl-oleoyl phosphatidylcholine from molecular dynamics simulations of hydrated fluid bilayer membranes. We compare results from the widely used lipid force field of Berger et al. with those from the most recent C36 release of the CHARMM force field for lipids. Only the CHARMM force field produces the chain inequivalence with sn-1 as leading chain that is characteristic of glycerolipid packing in fluid bilayers. The exposure and high partial charge of the backbone carbonyls in Berger lipids leads to artifactual binding of Na⁺ ions reported in the literature. Both force fields predict coupled, near-symmetrical distributions of headgroup dihedral angles, which is compatible with models of interconverting mirror-image conformations used originally to interpret NMR order parameters. The Berger force field produces rotamer populations that correspond to the headgroup conformation found in a phosphatidylcholine lipid bilayer crystal, whereas CHARMM36 rotamer populations are closer to the more relaxed crystal conformations of phosphatidylethanolamine and glycerophosphocholine. CHARMM36 alone predicts the correct relative signs of the time-average headgroup order parameters, and reasonably reproduces the full range of NMR data from the phosphate diester to the choline methyls. There is strong motivation to seek further experimental criteria for verifying predicted conformational distributions in the choline headgroup, including the ³¹P chemical shift anisotropy and ¹⁴N and CD₃ NMR quadrupole splittings.     

Cryo-Cooling Effect on DHFR Crystal Studied by Replica-Exchange Molecular Dynamics Simulations

Cryo-cooling is routinely performed before x-ray diffraction image collection to reduce the damage to crystals due to ionizing radiation. It has been suggested that although backbone structures are usually very similar between room temperature and cryo-temperature, cryo-cooling may hamper biologically relevant dynamics. In this study, the crystal of Escherichia coli dihydrofolate reductase is studied with replica-exchange molecular dynamics simulation, and the results are compared with the crystal structure determined at cryo-temperature and room temperature with the time-averaged ensemble method. Although temperature dependence of unit cell compaction and root mean-square fluctuation of Cα is found in accord with experiment, it is found that the protein structure at low temperature can be more heterogeneous than the ensemble of structures reported by using the time-averaged ensemble method, encouraging further development of the time-averaged ensemble method and indicating that data should be examined carefully to avoid overinterpretation of one average structure.     

Molecular Dynamics of the Association of L-Selectin and FERM Regulated by PIP2

Phosphatidylinositol 4,5-bisphosphate (PIP2) acts as a signaling lipid, mediating membrane trafficking and recruitment of proteins to membranes. A key example is the PIP2-dependent regulation of the adhesion of L-selectin to the cytoskeleton adaptors of the N-terminal subdomain of ezrin-radixin-moesin (FERM). The molecular details of the mediating behavior of multivalent anionic PIP2 lipids in this process, however, remain unclear. Here, we use coarse-grained molecular dynamics simulation to explore the mechanistic details of PIP2 in the transformation, translocation, and association of the FERM/L-selectin complex. We compare membranes of different compositions and find that anionic phospholipids are necessary for both FERM and the cytoplasmic domain of L-selectin to absorb on the membrane surface. The subsequent formation of the FERM/L-selectin complex is strongly favored by the presence of PIP2, which clusters around both proteins and triggers a conformational transition in the cytoplasmic domain of L-selectin. We are able to quantify the effect of PIP2 on the association free energy of the complex by means of a potential of mean force. We conclude that PIP2 behaves as an adhesive agent to enhance the stability of the FERM/L-selectin complex and identify key residues involved. The molecular information revealed in this study highlights the specific role of membrane lipids such as PIP2 in protein translocation and potential signaling.     

Decrypting the Heat Activation Mechanism of TRPV1 Channel by Molecular Dynamics Simulation

As a prototype cellular sensor, the TRPV1 cation channel undergoes a closed-to-open gating transition in response to various physical and chemical stimuli including noxious heat. Despite recent progress, the molecular mechanism of heat activation of TRPV1 gating remains enigmatic. Toward decrypting the structural basis of TRPV1 heat activation, we performed extensive molecular dynamics simulations (with cumulative simulation time of ～11 μs) for the wild-type channel and a constitutively active double mutant at different temperatures (30, 60, and 72°C), starting from a high-resolution closed-channel structure of TRPV1 solved by cryo-electron microscopy. In the wild-type simulations, we observed heat-activated conformational changes (e.g., expansion or contraction) in various key domains of TRPV1 (e.g., the S2-S3 and S4-S5 linkers) to prime the channel for gating. These conformational changes involve a number of dynamic hydrogen-bond interactions that were validated with previous mutational studies. Next, our mutant simulations observed channel opening after a series of conformational changes that propagate from the channel periphery to the channel pore via key intermediate domains (including the S2-S3 and S4-S5 linkers). The gating transition is accompanied by a large increase in the protein-water electrostatic interaction energy, which supports the contribution of desolvation of polar/charged residues to the temperature-sensitive TRPV1 gating. Taken together, our molecular dynamics simulations and analyses offered, to our knowledge, new structural, dynamic, and energetic information to guide future mutagenesis and functional studies of the TRPV1 channels and development of TRPV1-targeting drugs.     

Binding of Net-Fla,a Netropsin-Flavin Hybrid Molecule,to DNA: Molecular Mechanics and Dynamics Studies in vacuo and in Water Solution

Abstract

We have studied the binding of the hybrid netropsin-flavin (Net-Fla) molecule onto four sequences containing four A.T base pairs. Molecular mechanics minimizations in vacuo show numerous minimal conformations separated by one base pair. 400 ps molecular dynamics simulations in vacuo have been performed using the lowest minima as the starting conformations. During these simulations, the flavin moiety of the drug makes two hydrogen bonds with an amino group of a neighboring guanine. A 200 ps molecular dynamics simulation in explicit water solution suggests that the binding of Net-Fla upon the DNA substrate is enhanced by water bridges. A water molecule bridging the amidinium of Net-Fla to the N3 atom of an adenine seems to be stuck in the dmg-DNA complex during the whole simulation. The fluctuations of the DNA helical parameters and of the torsion angles of the sugar-phosphate backbone are very similar in the simulations in vacuo and in water. The time auto-correlation functions for the DNA helical parameters decrease rapidly in the picosecond range in vacuo. The same functions computed from the water solution molecular dynamics simulations seem to have two modes: the rapid mode is similar to the behavior in vacuo, and is followed by a slower mode in the 10 ps range.     

Dynamics of a Protein Chain Motor Driving Helical Bacteria under Stress

The wall-less, helical bacterial genus Spiroplasma has a unique propulsion system; it is not driven by propeller-like flagella but by a membrane-bound, cytoplasmic, linear motor that consists of a contractile chain of identical proteins spanning the entire cell length. By a coordinated spread of conformational changes of the proteins, kinks propagate in pairs along the cell body. However, the mechanisms for the initiation or delay of kinks and their coordinated spread remain unclear. Here, we show how we manipulate the initiation of kinks, their propagation velocities, and the time between two kinks for a single cell trapped in an optical line potential. By interferometric three-dimensional shape tracking, we measured the cells’ deformations in response to various external stress situations. We observed a significant dependency of force generation on the cells’ local ligand concentrations (likely ATP) and ligand hydrolysis, which we altered in different ways. We developed a mechanistic, mathematical model based on Kramer’s rates, describing the subsequent cooperative and conformational switching of the chain’s proteins. The model reproduces our experimental observations and can explain deformation characteristics even when the motor is driven to its extreme. Nature has invented a set of minimalistic mechanical driving concepts. To understand or even rebuild them, it is essential to reveal the molecular mechanisms of such protein chain motors, which need only two components—coupled proteins and ligands—to function.     

Speed of Conformational Change: Comparing Explicit and Implicit Solvent Molecular Dynamics Simulations

Adequate sampling of conformation space remains challenging in atomistic simulations, especially if the solvent is treated explicitly. Implicit-solvent simulations can speed up conformational sampling significantly. We compare the speed of conformational sampling between two commonly used methods of each class: the explicit-solvent particle mesh Ewald (PME) with TIP3P water model and a popular generalized Born (GB) implicit-solvent model, as implemented in the AMBER package. We systematically investigate small (dihedral angle flips in a protein), large (nucleosome tail collapse and DNA unwrapping), and mixed (folding of a miniprotein) conformational changes, with nominal simulation times ranging from nanoseconds to microseconds depending on system size. The speedups in conformational sampling for GB relative to PME simulations, are highly system- and problem-dependent. Where the simulation temperatures for PME and GB are the same, the corresponding speedups are approximately onefold (small conformational changes), between ∼1- and ∼100-fold (large changes), and approximately sevenfold (mixed case). The effects of temperature on speedup and free-energy landscapes, which may differ substantially between the solvent models, are discussed in detail for the case of miniprotein folding. In addition to speeding up conformational sampling, due to algorithmic differences, the implicit solvent model can be computationally faster for small systems or slower for large systems, depending on the number of solute and solvent atoms. For the conformational changes considered here, the combined speedups are approximately twofold, ∼1- to 60-fold, and ∼50-fold, respectively, in the low solvent viscosity regime afforded by the implicit solvent. For all the systems studied, 1) conformational sampling speedup increases as Langevin collision frequency (effective viscosity) decreases; and 2) conformational sampling speedup is mainly due to reduction in solvent viscosity rather than possible differences in free-energy landscapes between the solvent models.     

A Dynamic Hydrophobic Core and Surface Salt Bridges Thermostabilize a Designed Three-Helix Bundle

Thermostable proteins are advantageous in industrial applications, as pharmaceuticals or biosensors, and as templates for directed evolution. As protein-design methodologies improve, bioengineers are able to design proteins to perform a desired function. Although many rationally designed proteins end up being thermostable, how to intentionally design de novo, thermostable proteins is less clear. UVF is a de novo-designed protein based on the backbone structure of the Engrailed homeodomain (EnHD) and is highly thermostable (T_m > 99°C vs. 52°C for EnHD). Although most proteins generally have polar amino acids on their surfaces and hydrophobic amino acids buried in their cores, protein engineers followed this rule exactly when designing UVF. To investigate the contributions of the fully hydrophobic core versus the fully polar surface to UVF’s thermostability, we built two hybrid, chimeric proteins combining the sets of buried and surface residues from UVF and EnHD. Here, we determined a structural, dynamic, and thermodynamic explanation for UVF’s thermostability by performing 4 μs of all-atom, explicit-solvent molecular dynamics simulations at 25 and 100°C, Tanford-Kirkwood solvent accessibility Monte Carlo electrostatic calculations, and a thermodynamic analysis of 40 temperature runs by the weighted-histogram analysis method of heavy-atom, structure-based models of UVF, EnHD, and both chimeric proteins. Our models showed that UVF was highly dynamic because of its fully hydrophobic core, leading to a smaller loss of entropy upon folding. The charged residues on its surface made favorable electrostatic interactions that contributed enthalpically to its thermostability. In the chimeric proteins, both the hydrophobic core and charged surface independently imparted thermostability.     

PELDOR Spectroscopy Reveals Two Defined States of a Sialic Acid TRAP Transporter SBP in Solution

The tripartite ATP-independent periplasmic (TRAP) transporters are a widespread class of membrane transporters in bacteria and archaea. Typical substrates for TRAP transporters are organic acids including the sialic acid N-acetylneuraminic acid. The substrate binding proteins (SBP) of TRAP transporters are the best studied component and are responsible for initial high-affinity substrate binding. To better understand the dynamics of the ligand binding process, pulsed electron-electron double resonance (PELDOR, also known as DEER) spectroscopy was applied to study the conformational changes in the N-acetylneuraminic acid-specific SBP VcSiaP. The protein is the SBP of VcSiaPQM, a sialic acid TRAP transporter from Vibrio cholerae. Spin-labeled double-cysteine mutants of VcSiaP were analyzed in the substrate-bound and -free state and the measured distances were compared to available crystal structures. The data were compatible with two clear states only, which are consistent with the open and closed forms seen in TRAP SBP crystal structures. Substrate titration experiments demonstrated the transition of the population from one state to the other with no other observed forms. Mutants of key residues involved in ligand binding and/or proposed to be involved in domain closure were produced and the corresponding PELDOR experiments reveal important insights into the open-closed transition. The results are in excellent agreement with previous in vivo sialylation experiments. The structure of the spin-labeled Q54R1/L173R1 R125A mutant was solved at 2.1 Å resolution, revealing no significant changes in the protein structure. Thus, the loss of domain closure appears to be solely due to loss of binding. In conclusion, these data are consistent with TRAP SBPs undergoing a simple two-state transition from an open-unliganded to closed-liganded state during the transport cycle.     

Architectural Dynamics of CaMKII-Actin Networks

Calcium-calmodulin-dependent kinase II (CaMKII) has an important role in dendritic spine remodeling upon synaptic stimulation. Using fluorescence video microscopy and image analysis, we investigated the architectural dynamics of rhodamine-phalloidin stabilized filamentous actin (F-actin) networks cross-linked by CaMKII. We used automated image analysis to identify F-actin bundles and crossover junctions and developed a dimensionless metric to characterize network architecture. Similar networks were formed by three different CaMKII species with a 10-fold length difference in the linker region between the kinase domain and holoenzyme hub, implying linker length is not a primary determinant of F-actin cross-linking. Electron micrographs showed that at physiological molar ratios, single CaMKII holoenzymes cross-linked multiple F-actin filaments at random, whereas at higher CaMKII/F-actin ratios, filaments bundled. Light microscopy established that the random network architecture resisted macromolecular crowding with polyethylene glycol and blocked ATP-powered compaction by myosin-II miniature filaments. Importantly, the networks disassembled after the addition of calcium-calmodulin and were then spaced within 3 min into compacted foci by myosin motors or more slowly (30 min) aggregated by crowding. Single-molecule total internal reflection fluorescence microscopy showed CaMKII dissociation from surface-immobilized globular actin exhibited a monoexponential dwell-time distribution, whereas CaMKII bound to F-actin networks had a long-lived fraction, trapped at crossover junctions. Release of CaMKII from F-actin, triggered by calcium-calmodulin, was too rapid to measure with flow-cell exchange (<20 s). The residual bound fraction was reduced substantially upon addition of an N-methyl-D-aspartate receptor peptide analog but not ATP. These results provide mechanistic insights to CaMKII-actin interactions at the collective network and single-molecule level. Our findings argue that CaMKII-actin networks in dendritic spines maintain spine size against physical stress. Upon synaptic stimulation, CaMKII is disengaged by calcium-calmodulin, triggering network disassembly, expansion, and subsequent compaction by myosin motors with kinetics compatible with the times recorded for the poststimulus changes in spine volume.