Czech Academy of Sciences, Institute of Biophysics

In 2008 we published the first set of guidelines for standardizing research in autophagy. Since then, research on this topic has continued to accelerate, and many new scientists have entered the field. Our knowledge base and relevant new technologies have also been expanding. Accordingly, it is important to update these guidelines for monitoring autophagy in different organisms. Various reviews have described the range of assays that have been used for this purpose. Nevertheless, there continues to be confusion regarding acceptable methods to measure autophagy, especially in multicellular eukaryotes. For example, a key point that needs to be emphasized is thatthere is a difference between measurements that monitor the numbers or volume of autophagic elements (e.g., autophagosomes or autolysosomes) at any stage of the autophagic process versus those that measure flux through the autophagy pathway (i.e., the completeprocess including the amount and rate of cargo sequestered and degraded). In particular, a block in macroautophagy that results in autophagosome accumulation must be differentiated from stimuli that increase autophagic activity, defined as increasedautophagy induction coupled with increased delivery to, and degradation within, lysosomes (inmost higher eukaryotes and some protists such as Dictyostelium) or the vacuole (in plants and fungi). In other words, it is especially important that investigators new to the field understand that the appearance of more autophagosomes does not necessarily equate with more autophagy. In fact, in manycases, autophagosomes accumulate because of a block in trafficking to lysosomes without a concomitant change in autophagosome biogenesis, whereas an increase in autolysosomes may reflect a reduction in degradative activity. It is worth emphasizing here that lysosomal digestion is a stage of autophagy and evaluating its competence is a crucial part of the evaluation of autophagic flux, or complete autophagy. Here, we present a set of guidelines for the selection and interpretation of methods for use by investigators who aim to examine macroautophagy and related processes, as well as forreviewers who need to provide realistic and reasonable critiques of papers that are focused on these processes. These guidelines are not meant to be a formulaic set of rules, because the appropriate assays depend in part on the question being asked and the system being used. In addition, we emphasize that no individual assay is guaranteed to be the most appropriate one in every situation, and we strongly recommend the use of multipleassays to monitor autophagy. Along these lines, because of the potential for pleiotropic effects due to blocking autophagy through genetic manipulation, it is imperative to target by gene knockout or RNA interference more than one autophagyrelated protein. In addition, some individual Atg proteins, or groups of proteins, are involved in other cellular pathways implying that not all Atg proteins can be used as a specific marker for an autophagic process. In these guidelines, we consider these various methods of assessing autophagy and what information can, or cannot, be obtained from them. Finally, by discussing the merits and limits of particular assays, we hope to encourage technical innovation in the field.

Characterization of the molecular function of the human genome and its variation across individuals is essential for identifying the cellular mechanisms that underlie human genetic traits and diseases. The Genotype-Tissue Expression (GTEx) project aims to characterize variation in gene expression levels across individuals and diverse tissues of the human body, many of which are not easily accessible. Here we describe genetic effects on gene expression levels across 44 human tissues. We find that local genetic variation affects gene expression levels for the majority of genes, and we further identify inter-chromosomal genetic effects for 93 genes and 112 loci. On the basis of the identified genetic effects, we characterize patterns of tissue specificity, compare local and distal effects, and evaluate the functional properties of the genetic effects. We also demonstrate that multi-tissue, multi-individual data can be used to identify genes and pathways affected by human disease-associated variation, enabling a mechanistic interpretation of gene regulation and the genetic basis of disease.

BACKGROUND/AIMS: Anxious responses are evolutionarily adaptive, but excessive fear can become disabling and lead to anxiety disorders. Translational models of anxiety might be useful sources for understanding the neurobiology of fear and anxiety and can contribute to future proposals of therapeutic intervention for the disorders studied. Brain-derived neurotrophic factor (BDNF), which is known for its importance on neuroplasticity and contextual memory, has emerged as a relevant element for emotional memory. Recent studies show that the Val(66)Met BDNF polymorphism correlates with various psychiatric disorders, including anxiety, but there are several differences between experimental and clinical studies. METHODS: In this work, we review the literature focused on the BDNF Val(66)Met polymorphism and anxiety, and discuss biological findings from animal models to clinical studies. RESULTS: As occurs with other psychiatric disorders, anxiety correlates with anatomical, behavioral and physiological changes related to the BDNF polymorphism. In animal studies, it has been shown that a significant decrease in regulated secretion from both BDNFVal/Met and BDNFMet/Met neurons represented a significant decrease in available BDNF. CONCLUSION: These studies suggest that developing pharmacological strategies facilitating the release of BDNF from synapses or prolongation of the half-life of secreted BDNF may improve the therapeutic responses of humans expressing the BDNF polymorphism.

