Level Structure and Properties of Open f-Shell Elements

Open f-shell elements still constitute a great challenge for atomic theory owing to their (very) rich fine-structure and strong correlations among the valence-shell electrons. For these medium and heavy elements, many atomic properties are sensitive to the correlated motion of electrons and, hence, require large-scale computations in order to deal consistently with all relativistic, correlation and rearrangement contributions to the electron density. Often, different concepts and notations need to be combined for just classifying the low-lying level structure of these elements. With Jac, the Jena Atomic Calculator, we here provide a toolbox that helps to explore and deal with such elements with open d- and f-shell structures. Based on Dirac’s equation, Jac is suitable for almost all atoms and ions across the periodic table. As an example, we demonstrate how reasonably accurate computations can be performed for the low-lying level structure, transition probabilities and lifetimes for Th2+ ions with a 5f6d ground configuration. Other, and more complex, shell structures are supported as well, though often for a trade-off between the size and accuracy of the computations. Owing to its simple use, however, Jac supports both quick estimates and detailed case studies on open d- or f-shell elements.

Atomic resolution Protein Allostery from the multi-state Structure of a PDZ Domain

Recent methodological advances in solution NMR allow the determination of multi-state protein structures and provide insights into correlated motion at atomic resolution as demonstrated here for the well-studied PDZ2 domain of protein human tyrosine phosphatase 1E for which protein allostery was predicted. Two-state protein structures were calculated for both the free form and in complex with the RA-GEF2 peptide using the exact nuclear Overhauser effect (eNOE) method. In the apo protein states an allosteric conformational preselection step comprising almost 60% of the domain was detected with an "open" ligand welcoming state and a "closed" state that obstructs the binding site by the distance between the β-sheet, α-helix 2 and sidechains of residues Lys38 and Lys72. Observed apo-holo structural rearrangements of induced fit-type are in line with previously published evolution-based analysis covering ~25% of the domain with only a partial overlap with the protein allostery of the open form. These presented structural studies highlight the presence of a dedicated highly optimized dynamic interplay of the complexity of the PDZ2 domain owed by the structure-dynamics landscape.

PIP2-dependent coupling of voltage sensor and pore domains in Kv7.2 channel

Phosphatidylinositol-4,5-bisphosphate (PIP2) is a signaling lipid which regulates voltage-gated Kv7/KCNQ potassium channels. Altered PIP2 sensitivity of neuronal Kv7.2 channel is involved in KCNQ2 epileptic encephalopathy. However, the molecular action of PIP2 on Kv7.2 gating remains largely elusive. Here, we use molecular dynamics simulations and electrophysiology to characterize PIP2 binding sites in a human Kv7.2 channel. In the closed state, PIP2 localizes to the periphery of the voltage-sensing domain (VSD). In the open state, PIP2 binds to 4 distinct interfaces formed by the cytoplasmic ends of the VSD, the gate, intracellular helices A and B and their linkers. PIP2 binding induces bilayer-interacting conformation of helices A and B and the correlated motion of the VSD and the pore domain, whereas charge-neutralizing mutations block this coupling and reduce PIP2 sensitivity of Kv7.2 channels by disrupting PIP2 binding. These findings reveal the allosteric role of PIP2 in Kv7.2 channel activation.

Mutations in SARS-CoV-2 Variants Modulate the Microscopic Dynamics of Neutralizing Antibodies

Monoclonal antibodies have emerged as viable treatment for the COVID-19 disease caused by the SARS-CoV-2 virus. But the new viral variants can reduce the efficacy of the currently available antibodies, as well as diminish the vaccine induced immunity. Here, we demonstrate how the microscopic dynamics of the SARS-CoV-2 neutralizing monoclonal antibodies, can be modulated by the mutations present in the spike proteins of the variants currently circulating in the world population. We show that the dynamical perturbation in the antibody structure can be diverse, depending both on the nature of the antibody and on the location of the mutation. The correlated motion between the antibody and the receptor binding domain (RBD) can also be changed, altering the binding affinity. By constructing a protein graph connectivity network, we could delineate the mutant induced modifications in the allosteric information flow pathway through the antibody, and observed the presence of both localized and long distance effects. We identified a loop consisting of residues 470-490 in the RBD which works like an anchor preventing the detachment of the antibodies, and individual mutations in that region can significantly affect the antibody binding propensity. Our study provides fundamental and atomistically detailed insight on how virus neutralization by monoclonal antibody can be impacted by the mutations in the epitope, and can potentially facilitate the rational design of monoclonal antibodies, effective against the new variants of the novel coronavirus.

