Novel blended SNRPE-related spliceosomopathy phenotype characterized by microcephaly and congenital atrichia

Variants in genes encoding core components of the spliceosomes are associated with craniofacial syndromes, collectively called craniofacial spliceosomopathies. SNRPE encodes a core component of pre-mRNA processing U-rich small nuclear ribonuclear proteins (UsnRNPs). Heterozygous variants in SNRPE have been reported in six families with isolated hypotrichosis simplex in addition to one case of isolated non syndromic congenital microcephaly. Here, we report a patient with a novel blended phenotype of microcephaly and congenital atrichia with multiple congenital anomalies due to a de novo intronic SNRPE deletion, c.82-28_82-16del, which results in exon skipping. As discussed within, this phenotype, which we propose be named SNRPE-related syndromic microcephaly and hypotrichosis, overlaps other craniofacial splicesosomopathies.

Accelerating Hybrid Density Functional Theory Molecular Dynamics Simulations by Seminumerical Integration, Resolution-of-the-Identity Approximation, and Graphics Processing Units

The computationally very demanding evaluation of the 4-center-2-electron (4c2e) integrals and their respective integral derivatives typically represents the major bottleneck within hybrid Kohn-Sham density functional theory molecular dynamics simulations. Building upon our previous works on seminumerical exact-exchange (sn-LinK) [Laqua, H., Thompsons, T. H., Kussmann, J., Ochsenfeld, C., J. Chem. Theory Comput. 2020, 16, 1465] and resolution-of-the-identity Coulomb (RI-J) [Kussmann, J., Laqua, H., Ochsenfeld, C., J. Chem. Theory Comput. 2021, 17, 1512], the expensive 4c2e integral evaluation can be avoided entirely, resulting in a highly efficient electronic structure theory method, allowing for fast ab initio molecular dynamics (AIMD) simulations even with large basis sets. Moreover, we propose to combine the final self-consistent field (SCF) step with the subsequent nuclear forces evaluation, providing the forces at virtually no additional cost after a converged SCF calculation, reducing the total runtime of an AIMD simulation by about another 25%. In addition, multiple independent MD trajectories can be computed concurrently on a single node, leading to a greatly increased utilization of the available hardware─especially when combined with graphics processing unit acceleration─improving the overall throughput by up to another 5 times in this way. With all of those optimizations combined, our proposed method provides nearly 3 orders of magnitude faster execution times than traditional 4c2e integral-based methods. To demonstrate the practical utility of the approach, quantum-mechanical/molecular-mechanical dynamics simulations on double-stranded DNA were performed, investigating the relative hydrogen bond strength between adenine-thymine and guanine-cytosine base pairs. In addition, this illustrative application also contains a general accuracy assessment of the introduced approximations (integration grids, resolution-of-the-identity) within AIMD simulations, serving as a protocol on how to apply these new methods to practical problems.

An effective sub-quadratic scaling atomic-orbital reformulation of the scaled opposite-spin RI-CC2 ground-state model using Cholesky-decomposed densities and an attenuated Coulomb metric

An atomic-orbital reformulation of the Laplace-transformed scaled opposite-spin (SOS) coupled cluster singles and doubles (CC2) model within the resolution of the identity (RI) approximation (SOS-RI-CC2) is presented that extends its applicability to molecules with several hundreds of atoms and triple-zeta basis sets. We exploit sparse linear algebra and an attenuated Coulomb metric to decrease the disk space demands and the computational efforts. In this way, an effective sub-quadratic computational scaling is achieved with our ω-SOS-CDD-RI-CC2 model. Moreover, Cholesky decomposition of the ground-state one-electron density matrix reduces the prefactor, allowing for an early crossover with the molecular orbital formulation. The accuracy and performance of the presented method are investigated for various molecular systems.

Factors That Determine the Variation of Equilibrium and Kinetic Properties of QM/MM Enzyme Simulations: QM Region, Conformation, and Boundary Condition

