Discovery through the computational microscope

Illuminating ligand-induced dynamics in nuclear receptors through MD simulations

Nuclear receptors (NRs) regulate gene expression in critical physiological processes, with their functionality finely tuned by ligand-induced conformational changes. While NRs may sometimes undergo significant conformational motions in response to ligand-binding, these effects are more commonly subtle and challenging to study by traditional structural or biophysical methods. Molecular dynamics (MD) simulations are a powerful tool to bridge the gap between static protein-ligand structures and dynamical changes that govern NR function. Here, we summarize a handful of recent studies that apply MD simulations to study NRs. We present diverse methodologies for analyzing simulation data with a detailed examination of the information each method can yield. By delving into the strengths, limitations and unique contributions of these tools, this review provides guidance for extracting meaningful data from MD simulations to advance the goal of understanding the intricate mechanisms by which ligands orchestrate a range of functional outcomes in NRs.

Dynamics accelerate the kinetics of ion diffusion through channels: Continuous-time random walk models beyond the mean field approximation

Ion channels are proteins that play a significant role in physiological processes, including neuronal excitability and signal transduction. However, the precise mechanisms by which these proteins facilitate ion diffusion through cell membranes are not well understood. This is because experimental techniques to characterize ion channel activity operate on a time scale too large to understand the role of the various protein conformations on diffusion. Meanwhile, computational approaches operate on a time scale too short to rationalize the observed behavior at the microscopic scale. In this paper, we present a continuous-time random walk model that aims to bridge the scales between the atomistic models of ion channels and the experimental measurement of their conductance. We show how diffusion slows down in complex systems by using 3D lattices that map out the pore geometry of two channels: Nav1.7 and gramicidin. We also introduce spatial and dynamic site disorder to account for system heterogeneity beyond the mean field approximation. Computed diffusion coefficients show that an increase in spatial disorder slows down diffusion kinetics, while dynamic disorder has the opposite effect. Our results imply that microscopic or phenomenological models based on the potential of mean force data overlook the functional importance of protein dynamics on ion diffusion through channels.

Comparative Pore Structure and Dynamics for Bacterial Microcompartment Shell Protein Assemblies in Sheets or Shells

Bacterial microcompartments (BMCs) are protein-bound organelles found in some bacteria which encapsulate enzymes for enhanced catalytic activity. These compartments spatially sequester enzymes within semi-permeable shell proteins, analogous to many membrane-bound organelles. The shell proteins assemble into multimeric tiles; hexamers, trimers, and pentamers, and these tiles self-assemble into larger assemblies with icosahedral symmetry. While icosahedral shells are the predominant form in vivo, the tiles can also form nanoscale cylinders or sheets. The individual multimeric tiles feature central pores that are key to regulating transport across the protein shell. Our primary interest is to quantify pore shape changes in response to alternative component morphologies at the nanoscale. We use molecular modeling tools to develop atomically detailed models for both planar sheets of tiles and curved structures representative of the complete shells found in vivo. Subsequently, these models were animated using classical molecular dynamics simulations. From the resulting trajectories, we analyzed overall structural stability, water accessibility to individual residues, water residence time, and pore geometry for the hexameric and trimeric protein tiles from the Haliangium ochraceum model BMC shell. These exhaustive analyses suggest no substantial variation in pore structure or solvent accessibility between the flat and curved shell geometries. We additionally compare our analysis to hydroxyl radical footprinting data to serve as a check against our simulation results, highlighting specific residues where water molecules are bound for a long time. Although with little variation in morphology or water interaction, we propose that the planar and capsular morphology can be used interchangeably when studying permeability through BMC pores.

The energy landscape of the ribosome

The ribosome is a prototypical assembly that can be used to establish general principles and techniques for the study of biological molecular machines. Motivated by the fact that the dynamics of every biomolecule is governed by an underlying energy landscape, there has been great interest to understand and quantify ribosome energetics. In the present review, we will focus on theoretical and computational strategies for probing the interactions that shape the energy landscape of the ribosome, with an emphasis on more recent studies of the elongation cycle. These efforts include the application of quantum mechanical methods for describing chemical kinetics, as well as classical descriptions to characterize slower (microsecond to millisecond) large-scale (10-100 Å) rearrangements, where motion is described in terms of diffusion across an energy landscape. Together, these studies provide broad insights into the factors that control a diverse range of dynamical processes in this assembly.

AMBERff at Scale: Multimillion-Atom Simulations with AMBER Force Fields in NAMD

All-atom molecular dynamics (MD) simulations are an essential structural biology technique with increasing application to multimillion-atom systems, including viruses and cellular machinery. Classical MD simulations rely on parameter sets, such as the AMBER family of force fields (AMBERff), to accurately describe molecular motion. Here, we present an implementation of AMBERff for use in NAMD that overcomes previous limitations to enable high-performance, massively parallel simulations encompassing up to two billion atoms. Single-point potential energy comparisons and case studies on model systems demonstrate that the implementation produces results that are as accurate as running AMBERff in its native engine.

DNA supercoiling in bacteria: state of play and challenges from a viewpoint of physics based modeling

DNA supercoiling is central to many fundamental processes of living organisms. Its average level along the chromosome and over time reflects the dynamic equilibrium of opposite activities of topoisomerases, which are required to relax mechanical stresses that are inevitably produced during DNA replication and gene transcription. Supercoiling affects all scales of the spatio-temporal organization of bacterial DNA, from the base pair to the large scale chromosome conformation. Highlighted in vitro and in vivo in the 1960s and 1970s, respectively, the first physical models were proposed concomitantly in order to predict the deformation properties of the double helix. About fifteen years later, polymer physics models demonstrated on larger scales the plectonemic nature and the tree-like organization of supercoiled DNA. Since then, many works have tried to establish a better understanding of the multiple structuring and physiological properties of bacterial DNA in thermodynamic equilibrium and far from equilibrium. The purpose of this essay is to address upcoming challenges by thoroughly exploring the relevance, predictive capacity, and limitations of current physical models, with a specific focus on structural properties beyond the scale of the double helix. We discuss more particularly the problem of DNA conformations, the interplay between DNA supercoiling with gene transcription and DNA replication, its role on nucleoid formation and, finally, the problem of scaling up models. Our primary objective is to foster increased collaboration between physicists and biologists. To achieve this, we have reduced the respective jargon to a minimum and we provide some explanatory background material for the two communities.

