Atomic Detail Research Articles

Computation of conformational transitions and activation of Src tyrosine kinase.

A protein molecule can adapt different structural forms depending on solution conditions or ligation state. The stable forms can be defined experimentally, but the transition pathways and metastable states between them cannot. A computational method to explore equilibrium conformational transition pathways is the adaptive bias path optimization (ABPO) algorithm. ABPO, implemented in the CHARMM program, differs from other path methods by using enhanced sampling along a path between two known states without restraining the system to the reduced coordinates of the evolving path. ABPO converges to a physically reasonable path even in the case of a highly rugged energy landscape and is amenable to assessing the power of the reduced variables chosen to explore the transition. ABPO was applied to three cases where a conformational transition is localized to one part of the polypeptide chain. With proper selection of the reduced coordinate set, sampling without restraints converged to a physically reasonable path. We also applied ABPO to a more complex transition of Src tyrosine kinase catalytic domain that undergoes a more complex transition between enzymatically down‑regulated and activated conformations. The results using coarse-grained and all‑atom models identified a switched electrostatic network (SEN) that could confer an entropic advantage to reduce energetic barriers by limiting the search space of the activation loop. The richness in atomic detail provided by the ABPO path reveals how the SEN couples the dynamics of the key αC‑helix and the activation loop, which provided a new rationale for the highly conserved HRD sequence of kinases. Further, the results related to energetics of the rotation of the αC‑helix provide a rationale for the interactions observed in a number of kinase regulatory complexes, some for which the basis for regulation from inspection of the structure along is puzzling.

Structural basis underlying the synergism of NADase and SLO during group A Streptococcus infection

Group A Streptococcus (GAS) is a strict human pathogen possessing a unique pathogenic trait that utilizes the cooperative activity of NAD+-glycohydrolase (NADase) and Streptolysin O (SLO) to enhance its virulence. How NADase interacts with SLO to synergistically promote GAS cytotoxicity and intracellular survival is a long-standing question. Here, the structure and dynamic nature of the NADase/SLO complex are elucidated by X-ray crystallography and small-angle scattering, illustrating atomic details of the complex interface and functionally relevant conformations. Structure-guided studies reveal a salt-bridge interaction between NADase and SLO is important to cytotoxicity and resistance to phagocytic killing during GAS infection. Furthermore, the biological significance of the NADase/SLO complex in GAS virulence is demonstrated in a murine infection model. Overall, this work delivers the structure-functional relationship of the NADase/SLO complex and pinpoints the key interacting residues that are central to the coordinated actions of NADase and SLO in the pathogenesis of GAS infection.

Computational Approaches to Enzyme Inhibition by Marine Natural Products in the Search for New Drugs.

The exploration of biologically relevant chemical space for the discovery of small bioactive molecules present in marine organisms has led not only to important advances in certain therapeutic areas, but also to a better understanding of many life processes. The still largely untapped reservoir of countless metabolites that play biological roles in marine invertebrates and microorganisms opens new avenues and poses new challenges for research. Computational technologies provide the means to (i) organize chemical and biological information in easily searchable and hyperlinked databases and knowledgebases; (ii) carry out cheminformatic analyses on natural products; (iii) mine microbial genomes for known and cryptic biosynthetic pathways; (iv) explore global networks that connect active compounds to their targets (often including enzymes); (v) solve structures of ligands, targets, and their respective complexes using X-ray crystallography and NMR techniques, thus enabling virtual screening and structure-based drug design; and (vi) build molecular models to simulate ligand binding and understand mechanisms of action in atomic detail. Marine natural products are viewed today not only as potential drugs, but also as an invaluable source of chemical inspiration for the development of novel chemotypes to be used in chemical biology and medicinal chemistry research.

The Role of Transient Intermediate Structures in the Unfolding of the Trp-Cage Fast-Folding Protein: Generating Ensembles from Time-Resolved X-ray Solution Scattering with Genetic Algorithms.

The Trp-cage miniprotein is one of the smallest systems to exhibit a stable secondary structure and fast-folding dynamics, serving as an apt model system to study transient intermediates with both experimental and computational analyses. Previous spectroscopic characterizations that have been done on Trp-cage have inferred a single stable intermediate on a pathway from folded to unfolded basins. We aim to bridge the understanding of Trp-cage structural folding dynamics on microsecond-time scales, by utilizing time-resolved X-ray solution scattering to probe the temperature-induced unfolding pathway. Our results indicate the formation of a conformationally extended intermediate on the time scale of 1 μs, which undergoes complete unfolding within 5 μs. We further investigated the atomistic structural details of the unfolding pathway using a genetic algorithm to generate ensemble model fits to the scattering profiles. This analysis paves the way for direct benchmarking of theoretical models of protein folding ensembles produced with molecular dynamics simulations.