MP2 and CCSD(T) complete basis set (CBS) limit interaction energies and geometries for more than 100 DNA base pairs, amino acid pairs and model complexes are for the first time presented together. Extrapolation to the CBS limit is done by using two-point extrapolation methods and different basis sets (aug-cc-pVDZ - aug-cc-pVTZ, aug-cc-pVTZ - aug-cc-pVQZ, cc-pVTZ - cc-pVQZ) are utilized. The CCSD(T) correction term, determined as a difference between CCSD(T) and MP2 interaction energies, is evaluated with smaller basis sets (6-31G** and cc-pVDZ). Two sets of complex geometries were used, optimized or experimental ones. The JSCH-2005 benchmark set, which is now available to the chemical community, can be used for testing lower-level computational methods. For the first screening the smaller training set (S22) containing 22 model complexes can be recommended. In this case larger basis sets were used for extrapolation to the CBS limit and also CCSD(T) and counterpoise-corrected MP2 optimized geometries were sometimes adopted.

Nanozymes are nanomaterial-based artificial enzymes. By effectively mimicking catalytic sites of natural enzymes or harboring multivalent elements for reactions, nanozyme systems have successfully served as direct surrogates of traditional enzymes for catalysis. With the rapid development and ever-deepening understanding of nanotechnology, nanozymes offer higher catalytic stability, ease of modification and lower manufacturing cost than protein enzymes. Additionally, nanozymes possess inherent nanomaterial properties, providing not only a simple substitute of enzymes but also a multimodal platform interfacing complex biologic environments. Recent extensive research has focused on designing various nanozyme systems that are responsive to one or multiple substrates by tailored means. Catalytic activities of nanozymes can be regulated by pH, H2O2 and glutathione concentrations and levels of oxygenation in different microenvironments. Moreover, nanozymes can be remotely-controlled via different stimuli, including a magnetic field, light, ultrasound, and heat. Collectively, these factors can be adjusted to maximize the diagnostic and therapeutic efficacies of different diseases in biomedical settings. Therefore, by integrating the catalytic property and inherent nanomaterial nature of nanozyme systems, we anticipate that stimuli-responsive nanozymes will open up new horizons for diagnosis, treatment, and theranostics.

Here we review studies that provided important information about conformational properties of DNA using circular dichroic (CD) spectroscopy. The conformational properties include the B-family of structures, A-form, Z-form, guanine quadruplexes, cytosine quadruplexes, triplexes and other less characterized structures. CD spectroscopy is extremely sensitive and relatively inexpensive. This fast and simple method can be used at low- as well as high-DNA concentrations and with short- as well as long-DNA molecules. The samples can easily be titrated with various agents to cause conformational isomerizations of DNA. The course of detected CD spectral changes makes possible to distinguish between gradual changes within a single DNA conformation and cooperative isomerizations between discrete structural states. It enables measuring kinetics of the appearance of particular conformers and determination of their thermodynamic parameters. In careful hands, CD spectroscopy is a valuable tool for mapping conformational properties of particular DNA molecules. Due to its numerous advantages, CD spectroscopy significantly participated in all basic conformational findings on DNA.

This work studies the underlying nature of H‐bonds (HBs) of different types and strengths and tries to predict binding energies (BEs) based on the properties derived from wave function analysis. A total of 42 HB complexes constructed from 28 neutral and 14 charged monomers were considered. This set was designed to sample a wide range of HB strengths to obtain a complete view about HBs. BEs were derived with the accurate coupled cluster singles and doubles with perturbative triples correction (CCSD(T))(T) method and the physical components of the BE were investigated by symmetry‐adapted perturbation theory (SAPT). Quantum theory of atoms‐in‐molecules (QTAIM) descriptors and other HB indices were calculated based on high‐quality density functional theory wave functions. We propose a new and rigorous classification of H‐bonds (HBs) based on the SAPT decomposition. Neutral complexes are either classified as “very weak” HBs with a BE ≥ −2.5 kcal/mol that are mainly dominated by both dispersion and electrostatic interactions or as “weak‐to‐medium” HBs with a BE varying between −2.5 and −14.0 kcal/mol that are only dominated by electrostatic interactions. On the other hand, charged complexes are divided into “medium” HBs with a BE in the range of −11.0 to −15.0 kcal/mol, which are mainly dominated by electrostatic interactions, or into “strong” HBs whose BE is more negative than −15.0 kcal/mol, which are mainly dominated by electrostatic together with induction interactions. Among various explored correlations between BEs and wave function‐based HB descriptors, a fairly satisfactory correlation was found for the electron density at the bond critical point (BCP; ρ BCP ) of HBs. The fitted equation for neutral complexes is BE /kcal/mol = − 223.08 × ρ BCP /a. u. + 0.7423 with a mean absolute percentage error (MAPE) of 14.7%, while that for charged complexes is BE /kcal/mol = − 332.34 × ρ BCP /a. u. − 1.0661 with a MAPE of 10.0%. In practice, these equations may be used for a quick estimation of HB BEs, for example, for intramolecular HBs or large HB networks in biomolecules. © 2019 Wiley Periodicals, Inc.