Influence of Temperature-Guided SAM Growth on Wetting and Its Mass Transfer Models

The formation and growth of self-assembled monolayers (SAMs) composed of amphiphiles have garnered significant attention due to their diverse technical applications. This article reports the findings of molecular dynamics simulations aimed at elucidating the intricate relationship between the wetting behavior of amphiphiles, specifically n-alkanols, and the growth of their SAMs on a mica surface under varying temperature conditions. The investigation quantifies the structural characteristics of the formed SAMs, including density profiles, in-plane radial distribution functions, order parameters, and end-to-end length distributions of n-alkanol molecules within the SAM. Thermodynamic properties, such as the second virial coefficient and excess entropy, are examined in relation to temperature and time. The growth of the SAM is assessed by analyzing characteristic time scales at different temperatures and in-plane diffusion of n-alkanol molecules and utilizing classical theories of mass transfer to quantify the growth rate as a function of temperature. These results are then correlated with changes in the contact angle and spreading coefficient of n-alkanol droplets on the mica surface over time, providing insights into the impact of SAM growth on the wetting behavior and the mass transfer model of such systems.

The SIRAH force field: A suite for simulations of complex biological systems at the coarse-grained and multiscale levels

The different combinations of molecular dynamics simulations with coarse-grained representations have acquired considerable popularity among the scientific community. Especially in biocomputing, the significant speedup granted by simplified molecular models opened the possibility of increasing the diversity and complexity of macromolecular systems, providing realistic insights on large assemblies for more extended time windows. However, a holistic view of biological ensembles' structural and dynamic features requires a self-consistent force field, namely, a set of equations and parameters that describe the intra and intermolecular interactions among moieties of diverse chemical nature (i.e., nucleic and amino acids, lipids, solvent, ions, etc.). Nevertheless, examples of such force fields are scarce in the literature at the fully atomistic and coarse-grained levels. Moreover, the number of force fields capable of handling simultaneously different scales is restricted to a handful. Among those, the SIRAH force field, developed in our group, furnishes a set of topologies and tools that facilitate the setting up and running of molecular dynamics simulations at the coarse-grained and multiscale levels. SIRAH uses the same classical pairwise Hamiltonian function implemented in the most popular molecular dynamics software. In particular, it runs natively in AMBER and Gromacs engines, and porting it to other simulation packages is straightforward. This review describes the underlying philosophy behind the development of SIRAH over the years and across families of biological molecules, discussing current limitations and future implementations.

Unveiling the interdomain dynamics of type II DNA topoisomerase through all-atom simulations: Implications for understanding its catalytic cycle

Type IIA DNA topoisomerases are complex molecular nanomachines that manage topological states of the DNA molecule in the cell and play a crucial role in cellular processes such as cell division and transcription. They are also established targets of cancer chemotherapy. Starting from the available crystal structure of a fully catalytic topoisomerase IIA homodimer from Saccharomyces cerevisiae, we constructed three states of this molecular motor primarily changing the configurations of the DNA segment bound in the DNA gate and performed μs-long all-atom molecular simulations. A comprehensive analysis revealed a sliding motion within the DNA gate and a teamwork between the N-gate and DNA gate that may be associated with the necessary molecular events that allow passage of the T-segment of DNA. The observed movement of the ATPase dimer relative to the DNA domain was reflected in different interaction patterns between the K-loops of the transducer domain and the B-A-B form of the bound DNA. Based on the obtained results, we mapped simulated configurations to the structures in the proposed catalytic cycle through which type IIA topoisomerases exert their function and discussed the possible transition events. The results extend our understanding of the mechanism of action of type IIA topoisomerases and provide an atomistic interpretation of some of the observed features of these molecular motors.

Evolution, Stability, and Applicability of Surfactant Aggregates in Targeted Delivery

Self-assembly/self-aggregation of surfactant molecules in bulk and the vicinity of a surface has been a topic of interest for decades because of its utilization in numerous modern technical applications. In this article, the results of molecular dynamics simulations are reported to investigate the self-aggregation of sodium dodecyl sulfate (SDS) at an interface of mica and water. SDS molecules starting from lower to higher surface concentrations tend to create distinct aggregated structures in the vicinity of a mica surface. The structural properties, such as density profiles, radial distribution functions, and thermodynamic properties like excess entropy and second virial coefficient, are calculated to address the bits and pieces of the self-aggregation. The change in the free energy for aggregates of varied sizes approaching the surface from the bulk aqueous solution, along with the change in their shapes during the process in terms of change in the radius of gyration and its components, is reported respectively to model a generic pathway for a surfactant-based targeted delivery system.

Machines on Genes through the Computational Microscope

Macromolecular machines acting on genes are at the core of life's fundamental processes, including DNA replication and repair, gene transcription and regulation, chromatin packaging, RNA splicing, and genome editing. Here, we report the increasing role of computational biophysics in characterizing the mechanisms of "machines on genes", focusing on innovative applications of computational methods and their integration with structural and biophysical experiments. We showcase how state-of-the-art computational methods, including classical and ab initio molecular dynamics to enhanced sampling techniques, and coarse-grained approaches are used for understanding and exploring gene machines for real-world applications. As this review unfolds, advanced computational methods describe the biophysical function that is unseen through experimental techniques, accomplishing the power of the "computational microscope", an expression coined by Klaus Schulten to highlight the extraordinary capability of computer simulations. Pushing the frontiers of computational biophysics toward a pragmatic representation of large multimegadalton biomolecular complexes is instrumental in bridging the gap between experimentally obtained macroscopic observables and the molecular principles playing at the microscopic level. This understanding will help harness molecular machines for medical, pharmaceutical, and biotechnological purposes.

