Mapping structural diversity in networks sharing a given degree distribution and global clustering: Adaptive resolution grid search evolution with Diophantine equation-based mutations

Methods that generate networks sharing a given degree distribution and global clustering can induce changes in structural properties other than that controlled for. Diversity in structural properties, in turn, can affect the outcomes of dynamical processes operating on those networks. Since exhaustive sampling is not possible, we propose a novel evolutionary framework for mapping this structural diversity. The three main features of this framework are: (a) subgraph-based encoding of networks, (b) exact mutations based on solving systems of Diophantine equations, and (c) heuristic diversity-driven mechanism to drive resolution changes in the MapElite algorithm.We show that our framework can elicit networks with diversity in their higher-order structure and that this diversity affects the behaviour of the complex contagion model. Through a comparison with state of the art clustered network generation methods, we demonstrate that our approach can uncover a comparably diverse range of networks without needing computationally unfeasible mixing times. Further, we suggest that the subgraph-based encoding provides greater confidence in the diversity of higher-order network structure for low numbers of samples and is the basis for explaining our results with complex contagion model. We believe that this framework could be applied to other complex landscapes that cannot be practically mapped via exhaustive sampling

Crystallographic analyses of an active HIV-1 ribonuclease H domain show structural features that distinguish it from the inactive form

. An active recombinant preparation of the carboxy-terminal ribonuclease H (RNase H) domain of HIV-I reverse transcriptase has produced crystals of several different forms, including a trigonal prism form (P3(1); a = b = 52.03, c = 113.9 A with two molecules per asymmetric unit) and a hexagonal tablet form (P6(2)22 or P6(4)22; a = b = 93.5, c = 74.1 A with one molecule per asymmetric unit). The former appears to be isomorphous with crystals of a similar, but inactive, version of the enzyme that was used for a prior crystal structure determination [Davies, Hostomska, Hostomsky, Jordan & Matthews (1991). Science, 252, 88-95]. We have also obtained a structure solution for this crystal form and have refined it with 2.8 A resolution data (R = 0.216). We report here details of our crystallization studies and some initial structural results that verify that the preparation of active HIV-1 RNase H yields a protein that is not just enzymatically, but also structurally, distinguishable from the inactive form. Evidence suggests that region 538-542, which may be involved in the catalytic site and which is disordered in both molecules in the prior structure determination, is ordered in the crystal structure of the active enzyme, although the ordering may include more than one conformation for this loop. It should also be noted that, in the crystal structure of the trigonal form, RNase H monomers associate to form noncrystallographic twofold-symmetric dimers by fusing five-stranded mixed beta sheets into a single ten-stranded dimerwide sheet, an assembly that was not remarked upon by previous investigators

Electronic damage in S atoms in a native protein crystal induced by an intense X-ray free-electron laser pulse

Current hard X-ray free-electron laser (XFEL) sources can deliver doses to biological macromolecules well exceeding 1 GGy, in timescales of a few tens of femtoseconds. During the pulse, photoionization can reach the point of saturation in which certain atomic species in the sample lose most of their electrons. This electronic radiation damage causes the atomic scattering factors to change, affecting, in particular, the heavy atoms, due to their higher photoabsorption cross sections. Here, it is shown that experimental serial femtosecond crystallography data collected with an extremely bright XFEL source exhibit a reduction of the effective scattering power of the sulfur atoms in a native protein. Quantitative methods are developed to retrieve information on the effective ionization of the damaged atomic species from experimental data, and the implications of utilizing new phasing methods which can take advantage of this localized radiation damage are discussed

Catalytic cleavage of HEAT and subsequent covalent binding of the tetralone moiety by the SARS-CoV-2 main protease

Here we present the crystal structure of SARS-CoV-2 main protease (Mpro) covalently bound to 2-methyl-1-tetralone. This complex was obtained by co-crystallization of Mpro with HEAT (2-(((4-hydroxyphenethyl)amino)methyl)-3,4-dihydronaphthalen-1(2H)-one) in the framework of a large X-ray crystallographic screening project of Mpro against a drug repurposing library, consisting of 5632 approved drugs or compounds in clinical phase trials. Further investigations showed that HEAT is cleaved by Mpro in an E1cB-like reaction mechanism into 2-methylene-1-tetralone and tyramine. The catalytic Cys145 subsequently binds covalently in a Michael addition to the methylene carbon atom of 2-methylene-1-tetralone. According to this postulated model HEAT is acting in a pro-drug-like fashion. It is metabolized by Mpro, followed by covalent binding of one metabolite to the active site. The structure of the covalent adduct elucidated in this study opens up a new path for developing non-peptidic inhibitors