We report a reparameterization of the glycosidic torsion χ of the Cornell et al. AMBER force field for RNA, χ(OL). The parameters remove destabilization of the anti region found in the ff99 force field and thus prevent formation of spurious ladder-like structural distortions in RNA simulations. They also improve the description of the syn region and the syn-anti balance as well as enhance MD simulations of various RNA structures. Although χ(OL) can be combined with both ff99 and ff99bsc0, we recommend the latter. We do not recommend using χ(OL) for B-DNA because it does not improve upon ff99bsc0 for canonical structures. However, it might be useful in simulations of DNA molecules containing syn nucleotides. Our parametrization is based on high-level QM calculations and differs from conventional parametrization approaches in that it incorporates some previously neglected solvation-related effects (which appear to be essential for obtaining correct anti/high-anti balance). Our χ(OL) force field is compared with several previous glycosidic torsion parametrizations.

PSI4 is a free and open-source ab initio electronic structure program providing implementations of Hartree-Fock, density functional theory, many-body perturbation theory, configuration interaction, density cumulant theory, symmetry-adapted perturbation theory, and coupled-cluster theory. Most of the methods are quite efficient, thanks to density fitting and multi-core parallelism. The program is a hybrid of C++ and Python, and calculations may be run with very simple text files or using the Python API, facilitating post-processing and complex workflows; method developers also have access to most of PSI4's core functionalities via Python. Job specification may be passed using The Molecular Sciences Software Institute (MolSSI) QCSCHEMA data format, facilitating interoperability. A rewrite of our top-level computation driver, and concomitant adoption of the MolSSI QCARCHIVE INFRASTRUCTURE project, makes the latest version of PSI4 well suited to distributed computation of large numbers of independent tasks. The project has fostered the development of independent software components that may be reused in other quantum chemistry programs.

A growing body of literature suggests that there is a link between periodontitis and systemic diseases. These diseases include cardiovascular disease, gastrointestinal and colorectal cancer, diabetes and insulin resistance, and Alzheimer's disease, as well as respiratory tract infection and adverse pregnancy outcomes. The presence of periodontal pathogens and their metabolic by-products in the mouth may in fact modulate the immune response beyond the oral cavity, thus promoting the development of systemic conditions. A cause-and-effect relationship has not been established yet for most of the diseases, and the mediators of the association are still being identified. A better understanding of the systemic effects of oral microorganisms will contribute to the goal of using the oral cavity to diagnose and possibly treat non-oral systemic disease.

We report a comprehensive toolkit that enables targeted, specific modification of monocot and dicot genomes using a variety of genome engineering approaches. Our reagents, based on transcription activator-like effector nucleases (TALENs) and the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 system, are systematized for fast, modular cloning and accommodate diverse regulatory sequences to drive reagent expression. Vectors are optimized to create either single or multiple gene knockouts and large chromosomal deletions. Moreover, integration of geminivirus-based vectors enables precise gene editing through homologous recombination. Regulation of transcription is also possible. A Web-based tool streamlines vector selection and construction. One advantage of our platform is the use of the Csy-type (CRISPR system yersinia) ribonuclease 4 (Csy4) and tRNA processing enzymes to simultaneously express multiple guide RNAs (gRNAs). For example, we demonstrate targeted deletions in up to six genes by expressing 12 gRNAs from a single transcript. Csy4 and tRNA expression systems are almost twice as effective in inducing mutations as gRNAs expressed from individual RNA polymerase III promoters. Mutagenesis can be further enhanced 2.5-fold by incorporating the Trex2 exonuclease. Finally, we demonstrate that Cas9 nickases induce gene targeting at frequencies comparable to native Cas9 when they are delivered on geminivirus replicons. The reagents have been successfully validated in tomato (Solanum lycopersicum), tobacco (Nicotiana tabacum), Medicago truncatula, wheat (Triticum aestivum), and barley (Hordeum vulgare).