X-ray free-electron laser studies reveal correlated motion during isopenicillin N synthase catalysis

Isopenicillin N synthase (IPNS) catalyzes the unique reaction of l-δ-(α-aminoadipoyl)-l-cysteinyl-d-valine (ACV) with dioxygen giving isopenicillin N (IPN), the precursor of all natural penicillins and cephalosporins. X-ray free-electron laser studies including time-resolved crystallography and emission spectroscopy reveal how reaction of IPNS:Fe(II):ACV with dioxygen to yield an Fe(III) superoxide causes differences in active site volume and unexpected conformational changes that propagate to structurally remote regions. Combined with solution studies, the results reveal the importance of protein dynamics in regulating intermediate conformations during conversion of ACV to IPN. The results have implications for catalysis by multiple IPNS-related oxygenases, including those involved in the human hypoxic response, and highlight the power of serial femtosecond crystallography to provide insight into long-range enzyme dynamics during reactions presently impossible for nonprotein catalysts.

Getting the Nod: Characterizing pediatric head motion in movie- and resting-state fMRI.

Head motion continues to be a major problem in fMRI research, particularly in developmental studies where an inverse relationship exists between head motion and age. Despite multifaceted and costly efforts to mitigate motion and motion-related signal artifact, few studies have characterized in-scanner head motion itself. This study leverages a large public dataset (N=1388, age 5-21y, The Healthy Brain Network Biobank) to characterize pediatric head motion in space, frequency, and time. We focus on practical aspects of head motion that could impact future study design, including comparing motion across groups (low, medium, and high movers), across conditions (movie-watching and rest), and between males and females. Analyses showed that in all conditions, high movers exhibited a different pattern of motion than low and medium movers that was dominated by x-rotation, and z- and y-translation. High motion spikes (>0.3mm) from all participants also showed this pitch-z-y pattern. Problematic head motion is thus composed of a single type of biomechanical motion, which we infer to be a nodding movement, providing a focused target for motion reduction strategies. A second type of motion was evident via spectral analysis of raw displacement data. This was observed in low and medium movers and was consistent with respiration rates. We consider this to be a baseline of motion best targeted in data preprocessing. Further, we found that males moved more than, but not differently from, females. Significant cross-condition differences in head motion were found. Movies had lower mean motion, and especially in high movers, movie-watching reduced within-run linear increases in head motion (i.e., temporal drift). Finally, we used intersubject correlations of framewise displacement (FD-ISCs) to assess for stimulus-correlated motion trends. Subject motion was more correlated in movie than rest and stimulus-correlated stillness occurred more often than stimulus-correlated motion. Possible reasons and future implications of these findings are discussed.

PIP2-dependent coupling of voltage sensor and pore domains in Kv7.2

Phosphatidylinositol-4,5-bisphosphate (PIP2) is a signaling lipid which regulates voltage-gated Kv7/KCNQ potassium channels. Altered PIP2 sensitivity of neuronal Kv7.2 channel is involved in KCNQ2 epileptic encephalopathy. However, the molecular action of PIP2 on Kv7.2 gating remains largely elusive. Here, we use molecular dynamics simulations and electrophysiology to characterize PIP2 binding sites in a human Kv7.2 channel. In the closed state, PIP2 localizes to the periphery of the voltage-sensing domain (VSD). In the open state, PIP2 binds to 4 distinct interfaces formed by the cytoplasmic ends of the VSD, the gate, intracellular helices A and B and their linkers. PIP2 binding induces bilayer-interacting conformation of helices A and B and the correlated motion of the VSD and the pore domain, whereas charge-neutralizing mutations block this coupling and reduce PIP2 sensitivity of Kv7.2 channels by disrupting PIP2 binding. These findings reveal the allosteric role of PIP2 in Kv7.2 channel activation.