To analyze the impact of various technical details on the results of quantum mechanical (QM)/molecular mechanical (MM) enzyme simulations, including the QM region size, catechol-O-methyltransferase (COMT) is studied as a model system using an approximate QM/MM method (DFTB3/CHARMM). The results show that key equilibrium and kinetic properties for methyl transfer in COMT exhibit limited variations with respect to the size of the QM region, which ranges from ∼100 to ∼500 atoms in this study. With extensive sampling, local and global structural characteristics of the enzyme are largely conserved across the studied QM regions, while the nature of the transition state (e.g., secondary kinetic isotope effect) and reaction exergonicity are largely maintained. Deviations in the free energy profile with different QM region sizes are similar in magnitude to those observed with changes in other simulation protocols, such as different initial enzyme conformations and boundary conditions. Electronic structural properties, such as the covariance matrix of residual charge fluctuations, appear to exhibit rather long-range correlations, especially when the peptide backbone is included in the QM region; this observation holds when a range-separated DFT approach is used as the QM region, suggesting that delocalization error is unlikely the origin. Overall, the analyses suggest that multiple simulation details determine the results of QM/MM enzyme simulations with comparable contributions.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Accelerating seminumerical Fock-exchange calculations using mixed single- and double-precision arithmethic

We investigate the applicability of single-precision (fp32) floating point operations within our linear-scaling, seminumerical exchange method sn-LinK [Laqua et al., J. Chem. Theory Comput. 16, 1456 (2020)] and find that the vast majority of the three-center-one-electron (3c1e) integrals can be computed with reduced numerical precision with virtually no loss in overall accuracy. This leads to a near doubling in performance on central processing units (CPUs) compared to pure fp64 evaluation. Since the cost of evaluating the 3c1e integrals is less significant on graphic processing units (GPUs) compared to CPU, the performance gains from accelerating 3c1e integrals alone is less impressive on GPUs. Therefore, we also investigate the possibility of employing only fp32 operations to evaluate the exchange matrix within the self-consistent-field (SCF) followed by an accurate one-shot evaluation of the exchange energy using mixed fp32/fp64 precision. This still provides very accurate (1.8 µEh maximal error) results while providing a sevenfold speedup on a typical "gaming" GPU (GTX 1080Ti). We also propose the use of incremental exchange-builds to further reduce these errors. The proposed SCF scheme (i-sn-LinK) requires only one mixed-precision exchange matrix calculation, while all other exchange-matrix builds are performed with only fp32 operations. Compared to pure fp64 evaluation, this leads to 4-7× speedups for the whole SCF procedure without any significant deterioration of the results or the convergence behavior.

Combining Graphics Processing Units, Simplified Time-Dependent Density Functional Theory, and Finite-Difference Couplings to Accelerate Nonadiabatic Molecular Dynamics

Starting from our recently published implementation of nonadiabatic molecular dynamics (NAMD) on graphics processing units (GPUs), we explore further approaches to accelerate ab initio NAMD calculations at the time-dependent density functional theory (TDDFT) level of theory. We employ (1) the simplified TDDFT schemes of Grimme et al. and (2) the Hammes-Schiffer-Tully approach to obtain nonadiabatic couplings from finite-difference calculations. The resulting scheme delivers an accurate physical picture while virtually eliminating the two computationally most demanding steps of the algorithm. Combined with our GPU-based integral routines for SCF, TDDFT, and TDDFT derivative calculations, NAMD simulations of systems of a few hundreds of atoms at a reasonable time scale become accessible on a single compute node. To demonstrate this and to present a first, illustrative example, we perform TDDFT/MM-NAMD simulations of the rhodopsin protein.

Nonadiabatic Molecular Dynamics on Graphics Processing Units: Performance and Application to Rotary Molecular Motors

Nonadiabatic molecular dynamics (NAMD) simulations of molecular systems require the efficient evaluation of excited-state properties, such as energies, gradients, and nonadiabatic coupling vectors. Here, we investigate the use of graphics processing units (GPUs) in addition to central processing units (CPUs) to efficiently calculate these properties at the time-dependent density functional theory (TDDFT) level of theory. Our implementation in the FermiONs++ program package uses the J-engine and a preselective screening procedure for the calculation of Coulomb and exchange kernels, respectively. We observe good speed-ups for small and large molecular systems (comparable to those observed in ground-state calculations) and reduced (down to sublinear) scaling behavior with respect to the system size (depending on the spatial locality of the investigated excitation). As a first illustrative application, we present efficient NAMD simulations of a series of newly designed light-driven rotary molecular motors and compare their S1 lifetimes. Although all four rotors show different S1 excitation energies, their ability to rotate upon excitation is conserved, making the series an interesting starting point for rotary molecular motors with tunable excitation energies.