Understanding Self-Assembly of Silica-Precipitating Peptides to Control Silica Particle Morphology.

The most advanced materials are those found in nature. These evolutionary optimized substances provide highest efficiencies, e.g., in harvesting solar energy or providing extreme stability, and are intrinsically biocompatible. However, the mimicry of biological materials is limited to a few successful applications since there is still a lack of the tools to recreate natural materials. Herein, such means are provided based on a peptide library derived from the silaffin protein R5 that enables rational biomimetic materials design. It is now evident that biomaterials do not form via mechanisms observed in vitro. Instead, the material's function and morphology are predetermined by precursors that self-assemble in solution, often from a combination of protein and salts. These assemblies act as templates for biomaterials. The RRIL peptides used here are a small part of the silica-precipitation machinery in diatoms. By connecting RRIL motifs via varying central bi- or trifunctional residues, a library of stereoisomers is generated, which allows characterization of different template structures in the presence of phosphate ions by combining residue-resolved real-time NMR spectroscopy and molecular dynamics (MD) simulations. Understanding these templates in atomistic detail, the morphology of silica particles is controlled via manipulation of the template precursors.

Towards a structurally resolved human protein interaction network

Cellular functions are governed by molecular machines that assemble through protein-protein interactions. Their atomic details are critical to studying their molecular mechanisms. However, fewer than 5% of hundreds of thousands of human protein interactions have been structurally characterized. Here we test the potential and limitations of recent progress in deep-learning methods using AlphaFold2 to predict structures for 65,484 human protein interactions. We show that experiments can orthogonally confirm higher-confidence models. We identify 3,137 high-confidence models, of which 1,371 have no homology to a known structure. We identify interface residues harboring disease mutations, suggesting potential mechanisms for pathogenic variants. Groups of interface phosphorylation sites show patterns of co-regulation across conditions, suggestive of coordinated tuning of multiple protein interactions as signaling responses. Finally, we provide examples of how the predicted binary complexes can be used to build larger assemblies helping to expand our understanding of human cell biology.

Polarity Control by Inversion Domain Suppression in N-Polar III-Nitride Heterostructures

Nitrogen-polar III-nitride heterostructures offer advantages over metal-polar structures in high frequency and high power applications. However, polarity control in III-nitrides is difficult to achieve as a result of unintentional polarity inversion domains (IDs). Herein, we present a comprehensive structural investigation with both atomic detail and thermodynamic analysis of the polarity evolution in low- and high-temperature AlN layers on on-axis and 4° off-axis carbon-face 4H-SiC (0001̅) grown by hot-wall metal organic chemical vapor deposition. A polarity control strategy has been developed by variation of thermodynamic Al supersaturation and substrate misorientation angle in order to achieve the desired growth mode and polarity. We demonstrate that IDs are completely suppressed for high-temperature AlN nucleation layers when a step-flow growth mode is achieved on the off-axis substrates. We employ this approach to demonstrate high quality N-polar epitaxial AlGaN/GaN/AlN heterostructures.

Structural basis of bacterial effector protein azurin targeting tumor suppressor p53 and inhibiting its ubiquitination

Tumor suppressor p53 prevents tumorigenesis by promoting cell cycle arrest and apoptosis through transcriptional regulation. Dysfunction of p53 occurs frequently in human cancers. Thus, p53 becomes one of the most promising targets for anticancer treatment. A bacterial effector protein azurin triggers tumor suppression by stabilizing p53 and elevating its basal level. However, the structural and mechanistic basis of azurin-mediated tumor suppression remains elusive. Here we report the atomic details of azurin-mediated p53 stabilization by combining X-ray crystallography with nuclear magnetic resonance. Structural and mutagenic analysis reveals that the p28 region of azurin, which corresponds to a therapeutic peptide, significantly contributes to p53 binding. This binding stabilizes p53 by disrupting COP1-mediated p53 ubiquitination and degradation. Using the structure-based design, we obtain several affinity-enhancing mutants that enable amplifying the effect of azurin-induced apoptosis. Our findings highlight how the structure of the azurin-p53 complex can be leveraged to design azurin derivatives for cancer therapy.

Point-of-Care Portable 3D-Printed Multispectral Sensor for Real-Time Enzyme Activity Monitoring in Healthcare Applications

