Proximity to criticality predicts surface properties of biomolecular condensates

It has recently become appreciated that cells self-organize their interiors through the formation of biomolecular condensates. These condensates, typically formed through liquid-liquid phase separation of proteins, nucleic acids, and other biopolymers, exhibit reversible assembly/disassembly in response to changing conditions. Condensates play many functional roles, aiding in biochemical reactions, signal transduction, and sequestration of certain components. Ultimately, these functions depend on the physical properties of condensates, which are encoded in the microscopic features of the constituent biomolecules. In general, the mapping from microscopic features to macroscopic properties is complex, but it is known that near a critical point, macroscopic properties follow power laws with only a small number of parameters, making it easier to identify underlying principles. How far does this critical region extend for biomolecular condensates and what principles govern condensate properties in the critical regime? Using coarse-grained molecular-dynamics simulations of a representative class of biomolecular condensates, we found that the critical regime can be wide enough to cover the full physiological range of temperatures. Within this critical regime, we identified that polymer sequence influences surface tension predominately via shifting the critical temperature. Finally, we show that condensate surface tension over a wide range of temperatures can be calculated from the critical temperature and a single measurement of the interface width.

Modification of a Single Atom Affects the Physical Properties of Double Fluorinated Fmoc-Phe Derivatives

Supramolecular hydrogels formed by the self-assembly of amino-acid based gelators are receiving increasing attention from the fields of biomedicine and material science. Self-assembled systems exhibit well-ordered functional architectures and unique physicochemical properties. However, the control over the kinetics and mechanical properties of the end-products remains puzzling. A minimal alteration of the chemical environment could cause a significant impact. In this context, we report the effects of modifying the position of a single atom on the properties and kinetics of the self-assembly process. A combination of experimental and computational methods, used to investigate double-fluorinated Fmoc-Phe derivatives, Fmoc-3,4F-Phe and Fmoc-3,5F-Phe, reveals the unique effects of modifying the position of a single fluorine on the self-assembly process, and the physical properties of the product. The presence of significant physical and morphological differences between the two derivatives was verified by molecular-dynamics simulations. Analysis of the spontaneous phase-transition of both building blocks, as well as crystal X-ray diffraction to determine the molecular structure of Fmoc-3,4F-Phe, are in good agreement with known changes in the Phe fluorination pattern and highlight the effect of a single atom position on the self-assembly process. These findings prove that fluorination is an effective strategy to influence supramolecular organization on the nanoscale. Moreover, we believe that a deep understanding of the self-assembly process may provide fundamental insights that will facilitate the development of optimal amino-acid-based low-molecular-weight hydrogelators for a wide range of applications.

A Molecular Lid Mechanism of K+ Channel Blocker Action Revealed by a Cone Peptide

Many venomous organisms carry in their arsenal short polypeptides that block K⁺ channels in a highly selective manner. These toxins may compete with the permeating ions directly via a "plug" mechanism or indirectly via a "pore-collapse" mechanism. An alternative "lid" mechanism was proposed but remained poorly defined. Here we study the Drosophila Shaker channel block by Conkunitzin-S1 and Conkunitzin-C3, two highly similar toxins derived from cone venom. Despite their similarity, the two peptides exhibited differences in their binding poses and biophysical assays, implying discrete action modes. We show that while Conkunitzin-S1 binds tightly to the channel turret and acts via a "pore-collapse" mechanism, Conkunitzin-C3 does not contact this region. Instead, Conk-C3 uses a non-conserved Arg to divert the permeant ions and trap them in off-axis cryptic sites above the SF, a mechanism we term a "molecular-lid". Our study provides an atomic description of the "lid" K⁺ blocking mode and offers valuable insights for the design of therapeutics based on venom peptides.

Oncogenic potential of ATAD2 in stomach cancer and insights into the protein-protein interactions at its AAA + ATPase domain and bromodomain

ATAD2 has recently been shown to promote stomach cancer. However, nothing is known about the functional network of ATAD2 in stomach carcinogenesis. This study illustrates the oncogenic potential of ATAD2 and the participation of its ATPase and bromodomain in stomach malignancy. Expression of ATAD2 in stomach cancer is analyzed by in silico and in vitro techniques including western blot and immunofluorescence microscopy of stomach cancer cells (SCCs) and tissues. The oncogenic potential of ATAD2 is examined thoroughly using genetic alterations, driver gene prediction, survival analysis, identification of interacting partners, and analysis of canonical pathways. To understand the protein-protein interactions (PPI) at residue level, molecular docking and molecular dynamics simulations (1200 ns) are performed. Enhanced expression of ATAD2 is observed in H. pylori-infected SCCs, patient biopsy tissues, and all stages and grades of stomach cancer. High expression of ATAD2 is found to be negatively correlated with the survival of stomach cancer patients. ATAD2 is a cancer driver gene with 37 mutational sites and a predictable factor for stomach cancer prognosis with high accuracy. The top canonical pathways of ATAD2 indicate its participation in stomach malignancy. The ATAD2-PPI in stomach cancer identify top-ranked partners; ESR1, SUMO2, SPTN2, and MYC show preference for the bromodomain whereas NCOA3 and HDA11 have preference for the ATPase domain of ATAD2. The oncogenic characterization of ATAD2 provides strong evidence to consider ATAD2 as a stomach cancer biomarker. These studies offer an insight for the first time into the ATAD2-PPI interface presenting a novel target for cancer therapeutics. Communicated by Ramaswamy H. Sarma.

