Graphene BioFET sensors for SARS-CoV-2 detection: a multiscale simulation approach.

Biological Field-Effect Transistors (BioFETs) have already demonstrated enormous potential for detecting minute amounts of ions and molecules. The use of two-dimensional (2D) materials has been shown to boost their performance and to enable the design of new applications. This combination deserves special interest in the current pandemic caused by the SARS-CoV-2 virus which demands fast, reliable and cheap detection methods. However, in spite of the experimental advances, there is a lack of a comprehensive and in-depth computational approach to capture the mechanisms underlying the sensor behaviour. Here, we present a multiscale platform that combines detailed atomic models of the molecules with mesoscopic device-level simulations. The fine-level description exploited in this approach accounts for the charge distribution of the receptor, its reconfiguration when the target binds to it, and the consequences in terms of sensitivity on the transduction mechanism. The results encourage the further exploration of improved sensor designs and 2D materials combined with diverse receptors selected to achieve the desired specificity.

Protein conformational flexibility modulates kinetics and thermodynamics of drug binding.

Structure-based drug design has often been restricted by the rather static picture of protein-ligand complexes presented by crystal structures, despite the widely accepted importance of protein flexibility in biomolecular recognition. Here we report a detailed experimental and computational study of the drug target, human heat shock protein 90, to explore the contribution of protein dynamics to the binding thermodynamics and kinetics of drug-like compounds. We observe that their binding properties depend on whether the protein has a loop or a helical conformation in the binding site of the ligand-bound state. Compounds bound to the helical conformation display slow association and dissociation rates, high-affinity and high cellular efficacy, and predominantly entropically driven binding. An important entropic contribution comes from the greater flexibility of the helical relative to the loop conformation in the ligand-bound state. This unusual mechanism suggests increasing target flexibility in the bound state by ligand design as a new strategy for drug discovery.

Three steps to gold: mechanism of protein adsorption revealed by Brownian and molecular dynamics simulations.

The addition of three N-terminal histidines to ß-lactamase inhibitor protein was shown experimentally to increase its binding potency to an Au(111) surface substantially but the binding mechanism was not resolved. Here, we propose a complete adsorption mechanism for this fusion protein by means of a multi-scale simulation approach and free energy calculations. We find that adsorption is a three-step process: (i) recognition of the surface predominantly by the histidine fusion peptide and formation of an encounter complex facilitated by a reduced dielectric screening of water in the interfacial region, (ii) adsorption of the protein on the surface and adoption of a specific binding orientation, and (iii) adaptation of the protein structure on the metal surface accompanied by induced fit. We anticipate that the mechanistic features of protein adsorption to an Au(111) surface revealed here can be extended to other inorganic surfaces and proteins and will therefore aid the design of specific protein-surface interactions.

The vibrationally mediated photodissociation of Cl2.

The photodissociation of vibrationally excited Cl(2)(v = 1) has been investigated experimentally using the velocity mapped ion imaging technique. The experimental measurements presented here are compared with the results of time-dependent wavepacket calculations performed on a set of ab initio potential energy curves. The high level calculations allow prediction of all the dynamical information regarding the dissociation, including electronic polarization effects. Using a combination of theory and experiment it was found that there was negligible cooling of the vibrational degree of freedom of the parent molecule in the molecular beam. The results presented are compared with those following the photodissociation of Cl(2)(v = 0). Although the same electronic states are found to be important for Cl(2)(v = 1) as for Cl(2)(v = 0), significant differences were found regarding many of the observables. The overall level of agreement between theory and experiment was found to be reasonable and confirms previous assignments of the photodissociation mechanism.

A complete quantum mechanical study of chlorine photodissociation.

A fully quantum mechanical dynamical calculation on the photodissociation of molecular chlorine is presented. The magnitudes and phases of all the relevant photofragment T-matrices have been calculated, making this study the computational equivalent of a "complete experiment," where all the possible parameters defining an experiment have been determined. The results are used to simulate cross-sections and angular momentum polarization information which may be compared with experimental data. The calculations rigorously confirm the currently accepted mechanism for the UV photodissociation of Cl(2), in which the majority of the products exit on the C(1)Π(1u) state, with non-adiabatic couplings to the A(3)Π(1u) and several other Ω = 1 states, and a small contribution from the B(3)Π state present at longer wavelengths.

Electronic polarization effects in the photodissociation of Cl2.

Velocity mapped ion imaging and resonantly enhanced multiphoton ionization time-of-flight methods have been used to investigate the photodissociation dynamics of the diatomic molecule Cl(2) following excitation to the first UV absorption band. The experimental results presented here are compared with high level time dependent wavepacket calculations performed on a set of ab initio potential energy curves [D. B. Kokh, A. B. Alekseyev, and R. J. Buenker, J. Chem. Phys. 120, 11549 (2004)]. The theoretical calculations provide the first determination of all dynamical information regarding the dissociation of a system of this complexity, including angular momentum polarization. Both low rank K = 1, 2 and high rank K = 3 electronic polarization are predicted to be important for dissociation into both asymptotic product channels and, in general, good agreement is found between the recent theory and the measurements made here, which include the first experimental determination of high rank K = 3 orientation.

Use of exchange maximization to generate starting vectors for self-consistent field calculations on metal cluster/adsorbate systems.

Localized molecular orbitals (LMOs) derived from exchange maximization with respect to all atom-centered basis functions in the basis set are shown to generate a good starting electronic field for self-consistent field calculations on extended systems such as metal clusters, for which well-defined chemical bonds are not present. Examples studied are a cluster of 20 Ni atoms and the Pt(97)CO, Ag(43)/H(3)CNON, Ag(91)/H(2)CO, and vinylidene/Ni metal cluster plus adsorbate systems. It is also shown that improved starting vectors can be obtained by remixing a subset of the LMOs with the largest exchange eigenvalues through diagonalization of the Fock matrix computed with a null electronic field. Employing only a subset of the exchange-maximized LMOs in the first iterations, and then gradually expanding the space in which the diagonalizations are carried out in succeeding cycles, is shown to be an effective means of guiding the SCF procedure to the converged full-basis solution.