Absorbance spectroscopy finds many biomedical and physical applications ranging from studying the atomic and molecular details of the analyte to determination of unknown biological species and their concentration or activity in the samples. Commercially available laboratory-based spectrometers are usually bulky and require high power and laborious manual processing, making them unsuitable to be deployed in portable and space-constrained environments, thereby further limiting their utility for real-time on-site monitoring. To address these challenges, here we developed a portable 3D-printed multispectral spectrophotometer based on absorbance spectroscopy for real-time monitoring of enzyme molecular activity. Monitoring enzyme (such as tyrosinase) activity is critical, as it quantifies its reaction rate, which is dependent on many factors such as the enzyme and substrate concentrations, temperature, pH, and other regulators such as inhibitors and effectors. Tyrosinase is a critical enzyme responsible for melanin synthesis in living beings and exhibits enzymatic browning in fruits and vegetables. It finds various commercial applications in the fields of healthcare (skin pigmentation, wound healing, etc.), forensics, and food processing. Here, tyrosinase activity was monitored using a 3D-printed spectral sensor at different rates and compared against measurements obtained from laboratory instruments. The enzyme activity was also studied using kojic acid (i.e., a commonly employed commercial tyrosinase inhibitor) while varying its molar and volume concentrations to control the reaction rate at discrete activity levels. For tyrosinase activity monitoring, the fabricated device has shown significant correlation (R2 = 0.9999) compared to measurements from the standard table-top spectrophotometer. We also provide a performance comparison between the 3D-printed and the laboratory spectrophotometer instruments by studying tyrosinase enzyme activity with and without the influence of an inhibitor. Such a device can be translated into various absorbance spectroscopy-based point-of-care biomedical and healthcare applications.

Three-Electron Two-Centered Bond and Single-Electron Transfer Mechanism of Water Splitting via a Copper–Bipyridine Complex

We report the atomistic and electronic details of the mechanistic pathway of the oxygen-oxygen bond formation catalyzed by a copper-2,2'-bipyridine complex. Density functional theory-based molecular dynamics simulations and enhanced sampling methods were employed for this study. The thermodynamics and electronic structure of the oxygen-oxygen bond formation are presented in this study by considering the cis-bishydroxo, [CuIII(bpy)(OH)2]+, and cis-(hydroxo)oxo, [CuIV(bpy)(OH)(═O)]+, complexes as active catalysts. In the cis-bishydroxo complex, the hydroxide transfer requires a higher kinetic barrier than the proton transfer process. In the case of [CuIV(bpy)(OH)(═O)]+, the proton transfer requires a higher free energy than the hydroxide one. The peroxide bond formation is thermodynamically favorable for the [CuIV(bpy)(OH)(═O)]+ complex compared with the other. The hydroxide ion is transferred to one of the Cu-OH moieties, and the proton is transferred to the solvent. The free energy barrier for this migration is higher than that for the former transfer. From the analysis of molecular orbitals, it is found that the electron density is primarily present on the water molecules near the active sites in the highest occupied molecular orbital (HOMO) state and lowest unoccupied molecular orbital (LUMO) of the ligands. Natural bond orbital (NBO) analysis reveals the electron transfer process during the oxygen-oxygen bond formation. The σ*Cu(dxz)-O(p) orbitals are involved in the oxygen-oxygen bond formation. During the bond formation, three-electron two-centered (3e--2C) bonds are observed in [CuIII(bpy)(OH)2]+ during the transfer of the hydroxide before the formation of the oxygen-oxygen bond.

Recent advances on molecular dynamics-based techniques to address drug membrane permeability with atomistic detail

Several factors affect the passive membrane permeation of small molecules, including size, charge, pH, or the presence of specific chemical groups. Understanding these features is paramount to identifying or designing drug candidates with optimal ADMET properties and this can be achieved through experimental/knowledge-based methodologies or using computational approaches. Empirical methods often lack detailed information about the underlying molecular mechanism. In contrast, Molecular Dynamics-based approaches are a powerful strategy, providing an atomistic description of this process. This technique is continuously growing, featuring new related methodologies. In this work, the recent advances in this research area will be discussed.

Crystallographic fragment screening in academic cancer drug discovery.

