Out of the Randomness: Correlating Noise in Biological Systems

The study of the dynamics of biological systems requires one to follow relaxation processes in time with micron-size spatial resolution. This need has led to the development of different fluorescence correlation techniques with high spatial resolution and a tremendous (from nanoseconds to seconds) temporal dynamic range. Spatiotemporal information can be obtained even on complex dynamic processes whose time evolution is not forecast by simple Brownian diffusion. Our discussion of the most recent applications of image correlation spectroscopy to the study of anomalous sub- or superdiffusion suggests that this field still requires the development of multidimensional image analyses based on analytical models or numerical simulations. We focus in particular on the framework of spatiotemporal image correlation spectroscopy and examine the critical steps in getting information on anomalous diffusive processes from the correlation maps. We point out how a dual space-time correlative analysis, in both the direct and the Fourier space, can provide quantitative information on superdiffusional processes when these are analyzed through an empirical model based on intermittent active dynamics. We believe that this dual space-time analysis, potentially amenable to mathematical treatment and to the exact fit of experimental data, could be extended to include the rich phenomenology of subdiffusive processes, thereby quantifying relevant parameters for the various motivating biological problems of interest.     

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.     

The Competition between the Noise and Shear Motion Sensitivity of Cochlear Inner Hair Cell Stereocilia

Acoustical excitation of the organ of Corti induces radial fluid flow in the subtectorial space (STS) that excites the hair bundles (HBs) of the sensory inner hair cell of the mammalian cochlea. The inner hair cell HBs are bathed in endolymphatic fluid filling a thin gap in the STS between the tectorial membrane and the reticular lamina. According to the fluctuation dissipation theorem, the fluid viscosity gives rise to mechanical fluctuations that are transduced into current noise. Conversely, the stochastic fluctuations of the mechanically gated channels of the HBs also induce dissipation. We develop an analytic model of the STS complex in a cross section of the gerbil organ of Corti. We predict that the dominant noise at the apex is due to the channel stochasticity whereas viscous effects dominate at the base. The net root mean square fluctuation of the HB motion is estimated to be at least 1.18 nm at the base and 2.72 nm at the apex. By varying the HB height for a fixed STS gap, we find that taller HBs are better sensors with lower thresholds. An integrated active HB model is shown to reduce the hydrodynamic resistance through a cycle-by-cycle power addition through adaptation, reducing the thresholds of hearing, hinting at one potential role for HB activity in mammalian hearing. We determine that a Couette flow approximation in the STS underestimates the dissipation and that modeling the entire STS complex is necessary to correctly predict the low-frequency dissipation in the cochlea. Finally, the difference in the noise budget at the base and the apex of the cochlea indicate that a sensing modality other than the shear motion of the TM that may be used to achieve low-noise acoustic sensing at the apex.     

Quantifying Molecular Forces with Serially Connected Force Sensors

To understand the mechanical forces involved in cell adhesion, molecular force sensors have been developed to study tension through adhesion proteins. Recently, a class of molecular force sensors called tension gauge tethers (TGTs) have been developed that rely on irreversible force-dependent dissociation of a DNA duplex to study cell adhesion forces. Although the TGT offers a high signal-to-noise ratio and is ideal for studying fast/single-molecular adhesion processes, quantitative interpretation of experimental results has been challenging. Here, we use a computational approach to investigate how TGT fluorescence readout can be quantitatively interpreted. In particular, we studied force sensors made of a single TGT, multiplexed single TGTs, and two TGTs connected in series. Our results showed that fluorescence readout using a single TGT can result from drastically different combinations of force history and adhesion event density that span orders of magnitude. In addition, the apparent behavior of the TGT is influenced by the tethered receptor-ligand, making it necessary to calibrate the TGT with every new receptor-ligand. To solve this problem, we proposed a system of two serially connected TGTs. Our result shows that not only is the ratiometric readout of serial TGT independent of the choice of receptor-ligand, it is able to reconstruct force history with sub-pN force resolution. This is also not possible by simply multiplexing different types of TGTs together. Last, we systematically investigated how the sequence composition of the two serially connected TGTs can be tuned to achieve different dynamic range. This computational study demonstrated how serially connected irreversible molecular dissociation processes can accurately quantify molecular force and laid the foundation for subsequent experimental studies.     

Molecular Mechanisms of Tight Binding through Fuzzy Interactions

