Dynamic structure of retinylidene ligand of rhodopsin probed by molecular simulations

Rhodopsin is currently the only available atomic-resolution template for understanding biological functions of the G protein-coupled receptor (GPCR) family. The structural basis for the phenomenal dark state stability of 11-cis-retinal bound to rhodopsin and its ultrafast photoreaction are active topics of research. In particular, the beta-ionone ring of the retinylidene inverse agonist is crucial for the activation mechanism. We analyzed a total of 23 independent, 100 ns all-atom molecular dynamics simulations of rhodopsin embedded in a lipid bilayer in the microcanonical (N,V,E) ensemble. Analysis of intramolecular fluctuations predicts hydrogen-out-of-plane (HOOP) wagging modes of retinal consistent with those found in Raman vibrational spectroscopy. We show that sampling and ergodicity of the ensemble of simulations are crucial for determining the distribution of conformers of retinal bound to rhodopsin. The polyene chain is rigidly locked into a single, twisted conformation, consistent with the function of retinal as an inverse agonist in the dark state. Most surprisingly, the beta-ionone ring is mobile within its binding pocket; interactions are non-specific and the cavity is sufficiently large to enable structural heterogeneity. We find that retinal occupies two distinct conformations in the dark state, contrary to most previous assumptions. The beta-ionone ring can rotate relative to the polyene chain, thereby populating both positively and negatively twisted 6-s-cis enantiomers. This result, while unexpected, strongly agrees with experimental solid-state (2)H NMR spectra. Correlation analysis identifies the residues most critical to controlling mobility of retinal; we find that Trp265 moves away from the ionone ring prior to any conformational transition. Our findings reinforce how molecular dynamics simulations can challenge conventional assumptions for interpreting experimental data, especially where existing models neglect conformational fluctuations.     

Orientational relaxation time of bottom-heavy squirmers in a semi-dilute suspension

One of the important quantities to characterize unsteady behaviour of a cell suspension is the orientational relaxation time, which is the time scale for a micro-organism to re-orientate to its preferred direction from disorientated conditions. In this paper, a swimming micro-organism is modelled as a squirming sphere with prescribed tangential surface velocity, in which the centre of mass of the sphere is displaced from the geometric centre (bottom-heaviness). The orientational relaxation time of bottom-heavy squirmers in a suspension is investigated both analytically and numerically. The three-dimensional movement of 64 identical squirmers in a fluid otherwise at rest, contained in a cube with periodic boundary conditions, is dynamically computed, for random initial positions and orientations. The effects of volume fraction of squirmers, the bottom-heaviness and the squirming mode on the relaxation time are discussed. The results for a semi-dilute suspension show that both the mean stresslet strength and the orientational relaxation time decrease from those for a dilute suspension. We also observe a stress overshoot in some cases. The mechanism for this is different from that for a visco-elastic fluid, and is explained by the change with time of the orientation of squirmers.     

Visualization of High Speed Liquid Jet Impaction on a Moving Surface

Two apparatuses for examining liquid jet impingement on a high-speed moving surface are described: an air cannon device (for examining surface speeds between 0 and 25 m/sec) and a spinning disk device (for examining surface speeds between 15 and 100 m/sec). The air cannon linear traverse is a pneumatic energy-powered system that is designed to accelerate a metal rail surface mounted on top of a wooden projectile. A pressurized cylinder fitted with a solenoid valve rapidly releases pressurized air into the barrel, forcing the projectile down the cannon barrel. The projectile travels beneath a spray nozzle, which impinges a liquid jet onto its metal upper surface, and the projectile then hits a stopping mechanism. A camera records the jet impingement, and a pressure transducer records the spray nozzle backpressure. The spinning disk set-up consists of a steel disk that reaches speeds of 500 to 3,000 rpm via a variable frequency drive (VFD) motor. A spray system similar to that of the air cannon generates a liquid jet that impinges onto the spinning disc, and cameras placed at several optical access points record the jet impingement. Video recordings of jet impingement processes are recorded and examined to determine whether the outcome of impingement is splash, splatter, or deposition. The apparatuses are the first that involve the high speed impingement of low-Reynolds-number liquid jets on high speed moving surfaces. In addition to its rail industry applications, the described technique may be used for technical and industrial purposes such as steelmaking and may be relevant to high-speed 3D printing.     

Frequent side chain methyl carbon‐oxygen hydrogen bonding in proteins revealed by computational and stereochemical analysis of neutron structures

