Chapter MM-GB(PB)SA Calculations of Protein-Ligand Binding Free Energies

Optical physic

MM-GB(PB)SA Calculations of Protein-Ligand Binding Free Energies

Optical physic

CHARMM: The biomolecular simulation program

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983. © 2009 Wiley Periodicals, Inc.J Comput Chem, 2009.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/63074/1/21287_ftp.pd

The flexibility of myosin 7a

Myosin 7a is a molecular motor found in hair cells of the ear and the photoreceptor cells of the eye. Myosin 7a is comprised of an actin-binding motor domain, a lever; which is composed of 5 IQ motifs that can potentially bind 5 light chains followed by a single alpha helical (SAH) domain, and a tail composed of 2 MyTH4-FERM domains. The lever is an essential mechanical element in myosin 7a function, but an understanding of its mechanical properties and how these derive from its substructure is lacking. It has been observed in vitro that myosin 7a is able to regulate its activity through a head-tail interaction. How the flexibility of the sub-domains of the lever allows the molecule to fold up is not completely understood. To address this, the first aim of this study was to look for evidence of novel light chain binding partners in myosin 7a, which revealed calmodulin to be the preferred light chain. My second aim was to study the structure and flexibility of the lever of full-length myosin 7a using single-particle image processing of images from negative stain electron microscopy (EM). Image averaging revealed the lever to be much shorter than expected. Additionally, there was evidence of thermally-driven flexing at the motor-lever junction. A stiffness of 78 pN.nm.rad-2 for the flexing was inferred, which represents a significant compliance in the head. An investigation into lever bending analysis, by monitoring the decay of tangent-tangent correlations of the lever shapes, yielded a persistence length of 38 ± 3 nm. Finally, long time molecular dynamics (MD) simulations were compared with a novel coarse-grained (CG) simulation technique called Fluctuating Finite Element Analysis (FFEA), which treats proteins as visco-elastic continua subject to thermal noise to probe the flexibility of myosin 7a. FFEA allows sufficiently long time simulations that are computationally less expensive than corresponding all-atom MD simulations to allow myosin 7a to explore its full range of configurations. Extraction of flexibility data from all-atom MD simulations calculated the bending stiffness of the SAH domain to be 60.5 pN.nm2, with reasonable overlap of the major modes of motion between the all-atom and CG simulation types

Bioinformatics

This book is divided into different research areas relevant in Bioinformatics such as biological networks, next generation sequencing, high performance computing, molecular modeling, structural bioinformatics, molecular modeling and intelligent data analysis. Each book section introduces the basic concepts and then explains its application to problems of great relevance, so both novice and expert readers can benefit from the information and research works presented here

Navigating the Extremes of Biological Datasets for Reliable Structural Inference and Design

Structural biologists currently confront serious challenges in the effective interpretation of experimental data due to two contradictory situations: a severe lack of structural data for certain classes of proteins, and an incredible abundance of data for other classes. The challenge with small data sets is how to extract sufficient information to draw meaningful conclusions, while the challenge with large data sets is how to curate, categorize, and search the data to allow for its meaningful interpretation and application to scientific problems. Here, we develop computational strategies to address both sparse and abundant data sets. In the category of sparse data sets, we focus our attention on the problem of transmembrane (TM) protein structure determination. As X-ray crystallography and NMR data is notoriously difficult to obtain for TM proteins, we develop a novel algorithm which uses low-resolution data from protein cross-linking or scanning mutagenesis studies to produce models of TM helix oligomers and show that our method produces models with an accuracy on par with X-ray crystallography or NMR for a test set of known TM proteins. Turning to instances of data abundance, we examine how to mine the vast stores of protein structural data in the Protein Data Bank (PDB) to aid in the design of proteins with novel binding properties. We show how the identification of an anion binding motif in an antibody structure allowed us to develop a phosphate binding module that can be used to produce novel antibodies to phosphorylated peptides - creating antibodies to 7 novel phospho-peptides to illustrate the utility of our approach. We then describe a general strategy for designing binders to a target protein epitope based upon recapitulating protein interaction geometries which are over-represented in the PDB. We follow this by using data describing the transition probabilities of amino acids to develop a novel set of degenerate codons to create more efficient gene libraries. We conclude by describing a novel, real-time, all-atom structural search engine, giving researchers the ability to quickly search known protein structures for a motif of interest and providing a new interactive paradigm of protein design

Estudios computacionales de mecanismos moleculares de la inmunidad innata