Thermodynamic, structural and dynamic properties of selected non-associative neat liquids

Non-associative neat liquids benzene, acetone and carbon tetrachloride have been examined via molecular dynamics simulations. Several models of each neat liquid have been simulated, and selected thermodynamic and structural results are presented. However, the models have been compared mainly in terms of their dynamic quantities. Since models are rarely parametrized with the dynamic properties in mind, the principal goal of this work is to present quantities such as the power spectra, rotational correlation functions and relaxation times, diffusion coefficients and self and distinct parts of the van Hove functions in relation to available experimental data. The general trends of the calculated data provide a benchmark for the behavior of neat simple liquids which will be built upon in the cases of mixtures with associative liquids.

Self-Assembly, Adaptive Response, and in,out-Stereoisomerism of Large Orthoformate Cryptands

We report on triethylene glycol-based orthoformate cryptands, which adapt their bridgehead configurations in response to metal templates and intramolecular hydrogen bonding in a complex manner. In contrast to smaller 1.1.1-orthoformate cryptands, the inversion from out,out-2.2.2 to in,in-2.2.2 occurs spontaneously by thermal homeomorphic isomerization, i. e., without bond breakage. The global thermodynamic minimum of the entire network, which includes an unprecedented third isomer (in,out-2.2.2), could only be reached under conditions that facilitate dynamic covalent exchange. Both inversion processes were studied in detail, including DFT calculations and MD simulations, which were particularly helpful for explaining differences between equilibrium compositions in solvents chloroform and acetonitrile. Unexpectedly, the system could be driven to the in,out-2.2.2 state by using a metal template with a size mismatch with respect to the out,out-2.2.2 cage.

Unraveling the mechanism of selective ion transport in hydrophobic subnanometer channels

Recently reported synthetic organic nanopore (SONP) can mimic a key feature of natural ion channels, i.e., selective ion transport. However, the physical mechanism underlying the K(+)/Na(+) selectivity for the SONPs is dramatically different from that of natural ion channels. To achieve a better understanding of the selective ion transport in hydrophobic subnanometer channels in general and SONPs in particular, we perform a series of ab initio molecular dynamics simulations to investigate the diffusivity of aqua Na(+) and K(+) ions in two prototype hydrophobic nanochannels: (i) an SONP with radius of 3.2 Å, and (ii) single-walled carbon nanotubes (CNTs) with radii of 3-5 Å (these radii are comparable to those of the biological potassium K(+) channels). We find that the hydration shell of aqua Na(+) ion is smaller than that of aqua K(+) ion but notably more structured and less yielding. The aqua ions do not lower the diffusivity of water molecules in CNTs, but in SONP the diffusivity of aqua ions (Na(+) in particular) is strongly suppressed due to the rugged inner surface. Moreover, the aqua Na(+) ion requires higher formation energy than aqua K(+) ion in the hydrophobic nanochannels. As such, we find that the ion (K(+) vs. Na(+)) selectivity of the (8, 8) CNT is ∼20× higher than that of SONP. Hence, the (8, 8) CNT is likely the most efficient artificial K(+) channel due in part to its special interior environment in which Na(+) can be fully solvated, whereas K(+) cannot. This work provides deeper insights into the physical chemistry behind selective ion transport in nanochannels.

Computing membrane-AQP5-phosphatidylserine binding affinities with hybrid steered molecular dynamics approach

In order to elucidate how phosphatidylserine (PS6) interacts with AQP5 in a cell membrane, we developed a hybrid steered molecular dynamics (hSMD) method that involved: (1) Simultaneously steering two centers of mass of two selected segments of the ligand, and (2) equilibrating the ligand-protein complex with and without biasing the system. Validating hSMD, we first studied vascular endothelial growth factor receptor 1 (VEGFR1) in complex with N-(4-Chlorophenyl)-2-((pyridin-4-ylmethyl)amino)benzamide (8ST), for which the binding energy is known from in vitro experiments. In this study, our computed binding energy well agreed with the experimental value. Knowing the accuracy of this hSMD method, we applied it to the AQP5-lipid-bilayer system to answer an outstanding question relevant to AQP5's physiological function: Will the PS6, a lipid having a single long hydrocarbon tail that was found in the central pore of the AQP5 tetramer crystal, actually bind to and inhibit AQP5's central pore under near-physiological conditions, namely, when AQP5 tetramer is embedded in a lipid bilayer? We found, in silico, using the CHARMM 36 force field, that binding PS6 to AQP5 was a factor of 3 million weaker than "binding" it in the lipid bilayer. This suggests that AQP5's central pore will not be inhibited by PS6 or a similar lipid in a physiological environment.