[Update of the S2k guidelines : Surgical treatment of benign thyroid diseases]

Thyroid resections represent one of the most common operations with 76,140 interventions in the year 2016 in Germany (source Destatis). These are predominantly benign thyroid gland diseases. Recommendations for the operative treatment of benign thyroid diseases were last published by the CAEK in 2010 as S2k guidelines (Arbeitsgemeinschaft der Wissenschaftlichen Medizinischen Fachgesellschaften e.V. [AWMF] 003/002) against the background of increasingly more radical resection procedures. Hemithyroidectomy and thyroidectomy are routinely performed for benign thyroid disease in practice. The operation-specific risks show a clear increase with the extent of the resection. Therefore, weighing-up of the risk-indications ratio between unilateral lobectomy or thyroidectomy necessitates an independent evaluation of the indications for both sides. This principle in particular has been used to update the guidelines. In addition, the previously published recommendations of the CAEK for correct execution and consequences of intraoperative neuromonitoring were included into the guidelines, which in particular serve the aim to avoid bilateral recurrent laryngeal nerve paralysis. Moreover, the recommendations for the treatment of postoperative complications, such as hypoparathyroidism and postoperative infections were revised. The updated guidelines therefore represent the current state of the science as well as the resulting surgical practice.

Stationäre und ambulante Thromboembolieprophylaxe in der Chirurgie und der perioperativen Medizin

Präambel: Alle aufgeführten medizinischen Fachgesellschaften haben sich zusammengefunden, um die früheren verschiedenen Empfehlungen zur Thromboembolieprophylaxe zu aktualisieren und in einer gemeinsamen »Leitlinie zur stationären und ambulanten Thromboembolieprophylaxe in der perioperativen Medizin« zusammenzufassen. Dazu wurden in einem ersten Schritt durch eine Konsensuskonferenz mit nominalem Gruppenprozess die früheren Empfehlungen (Fassung vom Juni 2000) überarbeitet (S2-Leitlinie). Dies erschien notwendig, um neue Therapieprinzipien und neu zugelassene Medikamente zu berücksichtigen. Als nächster Schritt erfolgt die Weiterentwicklung der vorliegenden Leitlinie nach der Vorgaben der 3. Stufe der Leitlinienentwicklung der AWMF (S3-Leitlinie).

Der präoperative Gerinnungsstatus

Der Wert einer präoperativen gerinnungsanalytischen Screening-Untersuchung zur Erhöhung der Patientensicherheit wird kontrovers diskutiert. In einer prospektiven Beobachtungsstudie mit dem Endpunkt Blutungskomplikation wurde bei 638 Patienten aus der elektiven Viszeralchirurgie die Bedeutung der präoperativen Bestimmung von Quick-Wert, APTT und Thrombozytenzahl untersucht. Zusätzlich wurde bei allen Patienten eine standardisierte Anamnese bezüglich einer möglichen Hämostasestörung erhoben. Postoperative Blutungskomplikationen wurden bei 3,1% der Patienten beobachtet. Die Gerinnungsanalytik zeigte bei 5,8% der Patienten Auffälligkeiten. Die Inzidenz von Blutungskomplikationen bei Patienten mit normaler bzw. auffälliger Gerinnungsuntersuchung unterschied sich nicht (p = 0,4141). Da sich 90% der Blutungskomplikationen bei Patienten mit normalem Gerinnungsstatus ereigneten, war die Sensitivität der Gerinnungsanalytik für das Auftreten einer Blutungskomplikation mit 10% ausgesprochen gering. Durch die standardisierte Anamnese und eine sich daran anschließende Spezialdiagnostik wurden 3 Patienten mit einem v.-Willebrand-Typ I identifiziert. Bei diesen Patienten war der Gerinnungsstatus normal.