Engineering excitonically coupled dimers in an artificial protein for light harvesting via computational modeling

In photosynthesis, pigment-protein complexes achieve outstanding photoinduced charge separation efficiencies through a set of strategies in which excited states delocalization over multiple pigments ("excitons") and charge-transfer states play key roles. These concepts, and their implementation in bioinspired artificial systems, are attracting increasing attention due to the vast potential that could be tapped by realizing efficient photochemical reactions. In particular, de novo designed proteins provide a diverse structural toolbox that can be used to manipulate the geometric and electronic properties of bound chromophore molecules. However, achieving excitonic and charge-transfer states requires closely spaced chromophores, a non-trivial aspect since a strong binding with the protein matrix needs to be maintained. Here, we show how a general-purpose artificial protein can be optimized via molecular dynamics simulations to improve its binding capacity of a chlorophyll derivative, achieving complexes in which chromophores form two closely spaced and strongly interacting dimers. Based on spectroscopy results and computational modeling, we demonstrate each dimer is excitonically coupled, and propose they display signatures of charge-transfer state mixing. This work could open new avenues for the rational design of chromophore-protein complexes with advanced functionalities.

Mesoscale simulations: An indispensable approach to understand biomembranes

Computer simulation techniques form a versatile tool, a computational microscope, for exploring biological processes. This tool has been particularly effective in exploring different features of biological membranes. In recent years, thanks to elegant multiscale simulation schemes, some fundamental limitations of investigations by distinct simulation techniques have been resolved. As a result, we are now capable of exploring processes spanning multiple scales beyond the capacity of any single technique. In this perspective, we argue that mesoscale simulations require more attention and must be further developed to fill evident gaps in a quest toward simulating and modeling living cell membranes.

Molecular dynamics simulation of an entire cell

The ultimate microscope, directed at a cell, would reveal the dynamics of all the cell's components with atomic resolution. In contrast to their real-world counterparts, computational microscopes are currently on the brink of meeting this challenge. In this perspective, we show how an integrative approach can be employed to model an entire cell, the minimal cell, JCVI-syn3A, at full complexity. This step opens the way to interrogate the cell's spatio-temporal evolution with molecular dynamics simulations, an approach that can be extended to other cell types in the near future.

The dawn of interoperating spatial models in cell biology

Spatial simulations are becoming an increasingly ubiquitous component in the cycle of discovery, experimentation, and communication across the sciences. In cell biology, many researchers share a vision of developing multiscale models that recapitulate observable behaviors spanning from atoms to cells to tissues. For this dream to become a reality, however, simulation technologies must provide a means for integration and interoperability as they advance. Already, the field has developed numerous methods that span scales of length, time, and complexity to create an extensive body of effective simulation approaches, and although these approaches rarely interoperate, they collectively cover a large spectrum of knowledge that future models may handle in a more unified manner. Here, we discuss the importance of making the data, workflows, and outputs of spatial simulations shareable and interoperable; and how democratization could encourage diverse biologists to participate more easily in developing models to advance our understanding of biological systems.

Metamorphic proteins under a computational microscope: Lessons from a fold-switching RfaH protein

Metamorphic proteins constitute unexpected paradigms of the protein folding problem, as their sequences encode two alternative folds, which reversibly interconvert within biologically relevant timescales to trigger different cellular responses. Once considered a rare aberration, metamorphism may be common among proteins that must respond to rapidly changing environments, exemplified by NusG-like proteins, the only transcription factors present in every domain of life. RfaH, a specialized paralog of bacterial NusG, undergoes an all-α to all-β domain switch to activate expression of virulence and conjugation genes in many animal and plant pathogens and is the quintessential example of a metamorphic protein. The dramatic nature of RfaH structural transformation and the richness of its evolutionary history makes for an excellent model for studying how metamorphic proteins switch folds. Here, we summarize the structural and functional evidence that sparked the discovery of RfaH as a metamorphic protein, the experimental and computational approaches that enabled the description of the molecular mechanism and refolding pathways of its structural interconversion, and the ongoing efforts to find signatures and general properties to ultimately describe the protein metamorphome.

Protein structure prediction in the era of AI: Challenges and limitations when applying to in silico force spectroscopy

Mechanoactive proteins are essential for a myriad of physiological and pathological processes. Guided by the advances in single-molecule force spectroscopy (SMFS), we have reached a molecular-level understanding of how mechanoactive proteins sense and respond to mechanical forces. However, even SMFS has its limitations, including the lack of detailed structural information during force-loading experiments. That is where molecular dynamics (MD) methods shine, bringing atomistic details with femtosecond time-resolution. However, MD heavily relies on the availability of high-resolution structural data, which is not available for most proteins. For instance, the Protein Data Bank currently has 192K structures deposited, against 231M protein sequences available on Uniprot. But many are betting that this gap might become much smaller soon. Over the past year, the AI-based AlphaFold created a buzz on the structural biology field by being able to predict near-native protein folds from their sequences. For some, AlphaFold is causing the merge of structural biology with bioinformatics. Here, using an in silico SMFS approach pioneered by our group, we investigate how reliable AlphaFold structure predictions are to investigate mechanical properties of Staphylococcus bacteria adhesins proteins. Our results show that AlphaFold produce extremally reliable protein folds, but in many cases is unable to predict high-resolution protein complexes accurately. Nonetheless, the results show that AlphaFold can revolutionize the investigation of these proteins, particularly by allowing high-throughput scanning of protein structures. Meanwhile, we show that the AlphaFold results need to be validated and should not be employed blindly, with the risk of obtaining an erroneous protein mechanism.

