Assessing AF2's ability to predict structural ensembles of proteins.

Recent breakthroughs in protein structure prediction have enhanced the precision and speed at which protein configurations can be determined. Additionally, molecular dynamics (MD) simulations serve as a crucial tool for capturing the conformational space of proteins, providing valuable insights into their structural fluctuations. However, the scope of MD simulations is often limited by the accessible timescales and the computational resources available, posing challenges to comprehensively exploring protein behaviors.

Structure-Based Protein Assembly Simulations Including Various Binding Sites and Conformations.

Many biological functions are mediated by large complexes formed by multiple proteins and other cellular macromolecules. Recent progress in experimental structure determination, as well as in integrative modeling and protein structure prediction using deep learning approaches, has resulted in a rapid increase in the number of solved multiprotein assemblies. However, the assembly process of large complexes from their components is much less well-studied.

Liquid-Vapor Coexistence and Spontaneous Evaporation at Atmospheric Pressure of Common Rigid Three-Point Water Models in Molecular Simulations.

Molecular dynamics (MD) simulations are widely used to investigate molecular systems at atomic resolution including biomolecular structures, drug-receptor interactions, and novel materials. Frequently, MD simulations are performed in an aqueous solution with explicit models of water molecules. Commonly, such models are parameterized to reproduce the liquid phase of water under ambient conditions.

The Role of Force Fields and Water Models in Protein Folding and Unfolding Dynamics.

Protein folding is a fascinating, not fully understood phenomenon in biology. Molecular dynamics (MD) simulations are an invaluable tool to study conformational changes in atomistic detail, including folding and unfolding processes of proteins. However, the accuracy of the conformational ensembles derived from MD simulations inevitably relies on the quality of the underlying force field in combination with the respective water model.

Structural mechanism of Fab domain dissociation as a measure of interface stability.

Therapeutic antibodies should not only recognize antigens specifically, but also need to be free from developability issues, such as poor stability. Thus, the mechanistic understanding and characterization of stability are critical determinants for rational antibody design. In this study, we use molecular dynamics simulations to investigate the melting process of 16 antigen binding fragments (Fabs).

Challenges in antibody structure prediction.

Advances in structural biology and the exponential increase in the amount of high-quality experimental structural data available in the Protein Data Bank has motivated numerous studies to tackle the grand challenge of predicting protein structures. In 2020 AlphaFold2 revolutionized the field using a combination of artificial intelligence and the evolutionary information contained in multiple sequence alignments. Antibodies are one of the most important classes of biotherapeutic proteins.

Bispecific antibodies-effects of point mutations on CH3-CH3 interface stability.

A new format of therapeutic proteins is bispecific antibodies, in which two different heavy chains heterodimerize to obtain two different binding sites. Therefore, it is crucial to understand and optimize the third constant domain (CH3-CH3) interface to favor heterodimerization over homodimerization, and to preserve the physicochemical properties, as thermal stability. Here, we use molecular dynamics simulations to investigate the dissociation process of 19 CH3-CH3 crystal structures that differ from each other in few point mutations.

Complexes++: Efficient and versatile coarse-grained simulations of protein complexes and their dense solutions.

The interior of living cells is densely filled with proteins and their complexes, which perform multitudes of biological functions. We use coarse-grained simulations to reach the system sizes and time scales needed to study protein complexes and their dense solutions and to interpret experiments. To take full advantage of coarse-graining, the models have to be efficiently implemented in simulation engines that are easy to use, modify, and extend.

Increase of Radiative Forcing through Midinfrared Absorption by Stable CO Dimers?

We performed matrix-isolation infrared (MI-IR) spectroscopy of carbon dioxide monomers, CO, and dimers, (CO), trapped in neon and in air. On the basis of vibration configuration interaction (VCI) calculations accounting for mode coupling and anharmonicity, we identify additional infrared-active bands in the MI-IR spectra due to the (CO) dimer. These bands are satellite bands next to the established CO monomer bands, which appear in the infrared window of Earth's atmosphere at around 4 and 15 μm.

Corrigendum: Comparing Antibody Interfaces to Inform Rational Design of New Antibody Formats.