Efficient and Accurate Born-Oppenheimer Molecular Dynamics for Large Molecular Systems

An efficient scheme for the calculation of Born-Oppenheimer molecular dynamics (BOMD) simulations is introduced. It combines the corrected small basis set Hartree-Fock (HF-3c) method by Sure and Grimme [J. Comput. Chem. 2013, 43, 1672], extended Lagrangian BOMD (XL-BOMD) by Niklasson et al. [J. Chem. Phys. 2009, 130, 214109], and the calculation of the two electron integrals on graphics processing units (GPUs) [J. Chem. Phys. 2013, 138, 134114; J. Chem. Theory Comput. 2015, 11, 918]. To explore the parallel performance of our strong scaling implementation of the method, we present timings and extract, as its validation and first illustrative application, high-quality vibrational spectra from simulated trajectories of β-carotene, paclitaxel, and liquid water (up to 500 atoms). We conclude that the presented BOMD scheme may be used as a cost-efficient and reliable tool for computing vibrational spectra and thermodynamics of large molecular systems including explicit solvent molecules containing 500 atoms and more. Simulating 50 ps of maitotoxin (nearly 500 atoms) employing time steps of 0.5 fs requires ∼3 weeks on 12 CPUs (Intel Xeon E5 2620 v3) with 24 GPUs (AMD FirePro 3D W8100).

Monte Carlo simulations of single- and multistep enzyme-catalyzed reaction sequences: effects of diffusion, cell size, enzyme fluctuations, colocalization, and segregation

Following the discovery of slow fluctuations in the catalytic activity of an enzyme in single-molecule experiments, it has been shown that the classical Michaelis-Menten (MM) equation relating the average enzymatic velocity and the substrate concentration may hold even for slowly fluctuating enzymes. In many cases, the average velocity is that given by the MM equation with time-averaged values of the fluctuating rate constants and the effect of enzyme fluctuations is simply averaged out. The situation is quite different for a sequence of reactions. For colocalization of a pair of enzymes in a sequence to be effective in promoting reaction, the second must be active when the first is active or soon after. If the enzymes are slowly varying and only rarely active, the product of the first reaction may diffuse away before the second enzyme is active, and colocalization may have little value. Even for single-step reactions the interplay of reaction and diffusion with enzyme fluctuations leads to added complexities, but for multistep reactions the interplay of reaction and diffusion, cell size, compartmentalization, enzyme fluctuations, colocalization, and segregation is far more complex than for single-step reactions. In this paper, we report the use of stochastic simulations at the level of whole cells to explore, understand, and predict the behavior of single- and multistep enzyme-catalyzed reaction systems exhibiting some of these complexities. Results for single-step reactions confirm several earlier observations by others. The MM relationship, with altered constants, is found to hold for single-step reactions slowed by diffusion. For single-step reactions, the distribution of enzymes in a regular grid is slightly more effective than a random distribution. Fluctuations of enzyme activity, with average activity fixed, have no observed effects for simple single-step reactions slowed by diffusion. Two-step sequential reactions are seen to be slowed by segregation of the enzymes for each step, and results of the calculations suggest limits for cell size. Colocalization of enzymes for a two-step sequence is seen to promote reaction, and rates fall rapidly with increasing distance between enzymes. Low frequency fluctuations of the activities of colocalized enzymes, with average activities fixed, can greatly reduce reaction rates for sequential reactions.

Linear-scaling method for calculating nuclear magnetic resonance chemical shifts using gauge-including atomic orbitals within Hartree-Fock and density-functional theory

Details of a new density matrix-based formulation for calculating nuclear magnetic resonance chemical shifts at both Hartree-Fock and density functional theory levels are presented. For systems with a nonvanishing highest occupied molecular orbital-lowest unoccupied molecular orbital gap, the method allows us to reduce the asymptotic scaling order of the computational effort from cubic to linear, so that molecular systems with 1000 and more atoms can be tackled with today's computers. The key feature is a reformulation of the coupled-perturbed self-consistent field (CPSCF) theory in terms of the one-particle density matrix (D-CPSCF), which avoids entirely the use of canonical MOs. By means of a direct solution for the required perturbed density matrices and the adaptation of linear-scaling integral contraction schemes, the overall scaling of the computational effort is reduced to linear. A particular focus of our formulation is to ensure numerical stability when sparse-algebra routines are used to obtain an overall linear-scaling behavior.