Designing drugs when there is low data availability: one-shot learning and other approaches to face the issues of a long-term concern

INTRODUCTION

Modern drug discovery generally is accessed by useful information from previous large databases or uncovering novel data. The lack of biological and/or chemical data tends to slow the development of scientific research and innovation. Here, approaches that may help provide solutions to generate or obtain enough relevant data or improve/accelerate existing methods within the last five years were reviewed.

AREAS COVERED

One-shot learning (OSL) approaches, structural modeling, molecular docking, scoring function space (SFS), molecular dynamics (MD), and quantum mechanics (QM) may be used to amplify the amount of available data to drug design and discovery campaigns, presenting methods, their perspectives, and discussions to be employed in the near future.

EXPERT OPINION

Recent works have successfully used these techniques to solve a range of issues in the face of data scarcity, including complex problems such as the challenging scenario of drug design aimed at intrinsically disordered proteins and the evaluation of potential adverse effects in a clinical scenario. These examples show that it is possible to improve and kickstart research from scarce available data to design and discover new potential drugs.

py-MCMD: Python Software for Performing Hybrid Monte Carlo/Molecular Dynamics Simulations with GOMC and NAMD

py-MCMD, an open-source Python software, provides a robust workflow layer that manages communication of relevant system information between the simulation engines NAMD and GOMC and generates coherent thermodynamic properties and trajectories for analysis. To validate the workflow and highlight its capabilities, hybrid Monte Carlo/molecular dynamics (MC/MD) simulations are performed for SPC/E water in the isobaric-isothermal (NPT) and grand canonical (GC) ensembles as well as with Gibbs ensemble Monte Carlo (GEMC). The hybrid MC/MD approach shows close agreement with reference MC simulations and has a computational efficiency that is 2 to 136 times greater than traditional Monte Carlo simulations. MC/MD simulations performed for water in a graphene slit pore illustrate significant gains in sampling efficiency when the coupled-decoupled configurational-bias MC (CD-CBMC) algorithm is used compared with simulations using a single unbiased random trial position. Simulations using CD-CBMC reach equilibrium with 25 times fewer cycles than simulations using a single unbiased random trial position, with a small increase in computational cost. In a more challenging application, hybrid grand canonical Monte Carlo/molecular dynamics (GCMC/MD) simulations are used to hydrate a buried binding pocket in bovine pancreatic trypsin inhibitor. Water occupancies produced by GCMC/MD simulations are in close agreement with crystallographically identified positions, and GCMC/MD simulations have a computational efficiency that is 5 times better than MD simulations. py-MCMD is available on GitHub at https://github.com/GOMC-WSU/py-MCMD.

Computed Three-Dimensional Atomic Force Microscopy Images of Biopolymers Using the Jarzynski Equality

Three-dimensional atomic force microscopy (3D-AFM) has resolved three-dimensional distributions of solvent molecules at solid-liquid interfaces at the subnanometer scale. This method is now being extended to the imaging of biopolymer assemblies such as chromosomes or proteins in cells, with the expectation of being able to resolve their three-dimensional structures. Here, we have developed a computational method to simulate 3D-AFM images of biopolymers by using the Jarzynski equality. It is found that some parts of the fiber structure of biopolymers are indeed resolved in the 3D-AFM image. The dependency of 3D-AFM images on the vertical scanning velocity is investigated, and optimum scanning velocities are found. It is also clarified that forces in nonequilibrium processes are measured in 3D-AFM measurements when the dynamics of polymers are slower than the scanning of the probe.

Optimizing the Calculation of Free Energy Differences in Nonequilibrium Work SQM/MM Switching Simulations

A key step during indirect alchemical free energy simulations using quantum mechanical/molecular mechanical (QM/MM) hybrid potential energy functions is the calculation of the free energy difference ΔAlow→high between the low level (e.g., pure MM) and the high level of theory (QM/MM). A reliable approach uses nonequilibrium work (NEW) switching simulations in combination with Jarzynski's equation; however, it is computationally expensive. In this study, we investigate whether it is more efficient to use more shorter switches or fewer but longer switches. We compare results obtained with various protocols to reference free energy differences calculated with Crooks' equation. The central finding is that fewer longer switches give better converged results. As few as 200 sufficiently long switches lead to ΔAlow→high values in good agreement with the reference results. This optimized protocol reduces the computational cost by a factor of 40 compared to earlier work. We also describe two tools/ways of analyzing the raw data to detect sources of poor convergence. Specifically, we find it helpful to analyze the raw data (work values from the NEW switching simulations) in a quasi-time series-like manner. Principal component analysis helps to detect cases where one or more conformational degrees of freedom are different at the low and high level of theory.

Molecular insights into the resistance of phospholipid heads to the membrane penetration of graphene nanosheets

The interaction between nanomaterials and phospholipid membranes underlies many emerging biological applications. To what extent hydrophilic phospholipid heads shield the bilayer from the integration of hydrophobic nanomaterials remains unclear, and this open question contains important insights for understanding biological membrane physics. Here, we present molecular dynamics (MD) simulations to clarify the resistance of phospholipid heads to the membrane penetration of graphene nanosheets. With 130 simulation trials, we observed that ∼22% graphene nanosheets penetrate the POPC bilayer. Sharp corners of the nanosheets should have a lower energy barrier than nanosheet edges, but interestingly, the membrane penetration mainly starts from the edge-approaching orientation. We thoroughly analyzed the pentration pathway and propulsion, indicating that the membrane penetration of graphene nanosheets is dominated by the joint effects of nanosheet edges and corners. Furthermore, the molecular origin of the resistance is clarified by evaluating the bilayers of different phospholipids, which successfully correlates the penetration resistance of phospholipid heads with the correlated motions of neighboring phospholipids for the first time. These results are expected to inspire future studies on the dynamic behavior of phospholipids, bio-nano interfaces, and design of biological nanomaterials.