[This corrects the article DOI: 10.3389/fmolb.2022.

Comparing Antibody Interfaces to Inform Rational Design of New Antibody Formats.

As the current biotherapeutic market is dominated by antibodies, the design of different antibody formats, like bispecific antibodies and other new formats, represent a key component in advancing antibody therapy. When designing new formats, a targeted modulation of pairing preferences is key. Several existing approaches are successful, but expanding the repertoire of design possibilities would be desirable.

Germline-Dependent Antibody Paratope States and Pairing Specific V-V Interface Dynamics.

Antibodies have emerged as one of the fastest growing classes of biotherapeutic proteins. To improve the rational design of antibodies, we investigate the conformational diversity of 16 different germline combinations, which are composed of 4 different kappa light chains paired with 4 different heavy chains. In this study, we systematically show that different heavy and light chain pairings strongly influence the paratope, interdomain interaction patterns and the relative V-V interface orientations.

Shark Antibody Variable Domains Rigidify Upon Affinity Maturation-Understanding the Potential of Shark Immunoglobulins as Therapeutics.

Sharks and other cartilaginous fish are the phylogenetically oldest living organisms that have antibodies as part of their adaptive immune system. As part of their humoral adaptive immune response, they produce an immunoglobulin, the so-called immunoglobulin new antigen receptor (IgNAR), a heavy-chain only antibody. The variable domain of an IgNAR, also known as V , binds the antigen as an independent soluble domain.

Implementation of the Freely Jointed Chain Model to Assess Kinetics and Thermodynamics of Thermosensitive Coil-Globule Transition by Markov States.

We revived and implemented a method developed by Kuhn in 1934, originally only published in German, that is, the so-called "freely jointed chain" model. This approach turned out to be surprisingly useful for analyzing state-of-the-art computer simulations of the thermosensitive coil-globule transition of -Isopropylacrylamide 20-mer. Our atomistic computer simulations are orders of magnitude longer than those of previous studies and lead to a reliable description of thermodynamics and kinetics at many different temperatures.

Conformational Shifts of Stacked Heteroaromatics: Vacuum vs. Water Studied by Machine Learning.

Stacking interactions play a crucial role in drug design, as we can find aromatic cores or scaffolds in almost any available small molecule drug. To predict optimal binding geometries and enhance stacking interactions, usually high-level quantum mechanical calculations are performed. These calculations have two major drawbacks: they are very time consuming, and solvation can only be considered using implicit solvation.

pyPolyBuilder: Automated Preparation of Molecular Topologies and Initial Configurations for Molecular Dynamics Simulations of Arbitrary Supramolecules.

The construction of a molecular topology file is a prerequisite for any classical molecular dynamics simulation. However, the generation of such a file may be very challenging at times, especially for large supramolecules. While many tools are available to provide topologies for large proteins and other biomolecules, the scientific community researching nonbiological systems is not equally well equipped.

X-Entropy: A Parallelized Kernel Density Estimator with Automated Bandwidth Selection to Calculate Entropy.

is a Python package used to calculate the entropy of a given distribution, in this case, based on the distribution of dihedral angles. The dihedral entropy facilitates an alignment-independent measure of local protein flexibility. The key feature of our approach is a Gaussian kernel density estimation (KDE) using a plug-in bandwidth selection, which is fully implemented in a C++ backend and parallelized with OpenMP.

Surprisingly Fast Interface and Elbow Angle Dynamics of Antigen-Binding Fragments.

Fab consist of a heavy and light chain and can be subdivided into a variable (V and V ) and a constant region (C 1 and C ). The variable region contains the complementarity-determining region (CDR), which is formed by six hypervariable loops, shaping the antigen binding site, the paratope. Apart from the CDR loops, both the elbow angle and the relative interdomain orientations of the V -V and the C 1-C domains influence the shape of the paratope.

Thermosensitive Hydration of Four Acrylamide-Based Polymers in Coil and Globule Conformations.

To characterize the thermosensitive coil-globule transition in atomistic detail, the conformational dynamics of linear polymer chains of acrylamide-based polymers have been investigated at multiple temperatures. Therefore, molecular dynamic simulations of 30mers of polyacrylamide (AAm), poly--methylacrylamide (NMAAm), poly--ethylacrylamide (NEAAm), and poly--isopropylacrylamide (NIPAAm) have been performed at temperatures ranging from 250 to 360 K for 2 μs. While two of the polymers are known to exhibit thermosensitivity (NEAAm, NIPAAm), no thermosensitivity is observed for AAm and NMAAm in aqueous solution.

Macrocycle Cell Permeability Measured by Solvation Free Energies in Polar and Apolar Environments.

The relation of surface polarity and conformational preferences is decisive for cell permeability and thus bioavailability of macrocyclic drugs. Here, we employ grid inhomogeneous solvation theory (GIST) to calculate solvation free energies for a series of six macrocycles in water and chloroform as a measure of passive membrane permeability. We perform accelerated molecular dynamics simulations to capture a diverse structural ensemble in water and chloroform, allowing for a direct profiling of solvent-dependent conformational preferences.

Quantifying Protein-Protein Interactions in Molecular Simulations.

Interactions among proteins, nucleic acids, and other macromolecules are essential for their biological functions and shape the physicochemcial properties of the crowded environments inside living cells. Binding interactions are commonly quantified by dissociation constants , and both binding and nonbinding interactions are quantified by second osmotic virial coefficients . As a measure of nonspecific binding and stickiness, is receiving renewed attention in the context of so-called liquid-liquid phase separation in protein and nucleic acid solutions.

Coil-Globule Transition Thermodynamics of Poly(-isopropylacrylamide).

Thermosensitive polymers such as poly(-isopropylacrylamide) (PNIPAM) undergo a phase transition in aqueous solution from a random-coil structural ensemble to a globule structural ensemble at the lower critical solution temperature (LCST). Above this temperature, PNIPAM agglomerates and becomes insoluble, whereas it is soluble below the temperature. Thus, thermosensitive polymers represent essential targets for several applications, e.