This review in fact represents a small monograph with 924 references. All possible aspects of the nucleic acids electrochemistry are thoroughly discussed, including early history, advances in the development of DNA hybridization sensors, analysis of biologically relevant samples, etc.

BACKGROUND: The use of homologous recombination to precisely modify plant genomes has been challenging, due to the lack of efficient methods for delivering DNA repair templates to plant cells. Even with the advent of sequence-specific nucleases, which stimulate homologous recombination at predefined genomic sites by creating targeted DNA double-strand breaks, there are only a handful of studies that report precise editing of endogenous genes in crop plants. More efficient methods are needed to modify plant genomes through homologous recombination, ideally without randomly integrating foreign DNA. RESULTS: Here, we use geminivirus replicons to create heritable modifications to the tomato genome at frequencies tenfold higher than traditional methods of DNA delivery (i.e., Agrobacterium). A strong promoter was inserted upstream of a gene controlling anthocyanin biosynthesis, resulting in overexpression and ectopic accumulation of pigments in tomato tissues. More than two-thirds of the insertions were precise, and had no unanticipated sequence modifications. Both TALENs and CRISPR/Cas9 achieved gene targeting at similar efficiencies. Further, the targeted modification was transmitted to progeny in a Mendelian fashion. Even though donor molecules were replicated in the vectors, no evidence was found of persistent extra-chromosomal replicons or off-target integration of T-DNA or replicon sequences. CONCLUSIONS: High-frequency, precise modification of the tomato genome was achieved using geminivirus replicons, suggesting that these vectors can overcome the efficiency barrier that has made gene targeting in plants challenging. This work provides a foundation for efficient genome editing of crop genomes without the random integration of foreign DNA.

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Z-DNA duplexes are a particularly complicated test case for current force fields. We performed a set of explicit solvent molecular dynamics (MD) simulations with various AMBER force field parametrizations including our recent refinements of the ε/ζ and glycosidic torsions. None of these force fields described the ZI/ZII and other backbone substates correctly, and all of them underpredicted the population of the important ZI substate. We show that this underprediction can be attributed to an inaccurate potential for the sugar-phosphate backbone torsion angle β. We suggest a refinement of this potential, β(OL1), which was derived using our recently introduced methodology that includes conformation-dependent solvation effects. The new potential significantly increases the stability of the dominant ZI backbone substate and improves the overall description of the Z-DNA backbone. It also has a positive (albeit small) impact on another important DNA form, the antiparallel guanine quadruplex (G-DNA), and improves the description of the canonical B-DNA backbone by increasing the population of BII backbone substates, providing a better agreement with experiment. We recommend using β(OL1) in combination with our previously introduced corrections, εζ(OL1) and χ(OL4), (the combination being named OL15) as a possible alternative to the current β torsion potential for more accurate modeling of nucleic acids.

The utility of molecular dynamics (MD) simulations to model biomolecular structure, dynamics, and interactions has witnessed enormous advances in recent years due to the availability of optimized MD software and access to significant computational power, including GPU multicore computing engines and other specialized hardware. This has led researchers to routinely extend conformational sampling times to the microsecond level and beyond. The extended sampling time has allowed the community not only to converge conformational ensembles through complete sampling but also to discover deficiencies and overcome problems with the force fields. Accuracy of the force fields is a key component, along with sampling, toward being able to generate accurate and stable structures of biopolymers. The Amber force field for nucleic acids has been used extensively since the 1990s, and multiple artifacts have been discovered, corrected, and reassessed by different research groups. We present a direct comparison of two of the most recent and state-of-the-art Amber force field modifications, bsc1 and OL15, that focus on accurate modeling of double-stranded DNA. After extensive MD simulations with five test cases and two different water models, we conclude that both modifications are a remarkable improvement over the previous bsc0 force field. Both force field modifications show better agreement when compared to experimental structures. To ensure convergence, the Drew-Dickerson dodecamer (DDD) system was simulated using 100 independent MD simulations, each extended to at least 10 μs, and the independent MD simulations were concatenated into a single 1 ms long trajectory for each combination of force field and water model. This is significantly beyond the time scale needed to converge the conformational ensemble of the internal portions of a DNA helix absent internal base pair opening. Considering all of the simulations discussed in the current work, the MD simulations performed to assess and validate the current force fields and water models aggregate over 14 ms of simulation time. The results suggest that both the bsc1 and OL15 force fields render average structures that deviate significantly less than 1 Å from the average experimental structures. This can be compared to similar but less exhaustive simulations with the CHARMM 36 force field that aggregate to the ∼90 μs time scale and also perform well but do not produce structures as close to the DDD NMR average structures (with root-mean-square deviations of 1.3 Å) as the newer Amber force fields. On the basis of these analyses, any future research involving double-stranded DNA simulations using the Amber force fields should employ the bsc1 or OL15 modification.