Fragment-based drug discovery (FBDD) has brought several drugs to the clinic, notably to target proteins once considered to be challenging, or even undruggable. Screening in FBDD relies upon observing and/or measuring weak (millimolar-scale) binding events using biophysical techniques or crystallographic fragment screening. This latter structural approach provides no information about binding affinity but can reveal binding mode and atomic detail on protein-fragment interactions to accelerate hit-to-lead development. In recent years, high-throughput platforms have been developed at synchrotron facilities to screen thousands of fragment-soaked crystals. However, using accessible manual techniques it is possible to run informative, smaller-scale screens within an academic lab setting. This chapter describes general protocols for home laboratory-scale fragment screening, from fragment soaking through to structure solution and, where appropriate, signposts to background, protocols or alternatives elsewhere.

Molecular Dynamics Simulation of the Aggregation of Extracellular Polymeric Substance

Extracellular polymeric substance (EPS) play an important role in the transport and transformation of pollutants. At present, some scholars have conducted basic experimental characterization of EPS, but the molecular lev-el interactions among the various components of EPS still need to be fully characterized. Here, we use molecular dynamics (MD) simulation to explore the structural properties of EPS systems in atomic detail. The process and mechanism of the aggregation of different EPS components (polysaccharide, lipid, nucleic acid and protein) were revealed by the simulation results. EPS aggregation consist of a hydrophobic core and an amphiphilic exterior. Lipid tail, as a hydrophobic core, promoted the aggregation of EPS. But strong hy-drophilic nucleic acid and protein components inhibited the aggregation of EPS in water and were located outside the aggregation body. The details of the structure of EPS aggregation are revealed here, which provides a micro molecular mechanism for understanding the interaction between EPS.

Coverage-modulated halogen bond geometry transformation in supramolecular assemblies.

Halogen bonding (HB) has emerged as a promising route for designing supramolecular assemblies due to its directional nature and versatility in modifying interactions through the choice of halogens and molecular entities. Despite this, methods for tuning these interactions on surfaces, particularly in terms of directionality, are limited. In this study, we present a strategy for tuning the directionality of self-assembly processes in homomolecular organic compounds on inert metal surfaces. A variety of halogen-halogen geometries can promote highly-extended one-dimensional or two-dimensional self-assembly depending on the molecular coverage. Our results indicate that under lower molecular coverage conditions, robust one-dimensional (1D) structures promote the self-assembly of halogen-bonded molecules on Au(111). At certain coverage, a transformation from type-I to synthon halogen bonding is observed, leading to an extended hexagonal pattern of molecular assembly. The atomistic details of the structures are experimentally studied using high-resolution atomic force microscopy and supported by first-principle calculations. We employed DFT to evaluate the interplay between electrostatics and dispersion forces driving both type-I and synthon assemblies. The results reveal a halogen-bond geometry transformation induced by a subtle balance of molecule-molecule interaction. Finally, we investigate the capability of the halogen-bonded supramolecular assembly to periodically confine electronic quantum states and single atoms. Our findings demonstrate the versatility of sigma-bonding in regulating molecular assembly and provide new insights for tailoring functional molecular structures on an inert metal substrate.

How a single mutation alters the protein structure: a simulation investigation on protein tyrosine phosphatase SHP2.

Protein tyrosine phosphatase SHP2 is a key regulator modulating several signaling pathways. The oncogenic mutation E76K in SHP2 releases the enzyme from an autoinhibited, closed conformation into an active, open conformation. Here, we investigated the conformational dynamics of SHP2 and the effect of the E76K mutation on its conformational ensemble via extensive molecular dynamics (MD) and metadynamics (MetaD) simulations. Our simulations provide atomistic details on how the E76K mutated SHP2 prefers the open state and also reveal that the transition between the closed and the open states is highly collective. Several intermediate metastable states during the conformational transition between the closed and the open states were also investigated. Understanding how the single E76K mutation induces the conformational change in SHP2 could facilitate the further design of SHP2 inhibitors.

Analysis of Interaction Between Odorant Receptors and Flexible Spike of SARS CoV-2- Key to Loss of Smell.

The development of a vaccine for SARS-CoV-2 is primarily focused on the structure of the spike (S) protein. The heavy glycosylation of S with flexible hinges at the stalk shields from antibody attachment. This study deciphers the flexible nature of hinges responsible for binding the odorant receptor on neurons responsible for the loss of smell in COVID-19 patients. The 3D structure via EPIK in Maestro, protein docking with ligands via Maestro protein analysis tool, and molecular dynamic simulation at 30 ns run using DESMOND was prepared. The data of the study strongly suggest that strong and stable bond formation results from the reaction between R:14: Trp and Phe at the residue, targeting the flexible hinges of SARS-CoV-2. The difference in the conformational structure of the S protein and its binding with the odorant receptor in COVID-19 is the prime factor for the loss of smell and taste in patients, as supported by the concept of Antigen (epitope) Antibody interaction by the stable formation of a hydrogen bond among odorant receptor and the S protein. The flexibility of structural proteins determines the binding potential of antibodies or other defense proteins produced to participate in the antigen-antibody reaction. Molecular and atomic details potentiate the design and screening of small molecules that can inhibit the fusion at entry level or odorant receptors and potentially be used in the prevention and treatment of infection, particularly when formulated as nasal drops, paving a new approach for pharmacologists in the treatment of COVID-19 infection.