Tesis inédita de la Universidad Complutense de Madrid, Facultad de Farmacia, leída el 20-12-2022Antimicrobial Resistance (AMR) is a worldwide health emergency. ESKAPE pathogens include the most relevant AMR bacterial families. In particular, Gram-negative bacteria stand out due to their cell envelope complexity, which exhibits strong resistance to antimicrobials. A key element for AMR is the chemical structure of bacterial lipopolysaccharide (LPS), and the phospholipid composition of the membrane, inflecting the membrane permeability to antibiotics. We have applied coarse-grained molecular dynamics simulations to capture the role of the phospholipid composition and lipid A structure in the membrane properties and morphology of ESKAPE Gram-negative bacterial vesicles. Moreover, the reported antimicrobial peptides Cecropin B1, JB95, and PTCDA1-kf were used to unveil their implications for membrane disruption. This study opens a promising starting point for understanding the molecular keys of bacterial membranes and promoting the discovery of new antimicrobials to overcome AMR...La resistencia a los antimicrobianos (AMR) es una emergencia sanitaria mundial. Los patógenos ESKAPE incluyen las familias bacterianas más resistentes a antibióticos y son altamente virulentas. En particular, las bacterias Gram negativas destacan por la complejidad de su pared celular, que presenta una fuerte resistencia frente a los antibióticos. Un elemento clave para la AMR es la estructura química del lipopolisacárido bacteriano (LPS) y la composición de los fosfolípidos de la membrana bacteriana, que influyen en su permeabilidad a los antibióticos. Se han empleado simulaciones de dinámica molecular de grano grueso para captar el papel de la composición de los fosfolípidos y la estructura del LPS en las propiedades y morfología de modelos de vesículas bacterianas Gram negativas ESKAPE. Además, se han empleado los péptidos antimicrobianos Cecropin B1, JB95 y PTCDA1-kf para desvelar su mecanismo disrupción de la membrana bacteriana. Este estudio abre un prometedor punto de partida para comprender las claves moleculares de la resistencia en membranas bacterianas y acelerar el descubrimiento de nuevos antibióticos para hacer frente a la AMR...Fac. de FarmaciaTRUEunpu

A Structural Perspective of Antibody Neutralization of Dengue Virus

The four dengue viruses: DENV) are mosquito-borne flaviviruses and are considered the world\u27s most significant arboviruses in terms of worldwide disease burden. Symptoms of dengue disease are classified into dengue fever, a mild, febrile illness, and the potentially fatal severe dengue, which can include hemorrhaging and shock. Antibody protection against DENV correlates with the production of neutralizing antibodies against the envelope: E) glycoprotein. To understand the role of antibodies in DENV infection, we sought to dissect the relationship between epitope and function. Virologic studies had identified that the most potently neutralizing antibodies are against domain III: DIII) of the E protein. We have identified five epitopes within DENV DIII. Our data suggests that the most potently neutralizing antibodies are specific for a single serotype, while cross-reactive antibodies are relatively poorly neutralizing. Additionally, we were surprised to define neutralizing epitopes that were shown to be inaccessible on the surface of the virion in cryo-electron microscopy studies. Fine epitope mapping was used to define the epitopes of a panel of existing DENV-2 antibodies. Antibodies against the lateral ridge were the most potently neutralizing antibodies and reacted only with the DENV-2 serotype. The second epitope was centered on the DIII A-strand, and antibodies against this epitope reacted with several serotypes of DENV. Several poorly neutralizing antibodies reacted to all four DENV serotypes, as well as West Nile virus, a related flavivirus, mapped to the highly conserved AB loop of DIII. We expanded our studies of DIII-specific antibodies to the DENV-1 serotype. One antibody, E106, potently neutralized the five DENV-1 strains representing the five genotypes, and bound a composite epitope of the lateral ridge and A-strand epitopes. Despite the potency of E106-mediated neutralization, a combination of structural, biophysical, virologic data suggest that potent DENV-1 neutralization by E106 is coincident with bivalent engagement of the virus. Additionally, we determined the crystal structures of E111 bound to a novel fifth CC\u27 loop epitope on domain III: DIII) of the E protein from two different DENV-1 genotypes. The available atomic models of DENV virions revealed that the E111 epitope was inaccessible, suggesting that it recognizes an uncharacterized virus conformation. While the affinity of binding between E111 and DIII varied by genotype, we observed limited correlation with inhibitory activity. Instead, our results support the conclusion that potent neutralization depends on genotype-dependent exposure of the CC\u27 loop epitope. These findings establish new structural complexity of the DENV virion, which may be relevant for the choice of DENV strain for induction or analysis of neutralizing antibodies in the context of vaccine development