Programmed death-ligand 1 (PD-L1) expression on tumor cells is essential for T cell impairment, and PD-L1 blockade therapy has shown unprecedented durable responses in several clinical studies. Although higher expression of PD-L1 on tumor cells is associated with a better immune response after Ab blockade, some PD-L1-negative patients also respond to this therapy. In the current study, we explored whether PD-L1 on tumor or host cells was essential for anti-PD-L1-mediated therapy in 2 different murine tumor models. Using real-time imaging in whole tumor tissues, we found that anti-PD-L1 Ab accumulates in tumor tissues, regardless of the status of PD-L1 expression on tumor cells. We further observed that, while PD-L1 on tumor cells was largely dispensable for the response to checkpoint blockade, PD-L1 in host myeloid cells was essential for this response. Additionally, PD-L1 signaling in defined antigen-presenting cells (APCs) negatively regulated and inhibited T cell activation. PD-L1 blockade inside tumors was not sufficient to mediate regression, as limiting T cell trafficking reduced the efficacy of the blockade. Together, these findings demonstrate that PD-L1 expressed in APCs, rather than on tumor cells, plays an essential role in checkpoint blockade therapy, providing an insight into the mechanisms of this therapy.

Partial-nitritation anammox (PNA) is a novel wastewater treatment procedure for energy-efficient ammonium removal. Here we use genome-resolved metagenomics to build a genome-based ecological model of the microbial community in a full-scale PNA reactor. Sludge from the bioreactor examined here is used to seed reactors in wastewater treatment plants around the world; however, the role of most of its microbial community in ammonium removal remains unknown. Our analysis yielded 23 near-complete draft genomes that together represent the majority of the microbial community. We assign these genomes to distinct anaerobic and aerobic microbial communities. In the aerobic community, nitrifying organisms and heterotrophs predominate. In the anaerobic community, widespread potential for partial denitrification suggests a nitrite loop increases treatment efficiency. Of our genomes, 19 have no previously cultivated or sequenced close relatives and six belong to bacterial phyla without any cultivated members, including the most complete Omnitrophica (formerly OP3) genome to date.

Reflected-light microscopy of semitransparent material, such as unstained nervous tissue, is usually unsatisfactory because of low contrast and light scattering. In a new microscope both the object plane and the image plane were scanned in tandem so that only light reflected from the object plane was included in the image. The object was illuminated with nearly incoherent light passing through holes in one side of a rotating scanning disk (Nipkow wheel) which was imaged by the objective into the object plane. Reflected-light images of these spots were conducted to the opposite side of the same disk. Light could pass from the source to the object plane, and from the object to the image plane, only through optically congruent holes on opposite side of the rotating disk. The image obtained had better contrast and sharpness for some semitransparent material than possible in usual reflected-light microscopy.

Hydrogen-bonded nucleic acids base pairs substantially contribute to the structure and stability of nucleic acids. The study presents reference ab initio structures and interaction energies of selected base pairs with binding energies ranging from -5 to -47 kcal/mol. The molecular structures are obtained using the RI-MP2 (resolution of identity MP2) method with extended cc-pVTZ basis set of atomic orbitals. The RI-MP2 method provides results essentially identical with the standard MP2 method. The interaction energies are calculated using the Complete Basis Set (CBS) extrapolation at the RI-MP2 level. For some base pairs, Coupled-Cluster corrections with inclusion of noniterative triple contributions (CCSD(T)) are given. The calculations are compared with selected medium quality methods. The PW91 DFT functional with the 6-31G basis set matches well the RI-MP2/CBS absolute interaction energies and reproduces the relative values of base pairing energies with a maximum relative error of 2.6 kcal/mol when applied with Becke3LYP-optimized geometries. The Becke3LYP DFT functional underestimates the interaction energies by few kcal/mol with relative error of 2.2 kcal/mol. Very good performance of nonpolarizable Cornell et al. force field is confirmed and this indirectly supports the view that H-bonded base pairs are primarily stabilized by electrostatic interactions.