Contribution of Hydrophobic Interactions to Protein Mechanical Stability

The role of hydrophobic and polar interactions in providing thermodynamic stability to folded proteins has been intensively studied, but the relative contribution of these interactions to the mechanical stability is less explored. We used steered molecular dynamics simulations with constant-velocity pulling to generate force-extension curves of selected protein domains and monitor hydrophobic surface unravelling upon extension. Hydrophobic contribution was found to vary between one fifth and one third of the total force while the rest of the contribution is attributed primarily to hydrogen bonds. Moreover, hydrophobic force peaks were shifted towards larger protein extensions with respect to the force peaks attributed to hydrogen bonds. The higher importance of hydrogen bonds compared to hydrophobic interactions in providing mechanical resistance is in contrast with the relative importance of the hydrophobic interactions in providing thermodynamic stability of proteins. The different contributions of these interactions to the mechanical stability are explained by the steeper free energy dependence of hydrogen bonds compared to hydrophobic interactions on the relative positions of interacting atoms. Comparative analyses for several protein domains revealed that the variation of hydrophobic forces is modest, while the contribution of hydrogen bonds to the force peaks becomes increasingly important for mechanically resistant protein domains.

Molecular Modeling is an Enabling Approach to Complement and Enhance Channelopathy Research

Hundreds of human membrane proteins form channels that transport necessary ions and compounds, including drugs and metabolites, yet details of their normal function or how function is altered by genetic variants to cause diseases are often unknown. Without this knowledge, researchers are less equipped to develop approaches to diagnose and treat channelopathies. High-resolution computational approaches such as molecular modeling enable researchers to investigate channelopathy protein function, facilitate detailed hypothesis generation, and produce data that is difficult to gather experimentally. Molecular modeling can be tailored to each physiologic context that a protein may act within, some of which may currently be difficult or impossible to assay experimentally. Because many genomic variants are observed in channelopathy proteins from high-throughput sequencing studies, methods with mechanistic value are needed to interpret their effects. The eminent field of structural bioinformatics integrates techniques from multiple disciplines including molecular modeling, computational chemistry, biophysics, and biochemistry, to develop mechanistic hypotheses and enhance the information available for understanding function. Molecular modeling and simulation access 3D and time-dependent information, not currently predictable from sequence. Thus, molecular modeling is valuable for increasing the resolution with which the natural function of protein channels can be investigated, and for interpreting how genomic variants alter them to produce physiologic changes that manifest as channelopathies. © 2022 American Physiological Society. Compr Physiol 12:3141-3166, 2022.

Elastic versus brittle mechanical responses predicted for dimeric cadherin complexes

Cadherins are a superfamily of adhesion proteins involved in a variety of biological processes that include the formation of intercellular contacts, the maintenance of tissue integrity, and the development of neuronal circuits. These transmembrane proteins are characterized by ectodomains composed of a variable number of extracellular cadherin (EC) repeats that are similar but not identical in sequence and fold. E-cadherin, along with desmoglein and desmocollin proteins, are three classical-type cadherins that have slightly curved ectodomains and engage in homophilic and heterophilic interactions through an exchange of conserved tryptophan residues in their N-terminal EC1 repeat. In contrast, clustered protocadherins are straighter than classical cadherins and interact through an antiparallel homophilic binding interface that involves overlapped EC1 to EC4 repeats. Here we present molecular dynamics simulations that model the adhesive domains of these cadherins using available crystal structures, with systems encompassing up to 2.8 million atoms. Simulations of complete classical cadherin ectodomain dimers predict a two-phased elastic response to force in which these complexes first softly unbend and then stiffen to unbind without unfolding. Simulated α, β, and γ clustered protocadherin homodimers lack a two-phased elastic response, are brittle and stiffer than classical cadherins and exhibit complex unbinding pathways that in some cases involve transient intermediates. We propose that these distinct mechanical responses are important for function, with classical cadherin ectodomains acting as molecular shock absorbers and with stiffer clustered protocadherin ectodomains facilitating overlap that favors binding specificity over mechanical resilience. Overall, our simulations provide insights into the molecular mechanics of single cadherin dimers relevant in the formation of cellular junctions essential for tissue function.

Collective mechanical responses of cadherin-based adhesive junctions as predicted by simulations

Cadherin-based adherens junctions and desmosomes help stabilize cell-cell contacts with additional function in mechano-signaling, while clustered protocadherin junctions are responsible for directing neuronal circuits assembly. Structural models for adherens junctions formed by epithelial cadherin (CDH1) proteins indicate that their long, curved ectodomains arrange to form a periodic, two-dimensional lattice stabilized by tip-to-tip trans interactions (across junction) and lateral cis contacts. Less is known about the exact architecture of desmosomes, but desmoglein (DSG) and desmocollin (DSC) cadherin proteins are also thought to form ordered junctions. In contrast, clustered protocadherin (PCDH)-based cell-cell contacts in neuronal tissues are thought to be responsible for self-recognition and avoidance, and structural models for clustered PCDH junctions show a linear arrangement in which their long and straight ectodomains form antiparallel overlapped trans complexes. Here, we report all-atom molecular dynamics simulations testing the mechanics of minimalistic adhesive junctions formed by CDH1, DSG2 coupled to DSC1, and PCDHγB4, with systems encompassing up to 3.7 million atoms. Simulations generally predict a favored shearing pathway for the adherens junction model and a two-phased elastic response to tensile forces for the adhesive adherens junction and the desmosome models. Complexes within these junctions first unbend at low tensile force and then become stiff to unbind without unfolding. However, cis interactions in both the CDH1 and DSG2-DSC1 systems dictate varied mechanical responses of individual dimers within the junctions. Conversely, the clustered protocadherin PCDHγB4 junction lacks a distinct two-phased elastic response. Instead, applied tensile force strains trans interactions directly, as there is little unbending of monomers within the junction. Transient intermediates, influenced by new cis interactions, are observed after the main rupture event. We suggest that these collective, complex mechanical responses mediated by cis contacts facilitate distinct functions in robust cell-cell adhesion for classical cadherins and in self-avoidance signaling for clustered PCDHs.