All-Atom MD Simulations of the HBV Capsid Complexed with AT130 Reveal Secondary and Tertiary Structural Changes and Mechanisms of Allostery

The hepatitis B virus (HBV) capsid is an attractive drug target, relevant to combating viral hepatitis as a major public health concern. Among small molecules known to interfere with capsid assembly, the phenylpropenamides, including AT130, represent an important antiviral paradigm based on disrupting the timing of genome packaging. Here, all-atom molecular dynamics simulations of an intact AT130-bound HBV capsid reveal that the compound increases spike flexibility and improves recovery of helical secondary structure in the spike tips. Regions of the capsid-incorporated dimer that undergo correlated motion correspond to established sub-domains that pivot around the central chassis. AT130 alters patterns of correlated motion and other essential dynamics. A new conformational state of the dimer is identified, which can lead to dramatic opening of the intradimer interface and disruption of communication within the spike tip. A novel salt bridge is also discovered, which can mediate contact between the spike tip and fulcrum even in closed conformations, revealing a mechanism of direct communication across these sub-domains. Altogether, results describe a dynamical connection between the intra- and interdimer interfaces and enable mapping of allostery traversing the entire core protein dimer.

Chromatin Network Retards Droplet Coalescence

Nuclear bodies are membraneless condensates that may form via liquid-liquid phase separation. The viscoelastic chromatin network could impact their stability and may hold the key for understanding experimental observations that defy predictions of classical theories. However, quantitative studies on the role of the chromatin network in phase separation have remained challenging. Using a diploid human genome model parameterized with chromosome conformation capture (Hi-C) data, we studied the thermodynamics and kinetics of droplet formation inside the nucleus. Dynamical simulations predicted the formation of multiple droplets for protein particles that experience specific interactions with nucleolus-associated domains (NADs). Coarsening dynamics, surface tension, and coalescence kinetics of the simulated droplets are all in quantitative agreements with experimental measurements for nucleoli. Free energy calculations further supported that a two-droplet state, which is often observed for nucleoli seen in somatic cells, is metastable and separated from the single-droplet state with an entropic barrier. Our study suggests that protein-chromatin interactions facilitate the nucleation of droplets, but hinders their coarsening due to the correlated motion between droplets and the chromatin network: as droplets coalesce, the chromatin network becomes increasingly constrained. Therefore, protein-chromatin interactions arrest phase separation in multi-droplet states and may drive the variation of nuclear body numbers across cell types.

Exploring the Association Between Sialic Acid and SARS-CoV-2 Spike Protein Through a Molecular Dynamics-Based Approach

Recent experimental evidence demonstrated the capability of SARS-CoV-2 Spike protein to bind sialic acid molecules, which was a trait not present in SARS-CoV and could shed light on the molecular mechanism used by the virus for the cell invasion. This peculiar feature has been successfully predicted by in-silico studies comparing the sequence and structural characteristics that SARS-CoV-2 shares with other sialic acid-binding viruses, like MERS-CoV. Even if the region of the binding has been identified in the N-terminal domain of Spike protein, so far no comprehensive analyses have been carried out on the spike-sialic acid conformations once in the complex. Here, we addressed this aspect performing an extensive molecular dynamics simulation of a system composed of the N-terminal domain of the spike protein and a sialic acid molecule. We observed several short-lived binding events, reconnecting to the avidic nature of the binding, interestingly occurring in the surface Spike region where several insertions are present with respect to the SARS-CoV sequence. Characterizing the bound configurations via a clustering analysis on the Principal Component of the motion, we identified different possible binding conformations and discussed their dynamic and structural properties. In particular, we analyze the correlated motion between the binding residues and the binding effect on the stability of atomic fluctuation, thus proposing regions with high binding propensity with sialic acid.