Application of hs-CRP in Assessment of Tissue Inflammation Following Implantation of Biodegradable Polymer in Experiment.

Implants made of biodegradable polymers are replaced by regenerating tissues through inflammation. The changes occurring in tissues and the organism are of practical interest for studying the biocompatibility of the material and searching for systemic markers in the blood that reflect inflammation in peri-implantation tissues. The highly sensitive C-reactive protein (hs-CRP) measurements in blood and morphometric studies of tissue surrounding the implant were carried out in the experiment within three months of implantation of a biopolymer consisting of polycaprolactone (PCL) and polytrimethylene carbonate (PTMC).

High-throughput molecular simulations of SARS-CoV-2 receptor binding domain mutants quantify correlations between dynamic fluctuations and protein expression.

Prediction of protein fitness from computational modeling is an area of active research in rational protein design. Here, we investigated whether protein fluctuations computed from molecular dynamics simulations can be used to predict the expression levels of SARS-CoV-2 receptor binding domain (RBD) mutants determined in the deep mutational scanning experiment of Starr et al. [Science (New York, N.

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed.

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published.

Structure-Based Rational Design of Constrained Peptides as TIM-3 Inhibitors.

Blocking the immunosuppressive function of T-cell immunoglobulin mucin-3 (TIM-3) is an established therapeutic strategy to maximize the efficacy of immune checkpoint inhibitors for cancer immunotherapy. Currently, effective inhibition of TIM-3 interactions relies on monoclonal antibodies (mAbs), which come with drawbacks such as immunogenicity risk, limited tumor penetration, and high manufacturing costs. Guided by the X-ray cocrystal structures of TIM-3 with mAbs, we report an structure-based rational design of constrained peptides as potent TIM-3 inhibitors.

Free Energy Simulations of Receptor-Binding Domain Opening of the SARS-CoV-2 Spike Indicate a Barrierless Transition with Slow Conformational Motions.

Infection by sarbecoviruses begins with the attachment of the homotrimeric viral "spike" protein to the angiotensin-converting enzyme 2 receptor on the surface of mammalian cells. This requires one or more receptor-binding domains (RBDs) to be in the open (up) position. Here, we present the results of long molecular dynamics simulations with umbrella sampling (US) to compute a one-dimensional free energy profile of RBD opening/closing and the associated transition times.

Molecular Simulations of Conformational Transitions within the Insulin Receptor Kinase Reveal Consensus Features in a Multistep Activation Pathway.

Modulating the transitions between active and inactive conformations of protein kinases is the primary means of regulating their catalytic activity, achieved by phosphorylation of the activation loop (A-loop). To elucidate the mechanism of this conformational activation, we applied the string method to determine the conformational transition path of insulin receptor kinase between the active and inactive conformations and the corresponding free-energy profiles with and without A-loop phosphorylation. The conformational change was found to proceed in three sequential steps: first, the flipping of the DFG motif of the active site; second, rotation of the A-loop; finally, the inward movement of the αC helix.

On the Rapid Calculation of Binding Affinities for Antigen and Antibody Design and Affinity Maturation Simulations.

The accurate and efficient calculation of protein-protein binding affinities is an essential component in antibody and antigen design and optimization, and in computer modeling of antibody affinity maturation. Such calculations remain challenging despite advances in computer hardware and algorithms, primarily because proteins are flexible molecules, and thus, require explicit or implicit incorporation of multiple conformational states into the computational procedure. The astronomical size of the amino acid sequence space further compounds the challenge by requiring predictions to be computed within a short time so that many sequence variants can be tested.

ppdx: Automated modeling of protein-protein interaction descriptors for use with machine learning.

This paper describes ppdx, a python workflow tool that combines protein sequence alignment, homology modeling, and structural refinement, to compute a broad array of descriptors for characterizing protein-protein interactions. The descriptors can be used to predict various properties of interest, such as protein-protein binding affinities, or inhibitory concentrations (IC ), using approaches that range from simple regression to more complex machine learning models. The software is highly modular.

Multiscale affinity maturation simulations to elicit broadly neutralizing antibodies against HIV.

The design of vaccines against highly mutable pathogens, such as HIV and influenza, requires a detailed understanding of how the adaptive immune system responds to encountering multiple variant antigens (Ags). Here, we describe a multiscale model of B cell receptor (BCR) affinity maturation that employs actual BCR nucleotide sequences and treats BCR/Ag interactions in atomistic detail. We apply the model to simulate the maturation of a broadly neutralizing Ab (bnAb) against HIV.

A Coarse-Grained Model of Affinity Maturation Indicates the Importance of B-Cell Receptor Avidity in Epitope Subdominance.

The elicitation of broadly neutralizing antibodies (bnAbs) is a major goal in the design of vaccines against rapidly-mutating viruses. In the case of influenza, many bnAbs that target conserved epitopes on the stem of the hemagglutinin protein (HA) have been discovered. However, these antibodies are rare, are not boosted well upon reinfection, and often have low neutralization potency, compared to strain-specific antibodies directed to the HA head.

Molecular Simulation of Stapled Peptides.

Constrained peptides represent a relatively new class of biologic therapeutics, which have the potential to overcome several limitations of small-molecule drugs, and of designed antibodies. Because of their modest size, the rational design of such peptides is becoming increasingly amenable to computer simulation; multi-microsecond molecular dynamic (MD) simulations are now routinely possible on consumer-grade graphical processors (GPUs). Here, we describe the procedures for performing and analyzing MD simulations of hydrocarbon-stapled peptides using the CHARMM energy function, in isolation and in complex with a binding partner, to investigate their conformational properties and to compute changes in their binding affinity upon mutation.

Dynamic Connection between Enzymatic Catalysis and Collective Protein Motions.

Enzymes employ a wide range of protein motions to achieve efficient catalysis of chemical reactions. While the role of collective protein motions in substrate binding, product release, and regulation of enzymatic activity is generally understood, their roles in catalytic steps per se remain uncertain. Here, molecular dynamics simulations, enzyme kinetics, X-ray crystallography, and nuclear magnetic resonance spectroscopy are combined to elucidate the catalytic mechanism of adenylate kinase and to delineate the roles of catalytic residues in catalysis and the conformational change in the enzyme.

A restrained locally enhanced sampling method (RLES) for finding free energy minima in complex systems.

We present an extension of the locally enhanced sampling method. A restraint potential is introduced to drive the many-replica system to the canonical ensemble corresponding to the physical, single-replica system. Convergence properties are demonstrated using a model rugged two-dimensional potential, for which sampling by conventional equilibrium molecular dynamics is inefficient.

The trajectory of intrahelical lesion recognition and extrusion by the human 8-oxoguanine DNA glycosylase.

Efficient search for DNA damage embedded in vast expanses of the DNA genome presents one of the greatest challenges to DNA repair enzymes. We report here crystal structures of human 8-oxoguanine (oxoG) DNA glycosylase, hOGG1, that interact with the DNA containing the damaged base oxoG and the normal base G while they are nested in the DNA helical stack. The structures reveal that hOGG1 engages the DNA using different protein-DNA contacts from those observed in the previously determined lesion recognition complex and other hOGG1-DNA complexes.

Microsecond Molecular Dynamics Simulations of Proteins Using a Quasi-Equilibrium Solvation Shell Model.

We describe the development and implementation of a quasi-equilibrium hydration shell model of biomolecular solvation with adaptive boundaries. Applying the model to microsecond-long molecular dynamics simulations of several protein systems of varying complexity, we find that the model simulation results are of comparable quality to those obtained from simulations of fully solvated systems, but at a reduced computational cost. We discuss the dominant sources of error in the model and outline directions for future improvements.

Cytosolic expression, solution structures, and molecular dynamics simulation of genetically encodable disulfide-rich de novo designed peptides.

Disulfide-rich peptides represent an important protein family with broad pharmacological potential. Recent advances in computational methods have made it possible to design new peptides which adopt a stable conformation de novo. Here, we describe a system to produce disulfide-rich de novo peptides using Escherichia coli as the expression host.

Structure of the multidrug transporter and its use for inhibitor peptide design.

Small multidrug resistance (SMR) pumps represent a minimal paradigm of proton-coupled membrane transport in bacteria, yet no high-resolution structure of an SMR protein is available. Here, atomic-resolution structures of the efflux-multidrug resistance E () multidrug transporter in ligand-bound form are refined using microsecond molecular dynamics simulations biased using low-resolution data from X-ray crystallography. The structures are compatible with existing mutagenesis data as well as NMR and biochemical experiments, including pKas of the catalytic glutamate residues and the dissociation constant ([Formula: see text]) of the tetraphenylphosphonium cation.

Role of framework mutations and antibody flexibility in the evolution of broadly neutralizing antibodies.

Eliciting antibodies that are cross reactive with surface proteins of diverse strains of highly mutable pathogens (e.g., HIV, influenza) could be key for developing effective universal vaccines.

Coexisting origins of subdiffusion in internal dynamics of proteins.

Subdiffusion in conformational dynamics of proteins is observed both experimentally and in simulations. Although its origin has been attributed to multiple mechanisms, including trapping on a rugged energy landscape, fractional Brownian noise, or a fractal topology of the energy landscape, it is unclear which of these, if any, is most relevant. To obtain insights into the actual mechanism, we introduce an analytically tractable hierarchical trapping model and apply it to molecular dynamics simulation trajectories of three proteins in solution.

Regulation and Plasticity of Catalysis in Enzymes: Insights from Analysis of Mechanochemical Coupling in Myosin.

The mechanism of ATP hydrolysis in the myosin motor domain is analyzed using a combination of DFTB3/CHARMM simulations and enhanced sampling techniques. The motor domain is modeled in the pre-powerstroke state, in the post-rigor state, and as a hybrid based on the post-rigor state with a closed nucleotide-binding pocket. The ATP hydrolysis activity is found to depend on the positioning of nearby water molecules, and a network of polar residues facilitates proton transfer and charge redistribution during hydrolysis.

QM/MM free energy simulations: recent progress and challenges.

Due to the higher computational cost relative to pure molecular mechanical (MM) simulations, hybrid quantum mechanical/molecular mechanical (QM/MM) free energy simulations particularly require a careful consideration of balancing computational cost and accuracy. Here we review several recent developments in free energy methods most relevant to QM/MM simulations and discuss several topics motivated by these developments using simple but informative examples that involve processes in water. For chemical reactions, we highlight the value of invoking enhanced sampling technique (e.

A Simple and Accurate Method To Calculate Free Energy Profiles and Reaction Rates from Restrained Molecular Simulations of Diffusive Processes.

A method is developed to obtain simultaneously free energy profiles and diffusion constants from restrained molecular simulations in diffusive systems. The method is based on low-order expansions of the free energy and diffusivity as functions of the reaction coordinate. These expansions lead to simple analytical relationships between simulation statistics and model parameters.

Role of Protein Dynamics in Allosteric Control of the Catalytic Phosphoryl Transfer of Insulin Receptor Kinase.

The catalytic and allosteric mechanisms of insulin receptor kinase (IRK) are investigated by a combination of ab initio and semiempirical quantum mechanical and molecular mechanical (QM/MM) methods and classical molecular dynamics (MD) simulations. The simulations reveal that the catalytic reaction proceeds in two steps, starting with the transfer of a proton from substrate Tyr to the catalytic Asp1132, followed by the phosphoryl transfer from ATP to substrate Tyr. The enhancement of the catalytic rate of IRK upon phosphorylations in the enzyme's activation loop is found to occur mainly via changes to the free energy landscape of the proton transfer step, favoring the proton transfer in the fully phosphorylated enzyme.

Investigations of α-helix↔β-sheet transition pathways in a miniprotein using the finite-temperature string method.

A parallel implementation of the finite-temperature string method is described, which takes into account the invariance of coordinates with respect to rigid-body motions. The method is applied to the complex α-helix↔β-sheet transition in a β-hairpin miniprotein in implicit solvent, which exhibits much of the complexity of conformational changes in proteins. Two transition paths are considered, one derived from a linear interpolant between the endpoint structures and the other derived from a targeted dynamics simulation.

A simplified confinement method for calculating absolute free energies and free energy and entropy differences.

A simple and robust formulation of the path-independent confinement method for the calculation of free energies is presented. The simplified confinement method (SCM) does not require matrix diagonalization or switching off the molecular force field, and has a simple convergence criterion. The method can be readily implemented in molecular dynamics programs with minimal or no code modifications.