The combustion-supporting mechanism of fluoropolymers on aluminum particles studied using reactive dynamics simulations

In order to clarify the combustion-supporting mechanism of fluoropolymers, this work focuses on the atomic details of the oxidation of Al nanoparticles with the fluoropolymers and O2, using reactive dynamics simulations.

Hydrolysis Mechanism of Carbamate Methomyl by a Novel Esterase PestE: A QM/MM Approach.

Methomyl is one of the most important carbamates that has caused potential hazardous effects on both human beings and the environment. Here, we systematically investigated the hydrolysis mechanism of methomyl catalyzed by esterase PestE using molecular dynamics simulations (MD) and quantum mechanics/molecular mechanics (QM/MM) calculations. The hydrolysis mechanism involves two elementary steps: (Ⅰ) serine-initiated nucleophilic attack and (Ⅱ) C-O bond cleavage. Our work elicits the atomic level details of the hydrolysis mechanism and free energy profiles along the reaction pathway. The Boltzmann-weighted average potential barriers are 19.1 kcal/mol and 7.5 kcal/mol for steps Ⅰ and Ⅱ, respectively. We identified serine-initiated nucleophilic attack as the rate determining-step. The deep learning-based kcat prediction model indicated that the barrier of the rate-determining step is 15.4 kcal/mol, which is in good agreement with the calculated results using Boltzmann-weighted average method. We have elucidated the importance of the protein-substrate interactions and the roles of the key active site residues during the hydrolysis process through noncovalent interactions analysis and electrostatic potential (ESP) analysis. The results provide practical value for achieving efficient degradation of carbamates by hydrolases.

NARall: a novel tool for reconstruction of the all-atom structure of nucleic acids from heavily coarse-grained model

Nucleic acids are one of the most important cellular components. These molecules have been studied both experimentally and theoretically. As all-atom simulations are still limited to short time scales, coarse-grain modeling allows to study of those molecules on a longer time scale. Nucleic-Acid united RESidue (NARES-2P) is a low-resolution coarse-grained model with two centers of interaction per repeating unit. It has been successfully applied to study DNA self-assembly and telomeric properties. This force field enables the study of nucleic acids Behavior on a long time scale but lacks atomistic details. In this article, we present new software to reconstruct atomistic details from the NARES-2P model. It has been applied to RNA pseudoknot, nucleic acid four-way junction, G-quadruplex and DNA duplex converted to NARES-2P model and back. Moreover, it was applied to DNA structure folded and self-assembled with NARES-2P.

Early Steps in C-Type Inactivation of the hERG Potassium Channel.

Fast C-type inactivation confers distinctive functional properties to the hERG potassium channel, and its association to inherited and acquired cardiac arrythmias makes the study of the inactivation mechanism of hERG at the atomic detail of paramount importance. At present, two models have been proposed to describe C-type inactivation in K+-channels. Experimental data and computational work on the bacterial KcsA channel support the hypothesis that C-type inactivation results from a closure of the selectivity filter that sterically impedes ion conduction. Alternatively, recent experimental structures of a mutated Shaker channel revealed a widening of the extracellular portion of the selectivity filter, which might diminish conductance by interfering with the mechanism of ion permeation. Here, we performed molecular dynamics simulations of the wild-type hERG, a non-inactivating mutant (hERG-N629D), and a mutant that inactivates faster than the wild-type channel (hERG-F627Y) to find out which and if any of the two reported C-type inactivation mechanisms applies to hERG. Closure events of the selectivity filter were not observed in any of the simulated trajectories but instead, the extracellular section of the selectivity filter deviated from the canonical conductive structure of potassium channels. The degree of widening of the potassium binding sites at the extracellular entrance of the channel was directly related to the degree of inactivation with hERG-F627Y > wild-type hERG > hERG-N629D. These findings support the hypothesis that C-type inactivation in hERG entails a widening of the extracellular entrance of the channel rather than a closure of the selectivity filter.