Helical Packing of Discotic Hexaphenyl Hexa-peri-hexabenzocoronenes: Theory and Experiment

The arrangement of discotic hexa-peri-hexabenzocoronenes (HBCs) into columnar helical superstructures has been investigated in relation to their molecular architecture. The supramolecular structure of two hexaphenyl-substituted HBC derivatives, differing only in the chiral/achiral nature of the attached alkyl side chains, was studied by circular dichroism and temperature-dependent wide-angle X-ray diffraction on oriented filaments. A structural model in agreement with the experimental observations was developed on the basis of accompanying quantum-chemical calculations. The helical organization along the self-assembled columnar structures was induced by the steric requirements of the bulky phenyl rings near the aromatic core, i.e., by their rotation out-of-plane with respect to the aromatic core. On the other hand, a uniform handedness of the twist was generated by chiral alkyl substituents. At higher temperatures the degree of helical organization decreases due to lateral and longitudinal dynamics of the discotic molecules. Annealing at ambient conditions improved the long-range arrangement of the discs along the columnar structures. This reorganization indicated a self-healing of the plastic material which is desirable for application of discotics as active layers in electronic devices. The helical packing resulted in a considerable stability of the mesophase up to 500 degrees C, which has not been reported for a discotic so far.

Structure of Molecular Tweezer Complexes in the Solid State: NMR Experiments, X-ray Investigations, and Quantum Chemical Calculations

The structure of supramolecular complexes formed by a naphthalene-spaced tweezer molecule as host and 1,4-dicyanobenzene (DCNB), 1,2,4,5-tetracyanobenzene (TCNB), and 7,7,8,8-tetracyano-p-quinodimethane (TCNQ) as aromatic, electron-deficient guests is investigated by solid-state NMR and X-ray diffraction measurements. Quantum chemical calculations using linear scaling methods are applied to predict and to assign the 1H NMR chemical shifts of the complexes. By combining experiment and theory, insights into intra- and intermolecular effects influencing the proton chemical shifts of the host-guest system are provided in the solid state.

Molecular recognition in molecular tweezers systems: quantum-chemical calculation of NMR chemical shifts

Quantum-chemical calculations for molecular tweezers systems are presented, where the focus is not only on the recognition process in the host-guest systems, but on the self aggregation of the tweezers host as well. Such intermolecular interactions influence the corresponding NMR spectra strongly by up to 6 ppm for proton chemical shifts, since ring-current effects are particularly important. The quantum-chemical results allow one to reliably assign the spectra and to gain information both on the structure and on the importance of intra- and intermolecular interactions. In addition, we study the accuracy of a variety of density functionals for describing the present host-guest systems, where we observe a considerable underestimation of ring-current effects on (1)H NMR chemical shifts at the density functional theory (DFT) level using smaller basis sets such as 6-31G**, so that larger bases like TZP are required. This stands in contrast to the behavior of the Hartree-Fock scheme, where small basis sets, such as 6-31G**, provide reliable (1)H NMR shieldings for molecular tweezers systems.

Advances in methods and algorithms in a modern quantum chemistry program package

Advances in theory and algorithms for electronic structure calculations must be incorporated into program packages to enable them to become routinely used by the broader chemical community. This work reviews advances made over the past five years or so that constitute the major improvements contained in a new release of the Q-Chem quantum chemistry package, together with illustrative timings and applications. Specific developments discussed include fast methods for density functional theory calculations, linear scaling evaluation of energies, NMR chemical shifts and electric properties, fast auxiliary basis function methods for correlated energies and gradients, equation-of-motion coupled cluster methods for ground and excited states, geminal wavefunctions, embedding methods and techniques for exploring potential energy surfaces.