Modeling the Dynamics of Protein–Protein Interfaces, How and Why?

Protein–protein assemblies act as a key component in numerous cellular processes. Their accurate modeling at the atomic level remains a challenge for structural biology. To address this challenge, several docking and a handful of deep learning methodologies focus on modeling protein–protein interfaces. Although the outcome of these methods has been assessed using static reference structures, more and more data point to the fact that the interaction stability and specificity is encoded in the dynamics of these interfaces. Therefore, this dynamics information must be taken into account when modeling and assessing protein interactions at the atomistic scale. Expanding on this, our review initially focuses on the recent computational strategies aiming at investigating protein–protein interfaces in a dynamic fashion using enhanced sampling, multi-scale modeling, and experimental data integration. Then, we discuss how interface dynamics report on the function of protein assemblies in globular complexes, in fuzzy complexes containing intrinsically disordered proteins, as well as in active complexes, where chemical reactions take place across the protein–protein interface.

Between Two Walls: Modeling the Adsorption Behavior of β-Glucosidase A on Bare and SAM-Functionalized Gold Surfaces

The efficient immobilization of enzymes on surfaces remains a complex but central issue in the biomaterials field, which requires us to understand this process at the atomic level. Using a multiscale approach combining all-atom molecular dynamics and coarse-grain Brownian dynamics simulations, we investigated the adsorption behavior of β-glucosidase A (βGA) on bare and self-assembled monolayer (SAM)-functionalized gold surfaces. We monitored the enzyme position and orientation during the molecular dynamics (MD) trajectories and measured the contacts it forms with both surfaces. While the adsorption process has little impact on the protein conformation, it can nonetheless perturb its mechanical properties and catalytic activity. Our results show that compared to the SAM-functionalized surface, the adsorption of βGA on bare gold is more stable, but less specific, and more likely to disrupt the enzyme's function. This observation emphasizes the fact that the structural organization of proteins at the solid interface is a key point when designing devices based on enzyme immobilization, as one must find an acceptable stability-activity trade-off.

Energy dissipation of osteopontin at a HAp mineral interface: Implications for bone biomechanics

Osteopontin (OPN) is a one of the most abundant non-collagenous proteins in the bone's organic matrix. OPN is responsible for mediating bonding at mineral interfaces in the extrafibrillar space and recent evidence shows that it is a major contributor to bone's fracture resistance. While several experimental studies have identified an important role for calcium ions in mediating energy dissipation in OPN protein networks, the underlying molecular mechanisms remain largely unknown. In the current study, the role of calcium ions on energy dissipation at OPN interface with hydroxyapatite (HAp) as the main bone mineral was investigated. For the first time, the three-dimensional structure of OPN proteins were predicted, and it was found that calcium ions greatly influenced the final protein configuration and energy dissipation performance. Under small deformation, the compact cOPN structure, resulting from calcium ions presence, facilitated greater energy dissipation through sacrificial bond breaking and mechanisms mediated by the surface-bound calcium. At larger deformation, the compact structure also enabled cOPN to dissipate higher energy. Moreover, it was found that phosphorylation of OPN played an important role in energy dissipation. While previous studies have shown that OPN dissipated energy by forming aggregate networks, this study also showed that network formation is not necessary and that individual OPN proteins can dissipate large amounts of energy at HAp interfaces.

Computer Simulations of Deep Eutectic Solvents: Challenges, Solutions, and Perspectives

Deep eutectic solvents (DESs) are one of the most rapidly evolving types of solvents, appearing in a broad range of applications, such as nanotechnology, electrochemistry, biomass transformation, pharmaceuticals, membrane technology, biocomposite development, modern 3D-printing, and many others. The range of their applicability continues to expand, which demands the development of new DESs with improved properties. To do so requires an understanding of the fundamental relationship between the structure and properties of DESs. Computer simulation and machine learning techniques provide a fruitful approach as they can predict and reveal physical mechanisms and readily be linked to experiments. This review is devoted to the computational research of DESs and describes technical features of DES simulations and the corresponding perspectives on various DES applications. The aim is to demonstrate the current frontiers of computational research of DESs and discuss future perspectives.

A stepwise docking molecular dynamics approach for simulating antibody recognition with substantial conformational changes

Conformational changes or rearrangements are common events during inter-biomolecular recognition. Tracking these changes are essential for exploring the allosteric mechanism and it is usually achieved by molecular dynamics simulation in silico. We previously identified a broad-neutralizing antibody against H5 influenza virus, 13D4, and solved the crystal structures of the free 13D4 Fab and its complex with hemagglutinin (HA). Structural comparison of the unbound and bound 13D4 Fabs showed that the heavy chain complementarity-determining region 3 (HCDR3) undergoes a substantial conformational rearrangement when it recognizes the receptor-binding site (RBS). Here, we used molecular dynamics (MD) to simulate the conformational changes that occur during antibody recognition. We showed that neither conventional MD nor steered MD could recapitulate the loop fitting of the RBS structure contour. Consequently, to simulate these challenging conformational changes, we engaged a stepwise docking MD method that allowed for the gradual docking of the ligand to receptor. This new method recapitulates the bound shape of the HCDR3 and provides the best approximation of the shape rendered by the co-crystal structure, with an RMSD of 0.926 Å. This strategy affords a flexible MD approach for simulating complicated conformational changes that occur during molecular recognition, and helps to provide an understanding of the involved allosteric mechanism.