The propensity of backbone Cα atoms to engage in carbon‐oxygen (CH···O) hydrogen bonding is well‐appreciated in protein structure, but side chain CH···O hydrogen bonding remains largely uncharacterized. The extent to which side chain methyl groups in proteins participate in CH···O hydrogen bonding is examined through a survey of neutron crystal structures, quantum chemistry calculations, and molecular dynamics simulations. Using these approaches, methyl groups were observed to form stabilizing CH···O hydrogen bonds within protein structure that are maintained through protein dynamics and participate in correlated motion. Collectively, these findings illustrate that side chain methyl CH···O hydrogen bonding contributes to the energetics of protein structure and folding. Proteins 2015; 83:403–410. © 2014 Wiley Periodicals, Inc.     

Velocity of myosin-based actin sliding depends on attachment and detachment kinetics and reaches a maximum when myosin-binding sites on actin saturate

Molecular motors such as kinesin and myosin often work in groups to generate the directed movements and forces critical for many biological processes. Although much is known about how individual motors generate force and movement, surprisingly, little is known about the mechanisms underlying the macroscopic mechanics generated by multiple motors. For example, the observation that a saturating number, N, of myosin heads move an actin filament at a rate that is influenced by actin–myosin attachment and detachment kinetics is accounted for neither experimentally nor theoretically. To better understand the emergent mechanics of actin–myosin mechanochemistry, we use an in vitro motility assay to measure and correlate the N-dependence of actin sliding velocities, actin-activated ATPase activity, force generation against a mechanical load, and the calcium sensitivity of thin filament velocities. Our results show that both velocity and ATPase activity are strain dependent and that velocity becomes maximized with the saturation of myosin-binding sites on actin at a value that is 40% dependent on attachment kinetics and 60% dependent on detachment kinetics. These results support a chemical thermodynamic model for ensemble motor mechanochemistry and imply molecularly explicit mechanisms within this framework, challenging the assumption of independent force generation.     

Modeling of inhibitor–metalloenzyme interactions and selectivity using molecular mechanics grounded in quantum chemistry

We investigated the binding properties of the metalloprotease inhibitors hydroxamate, methanethiolate, and methylphosphoramidate to a model coordination site occurring in several Zn2+ metalloproteases, including thermolysin. This was carried out using both the SIBFA (sum of interactions between fragments ab initio-computed) molecular mechanics and the SCF/MP2 procedures for the purpose of evaluating SIBFA as a metalloenzyme modeling tool. The energy-minimized structures were closely similar to the X-ray crystallographic structures of related thermolysin-inhibitor complexes. We found that selectivity between alternative geometries and between inhibitors usually stemmed from multiple interaction components included in SIBFA. The binding strength sequence is hydroxamate > methanethiolate ≥ methylphosphoramidate from multiple interaction components included in SIBFA. The trends in interaction energy components, rankings, and preferences for mono- or bidentate binding were consistent in both computational procedures. We also compared the Zn2+ vs. Mg2+ selectivities in several other polycoordinated sites having various “hard” and “soft” qualities. This included a hexahydrate, a model representing Mg2+/Ca2+ binding sites, a chlorophyll-like structure, and a zinc finger model. The latter three favor Zn2+ over Mg2+ by a greater degree than the hydrated state, but the selectivity varies widely according to the ligand “softness.” SIBFA was able to match the ab initio binding energies by <2%, with the SIBFA terms representing dispersion and charge-transfer contributing the most to Zn2+/Mg2+ selectivity. These results showed this procedure to be a very capable modeling tool for metalloenzyme problems, in this case giving valuable information about details and limitations of “hard” and “soft” selectivity trends. Proteins 31:42–60, 1998. © 1998 Wiley-Liss, Inc.     

蓖麻毒素A链突变体(MRTA)的分子设计

利用同源模建的方法,借助分子力学优化、分子动力学模拟退火设计构建了删除部分氨基酸序列的蓖麻毒素Ａ链突变体（ＭＲＴＡ）。采用泊松—玻尔兹曼方程对比分析了蓖麻毒素Ａ链（ＲＴＡ）与ＭＲＴＡ表观静电势分布,研究ＲＴＡ与ＭＲＴＡ蛋白表面静电性质;通过半经验量子化学ＡＭ１与分子力学结合方法探讨ＲＴＡ与ＭＲＴＡ功能域氨基酸前线分子轨道性质、能级分布,从理论上预测ＭＲＴＡ功能活性     

活髓牙修复的实验室探讨

龋病修复体的长寿耐用不仅要求洞型兼顾抗力型和固位型,还要考虑修复体对牙体组织的力学影响。本实验模拟活髓牙下颌第一磨牙颌面Ｉ类洞采用三维有限元法建立数值模型,用八节点等参三维实体单元进行离散,利用ＳＡＰ８４（Ｖ４．２）程序。计算并作一系列力学分析寻求活髓牙修复时垫底材料与修复材料最佳厚度比的力学模型。实验结果表明应力曲线的峰值主要出现在不同材料交界面附近,因此该处的修复应尽量采用圆弧过度。另外基底材