SARS-CoV-2 M^pro responds to oxidation by forming disulfide and NOS/SONOS bonds

The main protease (Mpro) of SARS-CoV-2 is critical for viral function and a key drug target. Mpro is only active when reduced; turnover ceases upon oxidation but is restored by re-reduction. This suggests the system has evolved to survive periods in an oxidative environment, but the mechanism of this protection has not been confirmed. Here, we report a crystal structure of oxidized Mpro showing a disulfide bond between the active site cysteine, C145, and a distal cysteine, C117. Previous work proposed this disulfide provides the mechanism of protection from irreversible oxidation. Mpro forms an obligate homodimer, and the C117-C145 structure shows disruption of interactions bridging the dimer interface, implying a correlation between oxidation and dimerization. We confirm dimer stability is weakened in solution upon oxidation. Finally, we observe the protein’s crystallization behavior is linked to its redox state. Oxidized Mpro spontaneously forms a distinct, more loosely packed lattice. Seeding with crystals of this lattice yields a structure with an oxidation pattern incorporating one cysteine-lysine-cysteine (SONOS) and two lysine-cysteine (NOS) bridges. These structures further our understanding of the oxidative regulation of Mpro and the crystallization conditions necessary to study this structurally

QSIPrep: an integrative platform for preprocessing and reconstructing diffusion MRI data

Diffusion-weighted magnetic resonance imaging (dMRI) is the primary method for noninvasively studying the organization of white matter in the human brain. Here we introduce QSIPrep, an integrative software platform for the processing of diffusion images that is compatible with nearly all dMRI sampling schemes. Drawing on a diverse set of software suites to capitalize on their complementary strengths, QSIPrep facilitates the implementation of best practices for processing of diffusion images

Frustrated hierarchical synchronization and emergent complexity in the human connectome network

The spontaneous emergence of coherent behavior through synchronization plays a key role in neural function, and its anomalies often lie at the basis of pathologies. Here we employ a parsimonious (mesoscopic) approach to study analytically and computationally the synchronization (Kuramoto) dynamics on the actual human-brain connectome network. We elucidate the existence of a so-far-uncovered intermediate phase, placed between the standard synchronous and asynchronous phases, i.e. between order and disorder. This novel phase stems from the hierarchical modular organization of the connectome. Where one would expect a hierarchical synchronization process, we show that the interplay between structural bottlenecks and quenched intrinsic frequency heterogeneities at many different scales, gives rise to frustrated synchronization, metastability, and chimera-like states, resulting in a very rich and complex phenomenology. We uncover the origin of the dynamic freezing behind these features by using spectral graph theory and discuss how the emerging complex synchronization patterns relate to the need for the brain to access –in a robust though flexible way– a large variety of functional attractors and dynamical repertoires without ad hoc fine-tuning to a critical pointWe acknowledge financial support from J. de Andalucía, grant P09-FQM-4682 and we thank O. Sporns for providing us access to the human connectome data

Monitoring and Scoring Counter-Diffusion Protein Crystallization Experiments in Capillaries by in situ Dynamic Light Scattering

In this paper, we demonstrate the feasibility of using in situ Dynamic Light Scattering (DLS) to monitor counter-diffusion crystallization experiments in capillaries. Firstly, we have validated the quality of the DLS signal in thin capillaries, which is comparable to that obtained in standard quartz cuvettes. Then, we have carried out DLS measurements of a counter-diffusion crystallization experiment of glucose isomerase in capillaries of different diameters (0.1, 0.2 and 0.3 mm) in order to follow the temporal evolution of protein supersaturation. Finally, we have compared DLS data with optical recordings of the progression of the crystallization front and with a simulation model of counter-diffusion in 1D