Boosting the conformational sampling by combining replica exchange with solute tempering and well-sliced metadynamics

Methods that combine collective variable (CV) based enhanced sampling and global tempering approaches are used in speeding-up the conformational sampling and free energy calculation of large and soft systems with a plethora of energy minima. In this paper, a new method of this kind is proposed in which the well-sliced metadynamics approach (WSMTD) is united with replica exchange with solute tempering (REST2) method. WSMTD employs a divide-and-conquer strategy wherein high-dimensional slices of a free energy surface are independently sampled and combined. The method enables one to accomplish a controlled exploration of the CV-space with a restraining bias as in umbrella sampling, and enhance-sampling of one or more orthogonal CVs using a metadynamics like bias. The new hybrid method proposed here enables boosting the sampling of more slow degrees of freedom in WSMTD simulations, without the need to specify associated CVs, through a replica exchange scheme within the framework of REST2. The high-dimensional slices of the probability distributions of CVs computed from the united WSMTD and REST2 simulations are subsequently combined using the weighted histogram analysis method to obtain the free energy surface. We show that the new method proposed here is accurate, improves the conformational sampling, and achieves quick convergence in free energy estimates. We demonstrate this by computing the conformational free energy landscapes of solvated alanine tripeptide and Trp-cage mini protein in explicit water.

Photopharmacology of Ion Channels through the Light of the Computational Microscope

The optical control and investigation of neuronal activity can be achieved and carried out with photoswitchable ligands. Such compounds are designed in a modular fashion, combining a known ligand of the target protein and a photochromic group, as well as an additional electrophilic group for tethered ligands. Such a design strategy can be optimized by including structural data. In addition to experimental structures, computational methods (such as homology modeling, molecular docking, molecular dynamics and enhanced sampling techniques) can provide structural insights to guide photoswitch design and to understand the observed light-regulated effects. This review discusses the application of such structure-based computational methods to photoswitchable ligands targeting voltage- and ligand-gated ion channels. Structural mapping may help identify residues near the ligand binding pocket amenable for mutagenesis and covalent attachment. Modeling of the target protein in a complex with the photoswitchable ligand can shed light on the different activities of the two photoswitch isomers and the effect of site-directed mutations on photoswitch binding, as well as ion channel subtype selectivity. The examples presented here show how the integration of computational modeling with experimental data can greatly facilitate photoswitchable ligand design and optimization. Recent advances in structural biology, both experimental and computational, are expected to further strengthen this rational photopharmacology approach.

Interpretation of Single-Molecule Force Experiments on Proteins Using Normal Mode Analysis

Single-molecule force spectroscopy experiments allow protein folding and unfolding to be explored using mechanical force. Probably the most informative technique for interpreting the results of these experiments at the structural level makes use of steered molecular dynamics (MD) simulations, which can explicitly model the protein under load. Unfortunately, this technique is computationally expensive for many of the most interesting biological molecules. Here, we find that normal mode analysis (NMA), a significantly cheaper technique from a computational perspective, allows at least some of the insights provided by MD simulation to be gathered. We apply this technique to three non-homologous proteins that were previously studied by force spectroscopy: T4 lysozyme (T4L), Hsp70 and the glucocorticoid receptor domain (GCR). The NMA results for T4L and Hsp70 are compared with steered MD simulations conducted previously, and we find that we can recover the main results. For the GCR, which did not undergo MD simulation, our approach identifies substructures that correlate with experimentally identified unfolding intermediates. Overall, we find that NMA can make a valuable addition to the analysis toolkit for the structural analysis of single-molecule force experiments on proteins.

Plant Sterol Clustering Correlates with Membrane Microdomains as Revealed by Optical and Computational Microscopy

Local inhomogeneities in lipid composition play a crucial role in the regulation of signal transduction and membrane traffic. This is particularly the case for plant plasma membrane, which is enriched in specific lipids, such as free and conjugated forms of phytosterols and typical phytosphingolipids. Nevertheless, most evidence for microdomains in cells remains indirect, and the nature of membrane inhomogeneities has been difficult to characterize. We used a new push–pull pyrene probe and fluorescence lifetime imaging microscopy (FLIM) combined with all-atom multiscale molecular dynamics simulations to provide a detailed view on the interaction between phospholipids and phytosterol and the effect of modulating cellular phytosterols on membrane-associated microdomains and phase separation formation. Our understanding of the organization principles of biomembranes is limited mainly by the challenge to measure distributions and interactions of lipids and proteins within the complex environment of living cells. Comparing phospholipids/phytosterol compositions typical of liquid-disordered (Ld) and liquid-ordered (Lo) domains, we furthermore show that phytosterols play crucial roles in membrane homeostasis. The simulation work highlights how state-of-the-art modeling alleviates some of the prior concerns and how unrefuted discoveries can be made through a computational microscope. Altogether, our results support the role of phytosterols in the lateral structuring of the PM of plant cells and suggest that they are key compounds for the formation of plant PM microdomains and the lipid-ordered phase.

AI-driven multiscale simulations illuminate mechanisms of SARS-CoV-2 spike dynamics

We develop a generalizable AI-driven workflow that leverages heterogeneous HPC resources to explore the time-dependent dynamics of molecular systems. We use this workflow to investigate the mechanisms of infectivity of the SARS-CoV-2 spike protein, the main viral infection machinery. Our workflow enables more efficient investigation of spike dynamics in a variety of complex environments, including within a complete SARS-CoV-2 viral envelope simulation, which contains 305 million atoms and shows strong scaling on ORNL Summit using NAMD. We present several novel scientific discoveries, including the elucidation of the spike's full glycan shield, the role of spike glycans in modulating the infectivity of the virus, and the characterization of the flexible interactions between the spike and the human ACE2 receptor. We also demonstrate how AI can accelerate conformational sampling across different systems and pave the way for the future application of such methods to additional studies in SARS-CoV-2 and other molecular systems.