Many intrinsically disordered proteins (IDPs) form fuzzy complexes upon binding to their targets. Although many IDPs are weakly bound in fuzzy complexes, some IDPs form high-affinity complexes. One example is the nonstructural protein 1 (NS1) of the 1918 Spanish influenza A virus, which hijacks cellular CRKII through the strong binding affinity (K_d ～10 nM) of its proline-rich motif (PRM^NS1) to the N-terminal Src-homology 3 domain of CRKII. However, its molecular mechanism remains elusive. Here, we examine the interplay between structural disorder of a bound PRM^NS1 and its long-range electrostatic interactions. Using x-ray crystallography and NMR spectroscopy, we found that PRM^NS1 retains substantial conformational flexibility in the bound state. Moreover, molecular dynamics simulations showed that structural disorder of the bound PRM^NS1 increases the number of electrostatic interactions and decreases the mean distances between the positively charged residues in PRM^NS1 and the acidic residues in the N-terminal Src-homology 3 domain. These results are analyzed using a polyelectrostatic model. Our results provide an insight into the molecular recognition mechanism for a high-affinity fuzzy complex.     

Ultrasound Molecular Imaging with Targeted Microbubbles for Cancer Diagnostics: From Bench to Bedside

Background

Ultrasound plays an important role in cancer diagnosis. B-mode imaging and contrast-enhanced ultrasound are routinely used to detect cancerous lesions in breast and liver. The use of ultrasound contrast agents (UCAs) such as microbubbles (MBs), which can be functionalized with targeting ligands, has further enabled ultrasound molecular imaging (USMI) of specific molecular markers in pre-clinical and the first clinical studies. As targeted MBs have a diameter of 1–4 μm, they are limited to the blood vasculature upon intravenous injection, and can bind to markers of the vascular endothelium. USMI with targeted MBs was applied for imaging of markers of inflammation, angiogenesis, and the tumor endothelium.

Molecular Insights into Variable Electron Transfer in Amphibian Cryptochrome

Cryptochrome proteins are activated by the absorption of blue light, leading to the formation of radical pairs through electron transfer in the active site. Recent experimental studies have shown that once some of the amino acid residues in the active site of Xenopus laevis cryptochrome DASH are mutated, radical-pair formation is still observed. In this study, we computationally investigate electron-transfer pathways in the X. laevis cryptochrome DASH by extensively equilibrating a previously established homology model using molecular dynamics simulations and then mutating key amino acids involved in the electron transfer. The electron-transfer pathways are then probed by using tight-binding density-functional theory. We report the alternative electron-transfer pathways resolved at the molecular level and, through comparison of amino acid sequences for cryptochromes from different species, we demonstrate that one of these alternative electron-transfer pathways could be general for all cryptochrome DASH proteins.     

Mobility of Molecular Motors Regulates Contractile Behaviors of Actin Networks

Cells generate mechanical forces primarily from interactions between F-actin, cross-linking proteins, myosin motors, and other actin-binding proteins in the cytoskeleton. To understand how molecular interactions between the cytoskeletal elements generate forces, a number of in vitro experiments have been performed but are limited in their ability to accurately reproduce the diversity of motor mobility. In myosin motility assays, myosin heads are fixed on a surface and glide F-actin. By contrast, in reconstituted gels, the motion of both myosin and F-actin is unrestricted. Because only these two extreme conditions have been used, the importance of mobility of motors for network behaviors has remained unclear. In this study, to illuminate the impacts of motor mobility on the contractile behaviors of the actin cytoskeleton, we employed an agent-based computational model based on Brownian dynamics. We find that if motors can bind to only one F-actin like myosin I, networks are most contractile at intermediate mobility. In this case, less motor mobility helps motors stably pull F-actins to generate tensile forces, whereas higher motor mobility allows F-actins to aggregate into larger clustering structures. The optimal intermediate motor mobility depends on the stall force and affinity of motors that are regulated by mechanochemical rates. In addition, we find that the role of motor mobility can vary drastically if motors can bind to a pair of F-actins. A network can exhibit large contraction with high motor mobility because motors bound to antiparallel pairs of F-actins can exert similar forces regardless of their mobility. Results from this study imply that the mobility of molecular motors may critically regulate contractile behaviors of actin networks in cells.     

Molecular Insights on the Release of Avibactam from the Acyl-Enzyme Complex

Avibactam is a non-β-lactam β-lactamase inhibitor for treating complicated urinary tract and respiratory infections caused by multidrug-resistant bacterial pathogens, a serious public health threat. Despite its importance, the release mechanism of avibactam from the enzyme-inhibitor complex has been scarcely studied from first principles, considering the total protein environment. This information at the molecular level is essential for the rational design of new antibiotics and inhibitors. In this article, we addressed the release of avibactam from the complex CTX-M-15 by means of molecular dynamics simulations and quantum mechanics/molecular mechanics calculations. This study provides molecular information not available earlier, including exploration of the potential energy surfaces, characterization of the observed intermediate, and their critical points, as well. Our results show that unlike that observed in the acylation reaction, the residues Glu166 and Lys73 would be in their neutral forms. Release of avibactam follows a stepwise mechanism in which the first stage corresponds to the formation of a tetrahedral intermediate, whereas the second stage corresponds to the cleavage of the Ser70-C7 bond, mediated by Lys73, either directly or through Ser130.     

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.