From System Modeling to System Analysis: The Impact of Resolution Level and Resolution Distribution in the Computer-Aided Investigation of Biomolecules

The ever increasing computer power, together with the improved accuracy of atomistic force fields, enables researchers to investigate biological systems at the molecular level with remarkable detail. However, the relevant length and time scales of many processes of interest are still hardly within reach even for state-of-the-art hardware, thus leaving important questions often unanswered. The computer-aided investigation of many biological physics problems thus largely benefits from the usage of coarse-grained models, that is, simplified representations of a molecule at a level of resolution that is lower than atomistic. A plethora of coarse-grained models have been developed, which differ most notably in their granularity; this latter aspect determines one of the crucial open issues in the field, i.e. the identification of an optimal degree of coarsening, which enables the greatest simplification at the expenses of the smallest information loss. In this review, we present the problem of coarse-grained modeling in biophysics from the viewpoint of system representation and information content. In particular, we discuss two distinct yet complementary aspects of protein modeling: on the one hand, the relationship between the resolution of a model and its capacity of accurately reproducing the properties of interest; on the other hand, the possibility of employing a lower resolution description of a detailed model to extract simple, useful, and intelligible information from the latter.

TorchMD: A Deep Learning Framework for Molecular Simulations

Molecular dynamics simulations provide a mechanistic description of molecules by relying on empirical potentials. The quality and transferability of such potentials can be improved leveraging data-driven models derived with machine learning approaches. Here, we present TorchMD, a framework for molecular simulations with mixed classical and machine learning potentials. All force computations including bond, angle, dihedral, Lennard-Jones, and Coulomb interactions are expressed as PyTorch arrays and operations. Moreover, TorchMD enables learning and simulating neural network potentials. We validate it using standard Amber all-atom simulations, learning an ab initio potential, performing an end-to-end training, and finally learning and simulating a coarse-grained model for protein folding. We believe that TorchMD provides a useful tool set to support molecular simulations of machine learning potentials. Code and data are freely available at github.com/torchmd.

Atomistic dynamics of a viral infection process: Release of membrane lytic peptides from a non-enveloped virus

Molecular simulations have played an instrumental role in uncovering the structural dynamics and physical properties of virus capsids. In this work, we move beyond equilibrium physicochemical characterization of a virus system to study a stage of the infection process that is required for viral proliferation. Despite many biochemical and functional studies, the molecular mechanism of host cell entry by non-enveloped viruses remains largely unresolved. Flock House virus (FHV) is a model system for non-enveloped viruses and is the subject of the current study. FHV infects through the acid-dependent endocytic pathway, where low pH triggers externalization of membrane-disrupting (γ) peptides from the capsid interior. Using all-atom equilibrium and enhanced sampling simulations, the mechanism and energetics of γ peptide liberation and the effect of pH on this process are investigated. Our computations agree with experimental findings and reveal nanoscopic details regarding the pH control mechanism, which are not readily accessible in experiments.

Molecular Simulations Guidelines for Biological Nanomaterials: From Peptides to Membranes

In studying biological processes and focusing on the molecular mechanisms at the basis of these, molecular dynamics (MD) simulations have demonstrated to be a very useful tool for the past 50 years. This suite of computational methods calculates the time-dependent evolution of a molecular system using physics-based first principles. In this chapter, we give a brief introduction to the theory and practical use of molecular dynamics simulations, highlighting the different models and algorithms that have been developed to tackle specific problems, with a special focus on classical force fields. Some examples of how simulations have been used in the past will help the reader in discerning their power, limitations, and significance.

CHARMM-GUI supports the Amber force fields

As part of our ongoing efforts to support diverse force fields and simulation programs in CHARMM-GUI, this work presents the development of FF-Converter to prepare Amber simulation inputs with various Amber force fields within the current CHARMM-GUI workflow. The currently supported Amber force fields are ff14SB/ff19SB (protein), Bsc1 (DNA), OL3 (RNA), GLYCAM06 (carbohydrate), Lipid17 (lipid), GAFF/GAFF2 (small molecule), TIP3P/TIP4P-EW/OPC (water), and 12-6-4 ions, and more will be added if necessary. The robustness and usefulness of this new CHARMM-GUI extension are demonstrated by two exemplary systems: a protein/N-glycan/ligand/membrane system and a protein/DNA/RNA system. Currently, CHARMM-GUI supports the Amber force fields only for the Amber program, but we will expand the FF-Converter functionality to support other simulation programs that support the Amber force fields.

Moltemplate: A Tool for Coarse-Grained Modeling of Complex Biological Matter and Soft Condensed Matter Physics

Coarse-grained models have long been considered indispensable tools in the investigation of biomolecular dynamics and assembly. However, the process of simulating such models is arduous because unconventional force fields and particle attributes are often needed, and some systems are not in thermal equilibrium. Although modern molecular dynamics programs are highly adaptable, software designed for preparing all-atom simulations typically makes restrictive assumptions about the nature of the particles and the forces acting on them. Consequently, the use of coarse-grained models has remained challenging. Moltemplate is a file format for storing coarse-grained molecular models and the forces that act on them, as well as a program that converts moltemplate files into input files for LAMMPS, a popular molecular dynamics engine. Moltemplate has broad scope and an emphasis on generality. It accommodates new kinds of forces as they are developed for LAMMPS, making moltemplate a popular tool with thousands of users in computational chemistry, materials science, and structural biology. To demonstrate its wide functionality, we provide examples of using moltemplate to prepare simulations of fluids using many-body forces, coarse-grained organic semiconductors, and the motor-driven supercoiling and condensation of an entire bacterial chromosome.